Search

Water availability controls vegetation activity and the carbon balance of terrestrial ecosystems across a large portion of the global land surface. Although the influence of terrestrial water storage (TWS) on the land carbon balance is evident in globally aggregated measures, it remains unknown whether the large annual amplitudes in TWS are causally linked to water availability in the rooting zone of vegetation, or whether they reflect a correlation of plant water stress with water stored in other landscape elements that may not directly be connected to vegetation functioning (lakes, rivers, groundwater). Global models of the land surface typically ignore hillslope-scale variations in plant water availability, and water stores that are located beyond the soil, and beyond prescribed plant rooting depths. This simplification is partly owed to a lack of empirical information.Here, we approach this gap from two angles: from the site scale using eddy covariance observations, and from the global scale using earth observations. Water mass balance constraints derived from thermal infrared-based evapotranspiration (ET) estimates and precipitation reanalysis data indicate plant-available water stores that exceed the storage capacity of 2 m deep soils across 37% of the Earth’s vegetated surface. Large spatial variations of the rooting zone water storage capacity across topographic and hydro-climatic gradients are tightly linked to the sensitivity of vegetation activity (measured by sun-induced fluorescence and by the evaporative fraction) to water deficits. Similar patterns between ET and cumulative water deficits emerge from site-level flux measurements. We found large variations of the vegetation sensitivity to dry conditions across sites and at several sites a muted response of ET to dry conditions in spite of large (>300 mm) seasonal water deficits at some sites.Taken together, results we show here hint at a critical role of plant access to deep water stores and the need to extend the focus beyond moisture in the top 1-2 m of soil for understanding and simulating land-atmosphere exchange. Our results add to the emerging evidence that water stored in the weathered bedrock and plant access to groundwater may have a more important role in regulating land-atmosphere exchange and the carbon cycle than previously appreciated.

Robust information on consumptive water use (evapotranspiration, ET) derived from remote sensing can significantly benefit water decision-making in agriculture, informing irrigation schedules and water management plans over extended regions. To be of optimal utility for operational usage, these remote sensing ET data should be generated at the sub-field spatial resolution and daily-to-weekly timesteps commensurate with the scales of water management activities. However, current methods for field-scale ET retrieval based on thermal infrared (TIR) imaging, a valuable diagnostic of canopy stress and surface moisture status, are limited by the temporal revisit of available medium-resolution (100 m or finer) thermal satellite sensors. This study investigates the efficacy of a data fusion method for combining information from multiple medium-resolution sensors toward generating high spatiotemporal resolution ET products for water management. TIR data from Landsat and ECOSTRESS (both at ~ 100-m native resolution), and VIIRS (375-m native) are sharpened to a common 30-m grid using surface reflectance data from the Harmonized Landsat-Sentinel dataset. Periodic 30-m ET retrievals from these combined thermal data sources are fused with daily retrievals from unsharpened VIIRS to generate daily, 30-m ET image timeseries. The accuracy of this mapping method is tested over several irrigated cropping systems in the Central Valley of California in comparison with flux tower observations, including measurements over irrigated vineyards collected in the GRAPEX campaign. Results demonstrate the operational value added by the augmented TIR sensor suite compared to Landsat alone, in terms of capturing daily ET variability and reduced latency for real-time applications. The method also provides means for incorporating new sources of imaging from future planned thermal missions, further improving our ability to map rapid changes in crop water use at field scales.

Management zones (MZs) are efficient for applying site-specific management in agricultural fields. This study proposes an approach for generating MZs using time-series clustering (TSC) to also enable time-specific management. TSC was applied to daily remote sensing retrievals in a California vineyard during four growing seasons (2015–2018) using three datasets: evapotranspiration (ET), leaf area index (LAI), and normalized difference vegetation index (NDVI). Distinct MZs were delineated based on similarities in pixel-level temporal dynamics for each dataset, using dissimilarity index to determine the optimal number of clusters and compare TSC results. The differences between the cluster centers were calculated, along with the ratio between the centers’ differences and the range of each dataset, denoting the degree of difference between MZs centers. Similarity between MZs from each factor was quantified using Cramer’s V and Frechet distances. Finally, an aggregated (multi-factor) MZ map was generated using multivariate clustering. The resulting MZs were compared to a 2016 yield map to determine the significance of differences between means and distribution among MZs. The findings show that LAI TSC achieved the best cluster separation. The NDVI and LAI MZs maps were nearly identical (Cramer’s V of 0.97), while ET showed weaker similarities to NDVI and LAI (0.61 and 0.62, respectively). Similar findings were observed for the Frechet distances. The yield values were found to be significantly different among MZs for all TSC maps. TSC may be further utilized for defining within-field spatial variability and temporal dynamics for precision irrigation practices that account for spatial and temporal variability.

