ENHANCING CONSERVATION WITH HIGH RESOLUTION PRODUCTIVITY DATASETS FOR THE CONTERMINOUS UNITED STATES

Human driven alteration of the earth’s terrestrial surface is accelerating through land use changes, intensification of human activity, climate change, and other anthropogenic pressures. These changes occur at broad spatio-temporal scales, challenging our ability to effectively monitor and assess the impacts and subsequent conservation strategies. While satellite remote sensing (SRS) products enable monitoring of the earth’s terrestrial surface continuously across space and time, the practical applications for conservation and management of these products are limited. Often the processes driving ecological change occur at fine spatial resolutions and are undetectable given the resolution of available datasets. Additionally, the links between SRS data and ecologically meaningful metrics are weak. Recent advances in cloud computing technology along with the growing record of high resolution SRS data enable the development of SRS products that quantify ecologically meaningful variables at relevant scales applicable for conservation and management. The focus of my dissertation is to improve the applicability of terrestrial gross and net primary productivity (GPP/NPP) datasets for the conterminous United States (CONUS). In chapter one, I develop a framework for creating high resolution datasets of vegetation dynamics. I use the entire archive of Landsat 5, 7, and 8 surface reflectance data and a novel gap filling approach to create spatially continuous 30 m, 16-day composites of the normalized difference vegetation index (NDVI) from 1986 to 2016. In chapter two, I integrate this with other high resolution datasets and the MOD17 algorithm to create the first high resolution GPP and NPP datasets for CONUS. I demonstrate the applicability of these products for conservation and management, showing the improvements beyond currently available products. In chapter three, I utilize this dataset to evaluate the relationships between land ownership and terrestrial production across the CONUS domain. The main results of this work are three publically available datasets: 1) 30 m Landsat NDVI; 2) 250 m MODIS based GPP and NPP; and 3) 30 m Landsat based GPP and NPP. My goal is that these products prove useful for the wider scientific, conservation, and land management communities as we continue to strive for better conservation and management practices

Development and application of landsat-based wetland vegetation cover and unvegetated-vegetated marsh ratio (UVVR) for the conterminous United States

© The Author(s), 2022. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Ganju, N. K., Couvillion, B. R., Defne, Z., & Ackerman, K. Development and application of landsat-based wetland vegetation cover and unvegetated-vegetated marsh ratio (UVVR) for the conterminous United States. Estuaries and Coasts, (2022), https://doi.org/10.1007/s12237-022-01081-x.Effective management and restoration of salt marshes and other vegetated intertidal habitats require objective and spatially integrated metrics of geomorphic status and vulnerability. The unvegetated-vegetated marsh ratio (UVVR), a recently developed metric, can be used to establish present-day vegetative cover, identify stability thresholds, and quantify vulnerability to open-water conversion over a range of spatial scales. We developed a Landsat-based approach to quantify the within-pixel vegetated fraction and UVVR for coastal wetlands of the conterminous United States, at 30-m resolution for 2014–2018. Here we present the methodology used to generate the UVVR from spectral indices, along with calibration, validation, and spatial autocorrelation assessments. We then demonstrate multiple applications of the data across varying spatial scales: first, we aggregate the UVVR across individual states and estuaries to quantify total vegetated wetland area for the nation. On the state level, Louisiana and Florida account for over 50% of the nation’s total, while on the estuarine level, the Chesapeake Bay Estuary and selected Louisiana coastal areas each account for over 6% of the nation’s total vegetated wetland area. Second, we present cases where this dataset can be used to track wetland change (e.g., expansion due to restoration and loss due to stressors). Lastly, we propose a classification methodology that delineates areas vulnerable to open-water expansion based on the 5-year mean and standard deviation of the UVVR. Calculating the UVVR for the period-of-record back to 1985, as well as regular updating, will fill a critical gap for tracking national status of salt marshes and other vegetated habitats through time and space.This work was supported by the U.S. Geological Survey’s Coastal and Marine Hazards/Resources Program

Terrestrial primary production for the conterminous United States derived from Landsat 30 m and MODIS 250 m

