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 739 records · Page 41Linked to original sources

Machine learning predictions of mean ages of shallow well samples in the Great Lakes Basin, USA

The travel time or “age” of groundwater affects catchment responses to hydrologic changes , geochemical reactions, and time lags between management actions and responses at down-gradient streams and wells. Use of atmospheric tracers has facilitated the characterization of groundwater ages, but most wells lack such measurements. This study applied machine learning to predict ages in wells across a large region around the Great Lakes Basin using well, chemistry, and landscape characteristics. For a dataset of age tracers in 961 samples, the travel time from the land surface to the sample location was estimated for each sample using parametric functions. The mean travel times were then modeled using a gradient boosting machine (GBM) algorithm with cross validation tuning of model metaparameters. The GBM approach was able to closely match estimated ages for the training data (RMSE = 0.26 natural-log scale years) and provided a reasonable match to testing data (RMSE = 0.84). Of the variables tested, well characteristics (e.g. depth), land use, hydrologic indicators (e.g. topographic wetness index), and water chemistry (e.g. nitrate, fluoride, and pH), substantially affected the predictions of age. GBM prediction was applied to 14,335 groundwater samples with median sample depth of 5.4 m, indicating for the Great Lakes Basin a broad distribution of ages among wells with a median of 32.9 years. Lag times of decades are likely for these wells to respond to changing solute fluxes near land surface. While depth variables most strongly affected predicted mean ages, chemical constituents exhibited smooth trends with age, consistent with prevailing conceptual models of evolving sources and geochemistry flowpaths. The results provide proof of concept for use of readily available variables of well, landscape, and chemical characteristics to improve groundwater age estimates across large regions.

Great Lakes basin↗

Multi-lithofacies alluvial characterization via airborne electromagnetic-borehole fusion using ordinary interval kriging and geologic constraints

Airborne electromagnetic (AEM) survey provides extensive spatial coverage and detailed resolution at the near surface and can be used to develop hydrogeological models. However, utilization of AEM data is not straightforward because AEM resistivity is an indirect measurement for inferring sediment types. This study develops an ordinary interval kriging (OIK) algorithm and a resistivity-to-multi-lithofacies (R2ML) data fusion workflow for multi-lithofacies alluvial characterization. OIK utilizes irregular interval data to construct three-dimensional (3D) resistivity fields from one-dimensional inverted AEM resistivity models. The R2ML workflow maps the resistivity field generated from OIK into a multi-facies lithological model, incorporating geologic constraints derived from well logs and geological observations. The numerical and real-world cases demonstrate that OIK is computationally efficient, accounts for 3D anisotropy, and minimizes the smoothing effect, thereby preserving resistivity contrasts and reducing interpolation uncertainty. The methodology is applied to lithologic characterization of the Mississippi River Valley alluvial aquifer (MRVA) in the Shellmound area, Mississippi, U.S. A frequency-domain AEM survey was conducted to support groundwater studies for the managed aquifer recharge (MAR) to the MRVA. The resulting lithological model, including four types of lithofacies—clay, very fine sand, fine-medium sands, and graveliferous sands, illustrates the geomorphological processes of the MRVA and implies potential MAR. The alignment between the lithological model and existing geological and hydrogeological investigations demonstrates that OIK and R2ML workflow effectively capture the subsurface architecture of the MRVA. The methods have broad applicability for characterizing alluvial aquifers through AEM-borehole data fusion, supporting sustainable groundwater management.

Mississippi↗

Seismicity and tilt associated with the 2003 Anatahan eruption sequence

