Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “forecast error”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

Status of Cycle 23 Forecasts

Forecasts for the amplitude of cycle 23 that were reported prior to the start of, the cycle covered a full range of, values from very small to very large. A forecast reached by the consensus of a panel of forecasters convened at the time of minimum in 1996 [Joselyn et al., 1997] suggested that this cycle would be much larger than average with the smoothed International Sunspot Number reaching a maximum of 160 +/- 30 in the middle of the year 2000. A recent survey of solar cycle prediction techniques [Hathaway et al., 1999] found that the two most reliable techniques for forecasting the cycle prior to its start give similar predictions for this cycle's maximum -- 154 +/- 26 and 153 +/- 33. Curve-fitting and regression techniques can be used with some confidence now that cycle 23 is well underway. These techniques indicate a more modest sunspot cycle with a maximum of 112 + 10 -- only slightly larger than average. The current (May 2000) prediction using the combined predictions from both precursors and curve-fitting gives a cycle amplitude of about 135 +/- 20. This is within the errors given by the consensus and the precursor technique predictions but very close to their lower bounds. Inspection of the characteristics of cycle 23 thus far indicates that in some respects this cycle is an outlier but in other respects it is quite normal.

Hathaway, David H.↗

Predictability of quasi-geostrophic planetary waves in global and hemispheric domains

Factors which affect the long term predictability, and also the feasibility of defining an ultimate level of weather predictability, are discussed. It is noted that current models yield a one week forecast at the best. Long range predictability, however, is a restriction on the knowledge of evolutive processes of the atmosphere, and does not characterize the prediction of specific atmospheric quantities. The influences of local phenomena such as cyclones and other nonlinear processes can cause planetary waves which affect the predictability of weather everywhere. It is suggested that any successful prediction model will necessarily account for the evolution of ultralong waves. Although significant differences in predictions have been found between hemispheric and global predictions, the location of the errors in extrapolating the hemispheric to the global scale have been identified.

Roads, J. O.↗

A High-Resolution Merged Wind Dataset for DYNAMO: Progress and Future Plans

In order to support research on optimal data assimilation methods for the Cyclone Global Navigation Satellite System (CYGNSS), launching in 2016, work has been ongoing to produce a high‐resolution merged wind dataset for the Dynamics of the Madden Julian Oscillation (DYNAMO) field campaign, which took place during late 2011/early 2012. The winds are produced by assimilating DYNAMO observations into the Weather Research and Forecasting (WRF) three‐dimensional variational (3DVAR) system. Data sources from the DYNAMO campaign include the upper‐air sounding network, radial velocities from the radar network, vector winds from the Advanced Scatterometer (ASCAT) and Oceansat‐2 Scatterometer (OSCAT) satellite instruments, the NOAA High Resolution Doppler Lidar (HRDL), and several others. In order the prep them for 3DVAR, significant additional quality control work is being done for the currently available TOGA and SMART‐R radar datasets, including automatically dealiasing radial velocities and correcting for intermittent TOGA antenna azimuth angle errors. The assimilated winds are being made available as model output fields from WRF on two separate grids with different horizontal resolutions ‐ a 3‐km grid focusing on the main DYNAMO quadrilateral (i.e., Gan Island, the R/V Revelle, the R/V Mirai, and Diego Garcia), and a 1‐km grid focusing on the Revelle. The wind dataset is focused on three separate approximately 2‐week periods during the Madden Julian Oscillation (MJO) onsets that occurred in October, November, and December 2011. Work is ongoing to convert the 10‐m surface winds from these model fields to simulated CYGNSS observations using the CYGNSS End‐To‐End Simulator (E2ES), and these simulated satellite observations are being compared to radar observations of DYNAMO precipitation systems to document the anticipated ability of CYGNSS to provide information on the relationships between surface winds and oceanic precipitation at the mesoscale level. This research will improve our understanding of the future utility of CYGNSS for documenting key MJO processes.

Lang, Timothy J.↗

Fisher's Peak Ecological Forecasting: Mapping Biomass to Inform Conservation Planning of a Future State Park in Southern Colorado