The assimilation of L-band surface brightness temperature (Tb) into the land surface model (LSM) component of a numerical weather prediction (NWP) system is generally expected to improve the quality of summertime 2-m air temperature (T2m) forecasts during water-limited surface conditions. However, recent retrospective results from the European Centre for Medium-Range Weather Forecasts (ECMWF) suggest that the assimilation of L-band Tb from the European Space Agency’s (ESA) Soil Moisture Ocean Salinity (SMOS) mission may, under certain circumstances, degrade the accuracy of growing-season 24-h T2m forecasts within the central United States. To diagnose the source of this degradation, we evaluate ECMWF soil moisture (SM) and evapotranspiration (ET) forecasts using both in situ and remote sensing resources. Results demonstrate that the assimilation of SMOS Tb broadly improves the ECMWF SM analysis in the central United States while simultaneously degrading the quality of 24-h ET forecasts. Based on a recently derived map of true global SM–ET coupling and a synthetic fraternal twin data assimilation experiment, we argue that the spatial and temporal characteristics of ECMWF SM analyses and ET forecast errors are consistent with the hypothesis that the ECMWF LSM overcouples SM and ET and, as a result, is unable to effectively convert an improved SM analysis into enhanced ET and T2m forecasts. We demonstrate that this overcoupling is likely linked to the systematic underestimation of root-zone soil water storage capacity by LSMs within the U.S. Corn Belt region.

Highlights.Precision agriculture (PA) applications in irrigation are stymied by lack of decision support systems.Modern PA relies on sensor systems and near real-time feedback for irrigation decision support and control.Sophisticated understanding of biophysics and biological systems now guides site-specific irrigation.The internet of things (IOT) enables new ways to increase yield per unit of water used and nutrient use efficiency. Keywords: Crop water productivity, Decision support system, Internet of things, Remote sensing, SCADA, Soil water content.

Efficient use of available water resources is key to sustainable viticulture management in California (CA) and other regions with limited water availability in the western US and abroad. This requi...

Landsat-based monitoring of seasonal and near real-time evapotranspiration (ET) in California vineyards is currently challenged by its low temporal revisit period and missing data under cloudy cond...

Probabilistic forecasts of changes in soil moisture and an Evaporative Stress Index (ESI) on sub-seasonal time scales over the contiguous U.S. are developed. The forecasts use the current land surface conditions and numerical weather prediction forecasts from the Sub-seasonal to Seasonal (S2S) Prediction Project. Changes in soil moisture are quite predictable 8-14 days in advance with 50% or more of the variance explained over the majority of the contiguous U.S.; however, changes in ESI are significantly less predictable. A simple red noise model of predictability shows that the spatial variations in forecast skill are primarily a result of variations in the autocorrelation, or persistence, of the predicted variable, especially for the ESI. The difference in overall skill between soil moisture and ESI, on the other hand, is due to the greater soil moisture predictability by the numerical model forecasts. As the forecast lead time increases from 8-14 days to 15-28 days, however, the autocorrelation dominates the soil moisture and ESI differences as well. An analysis of modelled transpiration, and bare soil and canopy water evaporation contributions to total evaporation, suggests improvements to the ESI forecasts can be achieved by estimating the relative contributions of these components to the initial ESI state. The importance of probabilistic forecasts for reproducing the correct probability of anomaly intensification is also shown.

An operational GOES ET and Drought (GET-D) product system has been developed at NESDIS of NOAA for numerical weather prediction (NWP), national water model validation, and drought monitoring. The GET-D system was generating ET and Evaporative Stress Index (ESI) maps at 8 km spatial resolution using thermal infrared (TIR) observations from the Imagers on GOES-13 and GOES-15 satellites. With the primary operational GOES satellites transitioned to GOES-16 and GOES-17 with the Advanced Baseline Imagers (ABI), the GET-D system is being upgraded to seamlessly generate ET products using ABI observations from GOES-16/17 at high special resolution of up to 2 km. The new ET product will be comprehensively validated using ground ET measurements. Moreover, AMSR2 Ka band observations will be coupled with the GOES TIR channel using machine learning technique to fill in cloudy pixels where current ET retrievals cannot be obtained because of cloud contamination. The all-weather ET product will also be validated against in-situ observations. Preliminary validation results will be presented in this paper.

