Search

As the global demands for the use of freshwater resources continues to rise, it has become increasingly important to insure the sustainability of this resources. This is accomplished through the use of management strategies that often utilize monitoring and the use of hydrological models. However, monitoring at large scales is not feasible and therefore model applications are becoming challenging, especially when spatially distributed datasets, such as evapotranspiration, are needed to understand the model performances. Due to these limitations, most of the hydrological models are only calibrated for data obtained from site/point observations, such as streamflow. Therefore, the main focus of this paper is to examine whether the incorporation of remotely sensed and spatially distributed datasets can improve the overall performance of the model. In this study, actual evapotranspiration (ETa) data was obtained from the two different sets of satellite based remote sensing data. One dataset estimates ETa based on the Simplified Surface Energy Balance (SSEBop) model while the other one estimates ETa based on the Atmosphere-Land Exchange Inverse (ALEXI) model. The hydrological model used in this study is the Soil and Water Assessment Tool (SWAT), which was calibrated against spatially distributed ETa and single point streamflow records for the Honeyoey Creek-Pine Creek Watershed, located in Michigan, USA. Two different techniques, multi-variable and genetic algorithm, were used to calibrate the SWAT model. Using the aforementioned datasets, the performance of the hydrological model in estimating ETa was improved using both calibration techniques by achieving Nash-Sutcliffe efficiency (NSE) values >0.5 (0.73–0.85), percent bias (PBIAS) values within ±25% (±21.73%), and root mean squared error – observations standard deviation ratio (RSR) values

Indicators based on evapotranspiration (ET) provide useful information about surface water status, response of vegetation to drought stress, and potential growth limitations. The capability of ET-based indicators, including actual ET and the evaporative stress index (ESI), to predict crop yields of spring barley and winter wheat was analyzed for 33 districts of the Czech Republic. In this study, the ET-based indicators were computed using two different approaches: (i) a prognostic model, SoilClim, which computes the water balance based on ground weather observations and information about soil and land cover; (ii) the diagnostic Atmosphere–Land Exchange Inverse (ALEXI) model based primarily on remotely sensed land surface temperature data. The capability of both sets of indicators to predict yields of spring barley and winter wheat was tested using artificial neural networks (ANNs) applied to the adjusting number and timeframe of inputs during the growing season. Yield predictions based on ANNs were computed for both crops for all districts together, as well as for individual districts. The root mean square error (RMSE) and coefficient of determination (R2) between observed and predicted yields varied with date within the growing season and with the number of ANN inputs used for yield prediction. The period with the highest predictive capability started from early-June to mid-June. This optimal period for yield prediction was identifiable already at the lower number of ANN inputs, nevertheless, the accuracy of the prediction improved as more inputs were included within ANNs.The RMSE values for individual districts varied between 0.4 and 0.7 t ha–1 while R2 reached values of 0.5–0.8 during the optimal period. Results of the study demonstrated that ET-based indicators can be used for yield prediction in real time during the growing season and therefore have great potential for decision making at regional and district levels.

