Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Algorithms”

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 1,423 records · Page 79Linked to original sources

Potential underestimation of satellite fire radiative power retrievals over gas flares and wildland fires

Fire Radiative Power (FRP) is related to fire combustion rates and is used to quantify the atmospheric emissions of greenhouse gases and aerosols. FRP over gas flares and wildfires can be retrieved remotely using satellites that observe in shortwave infrared (SWIR) to middle infrared (MIR) wavelengths. Heritage techniques to retrieve FRP developed for wildland fires using the MIR 4 μm radiances have been adapted for the hotter burning gas flares using the SWIR 2 μm observations. Effects of atmosphere, including smoke and aerosols, are assumed to be minimal in these algorithms because of the use of longer than visual wavelengths. Here we use Moderate Resolution Imaging Spectroradiometer (MODIS), Visible Infrared Imaging Radiometer Suite (VIIRS) and Landsat 8 observations acquired before and during emergency oil and gas flaring in eastern Saudi Arabia to show that dark, sooty smoke affects both 4 μm and 2 μm observations. While the 2 μm observations used to retrieve gas FRP may be reliable during clear atmospheric conditions, performance is severely impacted by dark smoke. Global remote sensing-based inventories of wildfire and gas flaring need to consider the possibility that soot and dark smoke can potentially lead to an underestimation of FRP over fires.

Remote Sensing↗

Discrimination of biological scatterers in polarimetric weather radar data: Opportunities and challenges

For radar aeroecology studies, the identification of the type of scatterer is critically important. Here, we used a random forest (RF) algorithm to develop a variety of scatterer classification models based on the backscatter values in radar resolution volumes of six radar variables (reflectivity, radial velocity, spectrum width, differential reflectivity, correlation coefficient, and differential phase) from seven types of biological scatterers and one type of meteorological scatterer (rain). Models that discriminated among fewer classes and/or aggregated similar types into more inclusive classes classified with greater accuracy and higher probability. Bioscatterers that shared similarities in phenotype tended to misclassify against one another more frequently than against more dissimilar types, with the greatest degree of misclassification occurring among vertebrates. Polarimetric variables proved critical to classification performance and individual polarimetric variables played central roles in the discrimination of specific scatterers. Not surprisingly, purposely overfit RF models (in one case study) were our highest performing. Such models have a role to play in situations where the inclusion of natural history can play an outsized role in model performance. In the future, bioscatter classification will become more nuanced, pushing machine-learning model development to increasingly rely on independent validation of scatterer types and more precise knowledge of the physical and behavioral properties of the scatterer.

Remote Sensing↗

Isolating anthropogenic wetland loss by concurrently tracking inundation and land cover disturbance across the Mid-Atlantic Region, U.S.

Global trends in wetland degradation and loss have created an urgency to monitor wetland extent, as well as track the distribution and causes of wetland loss. Satellite imagery can be used to monitor wetlands over time, but few efforts have attempted to distinguish anthropogenic wetland loss from climate-driven variability in wetland extent. We present an approach to concurrently track land cover disturbance and inundation extent across the Mid-Atlantic region, United States, using the Landsat archive in Google Earth Engine. Disturbance was identified as a change in greenness, using a harmonic linear regression approach, or as a change in growing season brightness. Inundation extent was mapped using a modified version of the U.S. Geological Survey’s Dynamic Surface Water Extent (DSWE) algorithm. Annual (2015–2018) disturbance averaged 0.32% (1095 km 2 year -1 ) of the study area per year and was most common in forested areas. While inundation extent showed substantial interannual variability, the co-occurrence of disturbance and declines in inundation extent represented a minority of both change types, totaling 109 km 2 over the four-year period, and 186 km 2 , using the National Wetland Inventory dataset in place of the Landsat-derived inundation extent. When the annual products were evaluated with permitted wetland and stream fill points, 95% of the fill points were detected, with most found by the disturbance product (89%) and fewer found by the inundation decline product (25%). The results suggest that mapping inundation alone is unlikely to be adequate to find and track anthropogenic wetland loss. Alternatively, remotely tracking both disturbance and inundation can potentially focus efforts to protect, manage, and restore wetlands.