Formal planning and development of what became the first Landsat satellite commenced over 50 years ago in 1967. Now, having collected earth observation data for well over four decades since the 1972 launch of Landsat-1, the Landsat program is increasingly complex and vibrant. Critical programmatic elements are ensuring the continuity of high quality measurements for scientific and operational investigations, including ground systems, acquisition planning, data archiving and management, and provision of analysis ready data products. Free and open access to archival and new imagery has resulted in a myriad of innovative applications and novel scientific insights. The planning of future compatible satellites in the Landsat series, which maintain continuity while incorporating technological advancements, has resulted in an increased operational use of Landsat data. Governments and international agencies, among others, can now build an expectation of Landsat data into a given operational data stream. International programs and conventions (e.g., deforestation monitoring, climate change mitigation) are empowered by access to systematically collected and calibrated data with expected future continuity further contributing to the existing multi-decadal record. The increased breadth and depth of Landsat science and applications have accelerated following the launch of Landsat-8, with significant improvements in data quality. Herein, we describe the programmatic developments and institutional context for the Landsat program and the unique ability of Landsat to meet the needs of national and international programs. We then present the key trends in Landsat science that underpin many of the recent scientific and application developments and follow-up with more detailed thematically organized summaries. The historical context offered by archival imagery combined with new imagery allows for the development of time series algorithms that can produce information on trends and dynamics. Landsat-8 has figured prominently in these recent developments, as has the improved understanding and calibration of historical data. Following the communication of the state of Landsat science, an outlook for future launches and envisioned programmatic developments are presented. Increased linkages between satellite programs are also made possible through an expectation of future mission continuity, such as developing a virtual constellation with Sentinel-2. Successful science and applications developments create a positive feedback loop—justifying and encouraging current and future programmatic support for Landsat. © 2019 The United States Geological Survey (USGS) and the National Aeronautics and Space Administration (NASA) are gratefully acknowledged for support and encouragement of the 2012–2017 Landsat Science Team ( https://landsat.usgs.gov/landsat-science-teams ). The Editor and Reviewers are thanked for the valuable insights and constructive suggestions made to improve this manuscript.

This study examines the evolution of soil moisture, evapotranspiration, vegetation, and atmospheric conditions during an unusual flash drought–flash recovery sequence that occurred across the south-central United States during 2015. This event was characterized by a period of rapid drought intensification (flash drought) during late summer that was terminated by heavy rainfall at the end of October that eliminated the extreme drought conditions over a 2-week period (flash recovery). A detailed analysis was performed using time series of environmental variables derived from meteorological, remote sensing, and land surface modeling datasets. Though the analysis revealed a similar progression of cascading effects in each region, characteristics of the flash drought such as its onset time, rate of intensification, and vegetation impacts differed between regions due to variations in the antecedent conditions and the atmospheric anomalies during its growth. Overall, flash drought signals initially appeared in the near-surface soil moisture, followed closely by reductions in evapotranspiration. Total column soil moisture deficits took longer to develop, especially in the western part of the region where heavy rainfall during the spring and early summer led to large moisture surpluses. Large differences were noted in how land surface models in the North American Land Data Assimilation System depicted soil moisture evolution during the flash drought; however, the models were more similar in their assessment of conditions during the flash recovery period. This study illustrates the need to use multiple datasets to track the evolution and impacts of rapidly evolving flash drought and flash recovery events.

Global food production depends upon many factors that Earth observing satellites routinely measure about water, energy, weather, and ecosystems. Increasingly sophisticated, publicly-available satellite data products can improve efficiencies in resource management and provide earlier indication of environmental disruption. Satellite remote sensing provides a consistent, long-term record that can be used effectively to detect large-scale features over time, such as a developing drought. Accuracy and capabilities have increased along with the range of Earth observations and derived products that can support food security decisions with actionable information. This paper highlights major capabilities facilitated by satellite observations and physical models that have been developed and validated using remotely-sensed observations. Although we primarily focus on variables relevant to agriculture, we also include a brief description of the growing use of Earth observations in support of aquaculture and fisheries.

Relatively small fluctuations in the surface energy balance and evapotranspiration in semiarid and arid regions can be indicative of significant changes to ecosystem health. Therefore, it is imperative to have approaches for monitoring surface fluxes in these regions. The remote sensing-based two-source energy balance (TSEB) model is a suitable method for flux estimation over sparsely vegetated semiarid and arid landscapes since it explicitly considers surface energy flux contributions from soil and vegetation. However, previous studies indicate that TSEB generally underestimates sensible heat flux (H) and hence overestimates latent heat flux (LE) or evapotranspiration for these regions unless soil resistance coefficients are modified based on additional ground information. In this study, TSEB is applied over semiarid and arid regions on three continents using the original soil resistance formulation with modified coefficients and a recently developed physically based soil resistance formulation. Model sensitivity analysis demonstrates the high sensitivity of TSEB with original soil resistance formulation to soil resistance coefficients, while TSEB with the new soil resistance formulation has relatively low sensitivity to uncertainties in all coefficients. The performance of TSEB using different soil resistance formulations is evaluated by comparing modeled H against eddy covariance measurements in six semiarid and arid study sites and ranking the error statistics. Our results indicate that incorporating the new soil resistance formulation into TSEB would enhance its utility in flux estimation over heterogeneous landscapes by obviating its reliance on semiempirical coefficients and thus provide more robust fluxes over sparsely vegetated regions without model calibration and/or parameter tuning., Spanish Ministry of Science and Innovation. Grant Number: CGL2016‐78075‐P DINCOS. Grant Number: CGL2016‐78075‐P Swiss National Science Foundation. Grant Number: P2EZP2‐165244 Strategic Priority Research Program of Chinese Academy of Sciences. Grant Number: XDA19040504

Land surface temperature (LST) is a key diagnostic indicator of agricultural water use and crop stress. LST data retrieved from thermal infrared (TIR) band imagery, however, tend to have a coarser spatial resolution (e.g., 100 m for Landsat 8) than surface reflectance (SR) data collected from shortwave bands on the same instrument (e.g., 30 m for Landsat). Spatial sharpening of LST data using the higher resolution multi-band SR data provides an important path for improved agricultural monitoring at sub-field scales. A previously developed Data Mining Sharpener (DMS) approach has shown great potential in the sharpening of Landsat LST using Landsat SR data co-collected over various landscapes. This work evaluates DMS performance for sharpening ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS) LST (~70 m native resolution) and Visible Infrared Imaging Radiometer Suite (VIIRS) LST (375 m) data using Harmonized Landsat and Sentinel-2 (HLS) SR data, providing the basis for generating 30-m LST data at a higher temporal frequency than afforded by Landsat alone. To account for the misalignment between ECOSTRESS/VIIRS and Landsat/HLS caused by errors in registration and orthorectification, we propose a modified version of the DMS approach that employs a relaxed box size for energy conservation (EC). Sharpening experiments were conducted over three study sites in California, and results were evaluated visually and quantitatively against LST data from unmanned aerial vehicles (UAV) flights and from Landsat 8. Over the three sites, the modified DMS technique showed improved sharpening accuracy over the standard DMS for both ECOSTRESS and VIIRS, suggesting the effectiveness of relaxing EC box in relieving misalignment-induced errors. To achieve reasonable accuracy while minimizing loss of spatial detail due to the EC box size increase, an optimal EC box size of 180-270 m was identified for ECOSTRESS and about 780 m for VIIRS data based on experiments from the three sites. Results from this work will facilitate the development of a prototype system that generates high spatiotemporal resolution LST products for improved agricultural water use monitoring by synthesizing multi-source remote sensing data.