Terrestrial primary production is a fundamental ecological process and a crucial component in understanding the flow of energy through trophic levels. The global MODIS gross primary production (GPP) and net primary production (NPP) products (MOD17) are widely used for monitoring GPP and NPP at coarse resolutions across broad spatial extents. The coarse input datasets and global biome‐level parameters, however, are well‐known limitations to the applicability of the MOD17 product at finer scales. We addressed these limitations and created two improved products for the conterminous United States (CONUS) that capture the spatiotemporal variability in terrestrial production. The MOD17 algorithm was utilized with medium resolution land cover classifications and improved meteorological data specific to CONUS in order to produce: (a) Landsat derived 16‐day GPP and annual NPP at 30 m resolution from 1986 to 2016 (GPPL30 and NPPL30, respectively); and (b) MODIS derived 8‐day GPP and annual NPP at 250 m resolution from 2001 to 2016 (GPPM250 and NPPM250 respectively). Biome‐specific input parameters were optimized based on eddy covariance flux tower‐derived GPP data from the FLUXNET2015 database. We evaluated GPPL30 and GPPM250 products against the standard MODIS GPP product utilizing a select subset of representative flux tower sites, and found improvement across all land cover classes except croplands. We also found consistent interannual variability and trends across NPPL30, NPPM250, and the standard MODIS NPP product. We highlight the application potential of the production products, demonstrating their improved capacity for monitoring terrestrial production at higher levels of spatial detail across broad spatiotemporal scales

Analysis of Normalized Difference Vegetation Index change of the West Bank, Palestine, Using Multitemporal Satellite Remote Sensing Data

The West Bank is characterized by the diversity of its climate despite its small area. It includes four climatic regions:- a humid, semi-humid, arid, and semi-arid climate. This in turn affected the geographical distribution of vegetation cover seasonally and over the years. This study investigated changes in the West Bank, Palestine vegetation cover using multitemporal Landsat data. Four images were selected for this purpose – two corresponding to 2001 and the other two corresponding to 2021. Seasonal change of the Normalized Difference Vegetation Index (NDVI) was investigated for the acquired images. ArcGIS 10.8 software was used for image processing and analysis. Results showed that the negative change of the NDVI of autumn for both dates is much higher than the positive change. About 98.89 percent of the West Bank area scored a negative change in 2001, while about 80.79 percent of the West Bank scored a negative change in 2021. The negative change in autumn is due to the dryness of summer season, which in turn led to the drying of the green grasses and herbs that grow among the trees and on the eastern slopes of the highlands. Concerning the NDVI change between 2001 and 2021, results showed a considerable positive change in the NDVI. Around 90.91% of the surface area of the West Bank (5132.304 km²) has positive change, while only 9.09% (513.405 km²) which mainly represents urban centers has negative change. The main reason that the NDVI in 2021 was higher than the NDVI in 2001 outside the urban centers is that the amount of rain in that year was greater, especially in the northern areas. Thematic maps of NDVI for the two dates were produced, and changes in vegetation cover were extracted from the four maps. Based on the results of the study, it is recommended that decision-makers in the Ministry of Agriculture and Environmental organizations develop plans to increase permanent green spaces, especially in arid and semi-arid areas of the West Bank

Monitoring Landscape Dynamics in Central U.S. Grasslands with Harmonized Landsat-8 and Sentinel-2 Time Series Data