Delaware, MarylandPennsylvania, Virginia, West Vir↗

QCam: sUAS-based doppler radar for measuring river discharge

The U.S. Geological Survey is actively investigating remote sensing of surface velocity and river discharge (discharge) from satellite-, high altitude-, small, unmanned aircraft systems- (sUAS or drone), and permanent (fixed) deployments. This initiative is important in ungaged basins and river reaches that lack the infrastructure to deploy conventional streamgaging equipment. By coupling alternative discharge algorithms with sensors capable of measuring surface velocity, streamgage networks can be established in regions where data collection was previously impractical or impossible. To differentiate from satellite or high-altitude platforms, near-field remote sensing is conducted from sUAS or fixed platforms. QCam is a Doppler (velocity) radar mounted and integrated on a 3DR © Solo sUAS. It measures the along-track surface velocity by spot dwelling in a river cross section at a vertical where the maximum surface velocity is recorded. The surface velocity is translated to a mean-channel (mean) velocity using the probability concept (PC), and discharge is computed using the PC-derived mean velocity and cross-sectional area. Factors including surface-scatterer quality, flight altitude, propwash, wind drift, and sample duration may affect the radar-returns and the subsequent computation of mean velocity and river discharge. To evaluate the extensibility of the method, five science flights were conducted on four rivers of varying size and dynamics and included the Arkansas River, Colorado (CO), USA (two events); Salcha River near Salchaket, Alaska (AK), USA; South Platte River, CO, USA; and the Tanana River, AK, USA. QCam surface velocities and river discharges were compared to conventional streamgaging methods, which represented truth. QCam surface velocities for the Arkansas River, Salcha River, South Platte River, and Tanana River were 1.02 meters per second (m/s) and 1.43 m/s; 1.58 m/s; 0.90 m/s; and 2.17 m/s, respectively. QCam discharges (and percent differences) were 9.48 (0.3%) and 20.3 cubic meters per second (m 3 /s) (2.5%); 62.1 m 3 /s (−10.4%); 3.42 m 3 /s (7.3%), and 1579 m 3 /s (−18.8%). QCam results compare favorably with conventional streamgaging and are a viable near-field remote sensing technology that can be operationalized to deliver real-time surface velocity, mean velocity, and river discharge, if cross-sectional area is available.

Remote Sensing↗

The new Landsat Collection-2 Digital Elevation Model

The Landsat Collection-2 distribution introduces a new global Digital Elevation Model (DEM) for scene orthorectification. The new global DEM is a composite of the latest and most accurate freely available DEM sources and will include reprocessed Shuttle Radar Topographic Mission (SRTM) data (called NASADEM), high-resolution stereo optical data (ArcticDEM), a new National Elevation Dataset (NED) and various publicly available national datasets including the Canadian Digital Elevation Model (CDEM) and DEMs for Sweden, Norway and Finland (SNF). The new DEM will be available world-wide with few exceptions. It is anticipated that the transition from the Collection-1 DEM at 3 arcsecond to the new DEM will be seamless because processing methods to maintain a seamless transition were employed, void filling techniques were used, where persistent gaps were found, and the pixel spacing is the same between the two collections. Improvements to the vertical accuracy were realized by differencing accuracies of other elevation datasets to the new DEM. The greatest improvement occurred where ArcticDEM data were used, where an improvement of 35 m was measured. By using theses improved vertical values in a line of sight algorithm, horizontal improvements were noted in some of the most mountainous regions over multiple 30-m Landsat pixels. This new DEM will be used to process all of the scenes from Landsat 1-8 in Collection-2 processing and will be made available to the public by the end of 2020.

Remote Sensing↗

Geometric calibration updates to Landsat 7 ETM+ instrument for Landsat Collection 2 products

