Future studies are needed to assess the accuracy of SBG for such a classification and further evaluate its suitability to detect crop stress independent of a crop species map. The SBG mission will allow for further exploration of crop temperatures as they change temporally and as they differ across different climates and landscapes. In addition to SBG, with the ability to track diurnal changes in plant water globally, ECOSTRESS is well positioned to lead scientific understanding of crop water stress as it changes spatially and temporally . This study provides rationale for accounting for species and soil type when analyzing thermal data and enhances the importance of pairing the thermal data with ground data or VSWIR imagery from Landsat or MODIS that can enhance understanding of surface characteristics. Further, with paired VSWIR imagery, it suggests the potential for extrapolating the methodology in this paper for global study of crop tree stress. Atmospheric water vapor is a critical element of climate, an indicator of land surface hydrologic processes, and a potent greenhouse gas . As such, analysis of vapor patterns at a fine spatial scale can inform climate and plant water use studies . Imaging spectrometers such as NASA’s Airborne Visible Infrared Imaging Spectrometer measure reflected radiance at fine spatial and spectral resolution, and in so doing provide measurements of column water vapor as well as a highly detailed reflected signal from the land surface below . These two spatially corresponding products allow us to uniquely observe surface processes and characteristics as they relate to the atmospheric patterns above them in ways not previously possible. As such,blackberries in containers this research proposes to leverage water vapor and reflectance imagery to observe and assess spatial patterns of water vapor in the Central Valley of California to evaluate the assets and limitations of this dataset for evaluation of agricultural water use.
Observation and evaluation of water vapor over agricultural fields in California’s Central Valley have substantial value for water resource management, irrigation assessments, and regional climate patterns. The Central Valley contains one of the world’s largest contiguous areas of high irrigation density with more than 3.6 million irrigated hectares of farmland that use over 80% of the state’s managed water supply . Worldwide, arid and semi-arid regions such as the Central Valley see upwards of ninety percent of precipitated water returned back to the atmosphere via evapotranspiration . In the Central Valley, where precipitation is low and managed water inputs are extreme, annual ET exceeds precipitation by about 60% . These extreme irrigation inputs, therefore, significantly modify the spatial and temporal distribution of hydrologic flows across the region by transforming liquid water resources into transpired atmospheric water vapor that can be transported and distributed as rainfall elsewhere . Further, as local atmospheric water vapor is intensified by ET, it is indicative of agricultural water inputs and crop functioning throughout the region. Imaging spectrometers such as AVIRIS quantify column water vapor using several water absorption features across the infrared portion of the electromagnetic spectrum . Atmospheric water vapor absorption features occur at 0.94, 1.14, 1.38 and 1.88 µm, and the relative depth of these features can be used to derive atmospheric water content . Hyperspectral imagery is uniquely suited to estimate water vapor because its high spectral resolution captures water absorption features that multi-spectral sensors, such as Landsat, are designed to avoid. In addition, the fine spatial resolution of the retrievals enable observation of water vapor patterns that are obscured in spatially coarser imagery, such as MODIS and GPS.
While water vapor imagery is produced as a byproduct of most visible to shortwave infrared reflectance retrievals, few analyses have been conducted with this rich dataset, leaving many questions as to the utility of this data unanswered. A notable exception is the work of Ogunjemiyo et al. who studied water vapor over poplar plantations in Washington State to assess the feasibility of using AVIRIS-retrieved column water vapor as a tool to study plant ET. That study proposed a conceptual model of water vapor and its relationship to the surface , hypothesizing that plants with higher rates of transpiration will produce more water vapor, which will advect downwind and accumulate to a level detectable in the imagery, and that crops with higher water use rates will have steeper water vapor slopes, modified by wind speeds. Ogunjemiyo et al. found that the patterns and magnitude of retrieved water vapor in their study areas were consistent with wind direction and reasonable transpiration rates for poplars with unlimited access to water, concluding that AVIRIS water vapor is sensitive to ET under certain boundary layer conditions. However, while transpiration is the dominant source of water vapor in the atmospheric boundary layer , directly relating water vapor to the agricultural field that produced it will depend on multiple factors including atmospheric turbulence, time of day, wind magnitude and direction, and atmospheric stability. Here, we build upon the work of Ogunjemiyo et al. to further evaluate if AVIRIS is sensitive to field scale ET by testing a series of specific hypotheses to investigate how AVIRIS estimates of water vapor vary with the surface properties and atmospheric conditions that might be expected to influence water vapor in a complex agricultural environment . Relating water vapor patterns to crop water use will be more challenging to isolate in the Central Valley of California, which includes a complex arrangement of agricultural fields that vary in crop types, land management, irrigation practices, field sizes, and vegetated cover fractions.
Relative to the large, single crop, well-watered field studied in Ogunjemiyo et al. , the heterogeneity in crop type and water availability, smaller field sizes, and the semi-arid climate in our study area add complexity to the interaction between the atmosphere and underlying physical and plant physiological characteristics. Table 4.2 summarizes additional mechanisms that might influence the interactions summarized in the simple conceptual model . This study will also conduct analyses of water vapor at three scales – pixel, field, and scene – using three AVIRIS-derived water vapor images collected in early June of three different years. Our hypothesis are designed to refine understanding of whether consistent relationships between imagery, atmospheric conditions, and surface properties are derivable from hyperspectral imagery in a complex agricultural landscape and at which scales. Results will identify opportunities and limitations of using water vapor imagery to study ET at the ground surface in this important agricultural region. The Surface Biology and Geology mission has been identified in the 2017 Decadal Survey as a designated program element prioritized for development as a means to enhance our ability to monitor ecosystems,blackberry container natural hazards, and land use over time . The proposed sensor will capture a large number of spectral bands in the visible and shortwave infrared at a 30 m resolution, as well as multiple bands in the thermal infrared at a 60 m resolution. In order to simulate the capabilities of SBG, the HyspIRI Airborne Campaign flew the AVIRIS and MODIS-ASTER Simulator instruments throughout California seasonally from 2013 to 2017 at an altitude of 20 km. AVIRIS is a 224 band imaging spectrometer that captures wavelengths from 350 nm – 2500 nm at approximately 10 nm increments . MASTER captures 8 bands in the TIR region from 4-12 μm . NASA’s Jet Propulsion Lab preprocessed the imagery and produced orthorectified, high spectral resolution radiance and atmospherically-corrected reflectance imagery from AVIRIS at an 18 m spatial resolution and a land surface temperature product from MASTER at a 36 m spatial resolution. AVIRIS imagery was resampled to 36 m by pixel aggregation for consistency with LST imagery. Three dates of imagery from HAC were analyzed: June 6, 2013, June 3, 2014 and June 2, 2015 collected at 18:25, 21:41, and 18:59 UTC respectively. Pixel-level water vapor estimates were calculated from the radiance imagery using Atmospheric COR rection Now 6.0 atmospheric correction software . ACORN 6 models atmospheric gas absorptions and scattering effects using nonlinear least-squares spectral fitting with look-up tables of water column densities generated from MODTRAN 4 radiative transfer runs .ACORN 6 was run in mode 1.5 for atmospheric correction of hyperspectral data.
Chosen parameters include the use of both 940 nm and 1140 nm to derive water vapor, a midlatitude summer model, an average surface elevation of 100 meters, automated estimate of visibility and artifact suppression. The studied agricultural landscape was characterized using both the reflectance imagery and a geographic information system layer of field data. Multiple Endmember Spectra Mixture Analysis was run on the AVIRIS A reflectance imagery to estimate pixel-level fractions of green vegetation , soil, and non-photosynthetic vegetation . MESMA is a spectral mixture technique that models the spectral signature of each pixel using a linear combination of spectra from the classes within the scene. The chosen endmembers for each class are allowed to vary on a per-pixel basis to capture the diversity of spectra contained within GV, NPV, and soil. The resultant product is a sub-pixel map of vegetation fraction, both senesced and green, and soil throughout the study area. MESMA was run on each of the three image dates using one spectral library of 40 image-selected endmembers that were collected from all three dates. In addition to cover fractions, we obtained field size and crop species data. This information was obtained from GIS crop data layers provided by Tulare, Kings, Kern, and Fresno Counties. Field boundaries and crop information is gathered as part of a California’s required registration and permitting of agricultural fields that use pesticides. Field sizes ranged from <100 m2 to 2.7 km2 . Using the MESMA results and the crop map, which delineates field boundaries, a mean estimate of GV, Soil, and NPV was calculated for each field. An expected driving factor of water vapor patterns is wind, both its directionality and its magnitude. Therefore, we interpolated a map of wind over the study area during the studytimes and dates in order to create an estimate of wind activity against which to compare water vapor patterns. To calculate wind speed and wind direction, we relied on weather station data from 18 meteorological stations in or within close proximity to the study area . Weather stations used in the analysis are managed by various sources including government agencies, private firms, and educational institutions. Meteorological data were downloaded from the MesoWest and California Irrigation Management Information System networks. For each station, the wind speed and direction that most closely matched the flight time for each of the three dates were recorded. All observations were within 30 minutes of the flight time. Using the 18 data points for wind magnitude and direction, data were spatially interpolated across the study area using an inverse weighted distance formula. IDW relies on the idea that each estimated data point will be influenced by the known data surrounding it, and this influence will diminish with distance. IDW has been widely used in climatic studies for interpolating data such as temperature, rainfall, and wind and is computationally efficient , but has been found to have low accuracy when data is sparse or unevenly distributed . As our study area is relatively flat and stations are somewhat evenly disbursed, we felt IDW would be an appropriate and efficient method for wind interpolation. The result is a pixel-level estimation of wind speed and wind direction across the study area for each of the three study dates. Because small-scale trends of water vapor are observed in the scene and crops are not uniformly distributed throughout the area, direct comparisons of water vapor concentrations between fields or crop types would not be suitable for ET evaluation. Therefore, most of the analyses in this study area are conducted using field-scale, normalized water vapor values. In order to normalize comparisons, intra-field slope, trajectory, and intercept were proposed. To quantify the concentrations and gradients of water vapor as they vary over field of crops, a 3D linear trend surface was fit to the water vapor data over every field in the crop map for each year . Using the fitted surface, an intercept, slope, and trajectory were calculated for each field to explore the magnitude, rate of change, and directionality of the water vapor above it.