Search USGS⌕ Search

SEARCH · Search USGS

Results for “Remote Sensing of the Environment”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

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

At least 145 records · Page 8Linked to original sources
USGS2013

Hyperspectral versus multispectral crop-productivity modeling and type discrimination for the HyspIRI mission

Precise monitoring of agricultural crop biomass and yield quantities is critical for crop production management and prediction. The goal of this study was to compare hyperspectral narrowband (HNB) versus multispectral broadband (MBB) reflectance data in studying irrigated cropland characteristics of five leading world crops (cotton, wheat, maize, rice, and alfalfa) with the objectives of: 1. Modeling crop productivity, and 2. Discriminating crop types. HNB data were obtained from Hyperion hyperspectral imager and field ASD spectroradiometer, and MBB data were obtained from five broadband sensors: Landsat-7 Enhanced Thematic Mapper Plus (ETM +), Advanced Land Imager (ALI), Indian Remote Sensing (IRS), IKONOS, and QuickBird. A large collection of field spectral and biophysical variables were gathered for the 5 crops in Central Asia throughout the growing seasons of 2006 and 2007. Overall, the HNB and hyperspectral vegetation index (HVI) crop biophysical models explained about 25% greater variability when compared with corresponding MBB models. Typically, 3 to 7 HNBs, in multiple linear regression models of a given crop variable, explained more than 93% of variability in crop models. The evaluation of λ 1 (400–2500 nm) versus λ 2 (400–2500 nm) plots of various crop biophysical variables showed that the best two-band normalized difference HVIs involved HNBs centered at: (i) 742 nm and 1175 nm (HVI742-1175), (ii) 1296 nm and 1054 nm (HVI1296-1054), (iii) 1225 nm and 697 nm (HVI1225-697), and (iv) 702 nm and 1104 nm (HVI702-1104). Among the most frequently occurring HNBs in various crop biophysical models, 74% were located in the 1051–2331 nm spectral range, followed by 10% in the moisture sensitive 970 nm, 6% in the red and red-edge (630–752 nm), and the remaining 10% distributed between blue (400–500 nm), green (501–600 nm), and NIR (760–900 nm). Discriminant models, used for discriminating 3 or 4 or 5 crop types, showed significantly higher accuracies when using HNBs (> 90%) over MBBs data (varied between 45 and 84%). Finally, the study highlighted 29 HNBs of Hyperion that are optimal in the study of agricultural crops and potentially significant to the upcoming NASA HyspIRI mission. Determining optimal and redundant bands for a given application will help overcoming the Hughes' phenomenon (or curse of high dimensionality of data).

Syr Darya River Basin↗
USGS2014

Rapid dispersal of saltcedar (Tamarix spp.) biocontrol beetles (Diorhabda carinulata) on a desert river detected by phenocams, MODIS imagery and ground observations

We measured the rate of dispersal of saltcedar leaf beetles ( Diorhabda carinulata ), a defoliating insect released on western rivers to control saltcedar shrubs ( Tamarix spp.), on a 63 km reach of the Virgin River, U.S. Dispersal was measured by satellite imagery, ground surveys and phenocams. Pixels from the Moderate Resolution Imaging Spectrometer (MODIS) sensors on the Terra satellite showed a sharp drop in NDVI in midsummer followed by recovery, correlated with defoliation events as revealed in networked digital camera images and ground surveys. Ground surveys and MODIS imagery showed that beetle damage progressed downstream at a rate of about 25 km yr −1 in 2010 and 2011, producing a 50% reduction in saltcedar leaf area index and evapotranspiration by 2012, as estimated by algorithms based on MODIS Enhanced Vegetation Index values and local meteorological data for Mesquite, Nevada. This reduction is the equivalent of 10.4% of mean annual river flows on this river reach. Our results confirm other observations that saltcedar beetles are dispersing much faster than originally predicted in pre-release biological assessments, presenting new challenges and opportunities for land, water and wildlife managers on western rivers. Despite relatively coarse resolution (250 m) and gridding artifacts, single MODIS pixels can be useful in tracking the effects of defoliating insects in riparian corridors.