The Landsat 7 (L7) spacecraft and its instrument, the enhanced thematic mapper plus (ETM+), have been consistently characterized and calibrated since its launch in April of 1999. These performance metrics and calibration updates are determined through the U. S. Geological Survey (USGS) Landsat image assessment system (IAS), which has been performing this function since launch. Starting in 2016, the USGS adopted a tiered collection management structure for its Landsat data products that ensures a consistent method of processing for the Landsat archive within a given collection while allowing a set of calibration updates to be performed between any two given collections. The time frame between 2016 and the end of 2020 was part of the Landsat data Collection 1, in the middle of 2020 was the start of the Landsat Collection 2 data products. The start of a given collection initiates the reprocessing of the Landsat archive, which may involve one or more of a set of updated calibration parameters, improvements in the support data needed for product generation, and improved algorithms used in both the processing flow of products along with the characterization and calibration of the Landsat instruments and spacecraft. This paper discusses only the ETM+ geometric spacecraft and instrument calibration improvements for Collection 2. Three ETM+ calibration updates were made for the ETM+; updates to the thermal band odd-to-even detector alignment, sensor to attitude control system (ACS) alignment, and a cold-to-warm focal plane alignment adjustment. The sensor alignment updates impact only the accuracy of the systematic terrain products (L1GT), which are the products generated before applying any corrections based on the ground control used in registration. The band alignment changes impacted only bands 5, 6, and 7 within the focal plane. Other geometric calibration updates, such as scan mirror alignment, are done on a routine basis and are not part of the Collection 2 updates due to their more dynamic characteristics.

Remote Sensing↗

Mapping benthic algae and cyanobacteria in river channels from aerial photographs and satellite images: A proof-of-concept investigation on the Buffalo National River, AR, USA

Although rivers are of immense practical, aesthetic, and recreational value, these aquatic habitats are particularly sensitive to environmental changes. Increasingly, changes in streamflow and water quality are resulting in blooms of bottom-attached (benthic) algae, also known as periphyton, which have become widespread in many water bodies of US national parks. Because these blooms degrade visitor experiences and threaten human and ecosystem health, improved methods of characterizing benthic algae are needed. This study evaluated the potential utility of remote sensing techniques for mapping variations in algal density in shallow, clear-flowing rivers. As part of an initial proof-of-concept investigation, field measurements of water depth and percent cover of benthic algae were collected from two reaches of the Buffalo National River along with aerial photographs and multispectral satellite images. Applying a band ratio algorithm to these data yielded reliable depth estimates, although a shallow bias and moderate level of precision were observed. Spectral distinctions among algal percent cover values ranging from 0 to 100% were subtle and became only slightly more pronounced when the data were aggregated to four ordinal levels. A bagged trees machine learning model trained using the original spectral bands and image-derived depth estimates as predictor variables was used to produce classified maps of algal density. The spatial and temporal patterns depicted in these maps were reasonable but overall classification accuracies were modest, up to 64.6%, due to a lack of spectral detail. To further advance remote sensing of benthic algae and other periphyton, future studies could adopt hyperspectral approaches and more quantitative, continuous metrics such as biomass.

Arkansas↗

Absolute accuracy assessment of lidar point cloud using amorphous objects

The accuracy assessment of airborne lidar point cloud typically estimates vertical accuracy by computing RMSEz (root mean square error of the z coordinate) from ground check points (GCPs). Due to the low point density of the airborne lidar point cloud, there is often not enough accurate semantic context to find an accurate conjugate point. To advance the accuracy assessment in full three-dimensional (3D) context, geometric features, such as the three-plane intersection point or two-line intersection point, are often used. Although the point density is still low, geometric features are mathematically modeled from many points. Thus, geometric features provide a robust determination of the intersection point, and the point is considered as a GCP. When no regular built objects are available, we describe the process of utilizing features of irregular shape called amorphous natural objects, such as a tree or a rock. When scanned to a high-density point cloud, an amorphous natural object can be used as ground truth reference data to estimate 3D georeferencing errors of the airborne lidar point cloud. The algorithm to estimate 3D accuracy is the optimization that minimizes the sum of the distance between the airborne lidar points to the ground scanned data. The search volume partitioning was the most important procedure to improve the computational efficiency. We also performed an extensive study to address the external uncertainty associated with the amorphous object method. We describe an accuracy assessment using amorphous objects (108 trees) spread over the project area. The accuracy results for ∆ x , ∆ y , and ∆ z obtained using the amorphous object method were 3.1 cm, 3.6 cm, and 1.7 cm RMSE, along with a mean error of 0.1 cm, 0.1 cm, and 4.5 cm, respectively, satisfying the accuracy requirement of U.S. Geological Survey lidar base specification. This approach shows strong promise as an alternative to geometric feature methods when artificial targets are scarce. The relative convenience and advantages of using amorphous targets, along with its good performance shown here, make this amorphous object method a practical way to perform 3D accuracy assessment.