The term “flash drought” is frequently invoked to describe droughts that develop rapidly over a relatively short timescale. Despite extensive and growing research on flash drought processes, predictability, and trends, there is still no standard quantitative definition that encompasses all flash drought characteristics and pathways. Instead, diverse definitions have been proposed, supporting wide-ranging studies of flash drought but creating the potential for confusion as to what the term means and how to characterize it. Use of different definitions might also lead to different conclusions regarding flash drought frequency, predictability, and trends under climate change. In this study, we compared five previously published definitions, a newly proposed definition, and an operational satellite-based drought monitoring product to clarify conceptual differences and to investigate the sensitivity of flash drought inventories and trends to the choice of definition. Our analyses indicate that the newly introduced Soil Moisture Volatility Index definition effectively captures flash drought onset in both humid and semi-arid regions. Analyses also showed that estimates of flash drought frequency, spatial distribution, and seasonality vary across the contiguous United States depending upon which definition is used. Definitions differ in their representation of some of the largest and most widely studied flash droughts of recent years. Trend analysis indicates that definitions that include air temperature show significant increases in flash droughts over the past 40 years, but few trends are evident for definitions based on other surface conditions or fluxes. These results indicate that “flash drought” is a composite term that includes several types of events and that clarity in definition is critical when monitoring, forecasting, or projecting the drought phenomenon.

The thermal-based Two-Source Energy Balance (TSEB) model partitions the evapotranspiration (ET) and energy fluxes from vegetation and soil components providing the capability for estimating soil evaporation (E) and canopy transpiration (T). However, it is crucial for ET partitioning to retrieve reliable estimates of canopy and soil temperatures and net radiation, as the latter determines the available energy for water and heat exchange from soil and canopy sources. These two factors become especially relevant in row crops with wide spacing and strongly clumped vegetation such as vineyards and orchards. To better understand these effects, very high spatial resolution remote-sensing data from an unmanned aerial vehicle were collected over vineyards in California, as part of the Grape Remote sensing and Atmospheric Profile and Evapotranspiration eXperiment and used in four different TSEB approaches to estimate the component soil and canopy temperatures, and ET partitioning between soil and canopy. Two approaches rely on the use of composite $$T_\mathrm{rad}$$ , and assume initially that the canopy transpires at the Priestley–Taylor potential rate. The other two algorithms are based on the contextual relationship between optical and thermal imagery partition $$T_\mathrm{rad}$$ into soil and canopy component temperatures, which are then used to drive the TSEB without requiring a priori assumptions regarding initial canopy transpiration rate. The results showed that a simple contextual algorithm based on the inverse relationship of a vegetation index and $$T_\mathrm{rad}$$ to derive soil and canopy temperatures yielded the closest agreement with flux tower measurements. The utility in very high-resolution remote-sensing data for estimating ET and E and T partitioning at the canopy level is also discussed.

The unique vertical canopy structure and clumped plant distribution/row structure of vineyards and orchards creates an environment that is likely to cause the wind profile inside the canopy air space to deviate from how it is typically modelled for most crops. This in turn affects the efficiency of turbulent flux exchange and energy transport as well as their partitioning between the plant canopy and soil/substrate layers. The objective of this study was to evaluate a new wind profile formulation in the canopy air space that explicitly considers the unique vertical variation in plant biomass of vineyards. The validity of the new wind profile formulation was compared to a simpler wind attenuation profile that assumes attenuation through a homogeneous canopy. We evaluated both attenuation models using measurements of wind speed in a vineyard interrow, as well as turbulent flux estimates retrieved from a two-source energy balance model, which uses land surface temperature as the key boundary condition for flux estimation. This is relevant in developing a robust remote sensing-based energy balance modelling system for accurately monitoring vineyard water use or evapotranspiration that can be applied using satellite and airborne imagery for field-to-regional scale applications. These tools are needed in intensive agricultural production regions with arid climates such as the Central Valley of California, which experiences water shortages during extended drought periods requiring an effective water management policy based on robust water use estimates for allocating water resources. Results showed that the new wind profile model improved sensible heat flux estimates (RMSE reduction from 42 to 35 $$\text{W}\,\text{m}^{-2}$$ ) when the vine canopy is in early growth stage resulting in a strongly clumped canopy.

