Advances in space-based observations have provided the capacity to develop regional- to global-scale estimates of evaporation, offering insights into this key component of the hydrological cycle. However, the evaluation of large-scale evaporation retrievals is not a straightforward task. While a number of studies have intercompared a range of these evaporation products by examining the variance amongst them, or by comparison of pixel-scale retrievals against ground-based observations, there is a need to explore more appropriate techniques to comprehensively evaluate remote-sensing-based estimates. One possible approach is to establish the level of product agreement between related hydrological components: for instance, how well do evaporation patterns and response match with precipitation or water storage changes? To assess the suitability of this consistency-based approach for evaluating evaporation products, we focused our investigation on four globally distributed basins in arid and semi-arid environments, comprising the Colorado River basin, Niger River basin, Aral Sea basin, and Lake Eyre basin. In an effort to assess retrieval quality, three satellite-based global evaporation products based on different methodologies and input data, including CSIRO-PML, the MODIS Global Evapotranspiration product (MOD16), and Global Land Evaporation: the Amsterdam Methodology (GLEAM), were evaluated against rainfall data from the Global Precipitation Climatology Project (GPCP) along with Gravity Recovery and Climate Experiment (GRACE) water storage anomalies. To ensure a fair comparison, we evaluated consistency using a degree correlation approach after transforming both evaporation and precipitation data into spherical harmonics. Overall we found no persistent hydrological consistency in these dryland environments. Indeed, the degree correlation showed oscillating values between periods of low and high water storage changes, with a phase difference of about 2–3 months. Interestingly, after imposing a simple lag in GRACE data to account for delayed surface runoff or baseflow components, an improved match in terms of degree correlation was observed in the Niger River basin. Significant improvements to the degree correlations (from ∼ 0 to about 0.6) were also found in the Colorado River basin for both the CSIRO-PML and GLEAM products, while MOD16 showed only half of that improvement. In other basins, the variability in the temporal pattern of degree correlations remained considerable and hindered any clear differentiation between the evaporation products. Even so, it was found that a constant lag of 2 months provided a better fit compared to other alternatives, including a zero lag. From a product assessment perspective, no significant or persistent advantage could be discerned across any of the three evaporation products in terms of a sustained hydrological consistency with precipitation and water storage anomaly data. As a result, our analysis has implications in terms of the confidence that can be placed in independent retrievals of the hydrological cycle, raises questions on inter-product quality, and highlights the need for additional techniques to evaluate large-scale products.

The agricultural sector in Saudi Arabia has witnessed rapid growth in both production and area under cultivation over the last few decades. This has prompted some concern over the state and future availability of fossil groundwater resources, which have been used to drive this expansion. Large-scale studies using satellite gravimetric data show a declining trend over this region. However, water management agencies require much more detailed information on both the spatial distribution of agricultural fields and their varying levels of water exploitation through time than coarse gravimetric data can provide. Relying on self-reporting from farm operators or sporadic data collection campaigns to obtain needed information are not feasible options, nor do they allow for retrospective assessments. In this work, a water accounting framework that combines satellite data, meteorological output from weather prediction models, and a modified land surface hydrology model was developed to provide information on both irrigated crop water use and groundwater abstraction rates. Results from the local scale, comprising several thousand individual center-pivot fields, were then used to quantify the regional-scale response. To do this, a semi-automated approach for the delineation of center-pivot fields using a multi-temporal statistical analysis of Landsat 8 data was developed. Next, actual crop evaporation rates were estimated using a two-source energy balance (TSEB) model driven by leaf area index, land surface temperature, and albedo, all of which were derived from Landsat 8. The Community Atmosphere Biosphere Land Exchange (CABLE) model was then adapted to use satellite-based vegetation and related surface variables and forced with a 3 km reanalysis dataset from the Weather Research and Forecasting (WRF) model. Groundwater abstraction rates were then inferred by estimating the irrigation supplied to each individual center pivot, which was determined via an optimization approach that considered CABLE-based estimates of evaporation and TSEB-based satellite estimates. The framework was applied over two study regions in Saudi Arabia: a small-scale experimental facility of around 40 center pivots in Al Kharj that was used for an initial evaluation and a much larger agricultural region in Al Jawf province comprising more than 5000 individual fields across an area exceeding 2500 km². Total groundwater abstraction for the year 2015 in Al Jawf was estimated at approximately 5.5 billion cubic meters, far exceeding any recharge to the groundwater system and further highlighting the need for a comprehensive water management strategy. Overall, this novel data–model fusion approach facilitates the compilation of national-scale groundwater abstractions while also detailing field-scale information that allows both farmers and water management agencies to make informed water accounting decisions across multiple spatial and temporal scales.