Remote Sensing↗

Improving the operational simplified surface energy balance evapotranspiration model using the forcing and normalizing operation

Actual evapotranspiration modeling is providing useful information for researchers and resource managers in agriculture and water resources around the world. The performance of models depends on the accuracy of forcing inputs and model parameters. We developed an improved approach to the parameterization of the Operational Simplified Surface Energy Balance (SSEBop) model using the Forcing and Normalizing Operation (FANO). SSEBop has two key model parameters that define the model boundary conditions. The FANO algorithm computes the wet-bulb boundary condition using a linear FANO Equation relating surface temperature, surface psychrometric constant, and the Normalized Difference Vegetation Index (NDVI). The FANO parameterization was implemented on two computing platforms using Landsat and gridded meteorological datasets: (1) Google Earth Engine (GEE) and (2) Earth Resources Observation and Science (EROS) Center Science Processing Architecture (ESPA). Evaluation was conducted by comparing modeled actual evapotranspiration ( ETa ) estimates with AmeriFlux eddy covariance (EC) and water balance ETa from level-8 Hydrologic Unit Code sub-basins in the conterminous United States. FANO brought substantial improvements in model accuracy and operational implementation. Compared to the earlier version (v0.1.7), SSEBop FANO (v0.2.6) reduced grassland bias from 47% to −2% while maintaining comparable bias for croplands (11% versus −7%) against EC data. A water balance-based ETa bias evaluation showed an overall improvement from 7% to −1%. Climatology versus annual gridded reference evapotranspiration ( ETr ) produced comparable ETa results, justifying the use of climatology ETr for the global SSEBop Landsat ETa that is accessible through the ESPA website. Besides improvements in model accuracy, SSEBop FANO increases the spatiotemporal coverage of ET modeling due to the elimination of high NDVI requirements for model parameterization. Because of the existence of potential biases from forcing inputs and model parameters, continued evaluation and bias corrections are necessary to improve the absolute magnitude of ETa for localized water budget applications.

Remote Sensing↗

SaTSeaD: Satellite Triangulated Sea Depth open-source bathymetry module for NASA Ames Stereo Pipeline

We developed the first-ever bathymetric module for the NASA Ames Stereo Pipeline (ASP) open-source topographic software called Satellite Triangulated Sea Depth, or SaTSeaD, to derive nearshore bathymetry from stereo imagery. Correct bathymetry measurements depend on water surface elevation, and whereas previous methods considered the water surface horizontal, our bathymetric module accounts for the curvature of the Earth in the imagery. The process is semiautomatic, reliable, and repeatable, independent of any external bathymetry data eliminating user bias in selecting bathymetry calibration points, and it can generate a fully integrated and seamless topo-bathymetry digital elevation model (TBDEM) in the same coordinate system, comparable with the band-ratio method irrespective of the regression method used for the band-ratio algorithm. The ASP output can be improved by applying a camera bundle adjustment to minimize reprojection errors and by alignment to a more accurate topographic (above water) surface without any bathymetric input since the derived TBDEM is a rigid surface. These procedures can decrease bathymetry root mean square errors from 30 to 80 percent, depending on environmental conditions, the quality of satellite imagery, and the spectral band used (e.g., blue, green, or panchromatic).

Remote Sensing↗

Annual and inter-annual variability in the diffuse attenuation coefficient and turbidity in an urbanized Washington lake from 2013 to 2022 assessed using Landsat-8/9