On May 10, 2003, the first historical eruption of Anatahan volcano in the western Pacific Mariana Islands was fortuitously recorded by a broadband seismograph installed on the island only 4 days prior to the eruption. This station, located 7 km WNW of the active crater, together with another broadband seismograph on Sarigan Island 45 km to the north, continued to operate throughout the 2-month period of major eruptive activity in May and June and throughout the majority of the following year. In June 2003, the Saipan Emergency Management Office and the US Geological Survey installed two telemetered high-gain short-period seismic stations to monitor the activity in real-time. The only earthquakes detected in the 4-day period from the initial seismograph installation until 6 h prior to the eruption occurred approximately 20 km to the northeast of the island on May 8. The first volcano-tectonic (VT) event located near the volcano occurred at 01:53 GMT on May 10. The number of events per hour then increased dramatically and a period of about 80 discrete earthquakes per hour commenced at about 06:20 GMT, immediately prior to the estimated eruption time of 07:30 from the Volcanic Ash Advisory Center. A long-period tilt signal recorded on the horizontal components of the broadband seismograph, indicating upward movement of the crater region, also commenced at about 06:20. Inflation continued until 09:30, when the direction of tilt reversed. Deflation continued until 17:50, coinciding with a reduction in the number of VT events. The larger VT events were located with a linearized least-squares location algorithm. Magnitudes of located VT events on May 10 ranged from 2.0 to 3.2, but a period of larger VT events were recorded on May 11, with the largest M 4.2. After about 36 h of intense earthquake activity, the number of discrete VT events declined and was replaced by nearly continuous volcanic tremor for the next 6 weeks. Differing types of very long-period events may suggest complex non-destructive magmatic source mechanisms, persisting dominantly throughout the first 10 days of the initial eruption. Visual reports indicate that a small craggy dome extruded sometime between May 20 and June 5. From analogy with other volcanic dome extrusions, we believe the dome probably extruded just as the tremor amplitude decreased dramatically about May 24. This dome was then destroyed between June 13 and 14. Reduced displacement of the co-eruption tremor is estimated as approximately 40-80 cm2, suggesting an eruption with a Volcanic Explosivity Index of about 3. 2005 Elsevier B.V. All rights reserved.

Northern Mariana Islands↗

An automated SO2 camera system for continuous, real-time monitoring of gas emissions from Kīlauea Volcano's summit Overlook Crater

SO 2 camera systems allow rapid two-dimensional imaging of sulfur dioxide (SO 2 ) emitted from volcanic vents. Here, we describe the development of an SO 2 camera system specifically designed for semi-permanent field installation and continuous use. The integration of innovative but largely “off-the-shelf” components allowed us to assemble a robust and highly customizable instrument capable of continuous, long-term deployment at Kīlauea Volcano's summit Overlook Crater. Recorded imagery is telemetered to the USGS Hawaiian Volcano Observatory (HVO) where a novel automatic retrieval algorithm derives SO 2 column densities and emission rates in real-time. Imagery and corresponding emission rates displayed in the HVO operations center and on the internal observatory website provide HVO staff with useful information for assessing the volcano's current activity. The ever-growing archive of continuous imagery and high-resolution emission rates in combination with continuous data from other monitoring techniques provides insight into shallow volcanic processes occurring at the Overlook Crater. An exemplary dataset from September 2013 is discussed in which a variation in the efficiency of shallow circulation and convection, the processes that transport volatile-rich magma to the surface of the summit lava lake, appears to have caused two distinctly different phases of lake activity and degassing. This first successful deployment of an SO 2 camera for continuous, real-time volcano monitoring shows how this versatile technique might soon be adapted and applied to monitor SO 2 degassing at other volcanoes around the world.

Hawai'i↗

A retrospective look at the February 1993 east rift zone intrusion at Kīlauea volcano, Hawaii

The February 1993 dike intrusion in the East Rift Zone (ERZ) of Kīlauea Volcano, Hawai'i, was recognized from tilt and seismic data, but ground-based geodetic data were too sparse to constrain the characteristics of the intrusion. Analysis of Interferometric Synthetic Aperture Radar (InSAR) from the Japan Aerospace Exploration Agency (JAXA) JERS-1 satellite reveals a maximum of ~30 cm of line-of-sight (LOS) displacement occurring near Makaopuhi Crater in the middle ERZ of Kīlauea. We model this deformation signal as a subvertical dike using a 3D-Mixed Boundary Element Method (3D-MBEM) paired with a nonlinear inversion algorithm to find the best-fit model. The best-fit dike is located just to the west of Makaopuhi Crater striking N50°W, extends to within 100 m of the surface, is ~1.3 km in length by ~4.2 km in width along strike, and has a total volume of ~7.4 × 10 6 m 3 . In addition, a post-intrusion interferogram from JERS-1 spanning 1993–1997 was analyzed. Guided by previous results, our model for the 4-year period consists of opening of the deep rift zones by about 0.5 m at 3–8.5 km depth beneath the Southwest Rift Zone (SWRZ), ERZ and the summit. A sub-horizontal detachment fault is connected to the seaward side of the vertical dike-like source to mimic the décollement known to exist beneath the volcano. We classify the 1993 dike intrusion as a passive intrusion similar to those that occurred in 1997 and 1999. Passive intrusions lack precursory inflation at Kīlauea's summit, and the likely triggering mechanism is persistent deep rift opening combined with seaward motion of the south flank along the basal décollement. Passive intrusions make forecasting and hazard assessment difficult since they are not preceded by inflation nor by large increases in seismicity.