Arizona, Nevada↗
USGS2013

Characterizing regional soil mineral composition using spectroscopyand geostatistics

This work aims at improving the mapping of major mineral variability at regional scale using scale-dependent spatial variability observed in remote sensing data. Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) data and statistical methods were combined with laboratory-based mineral characterization of field samples to create maps of the distributions of clay, mica and carbonate minerals and their abundances. The Material Identification and Characterization Algorithm (MICA) was used to identify the spectrally-dominant minerals in field samples; these results were combined with ASTER data using multinomial logistic regression to map mineral distributions. X-ray diffraction (XRD)was used to quantify mineral composition in field samples. XRD results were combined with ASTER data using multiple linear regression to map mineral abundances. We testedwhether smoothing of the ASTER data to match the scale of variability of the target sample would improve model correlations. Smoothing was donewith Fixed Rank Kriging (FRK) to represent the mediumand long-range spatial variability in the ASTER data. Stronger correlations resulted using the smoothed data compared to results obtained with the original data. Highest model accuracies came from using both medium and long-range scaled ASTER data as input to the statistical models. High correlation coefficients were obtained for the abundances of calcite and mica (R2 = 0.71 and 0.70, respectively). Moderately-high correlation coefficients were found for smectite and kaolinite (R2 = 0.57 and 0.45, respectively). Maps of mineral distributions, obtained by relating ASTER data to MICA analysis of field samples, were found to characterize major soil mineral variability (overall accuracies for mica, smectite and kaolinite were 76%, 89% and 86% respectively). The results of this study suggest that the distributions of minerals and their abundances derived using FRK-smoothed ASTER data more closely match the spatial variability of soil and environmental properties at regional scale.

Remote Sensing of Environment↗
USGS2014

Landsat-8: Science and product vision for terrestrial global change research

Landsat 8, a NASA and USGS collaboration, acquires global moderate-resolution measurements of the Earth's terrestrial and polar regions in the visible, near-infrared, short wave, and thermal infrared. Landsat 8 extends the remarkable 40 year Landsat record and has enhanced capabilities including new spectral bands in the blue and cirrus cloud-detection portion of the spectrum, two thermal bands, improved sensor signal-to-noise performance and associated improvements in radiometric resolution, and an improved duty cycle that allows collection of a significantly greater number of images per day. This paper introduces the current (2012–2017) Landsat Science Team's efforts to establish an initial understanding of Landsat 8 capabilities and the steps ahead in support of priorities identified by the team. Preliminary evaluation of Landsat 8 capabilities and identification of new science and applications opportunities are described with respect to calibration and radiometric characterization; surface reflectance; surface albedo; surface temperature, evapotranspiration and drought; agriculture; land cover, condition, disturbance and change; fresh and coastal water; and snow and ice. Insights into the development of derived ‘higher-level’ Landsat products are provided in recognition of the growing need for consistently processed, moderate spatial resolution, large area, long-term terrestrial data records for resource management and for climate and global change studies. The paper concludes with future prospects, emphasizing the opportunities for land imaging constellations by combining Landsat data with data collected from other international sensing systems, and consideration of successor Landsat mission requirements.

Remote Sensing of Environment↗
USGS2014

Monitoring Everglades freshwater marsh water level using L-band synthetic aperture radar backscatter

The Florida Everglades plays a significant role in controlling floods, improving water quality, supporting ecosystems, and maintaining biodiversity in south Florida. Adaptive restoration and management of the Everglades requires the best information possible regarding wetland hydrology. We developed a new and innovative approach to quantify spatial and temporal variations in wetland water levels within the Everglades, Florida. We observed high correlations between water level measured at in situ gages and L-band SAR backscatter coefficients in the freshwater marsh, though C-band SAR backscatter has no close relationship with water level. Here we illustrate the complementarity of SAR backscatter coefficient differencing and interferometry (InSAR) for improved estimation of high spatial resolution water level variations in the Everglades. This technique has a certain limitation in applying to swamp forests with dense vegetation cover, but we conclude that this new method is promising in future applications to wetland hydrology research.