Remotely monitoring changes in central U.S. grasslands is challenging because these landscapes tend to respond quickly to disturbances and changes in weather. Such dynamic responses influence nutrient cycling, greenhouse gas contributions, habitat availability for wildlife, and other ecosystem processes and services. Traditionally, coarse-resolution satellite data acquired at daily intervals have been used for monitoring. Recently, the harmonized Landsat-8 and Sentinel-2 (HLS) data increased the temporal frequency of the data. Here we investigated if the increased data frequency provided adequate observations to characterize highly dynamic grassland processes. We evaluated HLS data available for 2016 to (1) determine if data from Sentinel-2 contributed to an improvement in characterizing landscape processes over Landsat-8 data alone, and (2) quantify how observation frequency impacted results. Specifically, we investigated into estimating annual vegetation phenology, detecting burn scars from fire, and modeling within-season wetland hydroperiod and growth of aquatic vegetation. We observed increased sensitivity to the start of the growing season (SOST) with the HLS data. Our estimates of the grassland SOST compared well with ground estimates collected at a phenological camera site. We used the Continuous Change Detection and Classification (CCDC) algorithm to assess if the HLS data improved our detection of burn scars following grassland fires and found that detection was considerably influenced by the seasonal timing of the fires. The grassland burned in early spring recovered too quickly to be detected as change events by CCDC; instead, the spectral characteristics following these fires were incorporated as part of the ongoing time-series models. In contrast, the spectral effects from late-season fires were detected both by Landsat-8 data and HLS data. For wetland-rich areas, we used a modified version of the CCDC algorithm to track within-season dynamics of water and aquatic vegetation. The addition of Sentinel-2 data provided the potential to build full time series models to better distinguish different wetland types, suggesting that the temporal density of data was sufficient for within-season characterization of wetland dynamics. Although the different data frequency, in both the spatial and temporal dimensions, could cause inconsistent model estimation or sensitivity sometimes; overall, the temporal frequency of the HLS data improved our ability to track within-season grassland dynamics and improved results for areas prone to cloud contamination. The results suggest a greater frequency of observations, such as from harmonizing data across all comparable Landsat and Sentinel sensors, is still needed. For our study areas, at least a 3-day revisit interval during the early growing season (weeks 14–17) is required to provide a \u3e50% probability of obtaining weekly clear observations

EFFECTS OF SPATIAL RESOLUTION AND LANDSCAPE STRUCTURE ON LAND COVER CHARACTERIZATION

This dissertation addressed problems in scaling, problems that are among the main challenges in remote sensing. The principal objective of the research was to investigate the effects of changing spatial scale on the representation of land cover. A second objective was to determine the relationship between such effects, characteristics of landscape structure and scaling procedures. Four research issues related to spatial scaling were examined. They included: 1) the upscaling of Normalized Difference Vegetation Index (NDVI); 2) the effects of spatial scale on indices of landscape structure; 3) the representation of land cover databases at different spatial scales; and 4) the relationships between landscape indices and land cover area estimations. The overall bias resulting from non-linearity of NDVI in relation to spatial resolution is generally insignificant as compared to other factors such as influences of aerosols and water vapor. The bias is, however, related to land surface characteristics. Significant errors may be introduced in heterogeneous areas where different land cover types exhibit strong spectral contrast. Spatially upscaled SPOT and TM NDVIs have information content comparable with the AVHRR-derived NDVI: Indices of landscape structure and spatial resolution are generally related, but the exact forms of the relationships are subject to changes in other factors including the basic patch unit constituting a landscape and the proportional area of foreground land cover under consideration. The extent of agreement between spatially aggregated coarse resolution land cover datasets and full resolution datasets changes with the properties of the original datasets, including the pixel size and class definition. There are close relationships between landscape structure and class areas estimated from spatially aggregated land cover databases. The relationships, however, do not permit extension from one area to another. Inversion calibration across different geographic/ecological areas is, therefore, not feasible. Different rules govern the land cover area changes across resolutions when different upscaling methods are used. Special attention should be given to comparison between land cover maps derived using different methods

Characterization of Landsat-7 to Landsat-8 Reflective Wavelength and Normalized Difference Vegetation Index Continuity