Hawaii↗

Rock alteration mapping in and around fossil shallow intrusions at Mt. Ruapehu New Zealand with laboratory and aerial hyperspectral imaging

Diagnostic absorption features in hyperspectral data can be used to identify a specific mineral or mineral associations. However, it is unknown how accurate hyperspectral mapping can be for identifying alteration mineral compositions at the resolution required to describe structures such as fossil intrusions, or whether it can accurately quantify the alteration present. This study compared petrographic observation with visible, near-infrared (VNIR), and shortwave infrared (SWIR) hyperspectral remote sensing at laboratory- (centimetre-scale) and aerial- (metre-scale) scales to characterise the abundance of surface hydrothermal rock alteration in and around a shallow fossil intrusion on Pinnacle Ridge, Mt. Ruapehu, New Zealand. We classified a high-resolution aerial hyperspectral image to develop a new surface alteration map using Spectral Angle Mapper (SAM) algorithm. The petrographic thin-section and the laboratory and aerial hyperspectral imaging revealed a spectrum of hydrous alteration phases as indicated by the presence of an absorption feature at 2207 nm. Moderate correlation exists between the depth of the absorption feature at 2207 nm and the point counting-derived alteration percent values, indicating reliability of laboratory-based hyperspectral analytical methods. In contrast, aerial hyperspectral data failed to provide any clear correlations to field-mapped alteration using a band-depth approach, and we interpret this due to ‘oversampling’ of surface (supergene) alteration, spectral mixing, and sensor limitations (e.g., bandwidth, signal-to-noise ratio). The hyperspectral image-derived alteration map, created using supervised image classification, can loosely be translated to a geotechnical map where porosity and permeability play a major role in localizing hydrothermal fluid flow and the formation of alteration mineral associations.

Mt. Ruapehu↗

Modeling surface wave dynamics in upper Delaware Bay with living shorelines

Living shorelines gain increasing attention because they stabilize shorelines and reduce erosion. This study leverages physics-based models and bagged regression tree (BRT) machine learning algorithm to simulate wave dynamics at a living shoreline composed of constructed oyster reefs (CORs) in upper Delaware Bay. The physics-based models consist of coupled Delft3D-FLOW and SWAN in four-level nested domains. The model accuracy converges with increasing mesh resolution. The simulated wave-induced current circulation substantiates the effectiveness of CORs in trapping sediments. The simulated yearly-averaged wave power correlates qualitatively with historical shoreline retreat rates. BRT is adopted to improve the model accuracy, identify key processes responsible for simulation errors in wave height ( H 8 ) and wave period ( T p ), and quantify their importance. In the CORs sheltered area, BRT reveals that simulation errors of wind seas mainly arise from wind forcing, wave breaking and wave triad interactions. Wave breaking is seven times more important than wind forcing for simulating H 8 , while wind forcing and triad interactions are of equal importance for simulating T p . Simulation errors of swells mostly stem from bottom friction and offshore wave boundary conditions. Results from this study can help the assessment and adaptive management of CORs-based living shoreline restoration projects under climate change.

New Jersey↗

Effect of roughness formulation on the performance of a coupled wave, hydrodynamic, and sediment transport model