We studied foraging patterns of Raramuri Criollo (heritage breed, “RC”) and Angus x Hereford (commercial beef breed, “AH”) cows in relation to ambient heat and body temperature during summer (2016 and 2017) and winter (2017 and 2018) in the Chihuahuan Desert. Cows of each breed grazed separately in two adjacent pastures (~1100 ha) in a crossover design for four weeks in each season/year. Animals were fitted with temperature loggers attached to blank CIDRs (Controlled Internal Drug Release device) devoid of hormones that recorded body temperature (BodyT), and GPS collars that recorded position and ambient temperature (CollarT). All sensor data were logged at 10 min intervals. A landscape thermal map (LandT) was developed from Landsat satellite imagery for habitat analysis using GPS locations chosen by individual collared cows, and air temperature (AirT) was recorded by a nearby weather station. Data were analyzed within four daytime segments: dawn (sunrise – 9AM); pre-noon (9AM – noon); post-noon (noon – 3PM); and dusk (3PM – sunset). ANOVA was used to determine whether BodyT, CollarT, LandT selection, or animal movement variables within each of the four daily segments differed (P

Daily continuous evapotranspiration (ET) estimates of 1 km spatial resolution can benefit agricultural water resources management at regional scales. Thermal infrared remote sensing-derived land surface temperature (LST) is a critical variable for ET estimation using energy balance-based models. However, missing LST information under cloudy conditions remains a long-standing barrier for spatiotemporally continuous monitoring of daily ET at regional scales. In this study, LST data of 1 km spatial resolution at 11:00 local solar time under all-weather conditions across the North China Plain (NCP) were first generated using a data fusion approach developed previously. Second, combined with the generated LST data, MODIS products, and meteorological forcing from the China Land Data Assimilation System, the Two-Source Energy Balance model (TSEB) and a temporal upscaling method were jointly used to estimate daily ET at 1 km spatial resolution across the NCP for a decade from 2008 to 2017. In particular, to better incorporate the impact of crop phenology on ET and improve the ET estimation, the fraction of greenness in TSEB was determined in terms of a leaf area index threshold during the crop growth period. Compared with observed instantaneous latent heat flux (LE) corrected for energy balance closure, the estimated LE reasonably captures inter- and intra-annual variations in LE measured at the Huailai, Daxing, Weishan, and Guantao flux towers, with R2 of 0.63–0.79. Estimated daily ET against in situ ET measurements with energy balance closure at the Huailai, Daxing, and Guantao sites showed good performance in terms of R2 greater than 0.70 and RMSE below 0.91 mm/d. These accuracies are comparable with published results, with our ET data set validated by many more observations than previous studies and featuring spatiotemporal continuity and high spatial resolution across the entire NCP for a decade. Furthermore, seasonal ET variations reflected by our product outperform two widely used global products in capturing water consumption characteristics in the winter wheat-summer maize rotation system. In terms of temporal trends, annual ET estimates across the NCP show a decreasing and then increasing trend over the past decade, which is attributed to the increased cropping intensity over the recent years reflected by an increase in leaf area index.

To effectively meet growing food demands, the global agronomic community will require a better understanding of factors that are currently limiting crop yields and where production can be viably expanded with minimal environmental consequences. Remote sensing can inform these analyses, providing valuable spatiotemporal information about yield-limiting moisture conditions and crop response under current climate conditions. In this paper we study correlations for the period 2003–2013 between yield estimates for major crops grown in Brazil and the Evaporative Stress Index (ESI) – an indicator of agricultural drought that describes anomalies in the actual/reference evapotranspiration (ET) ratio, retrieved using remotely sensed inputs of land surface temperature (LST) and leaf area index (LAI). The strength and timing of peak ESI-yield correlations are compared with results using remotely sensed anomalies in water supply (rainfall from the Tropical Rainfall Mapping Mission; TRMM) and biomass accumulation (LAI from the Moderate Resolution Imaging Spectroradiometer; MODIS). Correlation patterns were generally similar between all indices, both spatially and temporally, with the strongest correlations found in the south and northeast where severe flash droughts have occurred over the past decade, and where yield variability was the highest. Peak correlations tended to occur during sensitive crop growth stages. At the state scale, the ESI provided higher yield correlations for most crops and regions in comparison with TRMM and LAI anomalies. Using finer scale yield estimates reported at the municipality level, ESI correlations with soybean yields peaked higher and earlier by 10 to 25 days in comparison to TRMM and LAI, respectively. In most states, TRMM peak correlations were marginally higher on average with municipality-level annual corn yield estimates, although these estimates do not distinguish between primary and late season harvests. A notable exception occurred in the northeastern state of Bahia, where the ESI better captured effects of rapid cycling of moisture conditions on corn yields during a series of flash drought events. The results demonstrate that for monitoring agricultural drought in Brazil, value is added by combining LAI with LST indicators within a physically based model of crop water use.

Flash droughts are characterized by a period of unusually rapid drought intensification over sub-seasonal time scales that often take vulnerable stakeholders by surprise given their rapid onset. Various studies have shown that flash drought is more likely to develop when extreme weather conditions persist over the same region for several weeks or longer. Though precipitation deficits over some period of time are a prerequisite for drought, their presence alone is unlikely to lead to flash drought because a lack of precipitation is only one of several factors that contribute to rapid drought development. When below normal precipitation occurs alongside other extreme weather anomalies such as intense heat that enhance atmospheric evaporative demand, their co-occurrence can lead to a rapid depletion of root zone soil moisture content due to increased evapotranspiration. This in turn can lead to a rapid increase in vegetation moisture stress and the onset of flash drought conditions.Several recent studies have used quantitative definitions based on rapid changes in a given drought monitoring dataset to identify flash droughts in the climatological record. Here, we build upon these recent studies by developing a new flash drought intensity index that accounts not only for their rapid rate of intensification, but also for how severe the drought conditions become during and after the period of rapid intensification. The method includes two components that together capture the suddenness of flash drought development (faster intensification corresponds to a more severe flash drought) and the actual drought severity after the rapid intensification period ends (severe drought conditions lasting for a longer period correspond to a more severe flash drought). The motivation behind this method is the desire to account for both the “flash” and “drought” aspects of flash drought because both of these characteristics influence how people view flash droughts. Thus, a metric that considers both of these aspects provides a more comprehensive assessment of flash drought intensity and its impacts on the environment. In this talk, we will present the proposed flash drought intensity index methodology, along with results from individual case studies and a 40-year climatology to illustrate its use.

Leaf Area Index (LAI) is a fundamental vegetation biophysical variable serving as an essential input to many land surface and atmospheric models. Long-term LAI maps are typically generated with satellite images at moderate spatial resolution (0.25 to 1 km), such as those from the Moderate Resolution Imaging Spectroradiometer (MODIS). While useful for regional-scale land surface modeling, these moderate resolution products often cannot resolve spatial heterogeneity important for many agricultural and hydrological applications. This paper proposes an approach to map LAI at 30-m resolution based on Landsat images for the Contiguous US (CONUS) consistent with the MODIS product, aimed at multi-scale modeling applications. The algorithm was driven by 1.6 million spatially homogeneous samples derived from MODIS LAI and Landsat surface reflectance products from 2006 to 2018. Based on these samples, we trained separate random forest models to estimate LAI from Landsat surface reflectance for eight biomes of the National Land Cover Database (NLCD). A balanced sample design regarding the saturation status of MODIS LAI and a machine-learning-based noise detection technique were introduced to mitigate the trade-off in estimation accuracy between medium LAI (e.g., 3 to 4, unsaturated) and high LAI (e.g., 4–6, saturated). This approach was evaluated using ground measurements from 19 National Ecological Observatory Network (NEON) sites and eight independent sites from other sources. These sites comprise a representative sample of forests, grasslands, shrublands, and croplands across the US. For NEON sites, the LAI estimates show an overall Root Mean Squared Error (RMSE) of 0.8 with r2 of 0.88. For the eight independent sites, the Landsat LAI algorithm achieves RMSE between 0.52 and 0.91. The uncertainty in Landsat estimated LAI varies across biomes and locations. The proposed algorithm was implemented on the Google Earth Engine platform, allowing for the rapid generation of long-term high-resolution LAI records from the 1980s using Landsat images (code is available at https://github.com/yanghuikang/Landsat-LAI). Our findings also highlight the importance of sample balance on regression-based modeling in remote sensing applications.

Severe droughts in the Southern Great Plains (SGP: Kansas, Oklahoma, and Texas) in recent years have reduced the productivity of tallgrass prairie and resulted in substantial economic losses to the beef cattle industry in this region. Understanding spatial and temporal patterns of agricultural drought in the SGP can help ranchers to develop and implement drought mitigation strategies. In this study, the Land Surface Water Index (LSWI), calculated from the Moderate Resolution Imaging Spectroradiometer (MODIS) near infrared and shortwave infrared bands, was used to assess agricultural drought in the tallgrass prairie region of the SGP during 2000–2013. The number of consecutive days with LSWI

Observations of land surface temperature (LST) are crucial for the monitoring of surface energy fluxes from satellite. Methods that require high temporal resolution LST observations (e.g., from geostationary orbit) can be difficult to apply globally because several geostationary sensors are required to attain near-global coverage (60°N to 60°S). While these LST observations are available from polar-orbiting sensors, providing global coverage at higher spatial resolutions, the temporal sampling (twice daily observations) can pose significant limitations. For example, the Atmosphere Land Exchange Inverse (ALEXI) surface energy balance model, used for monitoring evapotranspiration and drought, requires an observation of the morning change in LST - a quantity not directly observable from polar-orbiting sensors. Therefore, we have developed and evaluated a data-mining approach to estimate the mid-morning rise in LST from a single sensor (2 observations per day) of LST from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor on the Aqua platform. In general, the data-mining approach produced estimates with low relative error (5 to 10%) and statistically significant correlations when compared against geostationary observations. This approach will facilitate global, near real-time applications of ALEXI at higher spatial and temporal coverage from a single sensor than currently achievable with current geostationary datasets.

Large-scale crop monitoring and yield estimation are important for both scientific research and practical applications. Satellite remote sensing provides an effective means for regional and global cropland monitoring, particularly in data-sparse regions that lack reliable ground observations and reporting. The conventional approach of using visible and near-infrared based vegetation index (VI) observations has prevailed for decades since the onset of the global satellite era. However, other satellite data encompass diverse spectral ranges that may contain complementary information on crop growth and yield, but have been largely understudied and underused. Here we conducted one of the first attempts at synergizing multiple satellite data spanning a diverse spectral range, including visible, near-infrared, thermal and microwave, into one framework to estimate crop yield for the U.S. Corn Belt, one of the world's most important food baskets. Specifically, we included MODIS Enhanced VI (EVI), estimated Gross Primary Production based on GOME-2 solar-induced fluorescence (SIF-GPP), thermal-based ALEXI Evapotranspiration (ET), QuikSCAT Ku-band radar backscatter, and AMSR-E X-band passive microwave Vegetation Optical Depth (VOD) in this study, benchmarked on USDA county-level crop yield statistics. We used Partial Least Square Regression (PLSR), an effective statistical model for dimension reduction, to distinguish commonly shared and unique individual information from the various satellite data and other ancillary climate information for crop yield estimation. In the PLSR model that includes all of the satellite data and climate variables from 2007 to 2009, we assessed the first two major PLSR components and found that the first component (an integrated proxy of crop aboveground biomass) explained 82% variability of modelled crop yield, and the second component (dominated by environmental stresses) explained 15% variability of modelled crop yield. We found that most of the satellite derived metrics (e.g. SIF-GPP, radar backscatter, EVI, VOD, ALEXI-ET) share common information related to aboveground crop biomass (i.e. the first component). For this shared information, the SIF-GPP and backscatter data contain almost the same amount of information as EVI at the county scale. When removing the above shared component from all of the satellite data, we found that EVI and SIF-GPP do not provide much extra information; instead, Ku-band backscatter, thermal-based ALEXI-ET, and X-band VOD provide unique information on environmental stresses that improves overall crop yield predictive skill. In particular, Ku-band backscatter and associated differences between morning and afternoon overpasses contribute unique information on crop growth and environmental stress. Overall, using satellite data from various spectral bands significantly improves regional crop yield predictions. The additional use of ancillary climate data (e.g. precipitation and temperature) further improves model skill, in part because the crop reproductive stage related to harvest index is highly sensitive to environmental stresses but they are not fully captured by the satellite data used in our study. We conclude that using satellite data across various spectral ranges can improve monitoring of large-scale crop growth and yield beyond what can be achieved from individual sensors. These results also inform the synergistic use and development of current and next generation satellite missions, including NASA ECOSTRESS, SMAP, and OCO-2, for agricultural applications.

Land surface temperature (LST) is a critical parameter in environmental studies and resource management. The MODIS LST data product has been widely used in various studies, such as drought monitoring, evapotranspiration mapping, soil moisture estimation and forest fire detection. However, cloud contamination affects thermal band observations and will lead to inconsistent LST results. In this study, we present a new Remotely Sensed DAily land Surface Temperature reconstruction (RSDAST) model that recovers clear sky LST for pixels covered by cloud using only clear-sky neighboring pixels from nearby dates. The reconstructed LST was validated using the original LST pixels. Model shows high accuracy for reconstructing one masked pixel with R 2 of 0.995, bias of −0.02 K and RMSE of 0.51 K. Extended spatial reconstruction results show a better accuracy for flat areas with R 2 of 0.72‒0.89, bias of −0.02–0.21 K, and RMSE of 0.92–1.16 K, and for mountain areas with R 2 of 0.81–0.89, bias of −0.35–−1.52 K, and RMSE of 1.42‒2.24 K. The reconstructed areas show spatial and temporal patterns that are consistent with the clear neighbor areas. In the reconstructed LST and NDVI triangle feature space which is controlled by soil moisture, LST values distributed reasonably and correspond well to the real soil moisture conditions. Our approach shows great potential for reconstructing clear sky LST under cloudy conditions and provides consistent daily LST which are critical for daily drought monitoring.

The U.S. Drought Monitor (USDM) classifies drought into five discrete dryness/drought categories based on expert synthesis of numerous data sources. In this study, an empirical methodology is presented for creating a nondiscrete USDM index that simultaneously 1) represents the dryness/wetness value on a continuum and 2) is most consistent with the time scales and processes of the actual USDM. A continuous USDM representation will facilitate USDM forecasting methods, which will benefit from knowledge of where, within a discrete drought class, the current drought state most probably lies. The continuous USDM is developed such that the actual discrete USDM can be reconstructed by discretizing the continuous USDM based on the 30th, 20th, 10th, 5th, and 2nd percentiles—corresponding with USDM definitions for the D4–D0 drought classes. Anomalies in precipitation, soil moisture, and evapotranspiration over a range of different time scales are used as predictors to estimate the continuous USDM. The methodology is fundamentally probabilistic, meaning that the probability density function (PDF) of the continuous USDM is estimated and therefore the degree of uncertainty in the fit is properly characterized. Goodness-of-fit metrics and direct comparisons between the actual and predicted USDM analyses during different seasons and years indicate that this objective drought classification method is well correlated with the current USDM analyses. In Part II, this continuous USDM index will be used to improve intraseasonal USDM intensification forecasts because it is capable of distinguishing between USDM states that are either far from or near to the next-higher drought category.

