Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Remote Sensing of 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 199 records · Page 11Linked to original sources
USGS2016

Landsat 8: The plans, the reality, and the legacy

Landsat 8, originally known as the Landsat Data Continuity Mission (LDCM), is a National Aeronautics and Space Administration (NASA)-U.S. Geological Survey (USGS) partnership that continues the legacy of continuous moderate resolution observations started in 1972. The conception of LDCM to the reality of Landsat 8 followed an arduous path extending over nearly 13 years, but the successful launch on February 11, 2013 ensures the continuity of the unparalleled Landsat record. The USGS took over mission operations on May 30, 2013 and renamed LCDM to Landsat 8. Access to Landsat 8 data was opened to users worldwide. Three years following launch we evaluate the science and applications impact of Landsat 8. With a mission objective to enable the detection and characterization of global land changes at a scale where differentiation between natural and human-induced causes of change is possible, LDCM promised incremental technical improvements in capabilities needed for Landsat scientific and applications investigations. Results show that with Landsat 8, we are acquiring more data than ever before, the radiometric and geometric quality of data are generally technically superior to data acquired by past Landsat missions, and the new measurements, e.g., the coastal aerosol and cirrus bands, are opening new opportunities. Collectively, these improvements are sparking the growth of science and applications opportunities. Equally important, with Landsat 7 still operational, we have returned to global imaging on an 8-day cycle, a capability that ended when Landsat 5 ceased operational Earth imaging in November 2011. As a result, the Landsat program is on secure footings and planning is underway to extend the record for another 20 or more years.

Remote Sensing of Environment↗
USGS2016

Marsh canopy structure changes and the Deepwater Horizon oil spill

Marsh canopy structure was mapped yearly from 2009 to 2012 in the Barataria Bay, Louisiana coastal region that was impacted by the 2010 Deepwater Horizon (DWH) oil spill. Based on the previously demonstrated capability of NASA's UAVSAR polarimetric synthetic aperture radar (PolSAR) image data to map Spartina alterniflora marsh canopy structure, structure maps combining the leaf area index (LAI) and leaf angle distribution (LAD, orientation) were constructed for yearly intervals that were directly relatable to the 2010 LAI-LAD classification. The yearly LAI-LAD and LAI difference maps were used to investigate causes for the previously revealed dramatic change in marsh structure from prespill (2009) to postspill (2010, spill cessation), and the occurrence of structure features that exhibited abnormal spatial and temporal patterns. Water level and salinity records showed that freshwater releases used to keep the oil offshore did not cause the rapid growth from 2009 to 2010 in marsh surrounding the inner Bay. Photointerpretation of optical image data determined that interior marsh patches exhibiting rapid change were caused by burns and burn recovery, and that the pattern of 2010 to 2011 LAI decreases in backshore marsh and extending along some tidal channels into the interior marsh were not associated with burns. Instead, the majority of 2010 to 2011 shoreline features aligned with vectors displaying the severity of 2010 shoreline oiling from the DWH spill. Although the association is not conclusive of a causal oil impact, the coexistent pattern is a significant discovery. PolSAR marsh structure mapping provided a unique perspective of marsh biophysical status that enhanced detection of change and monitoring of trends important to management effectiveness.

Louisiana↗
USGS2016

A note on the temporary misregistration of Landsat-8 Operational Land Imager (OLI) and Sentinel-2 Multi Spectral Instrument (MSI) imagery

The Landsat-8 and Sentinel-2 sensors provide multi-spectral image data with similar spectral and spatial characteristics that together provide improved temporal coverage globally. Both systems are designed to register Level 1 products to a reference image framework, however, the Landsat-8 framework, based upon the Global Land Survey images, contains residual geolocation errors leading to an expected sensor-to-sensor misregistration of 38 m (2σ). These misalignments vary geographically but should be stable for a given area. The Landsat framework will be readjusted for consistency with the Sentinel-2 Global Reference Image, with completion expected in 2018. In the interim, users can measure Landsat-to-Sentinel tie points to quantify the misalignment in their area of interest and if appropriate to reproject the data to better alignment.

Remote Sensing of Environment↗
USGS2016

Statistical correction of lidar-derived digital elevation models with multispectral airborne imagery in tidal marshes