Florida↗
USGS2014

Evaluation of sensor types and environmental controls on mapping biomass of coastal marsh emergent vegetation

There is a need to quantify large-scale plant productivity in coastal marshes to understand marsh resilience to sea level rise, to help define eligibility for carbon offset credits, and to monitor impacts from land use, eutrophication and contamination. Remote monitoring of aboveground biomass of emergent wetland vegetation will help address this need. Differences in sensor spatial resolution, bandwidth, temporal frequency and cost constrain the accuracy of biomass maps produced for management applications. In addition the use of vegetation indices to map biomass may not be effective in wetlands due to confounding effects of water inundation on spectral reflectance. To address these challenges, we used partial least squares regression to select optimal spectral features in situ and with satellite reflectance data to develop predictive models of aboveground biomass for common emergent freshwater marsh species, Typha spp. and Schoenoplectus acutus , at two restored marshes in the Sacramento–San Joaquin River Delta, California, USA. We used field spectrometer data to test model errors associated with hyperspectral narrowbands and multispectral broadbands, the influence of water inundation on prediction accuracy, and the ability to develop species specific models. We used Hyperion data, Digital Globe World View-2 (WV-2) data, and Landsat 7 data to scale up the best statistical models of biomass. Field spectrometer-based models of the full dataset showed that narrowband reflectance data predicted biomass somewhat, though not significantly better than broadband reflectance data [R 2 = 0.46 and percent normalized RMSE (%RMSE) = 16% for narrowband models]. However hyperspectral first derivative reflectance spectra best predicted biomass for plots where water levels were less than 15 cm (R 2 = 0.69, %RMSE = 12.6%). In species-specific models, error rates differed by species ( Typha spp.: %RMSE = 18.5%; S. acutus : %RMSE = 24.9%), likely due to the more vertical structure and deeper water habitat of S. acutus. The Landsat 7 dataset (7 images) predicted biomass slightly better than the WV-2 dataset (6 images) (R 2 = 0.56, %RMSE = 20.9%, compared to R 2 = 0.45, RMSE = 21.5%). The Hyperion dataset (one image) was least successful in predicting biomass (R 2 = 0.27, %RMSE = 33.5%). Shortwave infrared bands on 30 m-resolution Hyperion and Landsat 7 sensors aided biomass estimation; however managers need to weigh tradeoffs between cost, additional spectral information, and high spatial resolution that will identify variability in small, fragmented marshes common to the Sacramento–San Joaquin River Delta and elsewhere in the Western U.S.

California↗
USGS2014

Mapping irrigated areas in Afghanistan over the past decade using MODIS NDVI

Agricultural production capacity contributes to food security in Afghanistan and is largely dependent on irrigated farming, mostly utilizing surface water fed by snowmelt. Because of the high contribution of irrigated crops (> 80%) to total agricultural production, knowing the spatial distribution and year-to-year variability in irrigated areas is imperative to monitoring food security for the country. We used 16-day composites of the Normalized Difference Vegetation Index (NDVI) from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor to create 23-point time series for each year from 2000 through 2013. Seasonal peak values and time series were used in a threshold-dependent decision tree algorithm to map irrigated areas in Afghanistan for the last 14 years. In the absence of ground reference irrigated area information, we evaluated these maps with the irrigated areas classified from multiple snapshots of the landscape during the growing season from Landsat 5 optical and thermal sensor images. We were able to identify irrigated areas using Landsat imagery by selecting as irrigated those areas with Landsat-derived NDVI greater than 0.30–0.45, depending on the date of the Landsat image and surface temperature less than or equal to 310 Kelvin (36.9 ° C). Due to the availability of Landsat images, we were able to compare with the MODIS-derived maps for four years: 2000, 2009, 2010, and 2011. The irrigated areas derived from Landsat agreed well r 2 = 0.91 with the irrigated areas derived from MODIS, providing confidence in the MODIS NDVI threshold approach. The maps portrayed a highly dynamic irrigated agriculture practice in Afghanistan, where the amount of irrigated area was largely determined by the availability of surface water, especially snowmelt, and varied by as much as 30% between water surplus and water deficit years. During the past 14 years, 2001, 2004, and 2008 showed the lowest levels of irrigated area (~ 1.5 million hectares), attesting to the severe drought conditions in those years, whereas 2009, 2012 and 2013 registered the largest irrigated area (~ 2.5 million hectares) due to record snowpack and snowmelt in the region. The model holds promise the ability to provide near-real-time (by the end of the growing seasons) estimates of irrigated area, which are beneficial for food security monitoring as well as subsequent decision making for the country. While the model is developed for Afghanistan, it can be adopted with appropriate adjustments in the derived threshold values to map irrigated areas elsewhere.