The health of the Chesapeake Bay ecosystem has been declining for several decades due to high levels of nutrients and sediments largely tied to agricultural production systems. Therefore, monitoring of agricultural water use and hydrologic connections between crop lands and Bay tributaries has received increasing attention. Remote sensing retrievals of actual evapotranspiration (ET) can provide valuable information in support of these hydrologic modeling efforts, spatially and temporally describing consumptive water use by crops and natural vegetation and quantifying response to expansion of irrigated area occurring with Bay watershed. In this study, a multi-sensor satellite data fusion methodology, combined with a multi-scale ET retrieval algorithm, was applied over the Choptank River watershed located within the Lower Chesapeake Bay region on the Eastern Shore of Maryland, USA to produce daily 30-m resolution ET maps. ET estimates directly retrieved on Landsat satellite overpass dates have high accuracy with relative error (RE) of 9%, as evaluated using flux tower measurements. The fused daily ET time series have reasonable errors of 18% at the daily time step - an improvement from 27% errors using standard Landsat-only interpolation techniques. Annual water consumption by different land cover types was assessed, showing reasonable distributions of water use with cover class. Seasonal patterns in modeled crop transpiration and soil evaporation for dominant crop types were analyzed, and agree well with crop phenology at field scale. Additionally, effects of irrigation occurring during a period of rainfall shortage were captured by the fusion program. These results suggest that the ET fusion system will have utility for water management at field and regional scales over the Eastern Shore. Further efforts are underway to integrate these detailed water use datasets into watershed-scale hydrologic models to improve assessments of water quality and inform best management practices to reduce nutrient and sediment loads to the Chesapeake Bay.

Probabilistic forecasts of U.S. Drought Monitor (USDM) intensification over 2-, 4-, and 8-week time periods are developed based on recent anomalies in precipitation, evapotranspiration, and soil moisture. These statistical forecasts are computed using logistic regression with cross validation. While recent precipitation, evapotranspiration, and soil moisture do provide skillful forecasts, it is found that additional information on the current state of the USDM adds significant skill to the forecasts. The USDM state information takes the form of a metric that quantifies the “distance” from the next-higher drought category using a nondiscrete estimate of the current USDM state. This adds skill because USDM states that are close to the next-higher drought category are more likely to intensify than states that are farther from this threshold. The method shows skill over most of the United States but is most skillful over the north-central United States, where the cross-validated Brier skill score averages 0.20 for both 2- and 4-week forecasts. The 8-week forecasts are less skillful in most locations. The 2- and 4-week probabilities have very good reliability. The 8-week probabilities, on the other hand, are noticeably overconfident. For individual drought events, the method shows the most skill when forecasting high-amplitude flash droughts and when large regions of the United States are experiencing intensifying drought.

The annual flood cycle of the Sudd wetland in South Sudan plays an important role in the Nile River Basin water balance. The wetland, however, is extensive and sparsely instrumented, which has inhibited credible understanding of regional flooding across space and time. Here we explore the potential to apply low resolution C-band ENVISAT Advanced SAR imagery for remote estimation of Sudd flooded area. Over a five year study period (2007–2011) the time-averaged flooded Sudd area was found to be 18,033 km 2 with an average annual high of 29,702 km 2 in late September and a low of 10,128 km 2 in early May. Annual peak flood area ranges considerably from 19,259 km 2 in 2009 to 36,649 km 2 in 2007, but we found no systematic trend over the five year study period. Derived flood frequency maps identify areas of open water and permanent flooding (12% of total area), seasonal flooding (29%), and intermittent flooding (48%). To evaluate the certainty of our results, we consider their consistency with (1) prior studies, (2) evapotranspiration estimates from the Atmosphere-Land Exchange Inverse (ALEXI) surface energy balance algorithm, (3) watershed storage anomaly estimates from GRACE, (4) supervised classification of open water area using Landsat, and (5) a rough measure of water availability (antecedent precipitation). The analyses show reasonable temporal and spatial consistency with available lines of evidence. We conclude that low resolution C-band SAR imagery shows promise for study of Sudd wetland flood dynamics.

This paper presents an intercomparative study of 12 operationally produced large-scale datasets describing soil moisture, evapotranspiration (ET), and/or vegetation characteristics within agricultural regions of the contiguous United States (CONUS). These datasets have been developed using a variety of techniques, including, hydrologic modeling, satellite-based retrievals, data assimilation, and survey/in-field data collection. The objectives are to assess the relative utility of each dataset for monitoring crop yield variability, to quantitatively assess their capacity for predicting end-of-season corn and soybean yields, and to examine the evolution of the yield-index correlations during the growing season. This analysis is unique both with regards to the number and variety of examined yield predictor datasets and the detailed assessment of the water availability timing on the end-of-season crop production during the growing season. Correlation results indicate that over CONUS, at state-level soil moisture and ET indices can provide better information for forecasting corn and soybean yields than vegetation-based indices such as normalized difference vegetation index. The strength of correlation with corn and soybean yields strongly depends on the interannual variability in yield measured at a given location. In this case study, some of the remotely derived datasets examined provide skill comparable to that of in-situ field survey-based data—further demonstrating the utility of these remote sensing-based approaches for estimating crop yield.

The fate of the terrestrial biosphere is highly uncertain given recent and projected changes in climate. This is especially acute for impacts associated with changes in drought frequency and intensity on the distribution and timing of water availability. The development of effective adaptation strategies for these emerging threats to food and water security are compromised by limitations in our understanding of how natural and managed ecosystems are responding to changing hydrological and climatological regimes. This information gap is exacerbated by insufficient monitoring capabilities from local to global scales. Here, we describe how evapotranspiration (ET) represents the key variable in linking ecosystem functioning, carbon and climate feedbacks, agricultural management, and water resources, and highlight both the outstanding science and applications questions and the actions, especially from a space-based perspective, necessary to advance them.

Due to their shallow vertical support, remotely sensed surface soil moisture retrievals are commonly regarded as being of limited value for water budget applications requiring the characterization of temporal variations in total terrestrial water storage (dS ∕ dt). However, advances in our ability to estimate evapotranspiration remotely now allow for the direct evaluation of approaches for quantifying dS ∕ dt via water budget closure considerations. By applying an annual water budget analysis within a series of medium-scale (2000–10 000 km2) basins within the United States, we demonstrate that, despite their clear theoretical limitations, surface soil moisture retrievals derived from passive microwave remote sensing contain statistically significant information concerning dS ∕ dt. This suggests the possibility of using (relatively) higher-resolution microwave remote sensing products to enhance the spatial resolution of dS ∕ dt estimates acquired from gravity remote sensing.

Soil drainage is a widely used agricultural practice in the midwest USA to remove excess soil water to potentially improve the crop yield. Research shows an increasing trend in baseflow and streamflow in the midwest over the last 60 years, which may be related to artificial drainage. Subsurface drainage (i.e., tile) in particular may have strongly contributed to the increase in these flows, because of its extensive use and recent gain in the popularity as a yield-enhancement practice. However, how evapotranspiration (ET) is impacted by tile drainage on a regional level is not well-documented. To explore spatial and temporal ET patterns and their relationship to tile drainage, we applied an energy balance-based multisensor data fusion method to estimate daily 30-m ET over an intensively tile-drained area in South Dakota, USA, from 2005 to 2013. Results suggest that tile drainage slightly decreases the annual cumulative ET, particularly during the early growing season. However, higher mid-season crop water use suppresses the extent of the decrease of the annual cumulative ET that might be anticipated from widespread drainage. The regional water balance analysis during the growing season demonstrates good closure, with the average residual from 2005 to 2012 as low as -3 mm. As an independent check of the simulated ET at the regional scale, the water balance analysis lends additional confidence to the study. The results of this study improve our understanding of the influence of agricultural drainage practices on regional ET, and can affect future decision making regarding tile drainage systems.