Airborne light detection and ranging (lidar) is a valuable tool for collecting large amounts of elevation data across large areas; however, the limited ability to penetrate dense vegetation with lidar hinders its usefulness for measuring tidal marsh platforms. Methods to correct lidar elevation data are available, but a reliable method that requires limited field work and maintains spatial resolution is lacking. We present a novel method, the Lidar Elevation Adjustment with NDVI (LEAN), to correct lidar digital elevation models (DEMs) with vegetation indices from readily available multispectral airborne imagery (NAIP) and RTK-GPS surveys. Using 17 study sites along the Pacific coast of the U.S., we achieved an average root mean squared error (RMSE) of 0.072 m, with a 40–75% improvement in accuracy from the lidar bare earth DEM. Results from our method compared favorably with results from three other methods (minimum-bin gridding, mean error correction, and vegetation correction factors), and a power analysis applying our extensive RTK-GPS dataset showed that on average 118 points were necessary to calibrate a site-specific correction model for tidal marshes along the Pacific coast. By using available imagery and with minimal field surveys, we showed that lidar-derived DEMs can be adjusted for greater accuracy while maintaining high (1 m) resolution.

Remote Sensing of Environment↗
USGS2016

Detecting seasonal landslide movement within the Cascade landslide complex (Washington) using time-series SAR imagery

Detection of slow or limited landslide movement within broad areas of forested terrain has long been problematic, particularly for the Cascade landslide complex (Washington) located along the Columbia River Gorge. Although parts of the landslide complex have been found reactivated in recent years, the timing and magnitude of motion have not been systematically monitored or interpreted. Here we apply novel time-series strategies to study the spatial distribution and temporal behavior of the landslide movement between 2007 and 2011 using InSAR images from two overlapping L-band ALOS PALSAR-1 satellite tracks. Our results show that the reactivated part has moved approximately 700 mm downslope during the 4-year observation period, while other parts of the landslide complex have generally remained stable. However, we also detect about 300 mm of seasonal downslope creep in a terrain block upslope of the Cascade landslide complex—terrain previously thought to be stable. The temporal oscillation of the seasonal movement can be correlated with precipitation, implying that seasonal movement here is hydrology-driven. The seasonal movement also has a frequency similar to GPS-derived regional ground oscillations due to mass loading by stored rainfall and subsequent rebound but with much smaller magnitude, suggesting different hydrological loading effects. From the time-series amplitude information on terrain upslope of the headscarp, we also re-evaluate the incipient motion related to the 2008 Greenleaf Basin rock avalanche, not previously recognized by traditional SAR/InSAR methods. The approach used in this study can be used to identify active landslides in forested terrain, to track the seasonal movement of landslides, and to identify previously unknown landslide hazards.

Washington↗
USGS2017

Multi-temporal LiDAR and Landsat quantification of fire-induced changes to forest structure

Measuring post-fire effects at landscape scales is critical to an ecological understanding of wildfire effects. Predominantly this is accomplished with either multi-spectral remote sensing data or through ground-based field sampling plots. While these methods are important, field data is usually limited to opportunistic post-fire observations, and spectral data often lacks validation with specific variables of change. Additional uncertainty remains regarding how best to account for environmental variables influencing fire effects (e.g., weather) for which observational data cannot easily be acquired, and whether pre-fire agents of change such as bark beetle and timber harvest impact model accuracy. This study quantifies wildfire effects by correlating changes in forest structure derived from multi-temporal Light Detection and Ranging (LiDAR) acquisitions to multi-temporal spectral changes captured by the Landsat Thematic Mapper and Operational Land Imager for the 2012 Pole Creek Fire in central Oregon. Spatial regression modeling was assessed as a methodology to account for spatial autocorrelation, and model consistency was quantified across areas impacted by pre-fire mountain pine beetle and timber harvest. The strongest relationship (pseudo-r 2 = 0.86, p < 0.0001) was observed between the ratio of shortwave infrared and near infrared reflectance (d74) and LiDAR-derived estimate of canopy cover change. Relationships between percentage of LiDAR returns in forest strata and spectral indices generally increased in strength with strata height. Structural measurements made closer to the ground were not well correlated. The spatial regression approach improved all relationships, demonstrating its utility, but model performance declined across pre-fire agents of change, suggesting that such studies should stratify by pre-fire forest condition. This study establishes that spectral indices such as d74 and dNBR are most sensitive to wildfire-caused structural changes such as reduction in canopy cover and perform best when that structure has not been reduced pre-fire.