A variety of algorithms are available for parameterizing the hydrodynamic bottom roughness associated with grain size, saltation, bedforms, and wave–current interaction in coastal ocean models. These parameterizations give rise to spatially and temporally variable bottom-drag coefficients that ostensibly provide better representations of physical processes than uniform and constant coefficients. However, few studies have been performed to determine whether improved representation of these variable bottom roughness components translates into measurable improvements in model skill. We test the hypothesis that improved representation of variable bottom roughness improves performance with respect to near-bed circulation, bottom stresses, or turbulence dissipation. The inner shelf south of Martha’s Vineyard, Massachusetts, is the site of sorted grain-size features which exhibit sharp alongshore variations in grain size and ripple geometry over gentle bathymetric relief; this area provides a suitable testing ground for roughness parameterizations. We first establish the skill of a nested regional model for currents, waves, stresses, and turbulent quantities using a uniform and constant roughness; we then gauge model skill with various parameterization of roughness, which account for the influence of the wave-boundary layer, grain size, saltation, and rippled bedforms. We find that commonly used representations of ripple-induced roughness, when combined with a wave–current interaction routine, do not significantly improve skill for circulation, and significantly decrease skill with respect to stresses and turbulence dissipation. Ripple orientation with respect to dominant currents and ripple shape may be responsible for complicating a straightforward estimate of the roughness contribution from ripples. In addition, sediment-induced stratification may be responsible for lower stresses than predicted by the wave–current interaction model.

Ocean Modelling↗

HyWaves: Hybrid downscaling of multimodal wave spectra to nearshore areas

Long-term and accurate wave hindcast databases are often required in different coastal engineering projects. The assessment of the nearshore wave climate is often accomplished by using downscaling techniques to translate offshore waves to coastal areas. However, dynamical downscaling approaches may incur huge computational cost. Additionally, the common use of bulk parameterizations are often not accurate for multidimensional waves. To overcome these limitations, we present a hybrid downscaling approach that combines mathematical algorithms (statistical downscaling) and numerical modeling (dynamical downscaling) over the individual spectral partitions. Every wave partition is downscaled and aggregated afterward by using principles of wave linear theory. By assuming linearity in the propagation of the wave celerity, the application of the method is limited from offshore to intermediate water depths. In addition, the method proposed uses a technique to simplify the spectral boundary conditions in complex domains. The methodology has been applied and validated in the island states of Samoa, American Samoa, Majuro, and Kwajalein, showing good skill at reproducing the spectral hourly time series of significant wave height, peak period, and peak direction. Moreover, an accurate representation of the observed energy spectrum was achieved. This study provides insight into the numerical approximation of the combined sea-swell states while improving the quality of fast spectral forecasting and early warning systems.

Ocean Modeling↗

TreeMAC: Localized TDMA MAC protocol for real-time high-data-rate sensor networks

Earlier sensor network MAC protocols focus on energy conservation in low-duty cycle applications, while some recent applications involve real-time high-data-rate signals. This motivates us to design an innovative localized TDMA MAC protocol to achieve high throughput and low congestion in data collection sensor networks, besides energy conservation. TreeMAC divides a time cycle into frames and each frame into slots. A parent node determines the children's frame assignment based on their relative bandwidth demand, and each node calculates its own slot assignment based on its hop-count to the sink. This innovative 2-dimensional frame-slot assignment algorithm has the following nice theory properties. First, given any node, at any time slot, there is at most one active sender in its neighborhood (including itself). Second, the packet scheduling with TreeMAC is bufferless, which therefore minimizes the probability of network congestion. Third, the data throughput to the gateway is at least 1/3 of the optimum assuming reliable links. Our experiments on a 24-node testbed show that TreeMAC protocol significantly improves network throughput, fairness, and energy efficiency compared to TinyOS's default CSMA MAC protocol and a recent TDMA MAC protocol Funneling-MAC. Partial results of this paper were published in Song, Huang, Shirazi and Lahusen [W.-Z. Song, R. Huang, B. Shirazi, and R. Lahusen, TreeMAC: Localized TDMA MAC protocol for high-throughput and fairness in sensor networks, in: The 7th Annual IEEE International Conference on Pervasive Computing and Communications, PerCom, March 2009]. Our new contributions include analyses of the performance of TreeMAC from various aspects. We also present more implementation detail and evaluate TreeMAC from other aspects. ?? 2009 Elsevier B.V.

Pervasive and Mobile Computing↗

The imaging performance of the SRC on Mars Express