The ability to regionally monitor crop progress and condition through the growing season benefits both crop management and yield estimation. In the United States, these metrics are reported weekly at state or district (multiple counties) levels by the U.S. Department of Agriculture (USDA) National Agricultural Statistics Service (NASS) using field observations provided by trained local reporters. However, the ground data collection process supporting this effort is time consuming and subjective. Furthermore, operational crop management and yield estimation efforts require information with more granularity than at the state or district level. This paper evaluates remote sensing approaches for mapping crop phenology using vegetation index time-series generated by fusing Landsat and MODIS (Moderate Resolution Imaging Spectroradiometer) surface reflectance imagery to improve temporal sampling over that provided by Landsat alone. The case study focuses on an agricultural region in central Iowa from 2001 to 2014. Our objectives are 1) to assess Landsat-MODIS data fusion results over cropland; 2) to map crop phenology at 30 m resolution using fused surface reflectance data; and 3) to identify the relationships between remotely sensed crop phenology metrics and the crop progress stages reported by NASS. The results show that detailed spatial and temporal variability in vegetation development across this landscape can be identified using the fused Landsat-MODIS data. The mean difference (bias) in Normalized Difference Vegetation Index (NDVI) between actual Landsat observations and the fused Landsat-MODIS data, generated for Landsat overpass dates, is in the range of − 0.011 to 0.028 for every year. The derived phenological metrics show distinct features for different crops and natural vegetation at field scales. Strong correlations are observed between remotely sensed phenological stages, based on NDVI curve inflection points, and the observed crop physiological growth stages from the NASS Crop Progress (CP) reports. The green-up dates detected from remote sensing data typically occurred during crop vegetative stages when 2–4 leaves were developed for both corn and soybeans, or about 1–3 weeks after the reported emergence dates when the plant were first visible to ground-based observers. Despite being a lagging indicator, remotely sensed green-up can be used effectively to backcast emergence, e.g. as input to spatially distributed crop models. The differences in green-up date between corn and soybean were 8–10 days, consistent with the offset in emergence dates reported by NASS at district level. The reported harvest dates were typically about 2–3 weeks after the dormancy stage was detected via remote sensing for corn and about 1–2 weeks for soybeans. This suggests that probable harvest times for individual fields may be predicted 1–3 weeks ahead using remote sensing data. The results suggest that crop phenology and certain growth stages at field scales (30 m spatial resolution) can be linked and mapped by integrating imagery from multiple remote sensing platforms.

Evapotranspiration (ET) is a major component of the global and regional water cycle. An operational Geostationary Operational Environmental Satellite (GOES) ET and Drought (GET-D) product system has been developed by the National Environmental Satellite, Data and Information Service (NESDIS) in the National Oceanic and Atmospheric Administration (NOAA) for numerical weather prediction model validation, data assimilation, and drought monitoring. GET-D system was generating ET and Evaporative Stress Index (ESI) maps at 8 km spatial resolution using thermal observations of the Imagers on GOES-13 and GOES-15 before the primary operational GOES satellites transitioned to GOES-16 and GOES-17 with the Advanced Baseline Imagers (ABI). In this study, the GET-D product system is upgraded to ingest the thermal observations of ABI with the best spatial resolution of 2 km. The core of the GET-D system is the Atmosphere-Land Exchange Inversion (ALEXI) model, which exploits the mid-morning rise in the land surface temperature to deduce the land surface fluxes including ET. Satellite-based land surface temperature and solar insolation retrievals from ABI and meteorological forcing from NOAA NCEP Climate Forecast System (CFS) are the major inputs to the GET-D system. Ancillary data required in GET-D include land cover map, leaf area index, albedo and cloud mask. This paper presents preliminary results of ET from the upgraded GET-D system after a brief introduction of the ALEXI model and the architecture of GET-D system. Comparisons with in situ ET measurements showed that the accuracy of the GOES-16 ABI based ET is similar to the results from the legacy GET-D ET based on GOES-13/15 Imager data. The agreement with the in situ measurements is satisfactory with a correlation of 0.914 averaged from three Mead sites. Further evaluation of the ABI-based ET product, upgrade efforts of the GET-D system for ESI products, and conclusions for the ABI-based GET-D products are discussed.

In viticulture, deficit irrigation strategies are often implemented to control vine canopy growth and to impose stress at critical stages of vine growth to improve wine grape quality. To support deficit irrigation scheduling, remote sensing technologies can be employed in the mapping of evapotranspiration (ET) at the field to sub-field scales, quantifying time-varying vineyard water requirements and actual water use. In the current study, we investigate the utility of ET maps derived from thermal infrared satellite imagery over a vineyard in the Central Valley of California equipped with a variable rate drip irrigation (VRDI) system which enables differential water applications at the 30 × 30 m scale. To support irrigation management at that scale, we utilized a thermal-based multi-sensor data fusion approach to generate weekly total actual ET (ETa) estimates at 30 m spatial resolution, coinciding with the resolution of the Landsat reflectance bands. Crop water requirements (ETc) were defined with a vegetative index (VI)-based approach. To test capacity to capture stress signals, the vineyard was sub-divided into four blocks with different irrigation management strategies and goals, inducing varying degrees of stress during the growing season. Results indicate derived weekly total ET from the thermal-based data fusion approach match well with observations. The thermal-based method was also able to capture the spatial heterogeneity in ET over the vineyard due to a water stress event imposed on two of the four vineyard blocks. This transient stress event was not reflected in the VI-based ETc estimate, highlighting the value of thermal band imaging. While the data fusion system provided valuable information, latency in current satellite data availability, particularly from Landsat, impacts operational applications over the course of a growing season.

Forests provide many important service functions, including habitat for wildlife, timber production and watershed water regulation. Short-term drought can make forests more susceptible to wildfire and insect attack, while long-term drought can directly increase forest mortality. Evapotranspiration (ET) is a key parameter that links the hydrological and ecological processes. Studying the impact of drought on ET, especially at stand-scale, can provide useful information for forest management. In this study, we applied a multi-scale data fusion ET modeling method using remotely sensed data to estimate daily 30 m ET over a natural forest in central US from 2010 to 2012, with 2012 as an extreme drought year. The estimated ET agrees well with the observed ET. Drought impact on ET is further analyzed and demonstrates the value of remotely sensed ET in studying drought impact on forest water use.

Accurate estimation of crop yield before harvest is critical for sustaining agricultural markets and ensuring food security. Remote sensing data have been demonstrated as a useful tool for estimating crop yields. Previous studies show that high spatial and temporal resolution vegetation index from remote sensing data fusion better explained spatial variability of yield in central Iowa. This paper extends the study area over 10 major agricultural states in the United States and further examines the added value by using Landsat-8, Sentinel-2 and MODIS images from 2016 and 2017. The large area data processing with Landsat and Sentinel-2 has been enabled by using the Google Earth Engine (GEE) technology. Results show that the combination of Landsat-8 and Sentinel-2 provides spatial variability of yield at field scales and explains yield variability better than using single data source alone even at the state-level.

Agitation in persons with dementia (PWD) poses major health risks both for themselves and for their caregivers. Passive sensing based continuous behavior tracking can prevent the escalation of such episodes. But, predicting such behavior from sensor streams, especially in real-world residential settings, is still an active area of research. Major challenges include the sparsity, unpredictability, and variations in such behavior, as well as the "weak" annotations from real-world participants. This paper proposes a novel approach to overcome these issues in predicting agitation episodes from the PWD's wrist motion data. In a transdisciplinary study on dementia dyads residing in their homes, the PWD motion is continuously sensed from their smart watch inertial sensors, while agitation episodes are actively marked by the caregivers. The data from 10 residential deployments, each with 30 days duration, are analyzed in this paper, and multiple-instance learning (MIL) based models are implemented to learn from such sparse and weakly annotated data. These models are compared with single-instance models in predicting the agitated behavior. The results show the potential of MIL models in sparsely labeled behavior inference from wearables in-the-wild.

The Arctic has become generally a warmer place over the past decades leading to earlier snow melt, permafrost degradation and changing plant communities. Increases in precipitation and local evaporation in the Arctic, known as the acceleration components of the hydrologic cycle, coupled with land cover changes, have resulted in significant changes in the regional surface energy budget. An extensive local evaluation of the two-source energy balance model (TSEB) - a remote sensing-based model using thermal infrared retrievals of land-surface temperature - was performed using tower measurements collected over two flux towers in Greenland’s tundra types and two flux tower sites in Alaska’s boreal forest in all sky conditions in summer. Evaluation results showed a mean turbulent flux errors generally less than 50 W•m-2 at half-hourly timesteps, similar to errors typically reported in surface energy balance modelling studies conducted in more temperate climatic regimes.