Remote Sensing of Environment↗
USGS2014

Marsh dieback, loss, and recovery mapped with satellite optical, airborne polarimetric radar, and field data

Landsat Thematic Mapper and Satellite Pour l'Observation de la Terre (SPOT) satellite based optical sensors, NASA Uninhabited Aerial Vehicle synthetic aperture radar (UAVSAR) polarimetric SAR (PolSAR), and field data captured the occurrence and the recovery of an undetected dieback that occurred between the summers of 2010, 2011, and 2012 in the Spartina alterniflora marshes of coastal Louisiana. Field measurements recorded the dramatic biomass decrease from 2010 to 2011 and a biomass recovery in 2012 dominated by a decrease of live biomass, and the loss of marsh as part of the dieback event. Based on an established relationship, the near-infrared/red vegetation index (VI) and site-specific measurements delineated a contiguous expanse of marsh dieback encompassing 6649.9 ha of 18,292.3 ha of S. alterniflora marshes within the study region. PolSAR data were transformed to variables used in biophysical mapping, and of this variable suite, the cross-polarization HV (horizontal send and vertical receive) backscatter was the best single indicator of marsh dieback and recovery. HV backscatter exhibited substantial and significant changes over the dieback and recovery period, tracked measured biomass changes, and significantly correlated with the live/dead biomass ratio. Within the context of regional trends, both HV and VI indicators started higher in pre-dieback marshes and exhibited substantially and statistically higher variability from year to year than that exhibited in the non-dieback marshes. That distinct difference allowed the capturing of the S. alterniflora marsh dieback and recovery; however, these changes were incorporated in a regional trend exhibiting similar but more subtle biomass composition changes.

Louisiana↗
USGS2014

Modeling a historical mountain pine beetle outbreak using Landsat MSS and multiple lines of evidence

Mountain pine beetles are significant forest disturbance agents, capable of inducing widespread mortality in coniferous forests in western North America. Various remote sensing approaches have assessed the impacts of beetle outbreaks over the last two decades. However, few studies have addressed the impacts of historical mountain pine beetle outbreaks, including the 1970s event that impacted Glacier National Park. The lack of spatially explicit data on this disturbance represents both a major data gap and a critical research challenge in that wildfire has removed some of the evidence from the landscape. We utilized multiple lines of evidence to model forest canopy mortality as a proxy for outbreak severity. We incorporate historical aerial and landscape photos, aerial detection survey data, a nine-year collection of satellite imagery and abiotic data. This study presents a remote sensing based framework to (1) relate measurements of canopy mortality from fine-scale aerial photography to coarse-scale multispectral imagery and (2) classify the severity of mountain pine beetle affected areas using a temporal sequence of Landsat data and other landscape variables. We sampled canopy mortality in 261 plots from aerial photos and found that insect effects on mortality were evident in changes to the Normalized Difference Vegetation Index (NDVI) over time. We tested multiple spectral indices and found that a combination of NDVI and the green band resulted in the strongest model. We report a two-step process where we utilize a generalized least squares model to account for the large-scale variability in the data and a binary regression tree to describe the small-scale variability. The final model had a root mean square error estimate of 9.8% canopy mortality, a mean absolute error of 7.6% and an R 2 of 0.82. The results demonstrate that a model of percent canopy mortality as a continuous variable can be developed to identify a gradient of mountain pine beetle severity on the landscape.