Oregon↗
USGS2017

Thematic accuracy assessment of the 2011 National Land Cover Database (NLCD)

Accuracy assessment is a standard protocol of National Land Cover Database (NLCD) mapping. Here we report agreement statistics between map and reference labels for NLCD 2011, which includes land cover for ca. 2001, ca. 2006, and ca. 2011. The two main objectives were assessment of agreement between map and reference labels for the three, single-date NLCD land cover products at Level II and Level I of the classification hierarchy, and agreement for 17 land cover change reporting themes based on Level I classes (e.g., forest loss; forest gain; forest, no change) for three change periods (2001–2006, 2006–2011, and 2001–2011). The single-date overall accuracies were 82%, 83%, and 83% at Level II and 88%, 89%, and 89% at Level I for 2011, 2006, and 2001, respectively. Many class-specific user's accuracies met or exceeded a previously established nominal accuracy benchmark of 85%. Overall accuracies for 2006 and 2001 land cover components of NLCD 2011 were approximately 4% higher (at Level II and Level I) than the overall accuracies for the same components of NLCD 2006. The high Level I overall, user's, and producer's accuracies for the single-date eras in NLCD 2011 did not translate into high class-specific user's and producer's accuracies for many of the 17 change reporting themes. User's accuracies were high for the no change reporting themes, commonly exceeding 85%, but were typically much lower for the reporting themes that represented change. Only forest loss, forest gain, and urban gain had user's accuracies that exceeded 70%. Lower user's accuracies for the other change reporting themes may be attributable to the difficulty in determining the context of grass (e.g., open urban, grassland, agriculture) and between the components of the forest-shrubland-grassland gradient at either the mapping phase, reference label assignment phase, or both. NLCD 2011 user's accuracies for forest loss, forest gain, and urban gain compare favorably with results from other land cover change accuracy assessments.

Remote Sensing of Environment↗
USGS2017

A land cover change detection and classification protocol for updating Alaska NLCD 2001 to 2011

Monitoring and mapping land cover changes are important ways to support evaluation of the status and transition of ecosystems. The Alaska National Land Cover Database (NLCD) 2001 was the first 30-m resolution baseline land cover product of the entire state derived from circa 2001 Landsat imagery and geospatial ancillary data. We developed a comprehensive approach named AKUP11 to update Alaska NLCD from 2001 to 2011 and provide a 10-year cyclical update of the state's land cover and land cover changes. Our method is designed to characterize the main land cover changes associated with different drivers, including the conversion of forests to shrub and grassland primarily as a result of wildland fire and forest harvest, the vegetation successional processes after disturbance, and changes of surface water extent and glacier ice/snow associated with weather and climate changes. For natural vegetated areas, a component named AKUP11-VEG was developed for updating the land cover that involves four major steps: 1) identify the disturbed and successional areas using Landsat images and ancillary datasets; 2) update the land cover status for these areas using a SKILL model (System of Knowledge-based Integrated-trajectory Land cover Labeling); 3) perform decision tree classification; and 4) develop a final land cover and land cover change product through the postprocessing modeling. For water and ice/snow areas, another component named AKUP11-WIS was developed for initial land cover change detection, removal of the terrain shadow effects, and exclusion of ephemeral snow changes using a 3-year MODIS snow extent dataset from 2010 to 2012. The overall approach was tested in three pilot study areas in Alaska, with each area consisting of four Landsat image footprints. The results from the pilot study show that the overall accuracy in detecting change and no-change is 90% and the overall accuracy of the updated land cover label for 2011 is 86%. The method provided a robust, consistent, and efficient means for capturing major disturbance events and updating land cover for Alaska. The method has subsequently been applied to generate the land cover and land cover change products for the entire state of Alaska.

Alaska↗
USGS2017

Satellite-based water use dynamics using historical Landsat data (1984–2014) in the southwestern United States

