Large-scale crop monitoring and yield estimation are important for both scientific research and practical applications. Satellite remote sensing provides an effective means for regional and global ...cropland monitoring, particularly in data-sparse regions that lack reliable ground observations and reporting. The conventional approach of using visible and near-infrared based vegetation index (VI) observations has prevailed for decades since the onset of the global satellite era. However, other satellite data encompass diverse spectral ranges that may contain complementary information on crop growth and yield, but have been largely understudied and underused. Here we conducted one of the first attempts at synergizing multiple satellite data spanning a diverse spectral range, including visible, near-infrared, thermal and microwave, into one framework to estimate crop yield for the U.S. Corn Belt, one of the world's most important food baskets. Specifically, we included MODIS Enhanced VI (EVI), estimated Gross Primary Production based on GOME-2 solar-induced fluorescence (SIF-GPP), thermal-based ALEXI Evapotranspiration (ET), QuikSCAT Ku-band radar backscatter, and AMSR-E X-band passive microwave Vegetation Optical Depth (VOD) in this study, benchmarked on USDA county-level crop yield statistics. We used Partial Least Square Regression (PLSR), an effective statistical model for dimension reduction, to distinguish commonly shared and unique individual information from the various satellite data and other ancillary climate information for crop yield estimation. In the PLSR model that includes all of the satellite data and climate variables from 2007 to 2009, we assessed the first two major PLSR components and found that the first component (an integrated proxy of crop aboveground biomass) explained 82% variability of modelled crop yield, and the second component (dominated by environmental stresses) explained 15% variability of modelled crop yield. We found that most of the satellite derived metrics (e.g. SIF-GPP, radar backscatter, EVI, VOD, ALEXI-ET) share common information related to aboveground crop biomass (i.e. the first component). For this shared information, the SIF-GPP and backscatter data contain almost the same amount of information as EVI at the county scale. When removing the above shared component from all of the satellite data, we found that EVI and SIF-GPP do not provide much extra information; instead, Ku-band backscatter, thermal-based ALEXI-ET, and X-band VOD provide unique information on environmental stresses that improves overall crop yield predictive skill. In particular, Ku-band backscatter and associated differences between morning and afternoon overpasses contribute unique information on crop growth and environmental stress. Overall, using satellite data from various spectral bands significantly improves regional crop yield predictions. The additional use of ancillary climate data (e.g. precipitation and temperature) further improves model skill, in part because the crop reproductive stage related to harvest index is highly sensitive to environmental stresses but they are not fully captured by the satellite data used in our study. We conclude that using satellite data across various spectral ranges can improve monitoring of large-scale crop growth and yield beyond what can be achieved from individual sensors. These results also inform the synergistic use and development of current and next generation satellite missions, including NASA ECOSTRESS, SMAP, and OCO-2, for agricultural applications.
•A novel framework to estimate crop yield for US Corn Belt is provided.•Satellite optical, fluorescence, thermal and microwave data are synergized.•Share and unique information of different satellite data is identified using PLSR.•Climate information is necessary to be included to quantify harvest index.•A synergized use of various-spectral satellite data improves yield estimation.
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 midmorning rise in LST from a single sensor (two 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 data sets.
Key Points
The morning rise in LST from temporally sparse observations can be accurately estimated through data‐mining techniques
This methodology provides a mechanisms to rapidly mobilize energy balance models to monitor the global energy and water cycle
As a primary flux in the global water cycle, evapotranspiration (ET) connects hydrologic and biological processes and is directly affected by water and land management, land use change and climate ...variability. Satellite remote sensing provides an effective means for diagnosing ET patterns over heterogeneous landscapes; however, limitations on the spatial and temporal resolution of satellite data, combined with the effects of cloud contamination, constrain the amount of detail that a single satellite can provide. In this study, we describe an application of a multi-sensor ET data fusion system over a mixed forested/agricultural landscape in North Carolina, USA, during the growing season of 2013. The fusion system ingests ET estimates from the Two-Source Energy Balance Model (TSEB) applied to thermal infrared remote sensing retrievals of land surface temperature from multiple satellite platforms: hourly geostationary satellite data at 4 km resolution, daily 1 km imagery from the Moderate Resolution Imaging Spectroradiometer (MODIS) and biweekly Landsat thermal data sharpened to 30 m. These multiple ET data streams are combined using the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) to estimate daily ET at 30 m resolution to investigate seasonal water use behavior at the level of individual forest stands and land cover patches. A new method, also exploiting the STARFM algorithm, is used to fill gaps in the Landsat ET retrievals due to cloud cover and/or the scan-line corrector (SLC) failure on Landsat 7. The retrieved daily ET time series agree well with observations at two AmeriFlux eddy covariance flux tower sites in a managed pine plantation within the modeling domain: US-NC2 located in a mid-rotation (20-year-old) loblolly pine stand and US-NC3 located in a recently clear-cut and replanted field site. Root mean square errors (RMSEs) for NC2 and NC3 were 0.99 and 1.02 mm day−1, respectively, with mean absolute errors of approximately 29 % at the daily time step, 12 % at the monthly time step and 0.7 % over the full study period at the two flux tower sites. Analyses of water use patterns over the plantation indicate increasing seasonal ET with stand age for young to mid-rotation stands up to 20 years, but little dependence on age for older stands. An accounting of consumptive water use by major land cover classes representative of the modeling domain is presented, as well as relative partitioning of ET between evaporation (E) and transpiration (T) components obtained with the TSEB. The study provides new insights about the effects of management and land use change on water yield over forested landscapes.
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 25days 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.
•We investigate correlations between Brazilian crop yields and satellite indicators.•The strength and timing of peak correlations are determined for each index.•Correlations are the strongest in states with high interannual yield variability.•Evapotranspiration anomalies provide useful early signals of yield impacts.
California's Central Valley grows a significant fraction of grapes used for wine production in the United States. With increasing vineyard acreage, reduced water availability in much of California, ...and competing water use interests, it is critical to be able to monitor regional water use and evapotranspiration (ET) over large areas, but also in detail at individual field scales to improve water management within these viticulture production systems. This can be achieved by integrating remote sensing data from multiple satellite systems with different spatiotemporal characteristics. In this research, we evaluate the utility of a multi-scale system for monitoring ET as applied over two vineyard sites near Lodi, California during the 2013 growing season, leading into the drought in early 2014. The system employs a multi-sensor satellite data fusion methodology (STARFM: Spatial and Temporal Adaptive Reflective Fusion Model) combined with a multi-scale ET retrieval algorithm based on the Two-Source Energy Balance (TSEB) land-surface representation to compute daily ET at 30m resolution. In this system, TSEB is run using thermal band imagery from the Geostationary Environmental Operational Satellites (GOES; 4-km spatial resolution, hourly temporal sampling), the Moderate Resolution Imaging Spectroradiometer (MODIS) data (1km resolution, daily acquisition) and the new Landsat 8 satellite (sharpened to 30m resolution, ~16day acquisition). Estimates of daily ET generated in two neighboring fields of Pinot noir vines of different age agree with ground-based flux measurements acquired in-field during most of the 2013 season with relative mean absolute errors on the order of 19–23% (root mean square errors of approximately 1mmd−1), reducing to 14–20% at the weekly timestep relevant for irrigation management (~5mmwk−1). A model overestimation of ET in the early season was detected in the younger vineyard, perhaps relating to an inter-row grass cover crop. Spatial patterns of cumulative ET generally correspond to measured yield maps and indicate areas of variable crop moisture, soil condition, and yield within the vineyards that could require adaptive management. The results suggest that multi-sensor remote sensing observations provide a unique means for monitoring crop water use and soil moisture status at field-scales over extended growing regions, and may have value in supporting operational water management decisions in vineyards and other high value crops.
•Multi-scale, multi-sensor fusion methodology estimates vineyard evapotranspiration.•Combining Landsat 8, MODIS, and GOES data provides daily field scale ET estimates.•Modeled surface energy fluxes agree well with ground flux measurements.•Spatial distribution of ET corresponds with yield estimates.•ET overestimation in early season may be due to cover crop between vines.
•The differences between near real time and multi-year mean green vegetation fraction and albedo are examined.•The improved model performance using near real time products is quantified.•The results ...can be used as feedbacks to the model developers and users.
The current operational Noah land surface model (LSM) uses multi-year climatology of monthly green vegetation fraction (GVF) and the multi-year averages of land surface albedo data for several numerical weather predictions at National Centers for Environmental Predictions of National Oceanic and Atmospheric Administration. However, these static GVF and albedo data can only prescribe the multiannual means and lack the ability to capture near real time (NRT) vegetation status and land surface condition. In this study, the impact of NRT GVF and albedo on Noah LSM (version 3.2) performances are examined against in situ measurements of surface net long wave radiation and net short wave radiation from 7 U.S. Surface Radiation Budget Network stations, and soil temperature and soil moisture from 9 USDA Soil Climate Analysis Network sites. Large differences between the NRT GVF/surface albedo and their climatological averages are found over the global, which have significant influences on Noah LSM simulations. With respect to in situ measurements, the Noah LSM simulation improvements from using the weekly GVF data are 19.3% for surface soil moisture, 9.3% for surface soil temperature. The benefits from the weekly GVF and monthly albedo can reach to 2.7Wm−2 for surface net long wave radiation and 2.6Wm−2 for surface net short wave radiation. The results suggest to Noah model developers and users that the NRT GVF and albedo should be used for better model performance.
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.
•Multi-scale, multi-sensor fusion methodology estimates evapotranspiration.•Modeled surface energy fluxes agree well with ground flux measurements.•Evaporative Stress Index (ESI) is correlated to yield estimates at field scale.•Incorporating phenological information significantly improves yield correlations.•ESI at 30-m is able to resolve different crop types with varying phenology.
Twenty-six years of lightning data were paired with over 68,000 lightning-initiated wildfire (LIW) reports to understand lightning flash characteristics responsible for ignition in between 1995 and ...2020. Results indicate that 92% of LIW were started by negative cloud-to-ground (CG) lightning flashes and 57% were single stroke flashes. Moreover, 62% of LIW reports did not have a positive CG within 10 km of the start location, contrary to the science literature’s suggestion that positive CG flashes are a dominant fire-starting mechanism. Nearly 1/3rd of wildfire events were holdovers, meaning one or more days elapsed between lightning occurrence and fire report. However, fires that were reported less than a day after lightning occurrence statistically burned more acreage. Peak current was not found to be a statistically significant delineator between fire starters and non-fire starters for -CGs but was for positive CGs. Results highlighted the need for reassessing the role of positive CG lightning and subsequently long continuing current in wildfire ignition started by lightning. One potential outcome of this study’s results is the development of real-time tools to identify ignition potential during lightning events to aid in fire mitigation efforts.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
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.
Global food production in the developing world occurs within sub-hectare fields that are difficult to identify with moderate resolution satellite imagery. Knowledge about the distribution of these ...fields is critical in food security programs. We developed a semi-automated image segmentation approach using wall-to-wall sub-meter imagery with high-performance computing to map crop area (CA) throughout Tigray, Ethiopia that encompasses over 41,000 km2. Multiple processing streams were tested to minimize mapping error while applying five unique smoothing kernels to capture differences in land surface texture associated to CA. Typically, very-small fields (mean < 2 ha) have a smooth image roughness compared to natural scrub/shrub woody vegetation at the ~1 m scale and these features can be segmented in panchromatic imagery with multi-level histogram thresholding. Multi-temporal very-high resolution (VHR) panchromatic imagery with multi-spectral VHR are sufficient in extracting critical CA information needed in food security programs. A 2011 to 2015 CA map was produced, using over 3000 WorldView-1 panchromatic images wall-to-wall in 1/2° mosaics for Tigray, Ethiopia. CA was evaluated with nearly 3000 WorldView-2 2 m multispectral 250 × 250 m image subsets by seven expert interpretations, and with in-situ global positioning system photography. CA estimates ranged from 32 to 41% in sub regions of Tigray with median maximum per bin commission and omission errors of 11% and 1% respectively, with most of the error occurring in bins <15%. This empirical, simple, and low direct cost approach via U.S. government license agreement to access commercial VHR data, could be a viable big-data high-performance computing methodology to extract wall-to-wall CA for other regions of the world that have very-small agriculture fields with similar image texture.
•WorldView-1 was segmented wall-to-wall in Tigray Ethiopia to estimate crop area.•Subregion crop area estimates ranged from 32 to 41%.•Classification evaluation used nearly 3000 250 × 250 WorldView-2 subsets.•Median commission and omission error from all bins was 11 and 1% respectively.•This is a simple approach to map sub-hectare agriculture fields in Tigray Ethiopia.