Fisher’s Peak is a 77.5 km2 property southeast of Trinidad, Colorado that is planned to become Colorado’s newest state park. The area has experienced limited anthropogenic disturbance and is home to an abundance of unique habitats and species. A rapid, approximately 900 m change in elevation over the extent of the area nurtures a variety of plants and animals, including the endangered New Mexico meadow jumping mouse. In 2019, the State of Colorado obtained Fisher’s Peak with plans to make it Colorado’s second largest state park. A diverse group of collaborators, including the Colorado State Forest Service and The Nature Conservancy, worked closely to design the state park to maximize recreation opportunity while conserving the property’s rich habitats and biodiversity. The Fisher’s Peak Ecological Forecasting Team utilized Light Detection and Ranging (LiDAR) surveys, in situ forest inventory data, and Earth observations from Landsat8 Operational Land Imager (OLI), Sentinel-1 C-band Synthetic Aperture Radar (C-SAR), Sentinel-2 Multispectral Instrument (MSI), Advanced Land Observing Satellite 2 (ALOS-2) Phased Array type L-band Synthetic Aperture Radar (PALSAR-2) and the Shuttle Radar Topography Mission (SRTM) to quantify and map biomass over the extent of the study area. The results from modeling biomass had an out-of-bag root mean square error of 55 Mg/ha and an R2 of 12. The resulting map indicates areas where carbon storage on the property is high, informing decision-making processes for future park development. While more in situ training data may improve modeling capacity for biomass in the Fisher’s Peak area, this work represents a feasible attempt to better understand biomass distribution using earth observation.

DEVELOP Tech Paper↗

The 4th order GISS model of the global atmosphere

The new GISS 4th order model of the global atmosphere is described. It is based on 4th order quadratically conservative differences with the periodic application of a 16th order filter on the sea level pressure and potential temperature equations, a combination which is approximately enstrophy conserving. Several short range forecasts indicate a significant improvement over 2nd order forecasts with the same resolution (approximately 400 km). However the 4th order forecasts are somewhat inferior to 2nd order forecasts with double resolution. This is probably due to the presence of short waves in the range between 1000 km and 2000 km, which are computed more accurately by the 2nd order high resolution model. An operation count of the schemes indicates that with similar code optimization, the 4th order model will require approximately the same amount of computer time as the 2nd order model with the same resolution. It is estimated that the 4th order model with a grid size of 200 km provides enough accuracy to make horizontal truncation errors negligible over a period of a week for all synoptic scales (waves longer than 1000 km).

Kalnay-Rivas, E.↗

Fisher's Peak Ecological Forecasting - Mapping Biomass to Inform Conservation Planning of a Future State Park in Southern Colorado

Fisher’s Peak is a 77.5 km2 property southeast of Trinidad, Colorado that is planned to become Colorado’s newest state park. The area has experienced limited anthropogenic disturbance and is home to an abundance of unique habitats and species. A rapid, approximately 900 m change in elevation over the extent of the area nurtures a variety of plants and animals, including the endangered New Mexico meadow jumping mouse. In 2019, the State of Colorado obtained Fisher’s Peak with plans to make it Colorado’s second largest state park. A diverse group of collaborators, including the Colorado State Forest Service and The Nature Conservancy, worked closely to design the state park to maximize recreation opportunity while conserving the property’s rich habitats and biodiversity. The Fisher’s Peak Ecological Forecasting Team utilized Light Detection and Ranging (LiDAR) surveys, in situ forest inventory data, and Earth observations from Landsat 8 Operational Land Imager (OLI), Sentinel-1 C-band Synthetic Aperture Radar (C-SAR), Sentinel-2 Multispectral Instrument (MSI), Advanced Land Observing Satellite 2 (ALOS-2) Phased Array type L-band Synthetic Aperture Radar (PALSAR-2) and the Shuttle Radar Topography Mission (SRTM) to quantify and map biomass over the extent of the study area. The results from modeling biomass had an out-of-bag root mean square error of 55 Mg/ha and an R2 of 12. The resulting map indicates areas where carbon storage on the property is high, informing decision-making processes for future park development. While more in situ training data may improve modeling capacity for biomass in the Fisher’s Peak area, this work represents a feasible attempt to better understand biomass distribution using earth observation

Lauren Lad↗

Sensitivity of CONUS Summer Rainfall to the Selection of Cumulus Parameterization Schemes in NU-WRF Seasonal Simulations