The Mars Express spacecraft carries the pushbroom scanner high-resolution stereo camera (HRSC) and its added imaging subsystem super resolution channel (SRC). The SRC is equipped with its own optical system and a 1024×1024 framing sensor. SRC produces snapshots with 2.3 m ground pixel size from the nominal spacecraft pericenter height of 250 km, which are typically embedded in the central part of the large HRSC scenes. The salient features of the SRC are its light-weight optics, a reliable CCD detector, and high-speed read-out electronics. The quality and effective visibility of details in the SRC images unfortunately falls short of what has been expected. In cases where thermal balance cannot be reached, artifacts, such as blurring and “ghost features” are observed in the images. In addition, images show large numbers of blemish pixels and are plagued by electronic noise. As a consequence, we have developed various image improving algorithms, which are discussed in this paper. While results are encouraging, further studies of image restoration by dedicated processing appear worthwhile. The SRC has obtained more than 6940 images at the time of writing (1 September 2007), which often show fascinating details in surface morphology. SRC images are highly useful for a variety of applications in planetary geology, for studies of the Mars atmosphere, and for astrometric observations of the Martian satellites. This paper will give a full account of the design philosophy, technical concept, calibration, operation, integration with HRSC, and performance, as well as science accomplishments of the SRC.

Planetary and Space Science↗

Estimating abiotic thresholds for sagebrush condition class in the western United States

Sagebrush ecosystems of the western United States can transition from extended periods of relatively stable conditions to rapid ecological change if acute disturbances occur. Areas dominated by native sagebrush can transition from species-rich native systems to altered states where non-native annual grasses dominate, if resistance to annual grasses is low. The non-native annual grasses provide relatively little value to wildlife, livestock, and humans and function as fuel that increases fire frequency. The more land area covered by annual grasses, the higher the potential for fire, thus reducing the potential for native vegetation to reestablish, even when applying restoration treatments. Mapping areas of stability and areas of change using machine-learning algorithms allows both the identification of dominant abiotic variables that drive ecosystem dynamics and the variables’ important thresholds. We develop a decision-tree model with rulesets that estimate three classes of sagebrush condition (i.e., sagebrush recovery, tipping point [ecosystem degradation], and stable). We find rulesets that primarily drive development of the sagebrush recovery class indicate areas of midelevations (1 602 m), warm 30-yr July temperature maximums (tmax) (30.62°C), and 30-yr March precipitation (ppt) averages equal to 26.26 mm, about 10% of the 30-yr annual ppt values. Tipping point and stable classes occur at elevations that are lower (1 505 m) and higher (1 939 m), respectively, more mesic during March and annually, and experience lower 30-yr July tmax averages. These defined variable averages can be used to understand current dynamics of sagebrush condition and to predict where future transitions may occur under novel conditions.

Arizona, California, Colorado, Idaho, Montana, Neb↗

A comparative analysis of the Global Land Cover 2000 and MODIS land cover data sets

Accurate and up-to-date global land cover data sets are necessary for various global change research studies including climate change, biodiversity conservation, ecosystem assessment, and environmental modeling. In recent years, substantial advancement has been achieved in generating such data products. Yet, we are far from producing geospatially consistent high-quality data at an operational level. We compared the recently available Global Land Cover 2000 (GLC-2000) and MODerate resolution Imaging Spectrometer (MODIS) global land cover data to evaluate the similarities and differences in methodologies and results, and to identify areas of spatial agreement and disagreement. These two global land cover data sets were prepared using different data sources, classification systems, and methodologies, but using the same spatial resolution (i.e., 1 km) satellite data. Our analysis shows a general agreement at the class aggregate level except for savannas/shrublands, and wetlands. The disagreement, however, increases when comparing detailed land cover classes. Similarly, percent agreement between the two data sets was found to be highly variable among biomes. The identified areas of spatial agreement and disagreement will be useful for both data producers and users. Data producers may use the areas of spatial agreement for training area selection and pay special attention to areas of disagreement for further improvement in future land cover characterization and mapping. Users can conveniently use the findings in the areas of agreement, whereas users might need to verify the informaiton in the areas of disagreement with the help of secondary information. Learning from past experience and building on the existing infrastructure (e.g., regional networks), further research is necessary to (1) reduce ambiguity in land cover definitions, (2) increase availability of improved spatial, spectral, radiometric, and geometric resolution satellite data, and (3) develop advanced classification algorithms.

Remote Sensing of Environment↗

Mapping an invasive plant, Phragmites australis, in coastal wetlands using the EO-1 Hyperion hyperspectral sensor

