Physical environmental processes, such as the evolution of precipitation or the diffusion of chemical clouds in the atmosphere, can be approximated by numerical models based on the underlying physics, e.g., for the purpose of prediction. As the modeling process is often very complex and resource demanding, such models are sometimes replaced by those that use historic and current data for calibration. For atmospheric (e.g., precipitation) or oceanographic (e.g., sea surface temperature) fields, the data-driven methods often concern the horizontal displacement driven by transport processes (called advection). These methods rely on flow fields estimated from images of the phenomenon by computer vision techniques, such as optical flow (OF). In this work, an algorithm is proposed for estimating the motion of spatio-temporal fields with the nodes of a geosensor network (GSN) deployed in situ when images are not available. The approach adapts a well-known raster-based OF algorithm to the specifics of GSNs, especially to the spatial irregularity of data. In this paper, the previously introduced approach has been further developed by introducing an error model that derives probabilistic error measures based on spatial node configuration. Further, a more generic motion model is provided, as well as comprehensive simulations that illustrate the performance of the algorithm in different conditions (fields, motion behaviors, node densities and deployments) for the two error measures of motion direction and motion speed. Finally, the algorithm is applied to data sampled from weather radar images, and the algorithm performance is compared to that of a state-of-the-art OF algorithm applied to the weather radar images directly, as often done in nowcasting.