This study investigates the sensitivity of daily rainfall rates in regional seasonal simulations over the contiguous United States (CONUS) to different cumulus parameterization schemes. Daily rainfall fields were simulated at 24-km resolution using the NASA-Unified Weather Research and Forecasting (NU-WRF) Model for June-August 2000. Four cumulus parameterization schemes and two options for shallow cumulus components in a specific scheme were tested. The spread in the domain-mean rainfall rates across the parameterization schemes was generally consistent between the entire CONUS and most subregions. The selection of the shallow cumulus component in a specific scheme had more impact than that of the four cumulus parameterization schemes. Regional variability in the performance of each scheme was assessed by calculating optimally weighted ensembles that minimize full root-mean-square errors against reference datasets. The spatial pattern of the seasonally averaged rainfall was insensitive to the selection of cumulus parameterization over mountainous regions because of the topographical pattern constraint, so that the simulation errors were mostly attributed to the overall bias there. In contrast, the spatial patterns over the Great Plains regions as well as the temporal variation over most parts of the CONUS were relatively sensitive to cumulus parameterization selection. Overall, adopting a single simulation result was preferable to generating a better ensemble for the seasonally averaged daily rainfall simulation, as long as their overall biases had the same positive or negative sign. However, an ensemble of multiple simulation results was more effective in reducing errors in the case of also considering temporal variation.

Iguchi, Takamichi↗

Estimating Alpine Snow Depth by Combining Multi-Frequency Passive Radiance Observations with Ensemble Snowpack Modeling

This paper presents a physically-based snow depth retrieval algorithm adapted for deep mountainous snowpack and airborne multifrequency (10.7, 18.7, 37.0 and 89.0 GHz) passive microwave (PM) radiance observations from a single flight. The algorithm employs a single forecast-analysis cycle of a traditional sequential assimilation scheme. It uses an ensemble of multi-layer snowpack model runs to resolve snow microstructure and melt-refreeze crusts, andmicrowave radiative transfer models to relate snow properties to microwave measurements.Snow depth was retrieved at a 120 m spatial resolution over three 1 km2 Intensive Study Area (ISA) within the Rabbit Ears Meso-Cell Study Area (MSA) from the NASA Cold Land Processes Experiment (CLPX) in Colorado (United States) for one date in February 2003. When evaluated against in situ observations, root mean square error (RMSE) of the snow depth from the assimilation was 13.3 cm for areas with low (<5 %) forest cover, which was a reduction of 48 % in the RMSE compared with the modeled snow depth when the PM observations were not assimilated; indicating a ~5% relative error of the posterior snow depth with respect to the average snow depth (200 cm) measured at these pixels. For pixels with forest cover ranging from 5-15% and 15-30%, results were improved (R2 increased from 0.65 to 0.71 and from 0 to 0.38, respectively) by introducing a forest radiative transfer model during the assimilation.

Airborne Multifrequency↗

Version 8 of the SMAP Level-4 Soil Moisture Data Assimilation Product

The NASA Soil Moisture Active Passive (SMAP) mission Level-4 Soil Moisture (L4_SM) product provides global, 9-km resolution, 3-hourly surface (0-5 cm) and root-zone (0-100 cm) soil moisture from April 2015 to present with a mean latency of 2.5 days from the time of observation. The product is based on the assimilation of SMAP L-band (1.4 GHz) brightness temperature (Tb) observations into the NASA Catchment land surface model as the model is driven with observations-based precipitation forcing. This presentation discusses the improvements in the forthcoming Version 8 of L4_SM, including updates in the precipitation forcing, the Catchment model parameters, and the L-band microwave radiative transfer model (mwRTM). The precipitation observations used in L4_SM Version 8 outside of North America and the high latitudes are from the latest (Version 7) NASA Integrated Multi-satellitE Retrievals for the Global Precipitation Measurement mission (IMERG) products. Moreover, occasionally excessive precipitation rates in earlier versions of L4_SM along certain longitudes in North America were eliminated by a bug fix in the precipitation corrections algorithm. The Catchment model in L4_SM Version 8 uses climatological snow albedo values based on observations from the Moderate Resolution Imaging Spectroradiometer, replacing the look-up table parameterization of earlier versions. Additionally, corrected soil parameters were implemented for a small region in Argentina that had erroneously been classified as peat because of an error in the ancillary soil data. Finally, the mwRTM in L4_SM Version 8 uses the Mironov soil mixing approach and updated values of the L-band scattering albedo, soil roughness, and vegetation opacity climatology obtained from the latest (Release 19) SMAP Level-2 dual-channel soil moisture retrieval product. During the development of L4_SM Version 8, the change in the mwRTM parameterization resulted in a reduced unbiased RMSE of surface soil moisture when verified against in situ measurements. It also reduced the standard deviation of the Tb observation-minus-forecast residuals by ~0.15 K, highlighting the importance of the mwRTM for successful data assimilation. The bug fixes in the precipitation corrections algorithm and the Catchment model soil parameters in the Argentina region result in locally large improvements of the simulated land surface states. In summary, the ongoing refinements of the L4_SM product continue to improve its science quality and performance for global soil moisture monitoring