Remote sensing-based field-scale evapotranspiration (ET) maps are useful for characterizing water use patterns and assessing crop performance. The relative impact of climate variability and water management decisions are better studied and quantified using historical data that are derived using a set of consistent datasets and methodology. Historical (1984–2014) Landsat-based ET maps were generated for major irrigation districts in California, i.e., Palo Verde and eight other sub-basins in parts of the middle and lower Central Valley. A total of 3396 Landsat images were processed using the Operational Simplified Surface Energy Balance (SSEBop) model that integrates weather and remotely sensed images to estimate monthly and annual ET within the study sites over the 31 years. Model output evaluation and validation using gridded-flux data and water balance ET approaches indicated relatively good correspondence (R 2 up to 0.88, root mean square error as low as 14 mm/month) between SSEBop ET and validation datasets. In a pairwise comparison, annual variability of agro-hydrologic parameters of actual evapotranspiration ( ET a ), land surface temperature ( T s ), and runoff ( Q ) were found to be more variable than their corresponding climatic counterparts of atmospheric water demand ( ET o ), air temperature ( T a ), and precipitation ( P ), revealing process differences between regional climatic drivers and localized agro-hydrologic responses. However, only T a showed a consistent increase (up to 1.2 K) over study sites during the 31 years, whereas other climate variables such as ET o and P showed a generally neutral trend. This study demonstrates a useful application of “Big Data” science where large volumes of historical Landsat and weather datasets were used to quantify and understand the relative importance of water management and climate variability in crop water use dynamics in regards to the linkages among water management decisions, hydrologic processes and economic transactions. Irrigation district-wide ET a estimates were used to compute historical crop water use volumes and monetary equivalents of water savings for the Palo Verde Irrigation District (PVID). During the peak crop fallowing year in PVID, the water saved reached a maximum of ~ 107,200 acre-feet in 2011 with an estimated monetary payout value of $20.5 million. A significant decreasing trend in actual ET despite an increasing atmospheric demand in PVID highlights the role of management decisions in affecting local hydrologic processes. This study has importance for planning water resource allocation, managing water rights, sustaining agricultural production, and quantifying impacts of climate and land use/land cover changes on water resources. With increased computational efficiency, similar studies can be conducted in other parts of the world to help policy and decision makers understand and quantify various aspects of water resources management.

Arizona, California↗
USGS2017

Validation of the USGS Landsat Burned Area Essential Climate Variable (BAECV) across the conterminous United States

The Landsat Burned Area Essential Climate Variable (BAECV), developed by the U.S. Geological Survey (USGS), capitalizes on the long temporal availability of Landsat imagery to identify burned areas across the conterminous United States (CONUS) (1984–2015). Adequate validation of such products is critical for their proper usage and interpretation. Validation of coarse-resolution products often relies on independent data derived from moderate-resolution sensors (e.g., Landsat). Validation of Landsat products, in turn, is challenging because there is no corresponding source of high-resolution, multispectral imagery that has been systematically collected in space and time over the entire temporal extent of the Landsat archive. Because of this, comparison between high-resolution images and Landsat science products can help increase user's confidence in the Landsat science products, but may not, alone, be adequate. In this paper, we demonstrate an approach to systematically validate the Landsat-derived BAECV product. Burned area extent was mapped for Landsat image pairs using a manually trained semi-automated algorithm that was manually edited across 28 path/rows and five different years (1988, 1993, 1998, 2003, 2008). Three datasets were independently developed by three analysts and the datasets were integrated on a pixel by pixel basis in which at least one to all three analysts were required to agree a pixel was burned. We found that errors within our Landsat reference dataset could be minimized by using the rendition of the dataset in which pixels were mapped as burned if at least two of the three analysts agreed. BAECV errors of omission and commission for the detection of burned pixels averaged 42% and 33%, respectively for CONUS across all five validation years. Errors of omission and commission were lowest across the western CONUS, for example in the shrub and scrublands of the Arid West (31% and 24%, respectively), and highest in the grasslands and agricultural lands of the Great Plains in central CONUS (62% and 57%, respectively). The BAECV product detected most (> 65%) fire events > 10 ha across the western CONUS (Arid and Mountain West ecoregions). Our approach and results demonstrate that a thorough validation of Landsat science products can be completed with independent Landsat-derived reference data, but could be strengthened by the use of complementary sources of high-resolution data.

Remote Sensing of Environment↗
USGS2017

Mapping burned areas using dense time-series of Landsat data