Montana↗
USGS2015

The evolution of mapping habitat for northern spotted owls (Strix occidentalis caurina): A comparison of photo-interpreted, Landsat-based, and lidar-based habitat maps

Wildlife habitat mapping has evolved at a rapid pace over the last few decades. Beginning with simple, often subjective, hand-drawn maps, habitat mapping now involves complex species distribution models (SDMs) using mapped predictor variables derived from remotely sensed data. For species that inhabit large geographic areas, remote sensing technology is often essential for producing range wide maps. Habitat monitoring for northern spotted owls ( Strix occidentalis caurina ), whose geographic covers about 23 million ha, is based on SDMs that use Landsat Thematic Mapper imagery to create forest vegetation data layers using gradient nearest neighbor (GNN) methods. Vegetation data layers derived from GNN are modeled relationships between forest inventory plot data, climate and topographic data, and the spectral signatures acquired by the satellite. When used as predictor variables for SDMs, there is some transference of the GNN modeling error to the final habitat map. Recent increases in the use of light detection and ranging (lidar) data, coupled with the need to produce spatially accurate and detailed forest vegetation maps have spurred interest in its use for SDMs and habitat mapping. Instead of modeling predictor variables from remotely sensed spectral data, lidar provides direct measurements of vegetation height for use in SDMs. We expect a SDM habitat map produced from directly measured predictor variables to be more accurate than one produced from modeled predictors. We used maximum entropy (Maxent) SDM modeling software to compare predictive performance and estimates of habitat area between Landsat-based and lidar-based northern spotted owl SDMs and habitat maps. We explored the differences and similarities between these maps, and to a pre-existing aerial photo-interpreted habitat map produced by local wildlife biologists. The lidar-based map had the highest predictive performance based on 10 bootstrapped replicate models (AUC = 0.809 ± 0.011), but the performance of the Landsat-based map was within acceptable limits (AUC = 0.717 ± 0.021). As is common with photo-interpreted maps, there was no accuracy assessment available for comparison. The photo-interpreted map produced the highest and lowest estimates of habitat area, depending on which habitat classes were included (nesting, roosting, and foraging habitat = 9962 ha, nesting habitat only = 6036 ha). The Landsat-based map produced an estimate of habitat area that was within this range (95% CI: 6679–9592 ha), while the lidar-based map produced an area estimate similar to what was interpreted by local wildlife biologists as nesting (i.e., high quality) habitat using aerial imagery (95% CI: 5453–7216). Confidence intervals of habitat area estimates from the SDMs based on Landsat and lidar overlapped. We concluded that both Landsat- and lidar-based SDMs produced reasonable maps and area estimates for northern spotted owl habitat within the study area. The lidar-based map was more precise and spatially similar to what local wildlife biologists considered spotted owl nesting habitat. The Landsat-based map provided a less precise spatial representation of habitat within the relatively small geographic confines of the study area, but habitat area estimates were similar to both the photo-interpreted and lidar-based maps. Photo-interpreted maps are time consuming to produce, subjective in nature, and difficult to replicate. SDMs provide a framework for efficiently producing habitat maps that can be replicated as habitat conditions change over time, provided that comparable remotely sensed data are available. When the SDM uses predictor variables extracted from lidar data, it can produce a habitat map that is both accurate and useful at large and small spatial scales. In comparison, SDMs using Landsat-based data are more appropriate for large scale analyses of amounts and general spatial patterns of habitat at regional scales.

Remote Sensing of Environment↗
USGS2015

Examining the utility of satellite-based wind sheltering estimates for lake hydrodynamic modeling