Monitoring the effects of water availability on vegetation globally using satellites is important for applications such as drought early warning, precision agriculture, and food security as well as for more broadly understanding relationships between water and carbon cycles. In this global study, we examine how quickly several satellite-based indicators, assumed to have relationships with water availability, respond, on timescales of days to weeks, in comparison with variations in root-zone soil moisture (RZM) that extends to about 1 m depth. The satellite indicators considered are the normalized difference vegetation and infrared indices (NDVI and NDII, respectively) derived from reflectances obtained with moderately wide (20–40 nm) spectral bands in the visible and near-infrared (NIR) and evapotranspiration (ET) estimated from thermal infrared observations and normalized by a reference ET. NDVI is primarily sensitive to chlorophyll contributions and vegetation structure while NDII may contain additional information on water content in leaves and canopy. ET includes both the loss of root zone soil water through transpiration (modulated by stomatal conductance) as well as evaporation from bare soil. We find that variations of these satellite-based drought indicators on time scales of days to weeks have significant correlations with those of RZM in the same water-limited geographical locations that are dominated by grasslands, shrublands, and savannas whose root systems are generally contained within the 1 m RZM layer. Normalized ET interannual variations show generally a faster response to water deficits and enhancements as compared with those of NDVI and NDII, particularly in sparsely vegetated regions.

A newly developed microwave (MW) land surface temperature (LST) product is used to substitute thermal infrared (TIR) based LST in the two-source energy balance approach (TSEB) for estimating ET from space. This TSEB land surface scheme, the Atmosphere Land Exchange Inverse (ALEXI) model framework, is an approach that minimizes sensitivity to absolute biases in input records of LST through the analysis of the rate of temperature change in the morning. This experiment is therefore an important test of the ability to retrieve diurnal temperature information from a constellation of satellites with microwave radiometers that together provide 6–8 observations of Ka-band brightness temperature per location per day. This represents the first ever attempt at a global implementation of ALEXI with MW-based LST and is intended as the first step towards providing all-weather capability to the ALEXI framework. The leveraging of all sky capability of MW sensors is the main motivation of this work, as TIR-based ALEXI is limited to clear sky conditions. The analysis is based on a 9-year long record of ALEXI ET generated with MW-LST as an input, which is compared to an existing implementation of the same framework with thermal infrared based LST. In this study, the MW-LST sampling is restricted to the same clear sky days as in the IR-based implementation to be able to analyse the impact of changing the LST dataset separately from the impact of sampling all-sky conditions. The results show that long-term bulk ET estimates agree with a spatial correlation of 92 % for total ET in the Europe/Africa domain and agreement in seasonal (3-month) totals of 83–97 % depending on the time of year. Most importantly, the ALEXI-MW also matches ALEXI-IR very closely in terms of 3-month inter-annual anomalies, demonstrating its ability to capture the development and extent of drought conditions. The weekly ET output from the two parallel ALEXI implementations is further compared to a common ground measured reference provided by the FLUXNET consortium. Overall, they indicate a surprisingly close match in both performance metrics (correlation and RMSE) for all but the most challenging sites in terms of spatial heterogeneity and level of aridity. Moreover, merging MW- and IR-based ALEXI may provide estimates of ET with a reduced uncertainty, even during nominally clear sky days. It is concluded that a constellation of MW satellites can effectively be used to provide LST for estimating ET through TSEB, which is an important step towards all-sky satellite-based ET estimates.

Soil moisture deficiency is a major factor in determining crop yields in water-limited agricultural production regions. Evapotranspiration (ET), which consists of crop water use through transpiration and water loss through direct soil evaporation, is a good indicator of soil moisture availability and vegetation health. ET therefore has been an integral part of many yield estimation efforts. The Evaporative Stress Index (ESI) is an ET-based crop stress indicator that describes temporal anomalies in a normalized evapotranspiration metric as derived from satellite remote sensing. ESI has demonstrated the capacity to explain regional yield variability in water-limited regions. However, its performance in some regions where the vegetation cycle is intensively managed appears to be degraded due to interannual phenological variability. This investigation selected three study sites across the U.S. Corn Belt – Mead, NE, Ames, IA and Champaign, IL – to investigate the potential operational value of 30-m resolution, phenologically corrected ESI datasets for yield prediction. The analysis was conducted over an 8-year period from 2010 to 2017, which included both drought and pluvial conditions as well as a broad range in yield values. Detrended yield anomalies for corn and soybean were correlated with ESI computed using annual ET curves temporally aligned based on (1) calendar date, (2) crop emergence date, and (3) a growing degree day (GDD) scaled time axis. Results showed that ESI has good correlations with yield anomalies at the county scale and that phenological corrections to the annual temporal alignment of the ET timeseries improve the correlation, especially when the time axis is defined by GDD rather than the calendar date. Peak correlations occur in the silking stage for corn and the reproductive stage for soybean – phases when these crops are particularly sensitive to soil moisture deficiencies. Regression equations derived at the time of peak correlation were used to estimate yields at county scale using a leave-one-out cross-validation strategy. The ESI-based yield estimates agree well with the USDA National Agricultural Statistics Service (NASS) county-level crop yield data, with correlation coefficients ranging from 0.79 to 0.93 and percent root-mean-square errors of 5–8%. These results demonstrate that remotely sensed ET at high spatiotemporal resolution can convey valuable water stress information for forecasting crop yields across the Corn Belt if interannual phenological variability is considered.

Multi-sensor remote sensing data fusion technologies have been developed and widely applied in recent years, providing a feasible and economical solution to increase the availability of high spatial and temporal resolution data. These methods, however, have been challenging to apply in highly heterogeneous areas, especially in complex agricultural landscapes where there are rapid changes at small scales, while features at larger scales change more slowly. In this study, we developed a novel method to reconstruct daily 30 m Normalized Difference Vegetation Index (NDVI) using imagery from the Moderate Resolution Imaging Spectroradiometer (MODIS), Landsat and Landsat-like platforms, and the Cropland Data Layer (CDL). This method utilizes a crop reference curve (CRC) approach, in which a set of NDVI time series are extracted from pure MODIS pixels (250 m resolution) identified using the CDL, and then used to fit Landsat-like observations (30 m). The CRC based method was applied over a complex agricultural landscape in the Choptank River watershed on the Eastern Shore of Maryland. Landsat data from 2013 and 2014 and Harmonized Landsat and Sentinel-2 (HLS) data from 2018 were used to reconstruct 30 m daily NDVI maps for major crop types. Results show that the relative error (RE) in reconstructed NDVI is around 6–8% during periods of rapid crop growth, and 3–5% during peak periods when growth is slow. The accuracy of the CRC method outperforms a standard image pair-based data fusion algorithm (Spatial and Temporal Adaptive Reflectance Fusion Model; STARFM), which yields RE of 4–9% in slow-growth periods and 10–16% in fast-growth periods when clear Landsat images are scarce. The CRC method was also compared with time-series data fusion methods, including a harmonic fitting model and the SaTellite dAta IntegRation (STAIR) model. The results show that CRC gives similar results when the Landsat-like image availability is high (around 27 images per year), but outperforms other methods when availability is limited (less than 15 images per year). The reconstructed NDVI time series for corn, soybean, winter wheat/soybean and forest at 30-m resolution show clear phenological patterns at the sub-field scale. The resulting 30-m NDVI timeseries data provide useful information for mapping crop phenology and monitoring crop condition in complex agricultural landscapes, especially for complex double-cropping areas. However, the input requirement of an accurate 30-m crop classification map constrains its application to areas and periods where classifications are available.

