Vegetation has been effectively monitored using remote sensing time-series vegetation index (VI) data for several decades. Drought monitoring has been a common application with algorithms tuned to capturing anomalous temporal and spatial vegetation patterns. Drought stress models, such as the Vegetation Drought Response Index (VegDRI), often use VIs like the Normalized Difference Vegetation Index (NDVI). The EROS expedited Moderate Resolution Imaging Spectroradiometer (eMODIS)-based, 7-day NDVI composites are integral to the VegDRI. As MODIS satellite platforms (Terra and Aqua) approach mission end, the Visible Infrared Imaging Radiometer Suite (VIIRS) presents an alternate NDVI source, with daily collection, similar band passes, and moderate spatial resolution. This study provides a statistical comparison between EROS expedited VIIRS (eVIIRS) 375-m and eMODIS 250-m and tests the suitability of replacing MODIS NDVI with VIIRS NDVI for drought monitoring and vegetation anomaly detection. For continuity with MODIS NDVI, we calculated a geometric mean regression adjustment algorithm using 375-m resolution for an eMODIS-like NDVI (eVIIRS’) eVIIRS’ = 0.9887 × eVIIRS − 0.0398. The resulting statistical comparisons (eVIIRS’ vs. eMODIS NDVI) showed correlations consistently greater than 0.84 throughout the three years studied. The eVIIRS’ VegDRI results characterized similar drought patterns and hotspots to the eMODIS-based VegDRI, with near zero bias.
The United States has a geographically mature and stable land use and land cover system including land used as irrigated cropland; however, changes in irrigation land use frequently occur related to various drivers. We applied a consistent methodology at a 250 m spatial resolution across the lower 48 states to map and estimate irrigation dynamics for four map eras (2002, 2007, 2012, and 2017) and over four 5-year mapping intervals. The resulting geospatial maps (called the Moderate Resolution Imaging Spectroradiometer (MODIS) Irrigated Agriculture Dataset or MIrAD-US) involved inputs from county-level irrigated statistics from the U.S. Department of Agriculture, National Agricultural Statistics Service, agricultural land cover from the U.S. Geological Survey National Land Cover Database, and an annual peak vegetation index derived from expedited MODIS satellite imagery. This study investigated regional and periodic patterns in the amount of change in irrigated agriculture and linked gains and losses to proximal causes and consequences. While there was a 7% overall increase in irrigated area from 2002 to 2017, we found surprising variability by region and by 5-year map interval. Irrigation land use dynamics affect the environment, water use, and crop yields. Regionally, we found that the watersheds with the largest irrigation gains (based on percent of area) included the Missouri, Upper Mississippi, and Lower Mississippi watersheds. Conversely, the California and the Texas–Gulf watersheds experienced fairly consistent irrigation losses during these mapping intervals. Various drivers for irrigation dynamics included regional climate fluctuations and drought events, demand for certain crops, government land or water policies, and economic incentives like crop pricing and land values. The MIrAD-US (Version 4) was assessed for accuracy using a variety of existing regionally based reference data. Accuracy ranged between 70% and 95%, depending on the region.
Crop cover maps have become widely used in a range of research applications. Multiple crop cover maps have been developed to suite particular research interests. The National Agricultural Statistics Service (NASS) Cropland Data Layers (CDL) are a series of commonly used crop cover maps for the conterminous United States (CONUS) that span from 2008 to 2013. In this investigation, we sought to contribute to the availability of consistent CONUS crop cover maps by extending temporal coverage of the NASS CDL archive back eight additional years to 2000 by creating annual NASS CDL-like crop cover maps derived from a classification tree model algorithm. We used over 11 million records to train a classification tree algorithm and develop a crop classification model (CCM). The model was used to create crop cover maps for the CONUS for years 2000–2013 at 250 m spatial resolution. The CCM and the maps for years 2008–2013 were assessed for accuracy relative to resampled NASS CDLs. The CCM performed well against a withheld test data set with a model prediction accuracy of over 90%. The assessment of the crop cover maps indicated that the model performed well spatially, placing crop cover pixels within their known domains; however, the model did show a bias towards the ‘Other’ crop cover class, which caused frequent misclassifications of pixels around the periphery of large crop cover patch clusters and of pixels that form small, sparsely dispersed crop cover patches.
Integrating spatially explicit biogeophysical and remotely sensed data into regression-tree models enables the spatial extrapolation of training data over large geographic spaces, enhancing a more complete understanding of broad-scale ecosystem processes. This data release presents maps of estimates of annual gross primary production (GPP) and annual ecosystem respiration (RE) that were derived from weekly summaries of gross photosynthesis (Pg) and ecosytem respiration (Re). To conduct this study we used carbon data from flux towers that are scattered strategically across the conterminous United States (CONUS). We also calculate and present a map of average annual net ecosystem production (NEP). We present and analyze carbon flux dynamics in the CONUS to better understand how drought affects GPP, RE, and NEP. Our study experienced correlation coefficients (r) greater than or equal to 0.94 between training and estimated data for both GPP and RE. We conclude that this modeling method effectively measures carbon dynamics through time and allows the exploration of impacts of meteorological anomalies and vegetation types on carbon dynamics.
Soil microbial communities control critical ecosystem processes such as decomposition, nutrient cycling, and soil organic matter formation. Continental scale patterns in the composition and functioning of microbial communities are related to climatic, biotic, and edaphic factors such as temperature and precipitation, plant community composition, and soil carbon, nitrogen, and pH. Although these relationships have been well explored individually, the examination of the factors that may act directly on microbial communities vs. those that may act indirectly through other ecosystem properties has not been well developed. To further such understanding, we utilized structural equation modeling (SEM) to evaluate a set of hypotheses about the direct and indirect effects of climatic, biotic, and edaphic variables on microbial communities across the continental United States. The primary goals of this work were to test our current understanding of the interactions among climate, soils, and plants in affecting microbial community composition, and to examine whether variation in the composition of the microbial community affects potential rates of soil enzymatic activities. A model of interacting factors created through SEM shows several expected patterns. Distal factors such as climate had indirect effects on microbial communities by influencing plant productivity, soil mineralogy, and soil pH, but factors related to soil organic matter chemistry had the most direct influence on community composition. We observed that both plant productivity and soil mineral composition were important indirect influences on community composition at the continental scale, both interacting to affect organic matter content and microbial biomass and ultimately community composition. Although soil hydrolytic enzymes were related to the moisture regime and soil carbon, oxidative enzymes were also affected by community composition, reflected in the abundance of soil fungi. These results highlight that soil microbial communities can be modeled within the context of multiple interacting ecosystem properties acting both directly and indirectly on their composition and function, and this provides a rich and informative context with which to examine communities. This work also highlights that variation in climate, microbial biomass, and microbial community composition can affect maximum rates of soil enzyme activities, potentially influencing rates of decomposition and nutrient mineralization in soils.
The shortgrass steppe (SGS) occupies the southwestern part of the Great Plains. Half of the land is cultivated, but significant areas remain under natural vegetation. Despite previous studies of the SGS carbon cycle, not all aspects have been completely addressed, including gross productivity, ecosystem respiration, and ecophysiological parameters. Our analysis of 1998−2007 flux tower measurements at five Bowen ratio–energy balance (BREB) and three eddy covariance (EC) sites characterized seasonal and interannual variability of gross photosynthesis and ecosystem respiration. Identification of the nonrectangular hyperbolic equation for the diurnal CO2 exchange, with vapor pressure deficit (VPD) limitation and exponential temperature response, quantified quantum yield α, photosynthetic capacity Amax, and respiration rate rd with variation ranges (19 < α < 51 mmol mol−1, 0.48 < Amax < 2.1 mg CO2 m−2 s−1, 0.15 < rd < 0.49 mg CO2 m−2 s−1). Gross photosynthesis varied from 1 100 to 2 700 g CO2 m−2 yr−1, respiration from 900 to 3,000 g CO2 m−2 yr−1, and net ecosystem production from −900 to +700 g CO2 m−2 yr−1, indicating that SGS may switch from a sink to a source depending on weather. Comparison of the 2004−2006 measurements at two BREB and two parallel EC flux towers located at comparable SGS sites showed moderately higher photosynthesis, lower respiration, and higher net production at the BREB than EC sites. However, the difference was not related only to methodologies, as the normalized difference vegetation index at the BREB sites was higher than at the EC sites. Overall magnitudes and seasonal patterns at the BREB and the EC sites during the 3-yr period were similar, with trajectories within the ± 1.5 standard deviation around the mean of the four sites and mostly reflecting the effects of meteorology.
Integrating spatially explicit biogeophysical and remotely sensed data into regression-tree models enables the spatial extrapolation of training data over large geographic spaces, allowing a better understanding of broad-scale ecosystem processes. The current study presents annual gross primary production GPP and annual ecosystem respiration RE for 2000–2013 in several short-statured vegetation types using carbon flux data from towers that are located strategically across the conterminous United States CONUS. We calculate carbon fluxes annual net ecosystem production [NEP] for each year in our study period, which includes 2012 when drought and higher-than-normal temperatures influence vegetation productivity in large parts of the study area. We present and analyse carbon flux dynamics in the CONUS to better understand how drought affects GPP, RE, and NEP. Model accuracy metrics show strong correlation coefficients rr ≥ 94% between training and estimated data for both GPP and RE. Overall, average annual GPP, RE, and NEP are relatively constant throughout the study period except during 2012 when almost 60% less carbon is sequestered than normal. These results allow us to conclude that this modelling method effectively estimates carbon dynamics through time and allows the exploration of impacts of meteorological anomalies and vegetation types on carbon dynamics.
Regression tree models have been widely used for remote sensing-based ecosystem mapping. Improper use of the sample data (model training and testing data) may cause overfitting and underfitting effects in the model. The goal of this study is to develop an optimal sampling data usage strategy for any dataset and identify an appropriate number of rules in the regression tree model that will improve its accuracy and robustness. Landsat 8 data and Moderate-Resolution Imaging Spectroradiometer-scaled Normalized Difference Vegetation Index (NDVI) were used to develop regression tree models. A Python procedure was designed to generate random replications of model parameter options across a range of model development data sizes and rule number constraints. The mean absolute difference (MAD) between the predicted and actual NDVI (scaled NDVI, value from 0–200) and its variability across the different randomized replications were calculated to assess the accuracy and stability of the models. In our case study, a six-rule regression tree model developed from 80% of the sample data had the lowest MAD (MADtraining = 2.5 and MADtesting = 2.4), which was suggested as the optimal model. This study demonstrates how the training data and rule number selections impact model accuracy and provides important guidance for future remote-sensing-based ecosystem modeling.
This paper presents the methodology and results of two ecological-based net ecosystem production (NEP) regression tree models capable of up scaling measurements made at various flux tower sites throughout the U.S. Great Plains. Separate grassland and cropland NEP regression tree models were trained using various remote sensing data and other biogeophysical data, along with 15 flux towers contributing to the grassland model and 15 flux towers for the cropland model. The models yielded weekly mean daily grassland and cropland NEP maps of the U.S. Great Plains at 250 m resolution for 2000–2008. The grassland and cropland NEP maps were spatially summarized and statistically compared. The results of this study indicate that grassland and cropland ecosystems generally performed as weak net carbon (C) sinks, absorbing more C from the atmosphere than they released from 2000 to 2008. Grasslands demonstrated higher carbon sink potential (139 g C·m−2·year−1) than non-irrigated croplands. A closer look into the weekly time series reveals the C fluctuation through time and space for each land cover type.
Switchgrass is being evaluated as a potential feedstock source for cellulosic biofuels and is being cultivated in several regions of the United States. The recent availability of switchgrass land cover maps derived from the National Agricultural Statistics Service cropland data layer for the conterminous United States provides an opportunity to assess the environmental conditions of switchgrass over large areas and across different geographic locations. The main goal of this study is to develop a data-driven multiple regression switchgrass productivity model and identify the optimal climate and environment conditions for the highly productive switchgrass in the Great Plains (GP). Environmental and climate variables used in the study include elevation, soil organic carbon, available water capacity, climate, and seasonal weather. Satellite-derived growing season averaged Normalized Difference Vegetation Index (GSN) was used as a proxy for switchgrass productivity. Multiple regression analyses indicate that there are strong correlations between site environmental variables and switchgrass productivity (r = 0.95). Sufficient precipitation and suitable temperature during the growing season (i.e., not too hot or too cold) are favorable for switchgrass growth. Elevation and soil characteristics (e.g., soil available water capacity) are also an important factor impacting switchgrass productivity. An anticipated switchgrass biomass productivity map for the entire GP based on site environmental and climate conditions and switchgrass productivity model was generated. Highly productive switchgrass areas are mainly located in the eastern part of the GP. Results from this study can help land managers and biofuel plant investors better understand the general environmental and climate conditions influencing switchgrass growth and make optimal land use decisions regarding switchgrass development in the GP. (C) 2014 Elsevier Ltd. All rights reserved.
Monitoring systems benefit from high temporal frequency image data collected from the Moderate Resolution Imaging Spectroradiometer (MODIS) system. Because of near-daily global coverage, MODIS data are beneficial to applications that require timely information about vegetation condition related to drought, flooding, or fire danger. Rapid satellite data streams in operational applications have clear benefits for monitoring vegetation, especially when information can be delivered as fast as changing surface conditions. An expedited processing system called eMODIS operated by the U.S. Geological Survey provides rapid MODIS surface reflectance data to operational applications in less than 24 h offering tailored, consistently-processed information products that complement standard MODIS products. We assessed eMODIS quality and consistency by comparing to standard MODIS data. Only land data with known high quality were analyzed in a central U.S. study area. When compared to standard MODIS (MOD/MYD09Q1), the eMODIS Normalized Difference Vegetation Index (NDVI) maintained a strong, significant relationship to standard MODIS NDVI, whether from morning (Terra) or afternoon (Aqua) orbits. The Aqua eMODIS data were more prone to noise than the Terra data, likely due to differences in the internal cloud mask used in MOD/MYD09Q1 or compositing rules. Post-processing temporal smoothing decreased noise in eMODIS data.
Net CO 2 exchange data of legume crops at 17 flux tower sites in North America and three sites in Europe representing 29 site‐years of measurements were partitioned into gross photosynthesis and ecosystem respiration by using the nonrectangular hyperbolic light‐response function method. The analyses produced net CO 2 exchange data and new ecosystem‐scale ecophysiological parameter estimates for legume crops determined at diurnal and weekly time steps. Dynamics and annual totals of gross photosynthesis, ecosystem respiration, and net ecosystem production were calculated by gap filling with multivariate nonlinear regression. Comparison with the data from grain crops obtained with the same method demonstrated that CO 2 exchange rates and ecophysiological parameters of legumes were lower than those of maize ( Zea mays L.) but higher than for wheat ( Triticum aestivum L.) crops. Year‐round annual legume crops demonstrated a broad range of net ecosystem production, from sinks of 760 g CO 2 m −2 yr −1 to sources of –2100 g CO 2 m −2 yr −1 , with an average of –330 g CO 2 m −2 yr −1 , indicating overall moderate CO 2 –source activity related to a shorter period of photosynthetic uptake and metabolic costs of N 2 fixation. Perennial legumes (alfalfa, Medicago sativa L.) were strong sinks for atmospheric CO 2 , with an average net ecosystem production of 980 (range 550–1200) g CO 2 m −2 yr −1 .
The purpose of this study was to increase the spatial and temporal availability of crop classification data. In this study, nearly 16.2 million crop observation points were used in the training of the US Great Plains classification tree crop type model (CTM). Each observation point was further defined by weekly Normalized Difference Vegetation Index, annual climate, and a number of other biogeophysical environmental characteristics. This study accounted for the most prevalent crop types in the region, including, corn, soybeans, winter wheat, spring wheat, cotton, sorghum, and alfalfa. Annual CTM crop maps of the US Great Plains were created for 2000 to 2011 at a spatial resolution of 250 meters. The CTM achieved an 87 percent classification success rate on 1.8 million observation points that were withheld from model training. Product validation was performed on greater than 15,000 county records with a coefficient of determination of R-2 = 0.76.
Topobathymetric Digital Elevation Models (DEMs) are a merged rendering of both topography (land elevation) and bathymetry (water depth) that provides a seamless elevation product useful for inundation mapping, as well as for other earth science applications, such as the development of sediment-transport, sea-level rise, and storm-surge models. This 1/9-arc-second (approximately 3 meters) resolution model of Mobile Bay, Alabama was developed using multiple topographic and bathymetric datasets, collected on different dates. The topographic data were obtained primarily from the U.S. Geological Survey (USGS) National Elevation Dataset (NED) (http://ned.usgs.gov/) at 1/9-arc-second resolution; USGS Experimental Advanced Airborne Research Lidar (EAARL) data (2 meters) (http://pubs.usgs.gov/ds/400/); and topographic lidar data (2 meters) and Compact Hydrographic Airborne Rapid Total Survey (CHARTS) lidar data (2 meters) from the U.S. Army Corps of Engineers (USACE) (http://www.csc.noaa.gov/digitalcoast/data/coastallidar/). Bathymetry was derived from digital soundings obtained from the National Oceanic and Atmospheric Administration’s (NOAA) National Geophysical Data Center (NGDC) (http://www.ngdc.noaa.gov/mgg/geodas/geodas.html) and from water-penetrating lidar sources, such as EAARL and CHARTS. Mobile Bay is ecologically important as it is the fourth largest estuary in the United States. The Mobile and Tensaw Rivers drain into the bay at the northern end with the bay emptying into the Gulf of Mexico at the southern end. Dauphin Island (a barrier island) and the Fort Morgan Peninsula form the mouth of Mobile Bay. Mobile Bay is 31 miles (50 kilometers) long by a maximum width of 24 miles (39 kilometers) with a total area of 413 square miles (1,070 square kilometers). The vertical datum of the Mobile Bay topobathymetric model is the North American Vertical Datum of 1988 (NAVD 88). All the topographic datasets were originally referenced to NAVD 88 and no transformations were made to these input data. The NGDC hydrographic, multibeam, and trackline surveys were transformed from mean low water (MLW) or mean lower low water (MLLW) to NAVD 88 using VDatum (http://vdatum.noaa.gov). VDatum is a tool developed by the National Geodetic Survey (NGS) that performs transformations among tidal, ellipsoid-based, geoid-based, and orthometric datums using calibrated hydrodynamic models. The vertical accuracy of the input topographic data varied depending on the input source. Because the input elevation data were derived primarily from lidar, the vertical accuracy ranges from 6 to 20 centimeters in root mean square error (RMSE). he horizontal datum of the Mobile Bay topobathymetric model is the North American Datum of 1983 (NAD 83), geographic coordinates. All the topographic and bathymetric datasets were originally referenced to NAD 83, and no transformations were made to the input data. The bathymetric surveys were downloaded referenced to NAD 83 geographic, and therefore no horizontal transformations were required. The topbathymetric model of Mobile Bay and detailed metadata can be obtained from the USGS Web sites: http://nationalmap.gov/.
We analyzed net CO2 exchange data from 13 flux tower sites with 27 site-years of measurements over maize and wheat fields across midcontinent North America. A numerically robust "light-soil temperature-VPD"-based method was used to partition the data into photosynthetic assimilation and ecosystem respiration components. Year-round ecosystem-scale ecophysiological parameters of apparent quantum yield, photosynthetic capacity, convexity of the light response, respiration rate parameters, ecological light-use efficiency, and the curvature of the VPD-response of photosynthesis for maize and wheat crops were numerically identified and interpolated/extrapolated. This allowed us to gap-fill CO2 exchange components and calculate annual totals and budgets. VPD-limitation of photosynthesis was systematically observed in grain crops of the region (occurring from 20 to 120 days during the growing season, depending on site and year), determined by the VPD regime and the numerical value of the curvature parameter of the photosynthesis-VPD-response, sigma(VPD). In 78% of the 27 site-years of observations, annual gross photosynthesis in these crops significantly exceeded ecosystem respiration, resulting in a net ecosystem production of up to 2100 g CO2 m(-2) year(-1). The measurement-based photosynthesis, respiration, and net ecosystem production data, as well as the estimates of the ecophysiological parameters, provide an empirical basis for parameterization and validation of mechanistic models of grain crop production in this economically and ecologically important region of North America. (C) 2012 Elsevier B.V. All rights reserved.