Complete and accurate burned area data are needed to document patterns of fires, to quantify relationships between the patterns and drivers of fire occurrence, and to assess the impacts of fires on human and natural systems. Unfortunately, in many areas existing fire occurrence datasets are known to be incomplete. Consequently, the need to systematically collect burned area information has been recognized by the United Nations Framework Convention on Climate Change and the Intergovernmental Panel on Climate Change, which have both called for the production of essential climate variables (ECVs), including information about burned area. In this paper, we present an algorithm that identifies burned areas in dense time-series of Landsat data to produce the Landsat Burned Area Essential Climate Variable (BAECV) products. The algorithm uses gradient boosted regression models to generate burn probability surfaces using band values and spectral indices from individual Landsat scenes, lagged reference conditions, and change metrics between the scene and reference predictors. Burn classifications are generated from the burn probability surfaces using pixel-level thresholding in combination with a region growing process. The algorithm can be applied anywhere Landsat and training data are available. For this study, BAECV products were generated for the conterminous United States from 1984 through 2015. These products consist of pixel-level burn probabilities for each Landsat scene, in addition to, annual composites including: the maximum burn probability and a burn classification. We compared the BAECV burn classification products to the existing Global Fire Emissions Database (GFED; 1997–2015) and Monitoring Trends in Burn Severity (MTBS; 1984–2013) data. We found that the BAECV products mapped 36% more burned area than the GFED and 116% more burned area than MTBS. Differences between the BAECV products and the GFED were especially high in the West and East where the BAECV products mapped 32% and 88% more burned area, respectively. However, the BAECV products found less burned area than the GFED in regions with frequent agricultural fires. Compared to the MTBS data, the BAECV products identified 31% more burned area in the West, 312% more in the Great Plains, and 233% more in the East. Most pixels in the MTBS data were detected by the BAECV, regardless of burn severity. The BAECV products document patterns of fire similar to those in the GFED but also showed patterns of fire that are not well characterized by the existing MTBS data. We anticipate the BAECV products will be useful to studies that seek to understand past patterns of fire occurrence, the drivers that created them, and the impacts fires have on natural and human systems.

continental United States↗
USGS2019

MODIS phenology-derived, multi-year distribution of conterminous U.S. crop types

Innovative, open, and rapid methods to map crop types over large areas are needed for long-term cropland monitoring. We developed two novel and automated decision tree classification approaches to map crop types across the conterminous United States (U.S.) using MODIS 250 m resolution data: 1) generalized, and 2) year-specific classification. The classification approaches use similarities and dissimilarities in crop type phenologyderived from NDVI time-series data for the two approaches. The year-specific approach uses the training samples from one year and classifies crop types for that year only, whereas the generalized classification approach uses above-average, average, and below-average precipitation years for training to produce crop type maps for one or multiple years more robustly. We produced annual crop type maps using the generalized classification approach for 2001–2014 and the year-specific approach for 2008, 2010, 2011 and 2012. The year-specific classification had overall accuracies > 78%, while the generalized classifier had accuracies > 75% for the conterminous U.S. for 2008, 2010, 2011, and 2012. The generalized classifier enables automated and routine crop type mapping without repeated and expensive ground sample collection year after year. The resulting crop type maps for years prior to 2007 are new and especially important for long-term cropland monitoring and food security analysis because no other map products are currently available for 2001–2007.

Remote Sensing of Environment↗
USGS2018

Mapping forest change using stacked generalization: An ensemble approach

The ever-increasing volume and accessibility of remote sensing data has spawned many alternative approaches for mapping important environmental features and processes. For example, there are several viable but highly varied strategies for using time series of Landsat imagery to detect changes in forest cover. Performance among algorithms varies across complex natural systems, and it is reasonable to ask if aggregating the strengths of an ensemble of classifiers might result in increased overall accuracy. Relatively simple rules have been used in the past to aggregate classifications among remotely sensed maps (e.g. using majority predictions), and in other fields, empirical models have been used to create situationally specific algorithm weights. The latter process, called “stacked generalization” (or “stacking”), typically uses a parametric model for the fusion of algorithm outputs. We tested the performance of several leading forest disturbance detection algorithms against ensembles of the outputs of those same algorithms based upon stacking using both parametric and Random Forests-based fusion rules. Stacking using a Random Forests model cut omission and commission error rates in half in many cases in relation to individual change detection algorithms, and cut error rates by one quarter compared to more conventional parametric stacking. Stacking also offers two auxiliary benefits: alignment of outputs to the precise definitions built into a particular set of empirical calibration data; and, outputs which may be adjusted such that map class totals match independent estimates of change in each year. In general, ensemble predictions improve when new inputs are added that are both informative and uncorrelated with existing ensemble components. As increased use of cloud-based computing makes ensemble mapping methods more accessible, the most useful new algorithms may be those that specialize in providing spectral, temporal, or thematic information not already available through members of existing ensembles.