Satellite-based measurements of vegetation canopy structure have been in common use for the last decade but have never been used to estimate canopy's impact on wind sheltering of individual lakes. Wind sheltering is caused by slower winds in the wake of topography and shoreline obstacles (e.g. forest canopy) and influences heat loss and the flux of wind-driven mixing energy into lakes, which control lake temperatures and indirectly structure lake ecosystem processes, including carbon cycling and thermal habitat partitioning. Lakeshore wind sheltering has often been parameterized by lake surface area but such empirical relationships are only based on forested lakeshores and overlook the contributions of local land cover and terrain to wind sheltering. This study is the first to examine the utility of satellite imagery-derived broad-scale estimates of wind sheltering across a diversity of land covers. Using 30 m spatial resolution ASTER GDEM2 elevation data, the mean sheltering height, h s , being the combination of local topographic rise and canopy height above the lake surface, is calculated within 100 m-wide buffers surrounding 76,000 lakes in the U.S. state of Wisconsin. Uncertainty of GDEM2-derived h s was compared to SRTM-, high-resolution G-LiHT lidar-, and ICESat-derived estimates of h s , respective influences of land cover type and buffer width on h s are examined; and the effect of including satellite-based h s on the accuracy of a statewide lake hydrodynamic model was discussed. Though GDEM2 h s uncertainty was comparable to or better than other satellite-based measures of h s , its higher spatial resolution and broader spatial coverage allowed more lakes to be included in modeling efforts. GDEM2 was shown to offer superior utility for estimating h s compared to other satellite-derived data, but was limited by its consistent underestimation of h s , inability to detect within-buffer h s variability, and differing accuracy across land cover types. Nonetheless, considering a GDEM2 h s -derived wind sheltering potential improved the modeled lake temperature root mean square error for non-forested lakes by 0.72 °C compared to a commonly used wind sheltering model based on lake area alone. While results from this study show promise, the limitations of near-global GDEM2 data in timeliness, temporal and spatial resolution, and vertical accuracy were apparent. As hydrodynamic modeling and high-resolution topographic mapping efforts both expand, future remote sensing-derived vegetation structure data must be improved to meet wind sheltering accuracy requirements to expand our understanding of lake processes.

Remote Sensing of Environment↗
USGS2015

Oil detection in the coastal marshes of Louisiana using MESMA applied to band subsets of AVIRIS data

We mapped oil presence in the marshes of Barataria Bay, Louisiana following the Deepwater Horizon oil spill using Airborne Visible InfraRed Imaging Spectrometer (AVIRIS) data. Oil and non-photosynthetic vegetation (NPV) have very similar spectra, differing only in two narrow hydrocarbon absorption regions around 1700 and 2300 nm. Confusion between NPV and oil is expressed as an increase in oil fraction error with increasing NPV, as shown by Multiple Endmember Spectral Mixture Analysis (MESMA) applied to synthetic spectra generated with known endmember fractions. Significantly, the magnitude of error varied depending upon the type of NPV in the mixture. To reduce error, we used stable zone unmixing to identify a nine band subset that emphasized the hydrocarbon absorption regions, allowing for more accurate detection of oil presence using MESMA. When this band subset was applied to post-spill AVIRIS data acquired over Barataria Bay on several dates following the 2010 oil spill, accuracies ranged from 87.5% to 93.3%. Oil presence extended 10.5 m into the marsh for oiled shorelines, showing a reduced oil fraction with increasing distance from the shoreline.

Louisiana↗
USGS2015

A global reference database from very high resolution commercial satellite data and methodology for application to Landsat derived 30 m continuous field tree cover data

