Droughts impede water balance recovery from fires in the Western United States.

A steady rise in fires in the Western United States, coincident with intensifying droughts, imparts substantial modifications to the underlying vegetation, hydrology and overall ecosystem. Drought can compound the ecosystem disturbance caused by fire, although how these compound effects on hydrologic and ecosystem recovery vary among ecosystems is poorly understood. Here we use remote sensing-derived high-resolution evapotranspiration (ET) estimates from before and after 1,514 fires to show that ecoregions dominated by grasslands and shrublands are more susceptible to drought, which amplifies fire-induced ET decline and, subsequently, shifts water flux partitioning. In contrast, severely burned forests recover from fire slowly or incompletely, but are less sensitive to dry extremes. We conclude that moisture limitation caused by droughts influences the dynamics of water balance recovery in post-fire years. This finding explains why moderate to extreme droughts aggravate impacts on the water balance in non-forested vegetation, while moisture accessed by deeper roots in forests helps meet evaporative demands unless severe burns disrupt internal tree structure and deplete fuel load availability. Our results highlight the dominant control of drought on altering the resilience of vegetation to fires, with critical implications for terrestrial ecosystem stability in the face of anthropogenic climate change in the West.

Global patterns of water storage in the rooting zones of vegetation.

The rooting-zone water-storage capacity-the amount of water accessible to plants-controls the sensitivity of land-atmosphere exchange of water and carbon during dry periods. How the rooting-zone water-storage capacity varies spatially is largely unknown and not directly observable. Here we estimate rooting-zone water-storage capacity globally from the relationship between remotely sensed vegetation activity, measured by combining evapotranspiration, sun-induced fluorescence and radiation estimates, and the cumulative water deficit calculated from daily time series of precipitation and evapotranspiration. Our findings indicate plant-available water stores that exceed the storage capacity of 2-m-deep soils across 37% of Earth's vegetated surface. We find that biome-level variations of rooting-zone water-storage capacities correlate with observed rooting-zone depth distributions and reflect the influence of hydroclimate, as measured by the magnitude of annual cumulative water-deficit extremes. Smaller-scale variations are linked to topography and land use. Our findings document large spatial variations in the effective root-zone water-storage capacity and illustrate a tight link among the climatology of water deficits, rooting depth of vegetation and its sensitivity to water stress.

Improving the spatiotemporal resolution of remotely sensed ET information for water management through Landsat, Sentinel-2, ECOSTRESS and VIIRS data fusion.

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.

Changes in land use enhance the sensitivity of tropical ecosystems to fire-climate extremes.

The Pantanal, the largest contiguous wetland in the world with a high diversity of ecosystems and habitat for several endangered species, was impacted by record-breaking wildfires in 2020. In this study, we integrate satellite and modeling data that enable exploration of natural and human contributing factors to the unprecedented 2020 fires. We demonstrate that the fires were fueled by an exceptional multi-year drought, but dry conditions solely could not explain the spatial patterns of burning. Our analysis reveals how human-caused fires exacerbated drought effects on natural ecosystem within the Pantanal, with large burned fractions primarily over natural (52%), and low cattle density areas (44%) in 2020. The post-fire ecosystem and hydrology changes also had strong ecological effects, with vegetation productivity less than - 1.5 σ over more than 30% of the natural and conservation areas. In contrast to more managed areas, there was a clear decrease in evaporation (by ~ 9%) and an increase in runoff (by ~ 5%) over the natural areas, with long-term impacts on ecosystem recovery and fire risk. This study provides the first tropical evidence outside rainforests of the synergy between climate, land management and fires, and the associated impacts on the ecosystem and hydrology over the largest contiguous wetlands in the world.

Data Assimilation of High-Resolution Thermal and Radar Remote Sensing Retrievals for Soil Moisture Monitoring in a Drip-Irrigated Vineyard.

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.

Sharpening ECOSTRESS and VIIRS land surface temperature using harmonized Landsat-Sentinel surface reflectances.

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.

Spatial, Temporal, and Electrical Characteristics of Lightning in Reported Lightning-Initiated Wildfire Events.

Analysis was performed to determine if a lightning flash could be associated with every reported lightning-initiated wildfire that grew to at least 4 km2. In total, 905 lightning-initiated wildfires within CONUS between 2012 and 2015 were analyzed. Fixed and fire radius search methods showed that 81-88% of wildfires had a corresponding lightning flash within a 14 day period prior to the report date. The two methods showed that 52-60% of lightning-initiated wildfire were reported on the same day as the closest lightning flash. The fire radius method indicated the most promising spatial results, where the median distance between the closest lightning and the wildfire start location was 0.83 km, followed by a 75th percentile of 1.6 km, and a 95th percentile of 5.86 km. Ninety percent of the closest lightning flashes to wildfires were negative polarity. Maximum flash densities were less than 0.41 flashes km2 for the 24 hour period at the fire start location. The majority of lightning-initiated holdover events were observed in the Western CONUS, with a peak density in north-central Idaho. A twelve day holdover event from New Mexico was also discussed; outlining the opportunities and limitations of using lightning data to characterize wildfires.

Microwave implementation of two-source energy balance approach for estimating evapotranspiration.