Agricultural intensification has resulted in the depletion of groundwater resources in many regions of the world. A prime example is Saudi Arabia, which witnessed dramatic agricultural expansion since the 1970s. To explore the influence of policy interventions aimed to better manage water resources, accurate information on the changes in the number and acreage of center-pivot fields is required. To quantify these metrics, we apply a hybrid machine learning framework, consisting of Density-Based Spatial Clustering of Applications with Noise, Convolutional Neural Networks, and Spectral Clustering, to the annual maximum Normalized Differential Vegetation Index maps obtained from Landsat imagery collected between 1990 to 2021. When evaluated against more than 28,000 manually delineated fields, the approach demonstrated producer’s accuracies ranging from 83.7% to 94.8% and user’s accuracies ranging from 90.2% to 97.9%. The coefficient of determination ($R^2$) between framework-delineated and manually delineated fields was higher than 0.97. Nationally, we found that most fields pre-dated 1990 (covering 8841 km${}^2$ in that year) and were primarily located within the central regions covering Hail, Qassim, Riyadh, and Wadi ad-Dawasir. A small decreasing trend in field acreage was observed for the period 1990–2010. However, by 2015, the acreage had increased to approximately 33,000 fields covering 9310 km${}^2$. While a maximum extent was achieved in 2016, recent decreases have seen levels return to pre-1990 levels. The gradual decrease between 1990 to 2010 was related to policy initiatives designed to phase-out wheat, while increases between 2010 to 2015 were linked to fodder crop expansion. There is evidence of an agricultural uptick starting in 2021, which is likely in response to global influences such as the COVID-19 pandemic or the conflict in Ukraine. Overall, this work offers the first detailed assessment of long-term agricultural development in Saudi Arabia, and provides important insights related to production metrics such as crop types, crop water consumption, and crop phenology and the overarching impacts of agricultural policy interventions.

Planning for use of a dune field aquifer for managed aquifer recharge (MAR) requires that hydraulic properties need to be estimated over a large geographic area. Saturated hydraulic conductivity of dune sands is commonly estimated from grain size distribution data by employing some type of empirical equation. Over 50 samples from the Wadi Khulays dune field in Western Saudi Arabia were collected and the grain size distribution, porosity, and hydraulic conductivity were measured. An evaluation of 20 existing empirical equations showed a generally high degree of error in the predicted compared to the measured hydraulic conductivity values of these samples. Statistical analyses comparing estimated versus measured hydraulic conductivity demonstrated that there is a significant relationship between hydraulic conductivity and mud percentage (and skewness). The modified Beyer equation, which showed a generally low prediction error, was modified by adding a second term fitting parameter related to the mud concentration based on 25 of the 50 samples analyzed. An inverse optimization process was conducted to quantify the fitting parameter and a new empirical equation was developed. This equation was tested against the remaining 25 samples analyzed and produced an estimated saturated hydraulic conductivity with the lowest error of any empirical equation. This methodology can be used for large dune field hydraulic conductivity estimation and reduce planning costs for MAR systems.

The agricultural sector in Saudi Arabia has witnessed rapid growth in both production and area under cultivation over the last few decades. This has prompted some concern over the state and future availability of fossil groundwater resources, which have been used to drive this expansion. Large-scale studies using satellite gravimetric data show a declining trend over this region. However, water management agencies require much more detailed information on both the spatial distribution of agricultural fields, and their varying levels of water exploitation through time, than coarse gravimetric data can provide. Relying on self-reporting from farm operators or sporadic data collection campaigns to obtain needed information are not feasible options, nor do they allow for retrospective assessments. In this work, a water accounting framework that combines satellite data, meteorological output from weather prediction models, and a modified land surface hydrology model, was developed to provide information on both irrigated crop-water use and groundwater abstraction rates. Results from the local-scale, comprising several thousand individual center-pivot fields, were then used to quantify the regional-scale response. To do this, a semi-automated approach for the delineation of center-pivot fields using a multi-temporal statistical analysis of Landsat 8 data was developed. Next, actual crop evaporation rates were estimated using a two-source energy balance (TSEB) model driven by leaf area index, land surface temperature, and albedo inputs, all of which were derived from Landsat 8. The Community Atmosphere Biosphere Land Exchange (CABLE) model was then adapted to use satellite-based vegetation and related surface variables, and forced with a 3 km reanalysis dataset from the Weather Research and Forecasting (WRF) model. Groundwater abstraction rates were then inferred by estimating the irrigation supplied to each individual center-pivot, which was determined via an optimization approach that considered CABLE-based estimates of evaporation and TSEB-based satellite estimates. The framework was applied over two study regions in Saudi Arabia: a small-scale experimental facility of around 40 center-pivots in Al Kharj that was used for an initial evaluation, and a much larger agricultural region in Al Jawf province comprising more than 5,000 individual fields across an area exceeding 2,500 km². Total groundwater abstraction for the year 2015 in Al Jawf were estimated at approximately 5.5 billion cubic meters, far exceeding any recharge to the groundwater system and further highlighting the need for a comprehensive water management strategy. Overall, this novel data-model fusion approach facilitates the compilation of national-scale groundwater abstractions, while also detailing field-scale information that allows both farmers and water management agencies to make informed water accounting decisions across multiple spatial and temporal scales.