Mapping tools are needed to document the location and extent of Phragmites australis, a tall grass that invades coastal marshes throughout North America, displacing native plant species and degrading wetland habitat. Mapping Phragmites is particularly challenging in the freshwater Great Lakes coastal wetlands due to dynamic lake levels and vegetation diversity. We tested the applicability of Hyperion hyperspectral satellite imagery for mapping Phragmites in wetlands of the west coast of Green Bay in Wisconsin, U.S.A. A reference spectrum created using Hyperion data from several pure Phragmites stands within the image was used with a Spectral Correlation Mapper (SCM) algorithm to create a raster map with values ranging from 0 to 1, where 0 represented the greatest similarity between the reference spectrum and the image spectrum and 1 the least similarity. The final two-class thematic classification predicted monodominant Phragmites covering 3.4% of the study area. Most of this was concentrated in long linear features parallel to the Green Bay shoreline, particularly in areas that had been under water only six years earlier when lake levels were 66??cm higher. An error matrix using spring 2005 field validation points (n = 129) showed good overall accuracy-81.4%. The small size and linear arrangement of Phragmites stands was less than optimal relative to the sensor resolution, and Hyperion's 30??m resolution captured few if any pure pixels. Contemporary Phragmites maps prepared with Hyperion imagery would provide wetland managers with a tool that they currently lack, which could aid attempts to stem the spread of this invasive species. ?? 2006 Elsevier Inc. All rights reserved.

Remote Sensing of Environment↗

Using LiDAR and quickbird data to model plant production and quantify uncertainties associated with wetland detection and land cover generalizations

Spatiotemporal data from satellite remote sensing and surface meteorology networks have made it possible to continuously monitor global plant production, and to identify global trends associated with land cover/use and climate change. Gross primary production (GPP) and net primary production (NPP) are routinely derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard satellites Terra and Aqua, and estimates generally agree with independent measurements at validation sites across the globe. However, the accuracy of GPP and NPP estimates in some regions may be limited by the quality of model input variables and heterogeneity at fine spatial scales. We developed new methods for deriving model inputs (i.e., land cover, leaf area, and photosynthetically active radiation absorbed by plant canopies) from airborne laser altimetry (LiDAR) and Quickbird multispectral data at resolutions ranging from about 30??m to 1??km. In addition, LiDAR-derived biomass was used as a means for computing carbon-use efficiency. Spatial variables were used with temporal data from ground-based monitoring stations to compute a six-year GPP and NPP time series for a 3600??ha study site in the Great Lakes region of North America. Model results compared favorably with independent observations from a 400??m flux tower and a process-based ecosystem model (BIOME-BGC), but only after removing vapor pressure deficit as a constraint on photosynthesis from the MODIS global algorithm. Fine-resolution inputs captured more of the spatial variability, but estimates were similar to coarse-resolution data when integrated across the entire landscape. Failure to account for wetlands had little impact on landscape-scale estimates, because vegetation structure, composition, and conversion efficiencies were similar to upland plant communities. Plant productivity estimates were noticeably improved using LiDAR-derived variables, while uncertainties associated with land cover generalizations and wetlands in this largely forested landscape were considered less important.

Remote Sensing of Environment↗

Spectral assessment of new ASTER SWIR surface reflectance data products for spectroscopic mapping of rocks and minerals

ASTER reflectance spectra from Cuprite, Nevada, and Mountain Pass, California, were compared to spectra of field samples and to ASTER-resampled AVIRIS reflectance data to determine spectral accuracy and spectroscopic mapping potential of two new ASTER SWIR reflectance datasets: RefL1b and AST_07XT. RefL1b is a new reflectance dataset produced for this study using ASTER Level 1B data, crosstalk correction, radiance correction factors, and concurrently acquired level 2 MODIS water vapor data. The AST_07XT data product, available from EDC and ERSDAC, incorporates crosstalk correction and non-concurrently acquired MODIS water vapor data for atmospheric correction. Spectral accuracy was determined using difference values which were compiled from ASTER band 5/6 and 9/8 ratios of AST_07XT or RefL1b data subtracted from similar ratios calculated for field sample and AVIRIS reflectance data. In addition, Spectral Analyst, a statistical program that utilizes a Spectral Feature Fitting algorithm, was used to quantitatively assess spectral accuracy of AST_07XT and RefL1b data.Spectral Analyst matched more minerals correctly and had higher scores for the RefL1b data than for AST_07XT data. The radiance correction factors used in the RefL1b data corrected a low band 5 reflectance anomaly observed in the AST_07XT and AST_07 data but also produced anomalously high band 5 reflectance in RefL1b spectra with strong band 5 absorption for minerals, such as alunite. Thus, the band 5 anomaly seen in the RefL1b data cannot be corrected using additional gain adjustments. In addition, the use of concurrent MODIS water vapor data in the atmospheric correction of the RefL1b data produced datasets that had lower band 9 reflectance anomalies than the AST_07XT data. Although assessment of spectral data suggests that RefL1b data are more consistent and spectrally more correct than AST_07XT data, the Spectral Analyst results indicate that spectral discrimination between some minerals, such as alunite and kaolinite, are still not possible unless additional spectral calibration using site specific spectral data are performed. ?? 2010.