The goal of drought-related weather index insurance (WII) is to protect smallholder farmers against the risk of weather shocks and to increase their agricultural productivity. Estimates of precipitation and vegetation greenness are the two dominant satellite datasets. However, ignoring additional moisture- and energy-related processes that influence the response of vegetation to rainfall leads to an incomplete representation of the hydrologic cycle. This study evaluates the added value of considering multiple independent satellite-based variables to design, calibrate, and validate weather insurance indices on the African continent. The satellite data include two rainfall datasets, soil moisture, the evaporative stress index (ESI), and vegetation greenness. We limit artificial advantages by resampling all datasets to the same spatial (0.25°) and temporal (monthly) resolution, although datasets with a higher spatial resolution might have an added value, if considered as the single source of information for localized applications. A higher correlation coefficient between the moisture-focused variables and the normalized difference vegetation index (NDVI), an indicator for vegetation vigor, provides evidence for the datasets’ capability to capture agricultural drought conditions on the ground. The Climate Hazards Group Infrared Precipitation with Stations (CHIRPS) rainfall dataset, soil moisture, and ESI show higher correlations with the (lagged) NDVI in large parts of Africa, for different land covers and various climate zones, than the African Rainfall Climatology, version 2 (ARC2), rainfall dataset, which is often used in WII. A comparison to drought years as reported by farmers in Ethiopia, Senegal, and Zambia indicates a high “hit rate” of all satellite-derived anomalies regarding the detection of severe droughts but limitations regarding moderate drought events.

Water-limiting conditions in many California vineyards necessitate assessment of vine water stress to aid irrigation management strategies and decisions. This study was designed to evaluate the utility of a Crop Water Stress Index (CWSI) using multiple canopy temperature sensors and to study the diurnal signature in the stress index of an irrigated vineyard. A detailed instrumentation package comprised of eddy covariance instrumentation, ancillary surface energy balance components, soil water content sensors and a unique multi-canopy temperature sensor array were deployed in a production vineyard near Lodi, CA. The instrument package was designed to measure and monitor hourly growing season turbulent fluxes of heat and water vapor, radiation, air temperature, soil water content directly beneath a vine canopy, and vine canopy temperatures. April 30–May 02, June 10–12 and July 27–28, 2016 were selected for analysis as these periods represented key vine growth and production phases. Considerable variation in computed CWSI was observed between each of the hourly average individual canopy temperature sensors throughout the study; however, the diurnal trends remained similar: highest CWSI values in morning and lowest in the late afternoon. While meteorological conditions were favorable for plant stress to develop, soil water content near field capacity due to frequent irrigation allowed high evapotranspiration rates resulting in downward trending CWSI values during peak evaporative demand. While the CWSI is typically used to evaluate plant stress under the conditions of our study, the trend of the CWSI suggested a lowering of plant water stress as long as there was adequate soil water available to meet atmospheric demand.

As a key variable in the climate system, soil moisture (SM) plays a central role in the earth’s terrestrial water, energy, and biogeochemical cycles through its coupling with surface latent heat flux (LH). Despite the need to accurately represent SM/LH coupling in earth system models, we currently lack quantitative, observation-based, and unbiased estimates of its strength. Here, we utilize the triple collocation (TC) approach introduced in Crow et al. (2015) to SM and LH products obtained from multiple satellite remote sensing platforms and land surface models (LSMs) to obtain unbiased global maps of SM/LH coupling strength. Results demonstrate that, relative to coupling strength estimates acquired directly from remote sensing-based datasets, the application of TC generally enhances estimates of warm-season SM/LH coupling, especially in the western United States, the Sahel, Central Asia, and Australia. However, relative to triple collocation estimates, LSMs (still) over-predict SM/LH coupling strength along transitional climate regimes between wet and dry climates, such as the central Great Plains of North America, India, and coastal Australia. Specific climate zones with biased relations in LSMs are identified to geographically focus the re-examination of LSM parameterizations. TC-based coupling strength estimates are robust to our choice of LSM contributing SM and LH products to the TC analysis. Given their robustness, TC-based coupling strength estimates can serve as an objective benchmark for investigating model predicted SM/LH coupling.

The principle of maximum entropy (POME) can be used to develop vertical soil moisture (SM) profiles. The minimal inputs required by the POME model make it an excellent choice for remote sensing applications. Two of the major input requirements of the POME model are the surface boundary condition and profile-mean moisture content. Microwave-based SM estimates from the Advanced Microwave Scanning Radiometer (AMSR-E) can supply the surface boundary condition whereas thermal infrared-based moisture estimated from the Atmospheric Land EXchange Inverse (ALEXI) surface energy balance model can provide the mean moisture condition. A disaggregation approach was followed to downscale coarse-resolution (∼25 km) microwave SM estimates to match the finer resolution (∼5 km) thermal data. The study was conducted over multiple years (2006–2010) in the southeastern US. Disaggregated soil moisture estimates along with the developed profiles were compared with the Noah land surface model (LSM), as well as in situ measurements from 10 Natural Resource Conservation Services (NRCS) Soil Climate Analysis Network (SCAN) sites spatially distributed within the study region. The overall disaggregation results at the SCAN sites indicated that in most cases disaggregation improved the temporal correlations with unbiased root mean square differences (ubRMSD) in the range of 0.01–0.09 m3 m−3. The profile results at SCAN sites showed a mean bias of 0.03 and 0.05 (m3 m−3); ubRMSD of 0.05 and 0.06 (m3 m−3); and correlation coefficient of 0.44 and 0.48 against SCAN observations and Noah LSM, respectively. Correlations were generally highest in agricultural areas where values in the 0.6–0.7 range were achieved.