A newly developed microwave (MW) land surface temperature (LST) product is used to substitute thermal infrared (TIR) based LST in the Atmosphere Land Exchange Inverse (ALEXI) modelling framework for estimating ET from space. ALEXI implements a two-source energy balance (TSEB) land surface scheme in a time-differential approach, designed to minimize sensitivity to absolute biases in input records of LST through the analysis of the rate of temperature change in the morning. Thermal infrared (TIR) retrievals of the diurnal LST curve, traditionally from geostationary platforms, are hindered by cloud cover, reducing model coverage on any given day. This study tests the utility of diurnal temperature information retrieved 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 analysis is based on 9-year long, global records of ALEXI ET generated using both MW and TIR based diurnal LST information as input. 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 from both LST sources agree well, 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. Weekly ET output from the two parallel ALEXI implementations is further compared to a common ground measured reference provided by the FLUXNET consortium. Overall, the two model implementations generate similar performance metrics (correlation and RMSE) for all but the most challenging sites in terms of spatial heterogeneity and level of aridity. It is concluded that a constellation of MW satellites can effectively be used to provide LST for estimating ET through ALEXI, which is an important step towards all-sky satellite-based retrieval of ET using an energy balance framework.

Developing a remotely sensed drought monitoring indicator for Morocco.

Drought is one of the most serious climatic and natural disasters inflicting serious impacts on the socio-economy of Morocco, which is characterized both by low-average annual rainfall and high irregularity in the spatial distribution and timing of precipitation across the country. This work aims to develop a comprehensive and integrated method for drought monitoring based on remote sensing techniques. The main input parameters are derived monthly from satellite data at the national scale and are then combined to generate a composite drought index presenting different severity classes of drought. The input parameters are: Standardized Precipitation Index calculated from satellite based precipitation data since 1981 (CHIRPS), anomalies in the day-night difference of Land Surface Temperature as a proxy for soil moisture, Normalized Difference Vegetation Index anomalies from MODIS data and Evapotranspiration anomalies from surface energy balance modeling. All of these satellite-based indices are being used to monitor vegetation condition, rainfall and land surface temperature. The weighted combination of these input parameters into one composite indicator takes into account the importance of the rainfall based parameter (SPI). The composite drought index maps were generated during the growing seasons going back to 2003. These maps have been compared to both the historical, in situ precipitation data across Morocco and with the historical yield data across different provinces with information being available since 2000. The maps are disseminated monthly to several main stakeholders groups including the Ministry of Agriculture and Department of Water in Morocco.

A water balance based, spatiotemporal evaluation of terrestrial evapotranspiration products across the contiguous United States.

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 inter-annual variability in ET derived from a terrestrial water budget (ETWB) for 671 gaged basins across the CONUS. All three ETEB products have low spatial bias and accurately capture inter-annual variability of ETWB in the central US, 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 inter-annual timescales, ΔTWS and error in precipitation estimates are non-negligible 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.

Global Investigation of Soil Moisture and Latent Heat Flux Coupling Strength.

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.

Global relationships among traditional reflectance vegetation indices (NDVI and NDII), evapotranspiration (ET), and soil moisture variability on weekly timescales.

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.

An Assessment of the Impact of Land Thermal Infrared Observation on Regional Weather Forecasts Using Two Different Data Assimilation Approaches.

Recent studies have shown the unique value of satellite-observed land surface thermal infrared (TIR) information (e.g., skin temperature) and the feasibility of assimilating land surface temperature (LST) into land surface models (LSMs) to improve the simulation of land-atmosphere water and energy exchanges. In this study, two different types of LST assimilation techniques are implemented and the benefits from the techniques are compared. One of the techniques is to directly assimilate LST using ensemble Kalman filter (EnKF) data assimilation (DA) utilities. The other is to use the Atmosphere-Land Exchange Inversion model (ALEXI) as an "observation operator" that converts LST retrievals into the soil moisture (SM) proxy based on the ratio of actual to potential evapotranspiration (fPET), which is then assimilated into an LSM. While most current studies have shown some success in both directly the assimilating LST and assimilating ALEXI SM proxy into offline LSMs, the potential impact of the assimilation of TIR information through coupled numerical weather prediction (NWP) models is unclear. In this study, a semi-coupled Land Information System (LIS) and Weather Research and Forecast (WRF) system is employed to assess the impact of the two different techniques for assimilating the TIR observations from NOAA GOES satellites on WRF model forecasts. The NASA LIS, equipped with a variety of LSMs and advanced data assimilation tools (e.g., the ensemble Kalman Filter (EnKF)), takes atmospheric forcing data from the WRF model run, generates updated initial land surface conditions with the assimilation of either LST- or TIR-based SM and returns them to WRF for initializing the forecasts. The WRF forecasts using the daily updated initializations with the TIR data assimilation are evaluated against ground weather observations and re-analysis products. It is found that WRF forecasts with the LST-based SM assimilation have better agreement with the ground weather observations than those with the direct LST assimilation or without the land TIR data assimilation.

The role of local heating in the 2015 Indian Heat Wave.

India faced a major heat wave during the summer of 2015. Temperature anomalies peaked in the dry period before the onset of the summer monsoon, suggesting that local land-atmosphere feedbacks involving desiccated soils and vegetation might have played a role in driving the heat extreme. Upon examination of in situ data, reanalysis, satellite observations, and land surface models, we find that the heat wave included two distinct peaks: one in late May, and a second in early June. During the first peak we find that clear skies led to a positive net radiation anomaly at the surface, but there is no significant sensible heat flux anomaly within the core of the heat wave affected region. By the time of the second peak, however, soil moisture had dropped to anomalously low levels in the core heat wave region, net surface radiation was anomalously high, and a significant positive sensible heat flux anomaly developed. This led to a substantial local forcing on air temperature that contributed to the intensity of the event. The analysis indicates that the highly agricultural landscape of North and Central India can reinforce heat extremes under dry conditions.

Estimating Morning Change in Land Surface Temperature from MODIS Day/Night Observations: Applications for Surface Energy Balance Modeling.

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.