Remote Sensing of Environment↗
USGS2018

Detecting geothermal anomalies and evaluating LST geothermal component by combining thermal remote sensing time series and land surface model data

This paper explores for the first time the possibilities to use two land surface temperature (LST) time series of different origins (geostationary Meteosat Second Generation satellite data and Noah land surface modelling, LSM), to detect geothermal anomalies and extract the geothermal component of LST, the LST gt . We hypothesize that in geothermal areas the LSM time series will underestimate the LST as compared to the remote sensing data, since the former does not account for the geothermal component in its model. In order to extract LST gt , two approaches of different nature (physical based and data mining) were developed and tested in an area of about 560 × 560 km 2 centered at the Kenyan Rift. Pre-dawn data in the study area during the first 45 days of 2012 were analyzed. The results show consistent spatial and temporal LST gt patterns between the two approaches, and systematic differences of about 2 K. A geothermal area map from surface studies was used to assess LST gt inside and outside the geothermal boundaries. Spatial means were found to be higher inside the geothermal limits, as well as the relative frequency of occurrence of high LST gt . Results further show that areas with strong topography can result in anomalously high LST gt values (false positives), which suggests the need for a slope and aspect correction in the inputs to achieve realistic results in those areas. The uncertainty analysis indicates that large uncertainties of the input parameters may limit detection of LST gt anomalies. To validate the approaches, higher spatial resolution images from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) data over the Olkaria geothermal field were used. An established method to estimate radiant geothermal flux was applied providing values between 9 and 24 W/m 2 in the geothermal area, which coincides with the LST gt flux rates obtained with the proposed approaches. The proposed approaches are a first step in estimating LST gt at large spatial coverage from remote sensing and LSM data series, and provide an innovative framework for future improvements.

Remote Sensing of Environment↗
USGS2018

Analyzing floating and bedfast lake ice regimes across Arctic Alaska using 25 years of space-borne SAR imagery

Late-winter lake ice regimes are controlled by water depth relative to maximum ice thickness (MIT). When MIT exceeds maximum water depth, lakes freeze to the bottom with bedfast ice (BI) and when MIT is less than maximum water depth lakes have floating ice (FI). Both airborne radar and space-borne synthetic aperture radar (SAR) imagery (Ku-, X-, C-, and L-band) have been used previously to determine whether lakes have a BI or FI regime in a given year, across a number of years, or across large regions. In this study, we use a combination of ERS-1/2, RADARSAT-2, Envisat, and Sentinel-1 SAR imagery for seven lake-rich regions in Arctic Alaska to analyze lake ice regime extents and dynamics over a 25-year period (1992–2016). Our interactive threshold classification method determines a unique statistic-based intensity threshold for each SAR scene, allowing for the comparison of classification results from C-band SAR data acquired with different polarizations and incidence angles. Additionally, our novel method accommodates declining signal strength in aging extended-mission satellite SAR instruments. Comparison of SAR ice regime classifications with extensive field measurements from six years yielded a 93% accuracy. Significant declines in BI regimes were only observed in the Fish Creek area with 3% of lakes exhibiting transitional ice regimes—lakes that switch from BI to FI during this 25-year period. This analysis suggests that the potential conversion from BI to FI regimes is primarily a function of lake depth distributions in addition to regional differences in climate variability . Remote sensing of lake ice regimes with C-band SAR is a useful tool to monitor the associated thermal impacts on permafrost, since lake ice regimes can be used as a proxy for of sub-lake permafrost thaw, considered by the Global Climate Observing System as an Essential Climate Variable (ECV). Continued winter warming and variable snow conditions in the Arctic are expected and our long-term analysis provides a valuable baseline for predicting where potential future lake ice regimes shifts will be most pronounced.

Alaska↗
USGS2018

Optimizing a remote sensing production efficiency model for macro-scale GPP and yield estimation in agroecosystems