With satellite soil moisture (SM) retrievals becoming widely and continuously available, we aim to develop a method to objectively integrate the drought indices into one that is more accurate and consistently reliable. The datasets used in this paper include the Noah land surface model‐based SM estimations, Atmosphere‐Land‐Exchange‐Inverse model‐based Evaporative Stress Index, and the satellite SM products from the Advanced Scatterometer, WindSat, Soil Moisture and Ocean Salinity, and Soil Moisture Operational Product System. Using the Triple Collocation Error Model (TCEM) to quantify the uncertainties of these data, we developed an optically blended drought index (BDI_b) that objectively integrates drought estimations with the lowest TCEM‐derived root‐mean‐square‐errors in this paper. With respect to the reported drought records and the drought monitoring benchmarks including the U.S. Drought Monitor, the Palmer Drought Severity Index and the standardized precipitation evapotranspiration index products, the BDI_b was compared with the sample average blending drought index (BDI_s) and the RMSE‐weighted average blending drought indices (BDI_w). Relative to the BDI_s and the BDI_w, the BDI_b performs more consistently with the drought monitoring benchmarks. With respect to the official drought records, the developed BDI_b shows the best performance on tracking drought development in terms of time evolution and spatial patterns of 2010‐Russia, 2011‐USA, 2013‐New Zealand droughts and other reported agricultural drought occurrences over the 2009‐2014 period. These results suggest that model simulations and remotely sensed observations of SM can be objectively translated into useful information for drought monitoring and early warning, in turn can reduce drought risk and impacts.

Advanced Land Surface Models (LSM) offer a powerful tool for studying hydrological variability. Highly managed systems, however, present a challenge for these models, which typically have simplified or incomplete representations of human water use. Here we examine recent groundwater declines in the US High Plains Aquifer (HPA), a region that is heavily utilized for irrigation and that is also affected by episodic drought. To understand observed decline in groundwater and terrestrial water storage during a recent multi-year drought, we modify the Noah-MP LSM to include a groundwater irrigation scheme. To account for seasonal and interannual variability in active irrigated area, we apply a monthly time-varying greenness vegetation fraction (GVF) dataset within the model. A set of five experiments were performed to study the impact of groundwater irrigation on the simulated hydrological cycle of the HPA and to assess the importance of time-varying GVF when simulating drought conditions. The results show that including the groundwater irrigation scheme improves model agreement with ALEXI ET data, mascon-based GRACE TWS data and depth-to-groundwater measurements in the southern HPA, including Texas and Kansas, and that accounting for time-varying GVF is important for model realism under drought. Results for the HPA in Nebraska are mixed, likely due to the model's weaknesses in representing subsurface hydrology in this region. This study highlights the value of GRACE datasets for model evaluation and development and the potential to advance the dynamic representations of the interactions between human water use and the hydrological cycle.

Land degradation (LD) is one of the most catastrophic outcomes of long-lasting drought events and anthropogenic activities. Assessing climate and human-induced impacts on land can provide information for decision makers to mitigate the effects of these phenomena. The Northeastern region of Brazil (NEB) is the most populous dryland on the planet, making it a highly vulnerable ecosystem especially when considering the lingering drought that started in 2012. The present work consisted of detecting trends in biomass [leaf area index (LAI)] anomalies as indicators of LD in NEB. We also assessed how the loss of vegetation impacts the LD cycle, by measuring trends in albedo and evapotranspiration (ET). LAI, albedo and ET data were derived from MODIS sensors at 8-day temporal and 500 m spatial resolutions. For precipitation anomalies, we relied on CHIRPS-v2 10-day temporal at 5 km spatial resolution data. For detecting trends, we applied the Theil-Sen slope analysis on time series of MODIS LAI, albedo and ET images. Trend analysis was performed for the periods ranging from 2002–2012 (no severe droughts) to 2002–2016 (including the last drought). LAI trends were more pronounced and had a stronger signal than ET and albedo, therefore, LAI was our choice for mapping LD. The first analysis highlighted the human-induced LD prone areas whereas the last detected drought-induced LD prone areas. Considering only the trending areas, which was about 23.4% of the total, 4.5% of this area has undergone human-induced degradation whereas drought was responsible for 73%, although, not mutually exclusive. As reported in the literature and official data, grazing intensification might be a factor driving human-induced degradation. We noticed that the range of variation of LAI is narrow and even narrower for albedo, which demonstrates that land surface response is more influenced by soil reflectivity rather than the characteristic sparse vegetation coverage (LAI ranging from 0.04 to 0.4 in the Caatinga biome), which can barely alter albedo. Finally, the effects of LD on ET anomalies were assessed by Granger causality and impulse-response analyses as means to link land surface feature changes to the hydrological cycle. Albedo had a slightly weaker impulse than LAI on ET whereas precipitation played a major role. These relations are site-specific and, land surface features (biomass and albedo) showed to have a more substantial influence on ET in severely degraded areas. We concluded that drought led to trends indicating LD prone areas in NEB and the degradation cycle has positive feedback derived from ET reduction resulting in an increased net moisture deficit, although the latter statement has yet to be further investigated. The study warns of the desertification risk that NEB is facing and the need for the authorities to take action to mitigate degradation and drought effects on both traditionally surveyed (desertification nuclei) and newfound LD prone areas. We also highlight the limitation of confirming LD, as to date there is no post-drought data available and, lessons learned from the Sahel case make us cautious about claiming that an area is in fact degraded.