Land-surface temperature retrieved from thermal infrared (TIR) remote sensing has proven to be a valuable constraint in surface energy balance models for estimating evapotranspiration (ET). For optimal utility in agricultural water management applications, frequent thermal imaging (

California's Central Valley grows a significant fraction of grapes used for wine production in the United States. With increasing vineyard acreage, reduced water availability in much of California, and competing water use interests, it is critical to be able to monitor regional water use and evapotranspiration (ET) over large areas, but also in detail at individual field scales to improve water management within these viticulture production systems. This can be achieved by integrating remote sensing data from multiple satellite systems with different spatiotemporal characteristics. In this research, we evaluate the utility of a multi-scale system for monitoring ET as applied over two vineyard sites near Lodi, California during the 2013 growing season, leading into the drought in early 2014. The system employs a multi-sensor satellite data fusion methodology (STARFM: Spatial and Temporal Adaptive Reflective Fusion Model) combined with a multi-scale ET retrieval algorithm based on the Two-Source Energy Balance (TSEB) land-surface representation to compute daily ET at 30 m resolution. In this system, TSEB is run using thermal band imagery from the Geostationary Environmental Operational Satellites (GOES; 4-km spatial resolution, hourly temporal sampling), the Moderate Resolution Imaging Spectroradiometer (MODIS) data (1 km resolution, daily acquisition) and the new Landsat 8 satellite (sharpened to 30 m resolution, ~ 16 day acquisition). Estimates of daily ET generated in two neighboring fields of Pinot noir vines of different age agree with ground-based flux measurements acquired in-field during most of the 2013 season with relative mean absolute errors on the order of 19–23% (root mean square errors of approximately 1 mm d − 1 ), reducing to 14–20% at the weekly timestep relevant for irrigation management (~ 5 mm wk − 1 ). A model overestimation of ET in the early season was detected in the younger vineyard, perhaps relating to an inter-row grass cover crop. Spatial patterns of cumulative ET generally correspond to measured yield maps and indicate areas of variable crop moisture, soil condition, and yield within the vineyards that could require adaptive management. The results suggest that multi-sensor remote sensing observations provide a unique means for monitoring crop water use and soil moisture status at field-scales over extended growing regions, and may have value in supporting operational water management decisions in vineyards and other high value crops.

The recent paper by Morillas et al. [Morillas, L. et al. Using radiometric surface temperature for surface energy flux estimation in Mediterranean drylands from a two-source perspective, Remote Sens. Environ. 136, 234–246, 2013] evaluates the two-source model (TSM) of Norman et al. (1995) with revisions by Kustas and Norman (1999) over a semiarid tussock grassland site in southeastern Spain. The TSM - in its current incarnation, the two-source energy balance model (TSEB) - was applied to this landscape using ground-based infrared radiometer sensors to estimate both the composite surface radiometric temperature and component soil and canopy temperatures. Morillas et al. (2013) found the TSEB model substantially underestimated the sensible H (and overestimated the latent heat LE) fluxes. Using the same data set from Morillas et al. (2013), we were able to confirm their results. We also found energy transport and exchange behavior derived from primarily the observations themselves to differ significantly from a number of prior studies using land surface temperature for estimating heat fluxes with one-source modeling approaches in semi-arid landscapes. However, revisions to key vegetation inputs to TSEB and the soil resistance formulation resulted in a significant reduction in the bias and root mean square error (RMSE) between model output of H and LE and the measurements compared to the prior results from Morillas et al. (2013). These included more representative ground-based vegetation greenness and local leaf area index values as well as modifications to the coefficients of the soil resistance formulation to account for the very rough (rocky) soil surface conditions with a clumped canopy. This indicates that both limitations in remote estimates of biophysical indicators of the canopy at the site and the lack of adjustment in soil resistance formulation to account for site specific characteristics, contributed to the earlier findings of Morillas et al. (2013). This suggests further studies need to be conducted to reduce the uncertainties in the vegetation and land surface temperature input data in order to more accurately assess the effects of the transport exchange processes of this Mediterranean landscape on TSEB formulations.

Conventional methods to estimate land surface temperature (LST) from space rely on the thermal infrared (TIR) spectral window and is limited to cloud-free scenes. To also provide LST estimates during periods with clouds, a new method was developed to estimate LST based on passive-microwave (MW) observations. The MW-LST product is informed by six polar-orbiting satellites to create a global record with up to eight observations per day for each 0.25° resolution grid box. For days with sufficient observations, a continuous diurnal temperature cycle (DTC) was fitted. The main characteristics of the DTC were scaled to match those of a geostationary TIR-LST product.This paper tests the cloud tolerance of the MW-LST product. In particular, we demonstrate its stable performance with respect to flux tower observation sites (four in Europe and nine in the United States), over a range of cloudiness conditions up to heavily overcast skies. The results show that TIR-based LST has slightly better performance than MW-LST for clear-sky observations but suffers an increasing negative bias as cloud cover increases. This negative bias is caused by incomplete masking of cloud-covered areas within the TIR scene that affects many applications of TIR-LST. In contrast, for MW-LST we find no direct impact of clouds on its accuracy and bias. MW-LST can therefore be used to improve TIR cloud screening. Moreover, the ability to provide LST estimates for cloud-covered surfaces can help expand current clear-sky-only satellite retrieval products to all-weather applications.

Precipitation, soil moisture, and air temperature are the most commonly used climate variables to monitor drought; however, other climatic factors such as solar radiation, wind speed, and humidity can be important drivers in the depletion of soil moisture and evolution and persistence of drought. This work assesses the Evaporative Demand Drought Index (EDDI) at multiple time scales for several hydroclimates as the second part of a two-part study. EDDI and individual evaporative demand components were examined as they relate to the dynamic evolution of flash drought over the central United States, characterization of hydrologic drought over the western United States, and comparison to commonly used drought metrics of the U.S. Drought Monitor (USDM), Standardized Precipitation Index (SPI), Standardized Soil Moisture Index (SSI), and the evaporative stress index (ESI). Two main advantages of EDDI over other drought indices are that it is independent of precipitation (similar to ESI) and it can be decomposed to identify the role individual evaporative drivers have on drought onset and persistence. At short time scales, spatial distributions and time series results illustrate that EDDI often indicates drought onset well in advance of the USDM, SPI, and SSI. Results illustrate the benefits of physically based evaporative demand estimates and demonstrate EDDI’s utility and effectiveness in an easy-to-implement agricultural early warning and long-term hydrologic drought–monitoring tool with potential applications in seasonal forecasting and fire-weather monitoring.

Many operational drought indices focus primarily on precipitation and temperature when depicting hydroclimatic anomalies, and this perspective can be augmented by analyses and products that reflect the evaporative dynamics of drought. The linkage between atmospheric evaporative demand E0 and actual evapotranspiration (ET) is leveraged in a new drought index based solely on E0—the Evaporative Demand Drought Index (EDDI). EDDI measures the signal of drought through the response of E0 to surface drying anomalies that result from two distinct land surface–atmosphere interactions: 1) a complementary relationship between E0 and ET that develops under moisture limitations at the land surface, leading to ET declining and increasing E0, as in sustained droughts, and 2) parallel ET and E0 increases arising from increased energy availability that lead to surface moisture limitations, as in flash droughts. To calculate EDDI from E0, a long-term, daily reanalysis of reference ET estimated from the American Society of Civil Engineers (ASCE) standardized reference ET equation using radiation and meteorological variables from the North American Land Data Assimilation System phase 2 (NLDAS-2) is used. EDDI is obtained by deriving empirical probabilities of aggregated E0 depths relative to their climatologic means across a user-specific time period and normalizing these probabilities. Positive EDDI values then indicate drier-than-normal conditions and the potential for drought. EDDI is a physically based, multiscalar drought index that that can serve as an indicator of both flash and sustained droughts, in some hydroclimates offering early warning relative to current operational drought indices. The performance of EDDI is assessed against other commonly used drought metrics across CONUS in Part II.

Significant advances have been achieved in generating soil moisture (SM) products from satellite remote sensing and/or land surface modeling with reasonably good accuracy in recent years. However, the discrepancies among the different SM data products can be considerably large, which hampers their usage in various applications. The bias of one SM product from another is well recognized in the literature. Bias estimation and spatial correction methods have been documented for assimilating satellite SM product into land surface and hydrologic models. Nevertheless, understanding the characteristics of each of these SM data products is required for many applications where the most accurate data products are desirable. This study inter-compares five SM data products from three different sources with each other, and evaluates them against in situ SM measurements over 14-year period from 2000 to 2013. Specifically, three microwave (MW) satellite based data sets provided by ESA's Climate Change Initiative (CCI) (CCI-merged, -active and -passive products), one thermal infrared (TIR) satellite based product (ALEXI), and the Noah land surface model (LSM) simulations. The in-situ SM measurements are collected from the North American Soil Moisture Database (NASMD), which involves more than 600 ground sites from a variety of networks. They are used to evaluate the accuracies of these five SM data products. In general, each of the five SM products is capable of capturing the dry/wet patterns over the study period. However, the absolute SM values among the five products vary significantly. SM simulations from Noah LSM are more stable relative to the satellite-based products. All TIR and MW satellite based products are relatively noisier than the Noah LSM simulations. Even though MW satellite based SM retrievals have been predominantly used in the past years, SM retrievals of the ALEXI model based on TIR satellite observations demonstrate skills equivalent to all the MW satellite retrievals and even slightly better over certain regions. Compared to the individual active and passive MW products, the merged CCI product exhibits higher anomaly correlation with both Noah LSM simulations and in-situ SM measurements.