Water clarity, defined in this study using measurements of the downwelling diffuse light attenuation coefficient (K d ) and turbidity, is an important indicator of lake trophic status and ecosystem health. We used in-situ measurements to evaluate existing semi-analytical models for K d and turbidity, developed a regional turbidity model based on spectral shape, and evaluated the spatial and temporal trends in Lake Washington from 2013 to 2022 using Landsat-8/9 Operational Land Imager (OLI). We found no significant trends from 2013 to 2022 in K d or turbidity when both the annual and full datasets were considered. In addition to the spring peak lasting from April through June, autumn K d peaks were present at all sites, a pattern consistent with seasonal chlorophyll a and zooplankton concentrations. There existed no autumn peak in the monthly turbidity dataset, and the spring peak occurred two months before the K d peak, nearly mirroring seasonal variability in the Cedar River discharge rates over the same period. The K d and turbidity algorithms were thus each more sensitive to different sources of water clarity variability in Lake Washington.

Washington↗

A "Region-Specific Model Adaptation (RSMA)" based training data method in large-scale land cover mapping

An accurate and historical land cover monitoring dataset for Alaska could provide fundamental information for a range of studies, such as conservation habitats, biogeochemical cycles, and climate systems, in this distinctive region. This research addresses challenges associated with the extraction of training data for timely and accurate land cover classifications in Alaska over longer time periods (e.g., greater than 10 years). Specifically, we designed the “Region-Specific Model Adaptation (RSMA)” method for training data. The method integrates land cover information from the National Land Cover Database (NLCD), LANDFIRE’s Existing Vegetation Type (EVT), and the National Wetlands Inventory (NWI) and machine learning techniques to generate robust training samples based on the Anderson Level II classification legend. The assumption of the method is that spectral signatures vary across regions because of diverse land surface compositions; however, despite these variations, there are consistent, collective land cover characteristics that span the entire region. Building upon this assumption, this research utilized the classification power of deep learning algorithms and the generalization ability of RSMA to construct a model for the RSMA method. Additionally, we interpreted existing vegetation plot information for land cover labels as validation data to reduce inconsistency in the human interpretation. Our validation results indicate that the RSMA method improved the quality of the training data derived solely from the NLCD by approximately 30% for the overall accuracy. The validation assessment also demonstrates that the RSMA method can generate reliable training data on large scales in regions that lack sufficient reliable data.

Remote Sensing↗

Mapping river flow from thermal images in approximately real time: Proof of concept on the Sacramento River, California, USA

Image velocimetry has become an effective method of mapping flow conditions in rivers, but this analysis is typically performed in a post-processing mode after data collection is complete. In this study, we evaluated the potential to infer flow velocities in approximately real time as thermal images are being acquired from an uncrewed aircraft system (UAS). The sensitivity of thermal image velocimetry to environmental conditions was quantified by conducting 20 flights over four days and assessing the accuracy of image-derived velocity estimates via comparison to direct field measurements made with an acoustic Doppler current profiler (ADCP). This analysis indicated that velocity mapping was most reliable when the air was cooler than the water. We also introduced a workflow for River Velocity Measurement in Approximately Real Time (RiVMART) that involved transferring brief image sequences from the UAS to a ground station as distinct data packets. The resulting velocity fields were as accurate as those generated via post-processing. A new particle image velocimetry (PIV) algorithm based on staggered image sequences increased the number of image pairs available for a given image sequence duration and slightly improved accuracy relative to a standard PIV implementation. Direct, automated geo-referencing of image-derived velocity vectors based on information on the position and orientation of the UAS acquired during flight led to poor alignment with vectors that were geo-referenced manually by selecting ground control points from an orthophoto. This initial proof-of-concept investigation suggests that our workflow could enable highly efficient characterization of flow fields in rivers and might help support applications that require rapid response to changing conditions.

California↗

Understanding the influence of image enhancement on underwater object detection: A quantitative and qualitative study