At over 40 years, the Landsat satellites provide the longest temporal record of space-based land surface observations, and the successful 2013 launch of the Landsat-8 is continuing this legacy. Ideally, the Landsat data record should be consistent over the Landsat sensor series. The Landsat-8 Operational Land Imager (OLI) has improved calibration, signal to noise characteristics, higher 12-bit radiometric resolution, and spectrally narrower wavebands than the previous Landsat-7 Enhanced Thematic Mapper (ETM +). Reflective wavelength differences between the two Landsat sensors depend also on the surface reflectance and atmospheric state which are difficult to model comprehensively. The orbit and sensing geometries of the Landsat-8 OLI and Landsat-7 ETM + provide swath edge overlapping paths sensed only one day apart. The overlap regions are sensed in alternating backscatter and forward scattering orientations so Landsat bi-directional reflectance effects are evident but approximately balanced between the two sensors when large amounts of time series data are considered. Taking advantage of this configuration a total of 59 million 30 m corresponding sensor observations extracted from 6317 Landsat-7 ETM + and Landsat-8 OLI images acquired over three winter and three summer months for all the conterminous United States (CONUS) are compared. Results considering different stages of cloud and saturation filtering, and filtering to reduce one day surface state differences, demonstrate the importance of appropriate per-pixel data screening. Top of atmosphere (TOA) and atmospherically corrected surface reflectance for the spectrally corresponding visible, near infrared and shortwave infrared bands, and derived normalized difference vegetation index (NDVI), are compared and their differences quantified. On average the OLI TOA reflectance is greater than the ETM + TOA reflectance for all bands, with greatest differences in the near-infrared (NIR) and the shortwave infrared bands due to the quite different spectral response functions between the sensors. The atmospheric correction reduces the mean difference in the NIR and shortwave infrared but increases the mean difference in the visible bands. Regardless of whether TOA or surface reflectance are used to generate NDVI, on average, for vegetated soil and vegetation surfaces (0 ≤ NDVI ≤ 1), the OLI NDVI is greater than the ETM + NDVI. Statistical functions to transform between the comparable sensor bands and sensor NDVI values are presented so that the user community may apply them in their own research to improve temporal continuity between the Landsat-7 ETM + and Landsat-8 OLI sensor data. The transformation functions were developed using ordinary least squares (OLS) regression and were fit quite reliably (r2 values \u3e 0.7 for the reflectance data and \u3e 0.9 for the NDVI data, p-values \u3c 0.0001)

Intercomparison of phenological transition dates derived from the PhenoCam Dataset V1.0 and MODIS satellite remote sensing

Phenology is a valuable diagnostic of ecosystem health, and has applications to environmental monitoring and management. Here, we conduct an intercomparison analysis using phenological transition dates derived from near-surface PhenoCam imagery and MODIS satellite remote sensing. We used approximately 600 site-years of data, from 128 camera sites covering a wide range of vegetation types and climate zones. During both “greenness rising” and “greenness falling” transition phases, we found generally good agreement between PhenoCam and MODIS transition dates for agricultural, deciduous forest, and grassland sites, provided that the vegetation in the camera field of view was representative of the broader landscape. The correlation between PhenoCam and MODIS transition dates was poor for evergreen forest sites. We discuss potential reasons (including sub-pixel spatial heterogeneity, flexibility of the transition date extraction method, vegetation index sensitivity in evergreen systems, and PhenoCam geolocation uncertainty) for varying agreement between time series of vegetation indices derived from PhenoCam and MODIS imagery. This analysis increases our confidence in the ability of satellite remote sensing to accurately characterize seasonal dynamics in a range of ecosystems, and provides a basis for interpreting those dynamics in the context of tangible phenological changes occurring on the ground

Improved Time Series Land Cover Classification by Missing-Observation-Adaptive Nonlinear Dimensionality Reduction

Dimensionality reduction (DR) is a widely used technique to address the curse of dimensionality when high-dimensional remotely sensed data, such as multi-temporal or hyperspectral imagery, are analyzed. Nonlinear DR algorithms, also referred to as manifold learning algorithms, have been successfully applied to hyperspectral data and provide improved performance compared with linear DR algorithms. However, DR algorithms cannot handle missing data that are common in multi-temporal imagery. In this paper, the Laplacian Eigenmaps (LE) nonlinear DR algorithm was refined for application to multi-temporal satellite data with large proportions of missing data. Refined LE algorithms were applied to 52-week Landsat time series for three study areas in Texas, Kansas and South Dakota that have different amounts of missing data and land cover complexity. A series of random forest classifications were conducted on the refined LE DR bands using varying proportions of training data provided by the United States Department of Agriculture (USDA) National Agricultural Statistics Service (NASS) Cropland Data Layer (CDL); these classification results were compared with conventional metrics-based random forest classifications. Experimental results show that compared with the metrics approach, higher per-class and overall classification accuracies were obtained using the refined LE DR bands of multispectral reflectance time series, and the number of training samples required to achieve a given degree of classification accuracy was also reduced. The approach of applying the refined LE to multispectral reflectance time series is promising in that it is automated and provides dimensionality-reduced data with desirable classification properties. The implications of this research and possibilities for future algorithm development and application are discussed