Rolf Reichle↗

Coupled Stochastic Time-Inverted Lagrangian Transport/Weather Forecast and Research/Vegetation Photosynthesis and Respiration Model: Simulations of Tower-Based and Airborne CO2 Measurements - Part II

This paper evaluates simulations of atmospheric CO2 measured in 2004 at continental surface and airborne receptors, intended to test the capability to use data with high temporal and spatial resolution for analyses of carbon sources and sinks at regional and continental scales. The simulations were performed using the Stochastic Time-Inverted Lagrangian Transport (STILT) model driven by the Weather Forecast and Research (WRF) model, and linked to surface fluxes from the satellite-driven Vegetation Photosynthesis and Respiration Model (VPRM). The simulations provide detailed representations of hourly CO2 tower data and reproduce the shapes of airborne vertical profiles with high fidelity. WRF meteorology gives superior model performance compared with standard meteorological products, and the impact of including WRF convective mass fluxes in the STILT trajectory calculations is significant in individual cases. Important biases in the simulation are associated with the nighttime CO2 build-up and subsequent morning transition to convective conditions, and with errors in the advected lateral boundary condition. Comparison of STILT simulations driven by the WRF model against those driven by the Brazilian variant of the Regional Atmospheric Modeling System (BRAMS) shows that model-to-model differences are smaller than between an individual transport model and observations, pointing to systematic errors in the simulated transport. Future developments in the WRF model s data assimilation capabilities, basic research into the fundamental aspects of trajectory calculations, and intercomparison studies involving other transport models, are possible venues for reducing these errors. Overall, the STILT/WRF/VPRM offers a powerful tool for continental and regional scale carbon flux estimates.

Eluszkiewicz, Janusz↗

Assessing the Impact of SMAP Soil Moisture Data Assimilation on the Simulation and Prediction of Tropical Cyclone Idai

The role of soil moisture in the evolution of tropical cyclones (TCs) approaching land has long been recognized. Dry land surface conditions can lead to faster dissipation of a TC over land (often associated with changes in precipitation structure), whereas very wet conditions can help sustain or even re-intensify a TC. The ability to forecast post-landfall TC evolution and mitigate the associated socio-economic impact thus hinges on accurate knowledge of land surface conditions prior to landfall. The NASA Soil Moisture Active Passive (SMAP) mission provides accurate observations of soil moisture globally and at high revisit times of 2-3 days. It has been shown that the assimilation of SMAPbrightness temperatures (Tb) significantly improves modeled land surface states and thus has the potential to constrain land surface initial conditions in TC forecasts. In this presentation, we investigate this potential through an extensive set of Observing System Experiments that systematically assess the impact of assimilating SMAP Tbs on TC forecast skill in the Goddard Earth Observing System. Focusing on the case of TC Idai, we show that the assimilation of SMAP generally leads to a drier soil moisture analysis for the land surface underneath much of the storm’s circulation. An exception is the direct TC track, where SMAP assimilation increases soil moisture. These changes are reflected in the surface fluxes, which propagate the land surface state changes to the atmosphere. This results a TC that is overall weaker with lower winds, but with a better-defined eye and a slightly smaller track error than a simulation without SMAP DA. The changes in soil moisture resulting from the assimilation of SMAP also impact the total precipitation amounts as well as the precipitation structure associated with TC Idai.

SMAP↗

Advanced Observation Operators for GPS Radio Occultation: Validation of the 2D Ray Tracing Approach with CHAMP and SAC-C bending angle and refractivity data - Part 1