The methodology for selection, creation, and application of a global remote sensing validation dataset using high resolution commercial satellite data is presented. High resolution data are obtained for a stratified random sample of 500 primary sampling units (5 km × 5 km sample blocks), where the stratification based on Köppen climate classes is used to distribute the sample globally among biomes. The high resolution data are classified to categorical land cover maps using an analyst mediated classification workflow. Our initial application of these data is to evaluate a global 30 m Landsat-derived, continuous field tree cover product. For this application, the categorical reference classification produced at 2 m resolution is converted to percent tree cover per 30 m pixel (secondary sampling unit)for comparison to Landsat-derived estimates of tree cover. We provide example results (based on a subsample of 25 sample blocks in South America) illustrating basic analyses of agreement that can be produced from these reference data. Commercial high resolution data availability and data quality are shown to provide a viable means of validating continuous field tree cover. When completed, the reference classifications for the full sample of 500 blocks will be released for public use.

Remote Sensing of Environment↗
USGS2015

Characterization of shrubland ecosystem components as continuous fields in the northwest United States

Accurate and consistent estimates of shrubland ecosystem components are crucial to a better understanding of ecosystem conditions in arid and semiarid lands. An innovative approach was developed by integrating multiple sources of information to quantify shrubland components as continuous field products within the National Land Cover Database (NLCD). The approach consists of several procedures including field sample collections, high-resolution mapping of shrubland components using WorldView-2 imagery and regression tree models, Landsat 8 radiometric balancing and phenological mosaicking, medium resolution estimates of shrubland components following different climate zones using Landsat 8 phenological mosaics and regression tree models, and product validation. Fractional covers of nine shrubland components were estimated: annual herbaceous, bare ground, big sagebrush, herbaceous, litter, sagebrush, shrub, sagebrush height, and shrub height. Our study area included the footprint of six Landsat 8 scenes in the northwestern United States. Results show that most components have relatively significant correlations with validation data, have small normalized root mean square errors, and correspond well with expected ecological gradients. While some uncertainties remain with height estimates, the model formulated in this study provides a cross-validated, unbiased, and cost effective approach to quantify shrubland components at a regional scale and advances knowledge of horizontal and vertical variability of these components.

California, Idaho, Nevada, Oregon↗
USGS2015

Distribution of near-surface permafrost in Alaska: estimates of present and future conditions

High-latitude regions are experiencing rapid and extensive changes in ecosystem composition and function as the result of increases in average air temperature. Increasing air temperatures have led to widespread thawing and degradation of permafrost, which in turn has affected ecosystems, socioeconomics, and the carbon cycle of high latitudes. Here we overcome complex interactions among surface and subsurface conditions to map nearsurface permafrost through decision and regression tree approaches that statistically and spatially extend field observations using remotely sensed imagery, climatic data, and thematic maps of a wide range of surface and subsurface biophysical characteristics. The data fusion approach generated medium-resolution (30-m pixels) maps of near-surface (within 1 m) permafrost, active-layer thickness, and associated uncertainty estimates throughout mainland Alaska. Our calibrated models (overall test accuracy of ~85%) were used to quantify changes in permafrost distribution under varying future climate scenarios assuming no other changes in biophysical factors. Models indicate that near-surface permafrost underlies 38% of mainland Alaska and that near-surface permafrost will disappear on 16 to 24% of the landscape by the end of the 21st Century. Simulations suggest that near-surface permafrost degradation is more probable in central regions of Alaska than more northerly regions. Taken together, these results have obvious implications for potential remobilization of frozen soil carbon pools under warmer temperatures. Additionally, warmer and drier conditions may increase fire activity and severity, which may exacerbate rates of permafrost thaw and carbon remobilization relative to climate alone. The mapping of permafrost distribution across Alaska is important for land-use planning, environmental assessments, and a wide-array of geophysical studies.

Alaska↗
USGS2015

LiDAR based prediction of forest biomass using hierarchical models with spatially varying coefficients