The Evaporative Stress Index (ESI) quantifies temporal anomalies in a normalized evapotranspiration (ET) metric describing the ratio of actual-to-reference ET (fRET) as derived from satellite remote sensing. At regional scales (3–10 km pixel resolution), the ESI has demonstrated the capacity to capture developing crop stress and impacts on regional yield variability in water-limited agricultural regions. However, its performance in some regions where the vegetation cycle is intensively managed appears to be degraded due to spatial and temporal limitations in the standard ESI products. In this study, we investigated potential improvements to ESI by generating maps of ET, fRET, and fRET anomalies at high spatiotemporal resolution (30-m pixels, daily time steps) using a multi-sensor data fusion method, enabling separation of landcover types with different phenologies and resilience to drought. The study was conducted for the period 2010–2014 covering a region around Mead, Nebraska that includes both rainfed and irrigated crops. Correlations between ESI and measurements of maize yield were investigated at both the field and county level to assess the potential of ESI as a yield forecasting tool. To examine the role of crop phenology in yield-ESI correlations, annual input fRET time series were aligned by both calendar day and by biophysically relevant dates (e.g. days since planting or emergence). At the resolution of the operational U.S. ESI product (4 km), adjusting fRET alignment to a regionally reported emergence date prior to anomaly computation improves r2 correlations with county-level yield estimates from 0.28 to 0.80. At 30-m resolution, where pure maize pixels can be isolated from other crops and landcover types, county-level yield correlations improved from 0.47 to 0.93 when aligning fRET by emergence date rather than calendar date. Peak correlations occurred 68 days after emergence, corresponding to the silking stage for maize when grain development is particularly sensitive to soil moisture deficiencies. The results of this study demonstrate the utility of remotely sensed ET in conveying spatially and temporally explicit water stress information to yield prediction and crop simulation models.

The Soil Moisture Active Passive (SMAP) mission is dedicated toward global soil moisture mapping. Typically, an L-band microwave radiometer has spatial resolution on the order of 36-40 km, which is too coarse for many specific hydro-meteorological and agricultural applications. With the failure of the SMAP active radar within three months of becoming operational, an intermediate (9-km) and finer (3-km) scale soil moisture product solely from the SMAP mission is no longer possible. Therefore, the focus of this study is a disaggregation of the 36-km resolution SMAP passive-only surface soil moisture (SSM) using the Soil Evaporative Efficiency (SEE) approach to spatial scales of 3-km and 9-km. The SEE was computed using thermal-infrared (TIR) estimation of surface evaporation over Continental U.S. (CONUS). The disaggregation results were compared with the 3 months of SMAP-Active (SMAP-A) and Active/Passive (AP) products, while comparisons with SMAP-Enhanced (SMAP-E), SMAP-Passive (SMAP-P), as well as with more than 180 Soil Climate Analysis Network (SCAN) stations across CONUS were performed for a 19 month period. At the 9-km spatial scale, the TIR-Downscaled data correlated strongly with the SMAP-E SSM both spatially (r = 0.90) and temporally (r = 0.87). In comparison with SCAN observations, overall correlations of 0.49 and 0.47; bias of −0.022 and −0.019 and unbiased RMSD of 0.105 and 0.100 were found for SMAP-E and TIR-Downscaled SSM across the Continental U.S., respectively. At 3-km scale, TIR-Downscaled and SMAP-A had a mean temporal correlation of only 0.27. In terms of gain statistics, the highest percentage of SCAN sites with positive gains (>55%) was observed with the TIR-Downscaled SSM at 9-km. Overall, the TIR-based downscaled SSM showed strong correspondence with SMAP-E; compared to SCAN, and overall both SMAP-E and TIR-Downscaled performed similarly, however, gain statistics show that TIR-Downscaled SSM slightly outperformed SMAP-E.

Given the increasing use of the term “flash drought” by the media and scientific community, it is prudent to develop a consistent definition that can be used to identify these events and to understand their salient characteristics. It is generally accepted that flash droughts occur more often during the summer owing to increased evaporative demand; however, two distinct approaches have been used to identify them. The first approach focuses on their rate of intensification, whereas the second approach implicitly focuses on their duration. These conflicting notions for what constitutes a flash drought (i.e., unusually fast intensification vs short duration) introduce ambiguity that affects our ability to detect their onset, monitor their development, and understand the mechanisms that control their evolution. Here, we propose that the definition for “flash drought” should explicitly focus on its rate of intensification rather than its duration, with droughts that develop much more rapidly than normal identified as flash droughts. There are two primary reasons for favoring the intensification approach over the duration approach. First, longevity and impact are fundamental characteristics of drought. Thus, short-term events lasting only a few days and having minimal impacts are inconsistent with the general understanding of drought and therefore should not be considered flash droughts. Second, by focusing on their rapid rate of intensification, the proposed “flash drought” definition highlights the unique challenges faced by vulnerable stakeholders who have less time to prepare for its adverse effects.