Underwater image enhancement is often perceived as a disadvantageous process to object detection. We propose a novel analysis of the interactions between enhancement and detection, elaborating on the potential of enhancement to improve detection. In particular, we evaluate object detection performance for each individual image rather than across the entire set to allow a direct performance comparison of each image before and after enhancement. This approach enables the generation of unique queries to identify the outperforming and underperforming enhanced images compared to the original images. To accomplish this, we first produce enhanced image sets of the original images using recent image enhancement models. Each enhanced set is then divided into two groups: (1) images that outperform or match the performance of the original images and (2) images that underperform. Subsequently, we create mixed original-enhanced sets by replacing underperforming enhanced images with their corresponding original images. Next, we conduct a detailed analysis by evaluating all generated groups for quality and detection performance attributes. Finally, we perform an overlap analysis between the generated enhanced sets to identify cases where the enhanced images of different enhancement algorithms unanimously outperform, equally perform, or underperform the original images. Our analysis reveals that, when evaluated individually, most enhanced images achieve equal or superior performance compared to their original counterparts. The proposed method uncovers variations in detection performance that are not apparent in a whole set as opposed to a per-image evaluation because the latter reveals that only a small percentage of enhanced images cause an overall negative impact on detection. We also find that over-enhancement may lead to deteriorated object detection performance. Lastly, we note that enhanced images reveal hidden objects that were not annotated due to the low visibility of the original images.

Remote Sensing↗

Phenological classification of the United States: A geographic framework for extending multi-sensor time-series data

This study introduces a new geographic framework, phenological classification, for the conterminous United States based on Moderate Resolution Imaging Spectroradiometer (MODIS) Normalized Difference Vegetation Index (NDVI) time-series data and a digital elevation model. The resulting pheno-class map is comprised of 40 pheno-classes, each having unique phenological and topographic characteristics. Cross-comparison of the pheno-classes with the 2001 National Land Cover Database indicates that the new map contains additional phenological and climate information. The pheno-class framework may be a suitable basis for the development of an Advanced Very High Resolution Radiometer (AVHRR)-MODIS NDVI translation algorithm and for various biogeographic studies.

Remote Sensing↗

High spatial resolution WorldView-2 imagery for mapping NDVI and its relationship to temporal urban landscape evapotranspiration factors

Evapotranspiration estimation has benefitted from recent advances in remote sensing and GIS techniques particularly in agricultural applications rather than urban environments. This paper explores the relationship between urban vegetation evapotranspiration (ET) and vegetation indices derived from newly-developed high spatial resolution WorldView-2 imagery. The study site was Veale Gardens in Adelaide, Australia. Image processing was applied on five images captured from February 2012 to February 2013 using ERDAS Imagine. From 64 possible two band combinations of WorldView-2, the most reliable one (with the maximum median differences) was selected. Normalized Difference Vegetation Index (NDVI) values were derived for each category of landscape cover, namely trees, shrubs, turf grasses, impervious pavements, and water bodies. Urban landscape evapotranspiration rates for Veale Gardens were estimated through field monitoring using observational-based landscape coefficients. The relationships between remotely sensed NDVIs for the entire Veale Gardens and for individual NDVIs of different vegetation covers were compared with field measured urban landscape evapotranspiration rates. The water stress conditions experienced in January 2013 decreased the correlation between ET and NDVI with the highest relationship of ET-Landscape NDVI (Landscape Normalized Difference Vegetation Index) for shrubs (r 2 = 0.66) and trees (r 2 = 0.63). However, when the January data was excluded, there was a significant correlation between ET and NDVI. The highest correlation for ET-Landscape NDVI was found for the entire Veale Gardens regardless of vegetation type (r 2 = 0.95, p > 0.05) and the lowest one was for turf (r 2 = 0.88, p > 0.05). In support of the feasibility of ET estimation by WV2 over a longer period, an algorithm recently developed that estimates evapotranspiration rates based on the Enhanced Vegetation Index (EVI) from MODIS was employed. The results revealed a significant positive relationship between ETMODIS and ETWV2 (r 2 = 0.9857, p > 0.05). This indicates that the relationship between NDVI using high resolution WorldView-2 imagery and ground-based validation approaches could provide an effective predictive tool for determining ET rates from unstressed mixed urban landscape plantings.