Many studies and production inventory systems have shown the utility of coupling covariates derived from Light Detection and Ranging (LiDAR) data with forest variables measured on georeferenced inventory plots through regression models. The objective of this study was to propose and assess the use of a Bayesian hierarchical modeling framework that accommodates both residual spatial dependence and non-stationarity of model covariates through the introduction of spatial random effects. We explored this objective using four forest inventory datasets that are part of the North American Carbon Program, each comprising point-referenced measures of above-ground forest biomass and discrete LiDAR. For each dataset, we considered at least five regression model specifications of varying complexity. Models were assessed based on goodness of fit criteria and predictive performance using a 10-fold cross-validation procedure. Results showed that the addition of spatial random effects to the regression model intercept improved fit and predictive performance in the presence of substantial residual spatial dependence. Additionally, in some cases, allowing either some or all regression slope parameters to vary spatially, via the addition of spatial random effects, further improved model fit and predictive performance. In other instances, models showed improved fit but decreased predictive performance—indicating over-fitting and underscoring the need for cross-validation to assess predictive ability. The proposed Bayesian modeling framework provided access to pixel-level posterior predictive distributions that were useful for uncertainty mapping, diagnosing spatial extrapolation issues, revealing missing model covariates, and discovering locally significant parameters.

Colorado, Minnesota↗
USGS2015

Combining NLCD and MODIS to create a land cover-albedo database for the continental United States

Land surface albedo is an essential climate variable that is tightly linked to land cover, such that specific land cover classes (e.g., deciduous broadleaf forest, cropland) have characteristic albedos. Despite the normative of land-cover class specific albedos, there is considerable variability in albedo within a land cover class. The National Land Cover Database (NLCD) and the Moderate Resolution Imaging Spectroradiometer (MODIS) albedo product were combined to produce a long-term (14 years) integrated land cover-albedo database for the continental United States that can be used to examine the temporal behavior of albedo as a function of land cover. The integration identifies areas of homogeneous land cover at the nominal spatial resolution of the MODIS (MCD43A) albedo product (500 m × 500 m) from the NLCD product (30 m × 30 m), and provides an albedo data record per 500 m × 500 m pixel for 14 of the 16 NLCD land cover classes. Individual homogeneous land cover pixels have up to 605 albedo observations, and 75% of the pixels have at least 319 MODIS albedo observations (≥ 50% of the maximum possible number of observations) for the study period (2000–2013). We demonstrated the utility of the database by conducting a multivariate analysis of variance of albedo for each NLCD land cover class, showing that locational (pixel-to-pixel) and inter-annual variability were significant factors in addition to expected seasonal (intra-annual) and geographic (latitudinal) effects.

Remote Sensing of Environment↗
USGS2015

Quantifying soil carbon loss and uncertainty from a peatland wildfire using multi-temporal LiDAR

Peatlands are a major reservoir of global soil carbon, yet account for just 3% of global land cover. Human impacts like draining can hinder the ability of peatlands to sequester carbon and expose their soils to fire under dry conditions. Estimating soil carbon loss from peat fires can be challenging due to uncertainty about pre-fire surface elevations. This study uses multi-temporal LiDAR to obtain pre- and post-fire elevations and estimate soil carbon loss caused by the 2011 Lateral West fire in the Great Dismal Swamp National Wildlife Refuge, VA, USA. We also determine how LiDAR elevation error affects uncertainty in our carbon loss estimate by randomly perturbing the LiDAR point elevations and recalculating elevation change and carbon loss, iterating this process 1000 times. We calculated a total loss using LiDAR of 1.10 Tg C across the 25 km 2 burned area. The fire burned an average of 47 cm deep, equivalent to 44 kg C/m 2 , a value larger than the 1997 Indonesian peat fires (29 kg C/m 2 ). Carbon loss via the First-Order Fire Effects Model (FOFEM) was estimated to be 0.06 Tg C. Propagating the LiDAR elevation error to the carbon loss estimates, we calculated a standard deviation of 0.00009 Tg C, equivalent to 0.008% of total carbon loss. We conclude that LiDAR elevation error is not a significant contributor to uncertainty in soil carbon loss under severe fire conditions with substantial peat consumption. However, uncertainties may be more substantial when soil elevation loss is of a similar or smaller magnitude than the reported LiDAR error.

Virginia↗