Earth observation data are increasingly used to provide consistent eco-physiological information over large areas through time. Production efficiency models (PEMs) estimate Gross Primary Production (GPP) as a function of the fraction of photosynthetically active radiation absorbed by the canopy, which is derived from Earth observation. GPP can be summed over the growing season and adjusted by a crop-specific harvest index to estimate yield. Although PEMs have many advantages over other crop yield models, they are not widely used, because performance is relatively poor. Here, a new PEM is presented that addresses deficiencies for macro-scale application: Production Efficiency Model Optimized for Crops (PEMOC). It was developed by optimizing functions from the literature with GPP estimated by eddy covariance flux towers in the United States. The model was evaluated using newly developed Earth observation products and county-level yield statistics for major crops. PEMOC generally performed better at the field and county level than another commonly used PEM, the Moderate Resolution Imaging Spectroradiometer GPP (MOD17). PEMOC and MOD17 estimates of GPP had an R 2 and root mean squared error (RMSE) over the growing season of 0.71–0.89 (9.87–17.47 g CO 2 d −1 ) and 0.59–0.83 (6.86–22.20 g CO 2 d −1 ) with flux tower GPP. PEMOC produced R 2 s and RMSE of 0.70 (0.52), 0.60 (0.61), and 0.62 (0.59), while MOD17 produced R 2 s and RMSE of 0.65 (0.57), 0.53 (0.66), and 0.65 (0.57) with corn, soybean, and winter wheat crop yield anomalies. The sample size of rice was small, so yields were compared directly. PEMOC and MOD17 produced R 2 s and RMSE of 0.53 (3.42 t ha −1 ) and 0.40 (4.89 t ha −1 ). The most sizeable model improvements were seen for C 3 and C 4 crops during emergence/senescence and peak season, respectively. These improvements were attributed to C 3 and C 4 partitioning, optimized temperature and moisture constraints, and an evapotranspiration-based soil moisture index.

Remote Sensing of Environment↗
USGS2018

Induced earthquake and liquefaction hazards in Oklahoma, USA: Constraints from InSAR

Oklahoma experienced three earthquakes of M w 5.0 or greater in 2016: the 13-Feb. Fairview earthquake (M w 5.1), the 03-Sep. Pawnee earthquake (M w 5.8), and the 07-Nov. Cushing earthquake (M w 5.0). These events are the first earthquakes in the state exceeding M w 5.0 since the 2011 M w 5.7 Prague earthquake and likely result from wide-scale deep fluid-injection. We use interferometric synthetic aperture radar (InSAR) observations to quantify the magnitude and location of surface deformation associated with these three events, determine the depth ranges of fault slip, and assess the spatial relationship between fault slip and well-calibrated mainshock and aftershock locations. We also include newly reported, calibrated event locations for the Cushing earthquake. We find that the Pawnee earthquake ruptured within the crystalline basement with the shallowest slip occurring at depths of 3.1–4.3 km. We find a similar, though shallower, crystalline basement source for the Cushing earthquake with a minimum depth to slip of 1.6–2.3 km. Despite the smaller magnitude of the Cushing earthquake, it generated anomalously high ground motions and damage compared to the larger Pawnee and Fairview earthquakes. We postulate that the shallow source of the Cushing earthquakes provides one explanation for the higher than expected ground motions. The Fairview earthquake generated no detectable co-seismic displacements, which is consistent with a relatively deep earthquake source (~8.5 km). We do, however, identify a 16 km stretch of floodplain where widespread liquefaction occurred in response to the Fairview earthquake, and where 30 gas production wells were exposed to surface displacements exceeding 5 cm. Consequently, the depth to crystalline basement, which limits the depth of injection-induced earthquakes in Oklahoma, and the potential for liquefaction are important factors in assessing shaking risk in the central United States.

Oklahoma↗
USGS2019

Historical background and current developments for mapping burned area from satellite Earth observation

Fire has a diverse range of impacts on Earth's physical and social systems. Accurate and up to date information on areas affected by fire is critical to better understand drivers of fire activity, as well as its relevance for biogeochemical cycles , climate, air quality, and to aid fire management. Mapping burned areas was traditionally done from field sketches. With the launch of the first Earth observation satellites, remote sensing quickly became a more practical alternative to detect burned areas, as they provide timely regional and global coverage of fire occurrence. This review paper explores the physical basis to detect burned area from satellite observations, describes the historical trends of using satellite sensors to monitor burned areas, summarizes the most recent approaches to map burned areas and evaluates the existing burned area products (both at global and regional scales). Finally, it identifies potential future opportunities to further improve burned area detection from Earth observation satellites.

Remote Sensing of Environment↗