This work aims to assess a diagnostic approach which links evapotranspiration (ET) to land surface temperature (LST) measured by thermal remote sensing in the Alps. We estimated gridded ET, from field (30 m) to regional (1 km) scales, and we performed a specific study on grassland ecosystems in the Alps in South Tyrol (Italy), to evaluate the model sensitivity to different types of land management. The energy balance model TSEB ALEXI (Two Source Energy Balance Atmosphere Land EXchange Inverse) was first applied to Meteosat satellite data. Then ET was estimated by the flux disaggregation procedure DisALEXI, driven by MODIS and Landsat LST retrievals, which has never been applied before in a mountain region. We evaluated the model against eddy-covariance (EC) measurements from established stations in the Alps, and analyzed the main limitations which affect the model performance in mountainous regions. The TSEB model, applied in plot-scale mode using tower-based meteorological and LST input data, performed well with errors in daytime (6–18 UTC+1) latent heat flux around 30–60 W m−2 in comparison with flux measurements corrected for the lack of closure in the energy balance. For landscape ET retrievals, while Landsat resolution (30 m) is preferable for capturing small-scale heterogeneity in landscape moisture conditions, and for direct comparison with tower fluxes, persistent cloud cover resulted in no clear Landsat scenes during the study period. MODIS-based retrievals at 1 km resolution are too coarse to resolve the flux tower footprint in this complex landscape, yielding discrepancies of 100 W m−2 in model-measurement comparisons. Still, MODIS DisALEXI partitioning of the energy budget was reasonable and enabled to detect evaporative stress at regional scale expressed as the ratio between actual and potential ET, fPET. We evaluated fPET in comparison with a crop stress index based on cumulative air temperature and precipitation at different stations in the study area, and investigated ability to capture differential responses between managed and unmanaged grasslands. Results show that in the Alps i) moderate resolution thermal data can be used to monitor evaporative stress at the regional scale; ii) the spatial-temporal evolution of ET can be characterized from MODIS and Landsat thermal data with limitations which are due to the low availability of clear-sky scenes and to the small-scale (∼10 m) changes in soil moisture, topography and canopy density, which control ET patterns in mountainous regions; iii) solar radiation and leaf area index are critical variables which strongly affect the accuracy of the modeled energy fluxes.

Accurate gridded estimates of evapotranspiration (ET) are essential to the analysis of terrestrial water budgets. In this study, ET estimates from three gridded energy balance–based products (ETEB) with independent model formations and data forcings are evaluated for their ability to capture long-term climatology and interannual variability in ET derived from a terrestrial water budget (ETWB) for 671 gauged basins across the contiguous United States. All three ETEB products have low spatial bias and accurately capture interannual variability of ETWB in the central United States, where ETEB and ancillary estimates of change in total surface water storage (ΔTWS) from the GRACE satellite project appear to close terrestrial water budgets. In humid regions, ETEB products exhibit higher long-term bias, and the covariability of ETEB and ETWB decreases significantly. Several factors related to either failure of ETWB, such as errors in ΔTWS and precipitation, or failure of ETEB, such as treatment of snowfall and horizontal heat advection, explain some of these discrepancies. These results mirror and build on conclusions from other studies: on interannual time scales, ΔTWS and error in precipitation estimates are nonnegligible uncertainties in ET estimates based on a terrestrial water budget, and this confounds their comparison to energy balance ET models. However, there is also evidence that in at least some regions, climate and landscape features may also influence the accuracy and long-term bias of ET estimates from energy balance models, and these potential errors should be considered when using these gridded products in hydrologic applications.

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

Crop emergence is a critical stage for crop development modeling, crop condition monitoring, and biomass accumulation estimation. Green-up dates (or the start of the season) detected from remote sensing time series are related to, but generally lag, crop emergence dates. In this paper, we refine the within-season emergence (WISE) algorithm and extend application to five Corn Belt states (Iowa, Illinois, Indiana, Minnesota, and Nebraska) using routine harmonized Landsat and Sentinel-2 (HLS) data from 2018 to 2020. Green-up dates detected from the HLS time series were assessed using field observations and near-surface measurements from PhenoCams. Statistical descriptions of green-up dates for corn and soybeans were generated and compared to county-level planting dates and district- to state-level crop emergence dates reported by the National Agricultural Statistics Service (NASS). Results show that emergence dates for corn and soybean can be reliably detected within the season using the HLS time series acquired during the early growing season. Compared to observed crop emergence dates, green-up dates from HLS using WISE were ~3 days later at the field scale (30-m). The mean absolute difference (MAD) was ~7 days and the root mean square error (RMSE) was ~9 days. At the state level, the mean differences between median HLS green-up date and median crop emergence date were within 2 days for 2018–2020. At this scale, MAD was within 4 days, and RMSE was less than 5 days for both corn and soybeans. The R-squares were 0.73 and 0.87 for corn and soybean, respectively. The 2019 late emergence of crops in Corn Belt states (1–4 weeks to five-year average) was captured by HLS green-up date retrievals. This study demonstrates that routine within-season mapping of crop emergence/green-up at the field scale is practicable over large regions using operational satellite data. The green-up map derived from HLS during the growing season provides valuable information on spatial and temporal variability in crop emergence that can be used for crop monitoring and refining agricultural statistics used in broad-scale modeling efforts.

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.