Adelaide↗

An automated approach for mapping persistent ice and snow cover over high latitude regions

We developed an automated approach for mapping persistent ice and snow cover (glaciers and perennial snowfields) from Landsat TM and ETM+ data across a variety of topography, glacier types, and climatic conditions at high latitudes (above ~65°N). Our approach exploits all available Landsat scenes acquired during the late summer (1 August–15 September) over a multi-year period and employs an automated cloud masking algorithm optimized for snow and ice covered mountainous environments. Pixels from individual Landsat scenes were classified as snow/ice covered or snow/ice free based on the Normalized Difference Snow Index (NDSI), and pixels consistently identified as snow/ice covered over a five-year period were classified as persistent ice and snow cover. The same NDSI and ratio of snow/ice-covered days to total days thresholds applied consistently across eight study regions resulted in persistent ice and snow cover maps that agreed closely in most areas with glacier area mapped for the Randolph Glacier Inventory (RGI), with a mean accuracy (agreement with the RGI) of 0.96, a mean precision (user’s accuracy of the snow/ice cover class) of 0.92, a mean recall (producer’s accuracy of the snow/ice cover class) of 0.86, and a mean F-score (a measure that considers both precision and recall) of 0.88. We also compared results from our approach to glacier area mapped from high spatial resolution imagery at four study regions and found similar results. Accuracy was lowest in regions with substantial areas of debris-covered glacier ice, suggesting that manual editing would still be required in these regions to achieve reasonable results. The similarity of our results to those from the RGI as well as glacier area mapped from high spatial resolution imagery suggests it should be possible to apply this approach across large regions to produce updated 30-m resolution maps of persistent ice and snow cover. In the short term, automated PISC maps can be used to rapidly identify areas where substantial changes in glacier area have occurred since the most recent conventional glacier inventories, highlighting areas where updated inventories are most urgently needed. From a longer term perspective, the automated production of PISC maps represents an important step toward fully automated glacier extent monitoring using Landsat or similar sensors.

Remote Sensing↗

Mapping annual forest cover in sub-humid and semi-arid regions through analysis of landsat and PALSAR imagery

Accurately mapping the spatial distribution of forests in sub-humid to semi-arid regions over years is a challenging task and causes difficulty to forest management. Relatively large uncertainties still exist in the spatial distribution of forests and deforestation in the sub-humid and semi-arid regions. Numerous publications have used either optical or synthetic aperture radar (SAR) remote sensing imagery, but the resultant forest cover maps often have large errors. In this study, we proposed a pixel- and rule-based algorithm to identify and map annual forests from 2007 to 2010 in Oklahoma, USA, a transition region with various climate and landscapes, using the integration of the L-band ALOS PALSAR Fine Beam Dual Polarization (FBD) mosaic dataset and Landsat images. The overall accuracy and Kappa coefficient of the PALSAR/Landsat forest map were about 88.2% and 0.75 in 2010, with the user and producer accuracy about 93.4% and 75.7%, based on the 3,270 random ground plots collected in 2012 and 2013. Compared with the forest products from JAXA, NLCD, OKESM and OKFRA, the PALSAR/Landsat forest map showed great improvement. The area of the PALSAR/Landsat forest was about 40,149 km2 in 2010, which was close to the area from OKFRA (40,468 km2), but much larger than those from JAXA (32,403 km2) and NLCD (37,628 km2). We analyzed annual forest cover dynamics, and the results show extensive deforestation (2,761 km2, 6.9% of the total forest area in 2010) and reforestation (3,630 km2, 9.0%) in the southeast and central Oklahoma, and the total area of forests increased by 684 km2 from 2007 to 2010. This study clearly demonstrates the potential of data fusion between PALSAR and Landsat images for mapping annual forest cover dynamics in sub-humid to semi-arid regions, and the resultant forest maps would be helpful to forest management.

Oklahoma↗