Accurate crop yield assessments using satellite remote sensing-based methods are of interest for regional monitoring and the design of policies that promote agricultural resiliency and food security. However, the application of current vegetation productivity algorithms derived from global satellite observations is generally too coarse to capture cropland heterogeneity. The fusion of data from different sensors can provide enhanced information and overcome many of the limitations of individual sensors. In thitables study, we estimate annual crop yields for seven important crop types across Montana in the continental USA from 2008–2015, including alfalfa, barley, maize, peas, durum wheat, spring wheat and winter wheat. We used a satellite data-driven light use efficiency (LUE) model to estimate gross primary productivity (GPP) over croplands at 30-m spatial resolution and eight-day time steps using a fused NDVI dataset constructed by blending Landsat (5 or 7) and Terra MODIS reflectance data. The fused 30-m NDVI record showed good consistency with the original Landsat and MODIS data, but provides better spatiotemporal delineations of cropland vegetation growth. Crop yields were estimated at 30-m resolution as the product of estimated GPP accumulated over the growing season and a crop-specific harvest index (HI_GPP). The resulting GPP estimates capture characteristic cropland productivity patterns and seasonal variations, while the estimated annual crop production results correspond favorably with reported county-level crop production data (r = 0.96, relative RMSE = 37.0%, p < 0.05) from the U.S. Department of Agriculture (USDA). The performance of estimated crop yields at a finer (field) scale was generally lower, but still meaningful (r = 0.42, relative RMSE = 50.8%, p < 0.05). Our methods and results are suitable for operational applications of crop yield monitoring at regional scales, suggesting the potential of using global satellite observations to improve agricultural management, policy decisions and regional/global food security.

Satellite remote sensing observations show an increased greenness trend over land in recent decades. While greenness observations can indicate increased productivity, estimation of total annual productivity is highly dependent on vegetation response to climate and environmental conditions. Models have been struggling to determine how much carbon is taken up by plants as a result of increased atmospheric CO₂ fertilization. Current remote sensing light use efficiency (LUE) models contain considerable uncertainty due to the lack of spatial and temporal variability in maximum LUE parameter and climate sensitivity defined for global plant functional types (PFTs). We used the optimum LUE (LUE_opt) previously derived from the global FLUXNET network to improve estimation of global gross primary productivity (GPP) for the period 1982–2016. Our results indicate increasing GPP in northern latitudes owing to reduced cold temperature constraints on plant growth, thereby suggesting increasing negative carbon‐climate feedback in high latitudes. In the tropics, by contrast, our results indicate an emerging positive climate feedback, mainly due to increasing atmospheric vapor pressure deficit (VPD). Further pervasive VPD increase is likely to continue to reduce global GPP and amplify carbon emissions.

Studies seeking to understand the impacts of climate variability and change on the hydrology of a region need to take into account the dynamics of vegetation and its interaction with the hydrologic and energy cycles. Yet, most of the hydrologic models used for these kinds of studies assume that vegetation is static. This paper presents a dynamic, spatially explicit model that couples a vertical energy balance scheme (surface and canopy layer) to a hydrologic model and a forest growth component to capture the dynamic interactions between energy, vegetation, and hydrology at hourly to daily time scales. The model is designed to be forced with outputs from regional climate models. Lateral water transfers are simulated using a 1D kinematic wave model. Infiltration is simulated using the Green and Ampt approximation to Richard's equation. The dynamics of soil moisture and energy drives carbon assimilation and forest growth, which in turn affect the distribution of energy and water through leaf dynamics by altering light interception, shading, and enhanced transpiration. The model is demonstrated in two case studies simulating energy, water, and vegetation dynamics at two different spatial and temporal scales.

High frequency and spatially explicit irrigated land maps are important for understanding the patterns and impacts of consumptive water use by agriculture. We built annual, 30 m resolution irrigation maps using Google Earth Engine for the years 1986–2018 for 11 western states within the conterminous U.S. Our map classifies lands into four classes: irrigated agriculture, dryland agriculture, uncultivated land, and wetlands. We built an extensive geospatial database of land cover from each class, including over 50,000 human-verified irrigated fields, 38,000 dryland fields, and over 500,000 km ${}^2$ of uncultivated lands. We used 60,000 point samples from 28 years to extract Landsat satellite imagery, as well as climate, meteorology, and terrain data to train a Random Forest classifier. Using a spatially independent validation dataset of 40,000 points, we found our classifier has an overall binary classification (irrigated vs. unirrigated) accuracy of 97.8%, and a four-class overall accuracy of 90.8%. We compared our results to Census of Agriculture irrigation estimates over the seven years of available data and found good overall agreement between the 2832 county-level estimates (r ${}^2$ = 0.90), and high agreement when estimates are aggregated to the state level (r ${}^2$ = 0.94). We analyzed trends over the 33-year study period, finding an increase of 15% (15,000 km ${}^2$ ) in irrigated area in our study region. We found notable decreases in irrigated area in developing urban areas and in the southern Central Valley of California and increases in the plains of eastern Colorado, the Columbia River Basin, the Snake River Plain, and northern California.

In this paper a high‐resolution linked hydroeconomic model is demonstrated for drought conditions in a Brazilian river basin. The economic model of agriculture includes 13 decision variables that can be optimized to maximize farmers' yearly net revenues. The economic model uses a multi‐input multioutput nonlinear constant elasticity of substitution (CES) production function simulating agricultural production. The hydrologic component is a detailed physics‐based three‐dimensional hydrodynamic model that simulates changes in the hydrologic system derived from agricultural activity while in turn providing biophysical constraints to the economic system. The linked models capture the effects of the interactions between the hydrologic and the economic systems at high spatial and temporal resolutions, ensuring that the model converges to an optimal economic scenario that takes into account the spatial and temporal distribution of the water resources. The operation and usefulness of the models are demonstrated in a rural catchment area of about 10 km² within the São Francisco River Basin in Brazil. Two droughts of increasing intensity are simulated to investigate how farmers behave under rain shortfalls of different severity. The results show that farmers react to rainfall shortages to minimize their effects on farm profits, and that the impact on farmers depends, among other things, on their location in the watershed and on their access to groundwater.