Thermal and multispectral remote sensing data from low-altitude aircraft can provide high spatial resolution necessary for sub-field (≤ 10 m) and plant canopy (≤ 1 m) scale evapotranspiration (ET) monitoring. In this study, high-resolution (sub-meter-scale) thermal infrared and multispectral shortwave data from aircraft are used to map ET over vineyards in central California with the two-source energy balance (TSEB) model and with a simple model having operational immediate capabilities called DATTUTDUT (Deriving Atmosphere Turbulent Transport Useful To Dummies Using Temperature). The latter uses contextual information within the image to scale between radiometric land surface temperature (TR) values representing hydrologic limits of potential ET and a non-evaporative surface. Imagery from 5 days throughout the growing season is used for mapping ET at the sub-field scale. The performance of the two models is evaluated using tower-based measurements of sensible (H) and latent heat (LE) flux or ET. The comparison indicates that TSEB was able to derive reasonable ET estimates under varying conditions, likely due to the physically based treatment of the energy and the surface temperature partitioning between the soil/cover crop inter-row and vine canopy elements. On the other hand, DATTUTDUT performance was somewhat degraded presumably because the simple scaling scheme does not consider differences in the two sources (vine and inter-row) of heat and temperature contributions or the effect of surface roughness on the efficiency of heat exchange. Maps of the evaporative fraction (EF = LE/(H + LE)) from the two models had similar spatial patterns but different magnitudes in some areas within the fields on certain days. Large EF discrepancies between the models were found on 2 of the 5 days (DOY 162 and 219) when there were significant differences with the tower-based ET measurements, particularly using the DATTUTDUT model. These differences in EF between the models translate to significant variations in daily water use estimates for these 2 days for the vineyards. Model sensitivity analysis demonstrated the high degree of sensitivity of the TSEB model to the accuracy of the TR data, while the DATTUTDUT model was insensitive to systematic errors in TR as is the case with contextual-based models. However, it is shown that the study domain and spatial resolution will significantly influence the ET estimation from the DATTUTDUT model. Future work is planned for developing a hybrid approach that leverages the strengths of both modeling schemes and is simple enough to be used operationally with high-resolution imagery.

This study examines the evolution of several model-based and satellite-derived drought metrics sensitive to soil moisture and vegetation conditions during the extreme flash drought event that impacted major agricultural areas across the central U.S. during 2012. Standardized anomalies from the remote sensing based Evaporative Stress Index (ESI) and Vegetation Drought Response Index (VegDRI) and soil moisture anomalies from the North American Land Data Assimilation System (NLDAS) are compared to the United States Drought Monitor (USDM), surface meteorological conditions, and crop and soil moisture data compiled by the National Agricultural Statistics Service (NASS). Overall, the results show that rapid decreases in the ESI and NLDAS anomalies often preceded drought intensification in the USDM by up to 6 wk depending on the region. Decreases in the ESI tended to occur up to several weeks before deteriorations were observed in the crop condition datasets. The NLDAS soil moisture anomalies were similar to those depicted in the NASS soil moisture datasets; however, some differences were noted in how each model responded to the changing drought conditions. The VegDRI anomalies tracked the evolution of the USDM drought depiction in regions with slow drought development, but lagged the USDM and other drought indicators when conditions were changing rapidly. Comparison to the crop condition datasets revealed that soybean conditions were most similar to ESI anomalies computed over short time periods (2–4 wk), whereas corn conditions were more closely related to longer-range (8–12 wk) ESI anomalies. Crop yield departures were consistent with the drought severity depicted by the ESI and to a lesser extent by the NLDAS and VegDRI datasets.

Accurate partitioning of the evaporation (E) and transpiration (T) components of evapotranspiration (ET) in remote sensing models is important for evaluating strategies aimed at increasing crop water productivity. A two-source energy balance (TSEB) model designed for row crops solves the energy balance of the soil-canopy-atmosphere continuum using surface brightness temperature. By solving the energy balance of the soil and plant canopy separately, the TSEB model can calculate E and T, which cannot be done with single-source models. However, few studies have tested the TSEB model where E or T measurements were available, which until recently has impeded its advance. This article reviews recent physically based advances of the TSEB model. The advances were tested using measurements of E, T, and ET by microlysimeters, sap flow gauges, and weighing lysimeters, respectively, at Bushland, Texas, for irrigated cotton having a wide range of canopy cover. Root mean square error (RMSE) and mean bias error (MBE) were 0.54 and -0.19 mm d-1, respectively, between measured and calculated E. RMSE and MBE were 0.87 and 0.31 mm d-1, respectively, between measured and calculated T. This was deemed an improvement over previous TSEB model versions, which overestimated E and underestimate T, resulting in RMSE and MBE of up to 3.8 and -3.5 mm d-1, respectively. Ongoing research includes testing the TSEB model using different remote sensing platforms, from ground-based to satellite scales.

Mapping the spatial variability of actual evapotranspiration (ETa) across vineyards is useful for optimizing irrigation scheduling and efficiency, leading to conservation of water resources and more sustainable wine grape production. To support efficient irrigation strategies, we investigate the utility of thermal infrared-based ETa maps over a range of vineyards located throughout California, each representing a unique local climate, trellis design, grape variety, row orientation and management practice. ETa maps are derived by combining the Disaggregated Atmosphere Land Exchange Inverse (ALEXI/DisALEXI) surface energy balance model and the Spatial Temporal Adaptive Reflectance Fusion Model (STARFM) to generate ETa estimates at high spatial (30 m) and temporal (daily) resolution. Model output is evaluated for years 2017 and 2018 over vineyard sites located in Sonoma, Sacramento, and Madera counties in California that are being monitored as part of the Grape Remote sensing Atmospheric Profile and Evapotranspiration eXperiment (GRAPEX). Overall, modeled daily ET estimates compare well with flux tower observations, with average root mean square error (RMSE), mean absolute error (MAE) and mean bias error (MBE) of 0.88 mm day−1, 0.70 mm day−1, and 0.17 mm day−1 respectively, over all four individual vineyard locations, aligning with past GRAPEX studies. Despite general agreement, record wildfires in northern California during 2018 likely resulted in positive model bias, while misrepresentation of leaf area index within a double-trellis designed canopy at the southern-most vineyard resulted in negative model bias. Spatial analysis of monthly total ET highlights the advantages of utilizing a satellite-based approach to characterize the variability in water use within and surrounding the targeted vineyards. A reliable spatial ET product at scale has the potential to improve water allocation and conservation efforts by identifying areas of uneven water use due to variations in soil texture and composition and other environmental or anthropogenic factors.

Leaf area index (LAI) is an essential indicator of crop development and growth. For many agricultural applications, satellite-based LAI estimates at the farm-level often require near-daily imagery at medium to high spatial resolution. The combination of data from different ongoing satellite missions, Sentinel 2 (ESA) and Landsat 8 (NASA), provides this opportunity. In this study, we evaluated the leaf area index generated from three methods, namely, existing vegetation index (VI) relationships applied to Harmonized Landsat-8 and Sentinel-2 (HLS) surface reflectance produced by NASA, the SNAP biophysical model, and the THEIA L2A surface reflectance products from Sentinel-2. The intercomparison was conducted over the agricultural scheme in Bekaa (Lebanon) using a large set of in-field LAIs and other biophysical measurements collected in a wide variety of canopy structures during the 2018 and 2019 growing seasons. The major studied crops include herbs (e.g., cannabis: Cannabis sativa, mint: Mentha, and others), potato (Solanum tuberosum), and vegetables (e.g., bean: Phaseolus vulgaris, cabbage: Brassica oleracea, carrot: Daucus carota subsp. sativus, and others). Additionally, crop-specific height and above-ground biomass relationships with LAIs were investigated. Results show that of the empirical VI relationships tested, the EVI2-based HLS models statistically performed the best, specifically, the LAI models originally developed for wheat (RMSE:1.27), maize (RMSE:1.34), and row crops (RMSE:1.38). LAI derived through European Space Agency’s (ESA) Sentinel Application Platform (SNAP) biophysical processor underestimated LAI and provided less accurate estimates (RMSE of 1.72). Additionally, the S2 SeLI LAI algorithm (from SNAP biophysical processor) produced an acceptable accuracy level compared to HLS-EVI2 models (RMSE of 1.38) but with significant underestimation at high LAI values. Our findings show that the LAI-VI relationship, in general, is crop-specific with both linear and non-linear regression forms. Among the examined indices, EVI2 outperformed other vegetation indices when all crops were combined, and therefore it can be identified as an index that is best suited for a unified algorithm for crops in semi-arid irrigated regions with heterogeneous landscapes. Furthermore, our analysis shows that the observed height-LAI relationship is crop-specific and essentially linear with an R2 value of 0.82 for potato, 0.79 for wheat, and 0.50 for both cannabis and tobacco. The ability of the linear regression to estimate the fresh and dry above-ground biomass of potato from both observed height and LAI was reasonable, yielding R2: ~0.60.