Remote Sensing of Environment↗

A simple and effective method for filling gaps in Landsat ETM+ SLC-off images

The scan-line corrector (SLC) of the Landsat 7 Enhanced Thematic Mapper Plus (ETM+) sensor failed in 2003, resulting in about 22% of the pixels per scene not being scanned. The SLC failure has seriously limited the scientific applications of ETM+ data. While there have been a number of methods developed to fill in the data gaps, each method has shortcomings, especially for heterogeneous landscapes. Based on the assumption that the same-class neighboring pixels around the un-scanned pixels have similar spectral characteristics, and that these neighboring and un-scanned pixels exhibit similar patterns of spectral differences between dates, we developed a simple and effective method to interpolate the values of the pixels within the gaps. We refer to this method as the Neighborhood Similar Pixel Interpolator (NSPI). Simulated and actual SLC-off ETM+ images were used to assess the performance of the NSPI. Results indicate that NSPI can restore the value of un-scanned pixels very accurately, and that it works especially well in heterogeneous regions. In addition, it can work well even if there is a relatively long time interval or significant spectral changes between the input and target image. The filled images appear reasonably spatially continuous without obvious striping patterns. Supervised classification using the maximum likelihood algorithm was done on both gap-filled simulated SLC-off data and the original "gap free" data set, and it was found that classification results, including accuracies, were very comparable. This indicates that gap-filled products generated by NSPI will have relevance to the user community for various land cover applications. In addition, the simple principle and high computational efficiency of NSPI will enable processing large volumes of SLC-off ETM+ data.

Remote Sensing of Environment↗

A comprehensive change detection method for updating the National Land Cover Database to circa 2011

The importance of characterizing, quantifying, and monitoring land cover, land use, and their changes has been widely recognized by global and environmental change studies. Since the early 1990s, three U.S. National Land Cover Database (NLCD) products (circa 1992, 2001, and 2006) have been released as free downloads for users. The NLCD 2006 also provides land cover change products between 2001 and 2006. To continue providing updated national land cover and change datasets, a new initiative in developing NLCD 2011 is currently underway. We present a new Comprehensive Change Detection Method (CCDM) designed as a key component for the development of NLCD 2011 and the research results from two exemplar studies. The CCDM integrates spectral-based change detection algorithms including a Multi-Index Integrated Change Analysis (MIICA) model and a novel change model called Zone, which extracts change information from two Landsat image pairs. The MIICA model is the core module of the change detection strategy and uses four spectral indices (CV, RCVMAX, dNBR, and dNDVI) to obtain the changes that occurred between two image dates. The CCDM also includes a knowledge-based system, which uses critical information on historical and current land cover conditions and trends and the likelihood of land cover change, to combine the changes from MIICA and Zone. For NLCD 2011, the improved and enhanced change products obtained from the CCDM provide critical information on location, magnitude, and direction of potential change areas and serve as a basis for further characterizing land cover changes for the nation. An accuracy assessment from the two study areas show 100% agreement between CCDM mapped no-change class with reference dataset, and 18% and 82% disagreement for the change class for WRS path/row p22r39 and p33r33, respectively. The strength of the CCDM is that the method is simple, easy to operate, widely applicable, and capable of capturing a variety of natural and anthropogenic disturbances potentially associated with land cover changes on different landscapes.

Remote Sensing of Environment↗