Global Positioning System (GPS) Radio Occultations (RO) bending angles and refractivity data characterize mostly the vertical structure of the Earth's atmosphere. We answer the question whether proper simulation of GPS RO data for data assimilation can be obtained with one-dimensional vertical operators, or if accounting also for horizontal atmospheric structures via ray-tracing makes a positive difference when compared with real data. We present a detailed implementation of a geometrical optics multi-plane two-dimensional (2D) ray-tracing as an observation operator to simulate GPS RO bending angles and refractivities within the Finite Volume Data Assimilation System (FVDAS). Comparisons of the outputs of that 2D observation operator with those of simpler ID observation operators are used to generate estimates of errors induced by neglecting tangent point drift (TPD) and horizontal gradients (HG). These error estimates are then confronted with errors estimates derived using 6335 real CHAMP and SAC-C occultations. The agreement for TPD-induced (HG-induced) errors is remarkably positive at altitudes 10-30 km (below 10 km). Comparisons in bending angles O - B STD of the outputs of the multi-plane 2D ray-tracer with those of a vertical Abel transform show reductions of about 8% of the usual O - B bending angle STD due to TPD in the stratosphere (3% due to HG, in the troposphere only). In terms of refractivity, the O - B STD reductions are about 1520% for TPD and 3-5% for HG in the same regions. These reductions are obtained using either 6-hour forecasts or analyses as backgrounds, and using Geometrical Optics (GO) or Canonical Transform (CT) data.

Poli, P.↗

Application of parametric weight and cost estimating relationships to future transport aircraft

A model comprised of system level weight and cost estimating relationships for transport aircraft is presented. In order to determine the production cost of future aircraft its weight is first estimated based on performance parameters, and then the cost is estimated as a function of weight. For initial evaluation CERs were applied to actual system weights of six aircraft (3 military and 3 commercial) with mean empty weights ranging from 30,000 to 300,000 lb. The resulting cost estimates were compared with actual costs. The average absolute error was only 4.3%. Then the model was applied to five aircraft still in the design phase (Boeing 757, 767 and 777, and BAC HS146-100 and HS146-200). While the estimates for the 757 and 767 are within 2 to 3 percent of their assumed break-even costs, it is recognized that these are very sensitive to the validity of the estimated weights, inflation factor, the amount assumed for nonrecurring costs, etc., and it is suggested that the model may be used in conjunction with other information such as RDT&E cost estimates and market forecasts. The model will help NASA evaluate new technologies and production costs of future aircraft.

Beltramo, M. N.↗

Satellite sensing capabilities for surface temperature and meteorological parameters over the ocean

Remote meteorological measurements over the oceans have been used extensively for weather forecasting, ship safety, and global-scale studies of climate and sea conditions. A review is conducted regarding the satellite-sensing capabilities for sea-surface temperature and tropospheric marine meteorological mesurements. The tropospheric measurements are concerned with vertical temperature profiles, integrated water vapor and liquid water, rain, aerosols and surface winds. Attention is also given to several environmental-monitoring applications where satellite sensed synoptic data might be combined with ground and airborne chemical measurements to provide improved assessments of regional and global ennvironmental quality and climatology. Special emphasis is given to current satellite measurement accuracies (errors) and to programs involving U.S. technology. Aspects of sea-surface temperature are treated in detail in connection with extensive studies concentrated in this area.

Darnell, W. L.↗

Assessing the Usefulness of AIRS Radiance Observations in a 4D-Var Assimilation Scheme Using the Penn State/NCAR Mesoscale Model Version 5 (MM5) and a Stand Alone Radiative Transfer Algorithm (SARTA)

Short-term precipitation forecast skill remains poor despite further advances in numerical modeling and data analysis. Utilizing indirect observations, such as the radiance observations obtained from the Atmospheric Infrared Sounder (AIRS), in a 4D-Var assimilation scheme may lead to a rapid improvement in the skill of 12-h QPF. In order to investigate the usefulness of AIRS observations, the adjoint code of a radiative transfer model (SARTA) is developed and linked to MM5 4D-Var system. A statistical comparison of the AIRS radiances with those calculated using the linked MMS-SARTA model is carried out to assess errors in the radiative transfer model, MM5 and AIRS data, and to develop a quality control scheme. Outliers and any bias that may exist in radiative transfer model, MM5 and AIRS data are identified. An adjoint sensitivity study, utilizing the SARTA adjoint model, is then done to identify the relative sensitivities of each spectral channel to such parameters as atmospheric temperature, pressure, and mixing ratio. Another sensitivity study is done using the MM5 alone to determine the sensitivity of the 12-h QPF to model parameters such as atmospheric temperature, pressure, and mixing ratio under various weather regimes (i.e. summertime convection, winter snowfall). Results from these two types of adjoint sensitivity studies are used to determine which AIRS channels contribute most to 12-h QPF under differing weather regimes. This information is useful when conducting 4D-Var experiments using a subset of the 2378 AIRS spectral channels. Finally, 4D-Var experiments are done to analyze the impact of the AIRS observations on the quality of the 12-h QPF results. This project falls under the data assimilation general session topic as specified by the conference organizers.