The Central Valley of California relies on irrigation for crop production, but water resources, particularly groundwater, have reached a critical state due to extended drought periods and overuse by irrigated agriculture. The purpose of the Grape Remote sensing Atmospheric Profile and Evapotranspiration eXperiment (GRAPEX) project is to improve efficiency in vineyard irrigation through field measurements and modeling efforts using remote sensing. In this study, we analyze multi-year timeseries of ground-based soil moisture and evapotranspiration (ET) measurements collected during GRAPEX to identify patterns that may inform future irrigation recommendations based on remotely sensed ET. The study focuses on field data collected in two adjacent Pinot noir vineyards near Lodi, CA from 2013 to 2017. The data used for this analysis are reference evapotranspiration (ETref), actual evapotranspiration (ETa), mean daytime soil water content (SWC) measured to a depth of 90 cm, precipitation, and irrigation. The relationship between SWC and the ratio of ETa/ETref (used as a soil moisture proxy indicator) changes throughout the spring and summer months due to advancing phenological stages and management practices. In early spring, when the interrow cover crop is the primary source of ET, ETa is strongly correlated with SWC. As the vine canopy grows in, the strong correlation breaks down as the vines begin to access to water beyond the depth of the soil moisture sensors. As the soil profile dries out during the summer, the correlation between ETa and SWC once again emerges as the vines become strongly dependent on irrigation. The dynamic interaction between the upper and lower root zone profile and evaporative demand means that the key to understanding vineyard water status using relatively shallow SWC observations is to use them in conjunction with ETa data, so that when both SWC and ETa timeseries are highly correlated vine water status is well defined. This suggests, delaying the initiation of irrigation during the period of rapid vine growth until a clear relationship between SWC and ETa emerges would both reduce the total amount of water used and give the grower more control over vine growth and grape quality affected by water status.

Evapotranspiration (ET) is a measure of plant water use that is utilized regionally for drought detection and monitoring, and locally for agricultural water resource management. Understanding the uncertainty associated with this measurement is vital for science predictions and analysis and for water resource management decision making. In this manuscript, the uncertainty in disaggregated Atmosphere-Land Exchange (disALEXI) is quantified; disALEXI is an ET algorithm that utilizes land surface temperature (LST) derived from the ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS), as well ancillary inputs for landcover, elevation, vegetation parameters, and meteorological inputs. Since each of these inputs has an associated, and potentially unknown, uncertainty, in this study a Monte Carlo simulation based on a spatial statistical model is used to determine the algorithm's sensitivity to each of its inputs, and to quantify the probability distribution of algorithm outputs. Analysis shows that algorithm is most sensitive to LST (the input derived from ECOSTRESS). Significantly, the output uncertainty distribution is non-Gaussian, due to the non-linear nature of the algorithm. This means that ET uncertainty cannot be prescribed by accuracy and precision alone. Here, uncertainty was represented using five quantiles of the output distribution. The distribution was consistent across five different datasets (mean offset is 0.01 mm/day, and 95% of the data is contained within 0.3 mm/day). An additional two datasets with low ET, showed higher uncertainty (95% of the data is within 1 mm/day), and a positive bias (i.e., ET was overestimated by an average of 0.12 mm/day when ET was low).

Crop emergence date is a critical input to models of crop development and biomass accumulation. The ability to robustly detect and map emergence date using remote sensing would greatly benefit operational yield estimation and crop monitoring efforts; however, this has proven to be challenging. Previous remote-sensing phenology algorithms showed that crop stages can typically be detected starting only around the V3-V4 (3 to 4 established leaves) vegetative stage. Furthermore, traditional approaches have a strong assumption regarding the temporal evolution of plant growth and normally require a complete growth period of observations to define seasonal changes. Most approaches were not designed for within-season operational mapping, particularly in the early growing season. In the current paper, we describe a new within-season emergence (WISE) approach to mapping crop green-up date using satellite observations available during early growth stages. The approach was first optimized using high spatiotemporal resolution (10 m, 2-day revisit) imagery from the Vegetation and Environment monitoring New MicroSatellite (VENμS) research mission, and assessed using ground observations of early crop growth stages (emergence VE and one leaf V1 stages for corn, and emergence VE and unifoliolate VC stages for soybeans) collected over the Beltsville Agricultural Research Center (BARC) experimental fields in Beltsville, MD during the 2019 growing season. Results show that early crop growth stages can be reliably detected at sub-field scale about two weeks after crop emergence. The remote-sensing green-up dates were about 4–5 days after crop emergence on average. Coefficients of determination (R2) between green-up dates and the mid-point dates of the early growth stages were above 0.90. The mean absolute differences, standard deviations, and root mean square errors comparing to the early growth stage mid-point dates were within six days. The maximum differences were within ±10 days across all fields. The WISE approach was assessed using operational Sentinel-2 data (10 m, 5-day revisit) over BARC. The detected green-up dates from Sentinel-2 were consistent with those from VENμS. Some fields were not detected due to the lack of observations around the emergence dates. For independent evaluation, the WISE approach was applied over an agricultural watershed on the Maryland Eastern Shore using both VENμS and the Harmonized Landsat and Sentinel-2 (HLS) data (30 m, 3–4-day revisit). The detected green-up dates were compared with emergence dates in crop progress reports from the National Agricultural Statistics Service (NASS) at the state-level. The WISE-detected green-up dates at the regional scale are within VE stage ranges but slightly earlier than NASS crop progress reports at the state-level. The WISE approach uses remote-sensing observations during the early crop growth stages and has potential for operational application within the season using Sentinel-2 and HLS data.

Crop yield estimates over large areas are conventionally made using weather observations, but a comprehensive understanding of the effects of various environmental indicators, observation frequency, and the choice of prediction algorithm remains elusive. Here we present a thorough assessment of county-level maize yield prediction in U.S. Midwest using six statistical/machine learning algorithms (Lasso, Support Vector Regressor, Random Forest, XGBoost, Long-short term memory (LSTM), and Convolutional Neural Network (CNN)) and an extensive set of environmental variables derived from satellite observations, weather data, land surface model results, soil maps, and crop progress reports. Results show that seasonal crop yield forecasting benefits from both more advanced algorithms and a large composite of information associated with crop canopy, environmental stress, phenology, and soil properties (i.e. hundreds of features). The XGBoost algorithm outperforms other algorithms both in accuracy and stability, while deep neural networks such as LSTM and CNN are not advantageous. The compositing interval (8-day, 16-day or monthly) of time series variable does not have significant effects on the prediction. Combining the best algorithm and inputs improves the prediction accuracy by 5% when compared to a baseline statistical model (Lasso) using only basic climatic and satellite observations. Reasonable county-level yield foresting is achievable from early June, almost four months prior to harvest. At the national level, early-season (June and July) prediction from the best model outperforms that of the United States Department of Agriculture (USDA) World Agricultural Supply and Demand Estimates (WASDE). This study provides insights into practical crop yield forecasting and the understanding of yield response to climatic and environmental conditions.

Geographically isolated wetlands (GIWs) provide a portfolio of ecosystem services in low-gradient, Coastal Plain landscapes. Understanding how GIWs influence downstream waters is becoming increasingly important for conservation and management of these unique and important wetland ecosystems. Climatic conditions are known to be key drivers of water budgets at both individual GIW and landscape scales; however differences in hydrologic response across these scales may provide insights into how GIWs influence downstream waters. In this study, we use a combination of GIW water level, gaged streamflow, and climatic data to explore linkages between seasonal climatic drivers, GIW hydrology, and downstream discharges within the Coastal Plain of the Chesapeake Bay watershed. We first examine water balance components at the larger watershed scale, where climatic drivers result in an energy-limited wet season from December to May and a water-limited dry season from June to November. We compare long-term water levels of three GIWs with downstream discharges. GIW water level and downstream discharges are correlated at seasonal (R2: 0.52–0.60) and daily (R2: 0.52–0.76) time steps. However, during dry seasons, GIW water level receded at a faster rate than downstream discharges, highlighting the influence of evapotranspiration on surface and shallow subsurface water storage. Conversely during wet seasons, GIW water level receded slower than downstream discharges, highlighting a potential period for surface water connectivity between GIWs and downstream discharges. Cumulatively, these findings quantify the impact of seasonal climatic drivers on GIW hydrology and connectivity to downstream waters.

Efficient water use assessment and irrigation management is critical for the sustainability of irrigated agriculture, especially under changing climate conditions. Due to the impracticality of maintaining ground instrumentation over wide geographic areas, remote sensing and numerical model-based fine-scale mapping of soil water conditions have been applied for water resource applications at a range of spatial scales. Here, we present a prototype framework for integrating high-resolution thermal infrared (TIR) and synthetic aperture radar (SAR) remote sensing data into a soil-vegetation-atmosphere-transfer (SVAT) model with the aim of providing improved estimates of surface- and root-zone soil moisture that can support optimized irrigation management strategies. Specifically, remotely-sensed estimates of water stress (from TIR) and surface soil moisture retrievals (from SAR) are assimilated into a 30-m resolution SVAT model over a vineyard site in the Central Valley of California, U.S. The efficacy of our data assimilation algorithm is investigated via both the synthetic and real data experiments. Results demonstrate that a particle filtering approach is superior to an ensemble Kalman filter for handling the nonlinear relationship between model states and observations. In addition, biophysical conditions such as leaf area index are shown to impact the relationship between observations and states and must therefore be represented accurately in the assimilation model. Overall, both surface and root-zone soil moisture predicted via the SVAT model are enhanced through the assimilation of thermal and radar-based retrievals, suggesting the potential for improving irrigation management at the agricultural sub-field scale using a data assimilation strategy.