Carrier, M.↗

Using the Coronal Evolution to Successfully Forward Model CMEs' In Situ Magnetic Profiles

Predicting the effects of a coronal mass ejection (CME) impact requires knowing if impact will occur, which part of the CME impacts, and its magnetic properties. We explore the relation between CME deflections and rotations, which change the position and orientation of a CME, and the resulting magnetic profiles at 1 AU. For 45 STEREO-era, Earth-impacting CMEs, we determine the solar source of each CME, reconstruct its coronal position and orientation, and perform a ForeCAT (Forecasting a CME's Altered Trajectory) simulation of the coronal deflection and rotation. From the reconstructed and modeled CME deflections and rotations, we determine the solar cycle variation and correlations with CME properties. We assume no evolution between the outer corona and 1 AU and use the ForeCAT results to drive the ForeCAT In situ Data Observer (FIDO) in situ magnetic field model, allowing for comparisons with ACE and Wind observations. We do not attempt to reproduce the arrival time. On average FIDO reproduces the in situ magnetic field for each vector component with an error equivalent to 35 percent of the average total magnetic field strength when the total modeled magnetic field is scaled to match the average observed value. Random walk best fits distinguish between ForeCAT's ability to determine FIDO's input parameters and the limitations of the simple flux rope model. These best fits reduce the average error to 30 percent.The FIDO results are sensitive to changes of order a degree in the CME latitude, longitude, and tilt, suggesting that accurate space weather predictions require accurate measurements of a CME's position and orientation.

Kay, C.↗

Remote Sensing Based Yield Monitoring: Application to Winter Wheat in United States and Ukraine

Accurate and timely crop yield forecasts are critical for making informed agricultural policies and investments, as well as increasing market efficiency and stability. Earth observation data from space can contribute to agricultural monitoring, including crop yield assessment and forecasting. In this study, we present a new crop yield model based on the Difference Vegetation Index (DVI) extracted from Moderate Resolution Imaging Spectroradiometer (MODIS) data at 1 km resolution and the un-mixing of DVI at coarse resolution to a pure wheat signal (100% of wheat within the pixel). The model was applied to estimate the national and subnational winter wheat yield in the United States and Ukraine from 2001 to 2017. The model at the subnational level shows very good performance for both countries with a coefficient of determination higher than 0.7 and a root mean square error (RMSE) of lower than 0.6 t/ha (15–18%). At the national level for the United States (US) and Ukraine the model provides a strong coefficient of determination of 0.81 and 0.86, respectively, which demonstrates good performance at this scale. The model was also able to capture low winter wheat yields during years with extreme weather events, for example 2002 in US and 2003 in Ukraine. The RMSE of the model for the US at the national scale is 0.11 t/ha (3.7%) while for Ukraine it is 0.27 t/ha (8.4%).

remote sensing↗

Remote Sensing Based Yield Monitoring: Application to Winter Wheat in United States and Ukraine

Accurate and timely crop yield forecasts are critical for making informed agricultural policies and investments, as well as increasing market efficiency and stability. Earth observation data from space can contribute to agricultural monitoring, including crop yield assessment and forecasting. In this study, we present a new crop yield model based on the Difference Vegetation Index (DVI) extracted from Moderate Resolution Imaging Spectroradiometer (MODIS) data at 1 km resolution and the un-mixing of DVI at coarse resolution to a pure wheat signal (100 percent of wheat within the pixel). The model was applied to estimate the national and subnational winter wheat yield in the United States and Ukraine from 2001 to 2017. The model at the subnational level shows very good performance for both countries with a coefficient of determination higher than 0.7 and a root mean square error (RMSE) of lower than 0.6 t/ha (tonnes per hectare) (15-18 percent). At the national level for the United States (US) and Ukraine the model provides a strong coefficient of determination of 0.81 and 0.86, respectively, which demonstrates good performance at this scale. The model was also able to capture low winter wheat yields during years with extreme weather events, for example 2002 in US and 2003 in Ukraine. The RMSE of the model for the US at the national scale is 0.11 t/ha (3.7 percent) while for Ukraine it is 0.27 t/ha (8.4 percent).

DVI↗