Search USGS⌕ Search

SEARCH · Search USGS

Results for “Stacks Journal”

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 37 records · Page 2Linked to original sources
USGS2012

Monitoring glacier surface seismicity in time and space using Rayleigh waves

Sliding glaciers and brittle ice failure generate seismic body and surface wave energy characteristic to the source mechanism. Here we analyze continuous seismic recordings from an array of nine short-period passive seismometers located on Bench Glacier, Alaska (USA) (61.033°N, 145.687°W). We focus on the arrival-time and amplitude information of the dominant Rayleigh wave phase. Over a 46-hour period we detect thousands of events using a cross-correlation based event identification method. Travel-time inversion of a subset of events (7% of the total) defines an active crevasse, propagating more than 200 meters in three hours. From the Rayleigh wave amplitudes, we estimate the amount of volumetric opening along the crevasse as well as an average bulk attenuation ( = 42) for the ice in this part of the glacier. With the remaining icequake signals we establish a diurnal periodicity in seismicity, indicating that surface run-off and subglacial water pressure changes likely control the triggering of these surface events. Furthermore, we find that these events are too weak (i.e., too noisy) to locate individually. However, stacking individual events increases the signal-to-noise ratio of the waveforms, implying that these periodic sources are effectively stationary during the recording period.

Journal of Geophysical Research F: Earth Surface↗
USGS1996

Depositional controls on coal distribution and quality in the Eocene Brunner Coal Measures, Buller Coalfield, South Island, New Zealand

The Buller Coalfield on the West Coast of the South Island, New Zealand, contains the Eocene Brunner Coal Measures. The coal measures unconformably overlie Paleozoic-Cretaceous basement rocks and are conformably overlain by, and laterally interfinger with, the Eocene marine Kaiata Formation. This study examines the lithofacies frameworks of the coal measures in order to interpret their depositional environments. The lower part of the coal measures is dominated by conglomeratic lithofacies that rest on a basal erosional surface and thicken in paleovalleys incised into an undulating peneplain surface. These lithofacies are overlain by sandstone, mudstone and organic-rich lithofacies of the upper part of the coal measures. The main coal seam of the organic-rich lithofacies is thick (10-20 m), extensive, locally split, and locally absent. This seam and associated coal seams in the Buller Coalfield are of low- to high-volatile bituminous rank (vitrinite reflectance between 0.65% and 1.75%). The main seam contains a variable percentage of ash and sulphur. These values are related to the thickening and areal distribution of the seam, which in turn, were controlled by the nature of clastic deposition and peat-forming mire systems, marine transgression and local tidal incursion. The conglomeratic lithofacies represent deposits of trunk and tributary braided streams that rapidly aggraded incised paleovalleys during sea-level stillstands. The main seam represents a deposit of raised mires that initially developed as topogenous mires on abandoned margins of inactive braidbelts. Peat accumulated in mires as a response to a rise in the water table, probably initially due to gradual sea-level rise and climate, and the resulting raised topography served as protection from floods. The upper part of the coal measures consists of sandstone lithofacies of flu vial origin and bioturbated sandstone, mudstone and organic-rich lithofacies, which represent deposits of paralic (deltaic, barrier shoreface, tidal and mire) and marine environments. The fluvial sandstone lithofacies accumulated in channels during a sea-level stillstand. The channels were infilled by coeval braided and meandering streams prior to transgression. Continued transgression, ranging from tidal channel-estuarine incursions to widespread but uneven paleoshoreline encroachment, accompanied by moderate basin subsidence, is marked by a stacked, back-stepping geometry of bioturbated sandstone and marine mudstone lithofacies. Final retrogradation (sea-level highstand) is marked by backfilling of estuaries and by rapid landward deposition of the marine Kaiata Formation in the late Eocene.

International Journal of Coal Geology↗
USGS2026

Bulk location and source mechanisms of repetitive long-period seismicity at Mount St. Helens in July 2005

The 2004–2008 eruption of Mount St. Helens exemplified the practical challenges of processing high seismic event rates and data volumes. This protracted multi-year eruptive sequence was accompanied by a near-continuous occurrence of long-period (LP, 0.5–5 Hz) seismicity, but network-based picks and locations have previously been performed for only a relatively small subset of the overall seismicity, and full-waveform source–mechanism inversions have previously been performed for just two event types during the eruption. We present a new workflow for the bulk detection, cluster picking, 3-D location, and full-waveform source–mechanism inversion of repetitive LP volcanic seismicity recorded by a local broadband seismic network. Our emphasis is on methods that: exploit waveform similarity (multiplets) commonly observed for volcanic LP seismicity swarms; can take into account 3-D topographic and velocity structure heterogeneity typical of volcanoes; can process high event rates with relatively low analyst burden; and do not require a priori catalogue picks or known waveform templates. We develop and illustrate the workflow using data from Mount St. Helens from 2005 June 29 to July 31, a time-window characterized by high-rate LP seismicity, including numerous lower amplitude small LPs, and for which relatively dense seismic network coverage is available. Approximately 40 per cent of seismicity (20 561 events out of 52 141 triggers) during 2005 June 29 to July 31 is repetitive enough to be grouped into 17 multiplets of at least 300 events. We present 3-D locations of the 17 multiplets, representing 20 076 events having picks from at least eight stations, along with full-waveform source mechanism inversions of the 17 separate multiplet stacks. Collectively, our seismic 3-D hypocentre distributions for 2005 June 29 to July 31 are compatible with sources around the margins of the upper conduit section that is associated with plug solidification. All location results indicate compact seismicity clusters, with the 17 multiplets spatially connected and collectively located at shallow depth, approximately ~ 775 m depth beneath the 1980s crater floor and directly beneath, or immediately adjacent to, the extruding dacite spine. The full-waveform source–mechanism inversions for the 17 multiplets generally indicate mechanisms dominated by volumetric components consistent with subhorizontal cracks along with minor contributions from subvertical cracks. The overall picture is that of a stack or plexus of subhorizontal cracks located at shallow depth beneath the extruding dacite spines, consistent with subhorizontal hydrofractures extending from the conduit margins into the surrounding wall rock; minor subvertical crack components of the mechanisms could possibly arise from surfaces parallel to the conduit wall, such as within the damage halo associated with degassing pathways.

Washington↗
USGS2000

High-resolution seismic reflection to delineate shallow gas in Eastern Kansas

Unique amplitude characteristics of shallow gas sands within Pennsylvanian clastic-carbonate dominated sequences are discernible on high-resolution seismic reflection data in eastern Kansas. Upward grading sequences of sand into shale represent a potential gas reservoir with a low-impedence acoustic contrast at the base of the encasing layer. The gas sand and encasing shale, which define the gas reservoir studied here, are part of an erosional incised valley where about 30 m of carbonates and shale have been replaced by sandstone and shale confined to the incised valley. These consolidated geologic settings would normally possess high impedence gas sand reservoirs as defined by abrupt contacts between the gas sand and encasing shale. Based orr core and borehole logs, the gas sand studied here grades from sand into shale in a fashion analogous to that observed in classic low-impedance environments. Amplitude and phase characteristics of high-resolution seismic data across this approximately 400-m wide gas sand are indicative of a low-impedance reservoir. Shot gathers possess classic amplitude with offsett-dependent characteristics which are manifeted on the stacked section as "bright spots." Dominant Frequencies of around 120 Hz allow detection of several reflectors within the 30+ meters of sand/shale that make up this localized gas-rich incised valley fill. The gradational nature of the trapping mechanism observed in this gas reservoir would make detection with conventional seismic reflection methods unlikely.

The Compass: Earth Science Journal of Sigma Gamma ↗
USGS2016

The potassic sedimentary rocks in Gale Crater, Mars, as seen by ChemCam Onboard Curiosity

The Mars Science Laboratory rover Curiosity encountered potassium-rich clastic sedimentary rocks at two sites in Gale Crater, the waypoints Cooperstown and Kimberley. These rocks include several distinct meters thick sedimentary outcrops ranging from fine sandstone to conglomerate, interpreted to record an ancient fluvial or fluvio-deltaic depositional system. From ChemCam Laser-Induced Breakdown Spectroscopy (LIBS) chemical analyses, this suite of sedimentary rocks has an overall mean K 2 O abundance that is more than 5 times higher than that of the average Martian crust. The combined analysis of ChemCam data with stratigraphic and geographic locations reveals that the mean K 2 O abundance increases upward through the stratigraphic section. Chemical analyses across each unit can be represented as mixtures of several distinct chemical components, i.e., mineral phases, including K-bearing minerals, mafic silicates, Fe-oxides, and Fe-hydroxide/oxyhydroxides. Possible K-bearing minerals include alkali feldspar (including anorthoclase and sanidine) and K-bearing phyllosilicate such as illite. Mixtures of different source rocks, including a potassium-rich rock located on the rim and walls of Gale Crater, are the likely origin of observed chemical variations within each unit. Physical sorting may have also played a role in the enrichment in K in the Kimberley formation. The occurrence of these potassic sedimentary rocks provides additional evidence for the chemical diversity of the crust exposed at Gale Crater.

Journal of Geophysical Research↗
USGS2008

Multi-residue method for the analysis of 85 current-use and legacy pesticides in bed and suspended sediments

A multi-residue method was developed for the simultaneous determination of 85 current-use and legacy organochlorine pesticides in a single sediment sample. After microwave-assisted extraction, clean-up of samples was optimized using gel permeation chromatography and either stacked carbon and alumina solid-phase extraction cartridges or a deactivated Florisil column. Analytes were determined by gas chromatography with ion-trap mass spectrometry and electron capture detection. Method detection limits ranged from 0.6 to 8.9 ??g/kg dry weight. Bed and suspended sediments from a variety of locations were analyzed to validate the method and 29 pesticides, including at least 1 from every class, were detected.

Journal of Chromatography A↗
USGS2020

Geodetic measurements of slow slip events southeast of Parkfield, CA

Tremor and low-frequency earthquakes are presumed to be indicative of surrounding slow, aseismic slip that is often below geodetic detection thresholds. This study uses data from borehole seismometers and long-baseline laser strainmeters to observe both the seismic and geodetic signatures of episodic tremor and slip on the Parkfield region of the San Andreas Fault near Cholame, CA. The observed occurrence rates of both the tremors and co-located families of low-frequency earthquakes are not steady but instead exhibit quasiperiodic bursts of increased activity. We show that these periods of elevated seismic activity correlate with statistically significant stacked strain signals consisting of 44 slow-slip events. Modeled individual slow-slip events and their total summed moment, which are constrained by seismic signals and stacked strain, respectively, indicate that the individual moment magnitudes of these events range from 4.6–5.2. We find that the measured geodetic signal likely precedes the seismic signal by several hours, consistent with the aseismic slip preceding and driving the observed seismic tremor activity. We confirm that strike-slip faults, in addition to subduction zones, are capable of producing episodic tremor and slip.

California↗
USGS1999

Sequential filling of a late paleozoic foreland basin

Through the use of an extensive data base of geophysical well logs, parasequence-scale subdivisions within a late Paleozoic synorogenic clastic wedge resolve cycles of sequential subsidence of a foreland basin, sediment progradation, subsidence of a carbonate shelf edge, diachronously subsiding discrete depositional centers, and basinwide transgression. Although temporal resolution of biostratigraphic markers is less precise in Paleozoic successions than in younger basins, parasequence-scale subdivisions provide more detailed resolution within marker-defined units in Paleozoic strata. As an example, the late Paleozoic Black Warrior basin in the foreland of the Ouachita thrust belt is filled with a synorogenic clastic wedge, the lower part of which intertongues with the fringe of a cratonic carbonate facies in the distal part of the basin. The stratal geometry of one tongue of the carbonate facies (lower tongue of Bangor Limestone) defines a ramp that grades basinward into a thin black shale. An overlying tongue of the synorogenic clastic wedge (lower tongue of Parkwood Formation) consists of cyclic delta and delta-front deposits, in which parasequences are defined by marine-flooding surfaces above coarsening- and shallowing-upward successions of mudstone and sandstone. Within the lower Parkwood tongue, two genetic stratigraphic sequences (A and B) are defined by parasequence offlap and downlap patterns and are bounded at the tops by basinwide maximum-flooding surfaces. The distribution of parasequences within sequences A and B indicates two cycles of sequential subsidence (deepening) and progradation, suggesting subsidence during thrust advance and progradation during thrust quiescence. Parasequence stacking in sequences A and B also indicates diachronous differential tectonic subsidence of two discrete depositional centers within the basin. The uppermost sequence (C) includes reworked sandstones and an overlying shallow-marine limestone, a vertical succession that reflects no tectonic subsidence, a very minor or null sediment supply, and basinwide transgression. The temporal resolution at parasequence scale significantly improves the resolution of the tectonic history of the thrust belt-foreland basin system.

Journal of Sedimentary Research↗
USGS2009

Rocky Mountain evolution: Tying Continental Dynamics of the Rocky Mountains and Deep Probe seismic experiments with receiver functions

In this study, we have determined the crustal structure using three different receiver function methods using data collected from the northern transect of the Continental Dynamics of the Rocky Mountains (CD-ROM) experiment. The resulting migrated image and crustal thickness determinations confirm and refine prior crustal thickness measurements based on the CD-ROM and Deep Probe experiment data sets. The new results show a very distinct and thick lower crustal layer beneath the Archean Wyoming province. In addition, we are able to show its termination at 42??N latitude, which provides a seismic tie between the CD-ROM and Deep Probe seismic experiments and thus completes a continuous north-south transect extending from New Mexico into Alberta, Canada. This new tie is particularly important because it occurs close to a major tectonic boundary, the Cheyenne belt, between an Archean craton and a Proterozoic terrane. We used two different stacking techniques, based on a similar concept but using two different ways to estimate uncertainties. Furthermore, we used receiver function migration and common conversion point (CCP) stacking techniques. The combined interpretation of all our results shows (1) crustal thinning in southern Wyoming, (2) strong northward crustal thickening beginning in central Wyoming, (3) the presence of an unusually thick and high-velocity lower crust beneath the Wyoming province, and (4) the abrupt termination of this lower crustal layer north of the Cheyenne belt at 42??N latitude. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2007

Multi-interferogram method for measuring interseismic deformation: Denali Fault, Alaska

Studies of interseismic strain accumulation are crucial to our understanding of continental deformation, the earthquake cycle and seismic hazard. By mapping small amounts of ground deformation over large spatial areas, InSAR has the potential to produce continental-scale maps of strain accumulation on active faults. However, most InSAR studies to date have focused on areas where the coherence is relatively good (e.g. California, Tibet and Turkey) and most analysis techniques (stacking, small baseline subset algorithm, permanent scatterers, etc.) only include information from pixels which are coherent throughout the time-span of the study. In some areas, such as Alaska, where the deformation rate is small and coherence very variable, it is necessary to include information from pixels which are coherent in some but not all interferograms. We use a three-stage iterative algorithm based on distributed scatterer interferometry. We validate our method using synthetic data created using realistic parameters from a test site on the Denali Fault, Alaska, and present a preliminary result of 10.5 ?? 5.0 mm yr-1 for the slip rate on the Denali Fault based on a single track of radar data from ERS1/2. ?? 2007 The Authors Journal compilation ?? 2007 RAS.

Alaska↗
USGS2012

Cross-sensor comparisons between Landsat 5 TM and IRS-P6 AWiFS and disturbance detection using integrated Landsat and AWiFS time-series images

Routine acquisition of Landsat 5 Thematic Mapper (TM) data was discontinued recently and Landsat 7 Enhanced Thematic Mapper Plus (ETM+) has an ongoing problem with the scan line corrector (SLC), thereby creating spatial gaps when covering images obtained during the process. Since temporal and spatial discontinuities of Landsat data are now imminent, it is therefore important to investigate other potential satellite data that can be used to replace Landsat data. We thus cross-compared two near-simultaneous images obtained from Landsat 5 TM and the Indian Remote Sensing (IRS)-P6 Advanced Wide Field Sensor (AWiFS), both captured on 29 May 2007 over Los Angeles, CA. TM and AWiFS reflectances were compared for the green, red, near-infrared (NIR), and shortwave infrared (SWIR) bands, as well as the normalized difference vegetation index (NDVI) based on manually selected polygons in homogeneous areas. All R 2 values of linear regressions were found to be higher than 0.99. The temporally invariant cluster (TIC) method was used to calculate the NDVI correlation between the TM and AWiFS images. The NDVI regression line derived from selected polygons passed through several invariant cluster centres of the TIC density maps and demonstrated that both the scene-dependent polygon regression method and TIC method can generate accurate radiometric normalization. A scene-independent normalization method was also used to normalize the AWiFS data. Image agreement assessment demonstrated that the scene-dependent normalization using homogeneous polygons provided slightly higher accuracy values than those obtained by the scene-independent method. Finally, the non-normalized and relatively normalized ‘Landsat-like’ AWiFS 2007 images were integrated into 1984 to 2010 Landsat time-series stacks (LTSS) for disturbance detection using the Vegetation Change Tracker (VCT) model. Both scene-dependent and scene-independent normalized AWiFS data sets could generate disturbance maps similar to what were generated using the LTSS data set, and their kappa coefficients were higher than 0.97. These results indicate that AWiFS can be used instead of Landsat data to detect multitemporal disturbance in the event of Landsat data discontinuity.

International Journal of Remote Sensing↗
USGS2001

Ancient drainage basin of the Tharsis region, Mars: Potential source for outflow channel systems and putative oceans or paleolakes

Paleotopographic reconstructions based on a synthesis of published geologic information and high-resolution topography, including topographic profiles, reveal the potential existence of an enormous drainage basin/aquifer system in the eastern part of the Tharsis region during the Noachian Period. Large topographic highs formed the margin of the gigantic drainage basin. Subsequently, lavas, sediments, and volatiles partly infilled the basin, resulting in an enormous and productive regional aquifer. The stacked sequences of water-bearing strata were then deformed locally and, in places, exposed by magmatic-driven uplifts, tectonic deformation, and erosion. This basin model provides a potential source of water necessary to carve the large outflow channel systems of the Tharsis and surrounding regions and to contribute to the formation of putative northern-plains ocean(s) and/or paleolakes.

Journal of Geophysical Research E: Planets↗
USGS2007

Spectral element modelling of fault-plane reflections arising from fluid pressure distributions

The presence of fault-plane reflections in seismic images, besides indicating the locations of faults, offers a possible source of information on the properties of these poorly understood zones. To better understand the physical mechanism giving rise to fault-plane reflections in compacting sedimentary basins, we numerically model the full elastic wavefield via the spectral element method (SEM) for several different fault models. Using well log data from the South Eugene Island field, offshore Louisiana, we derive empirical relationships between the elastic parameters (e.g. P-wave velocity and density) and the effective-stress along both normal compaction and unloading paths. These empirical relationships guide the numerical modelling and allow the investigation of how differences in fluid pressure modify the elastic wavefield. We choose to simulate the elastic wave equation via SEM since irregular model geometries can be accommodated and slip boundary conditions at an interface, such as a fault or fracture, are implemented naturally. The method we employ for including a slip interface retains the desirable qualities of SEM in that it is explicit in time and, therefore, does not require the inversion of a large matrix. We performa complete numerical study by forward modelling seismic shot gathers over a faulted earth model using SEM followed by seismic processing of the simulated data. With this procedure, we construct post-stack time-migrated images of the kind that are routinely interpreted in the seismic exploration industry. We dip filter the seismic images to highlight the fault-plane reflections prior to making amplitude maps along the fault plane. With these amplitude maps, we compare the reflectivity from the different fault models to diagnose which physical mechanism contributes most to observed fault reflectivity. To lend physical meaning to the properties of a locally weak fault zone characterized as a slip interface, we propose an equivalent-layer model under the assumption of weak scattering. This allows us to use the empirical relationships between density, velocity and effective stress from the South Eugene Island field to relate a slip interface to an amount of excess pore-pressure in a fault zone. ?? 2007 The Authors Journal compilation ?? 2007 RAS.

Geophysical Journal International↗
USGS2022

Orbital and in-situ investigation of periodic bedrock ridges in Glen Torridon, Gale Crater, Mars

Wind has been the dominant agent of landscape modification on Mars for the past ~3 billion years. Among the diversity of features formed by aeolian abrasion on the surface of Mars are periodic bedrock ridges (PBRs), a relatively recently recognized class of erosional bedforms on Mars for which Earth analogues are rare. Gale crater, the field site for NASA’s Mars Science Laboratory Curiosity rover since it landed there in 2012, contains a diverse and extensive record of aeolian deposition and erosion. This study focuses on a series of periodic, linear bedrock ridges that occur within the Fe/Mg-smectite clay-bearing Glen Torridon region of Aeolis Mons (informally Mount Sharp). During Curiosity’s exploration of the Glen Torridon region between sols ~2300-3080, the rover drove through this field of ridges, providing the first opportunity for the in situ observation of these enigmatic erosional features. This study characterizes the Glen Torridon ridges using orbiter and rover data to determine their morphology, spatial distribution, compositional and material properties, and association with other aeolian features in the area. Based on these observations, the Glen Torridon ridges are interpreted to be consistent with an origin as wind-eroded periodic bedrock ridges carved during the most recent exhumation of Mount Sharp into the present-day mound. Although there is evidence for multidirectional winds in the Glen Torridon region based on the orientation of modern ripples, megaripples, TARs and other bedrock indicators, the consistent orientation of the Glen Torridon ridges, coupled with morphologic asymmetries within the ridges, support formation and elongation of the Glen Torridon PBRs forms parallel to a net regional northerly wind direction in and around Gale crater.

Journal of Geophysical Research E: Planets↗
USGS2020

Evidence for a diagenetic origin of Vera Rubin Ridge, Gale Crater, Mars: Summary and synthesis of Curiosity's exploration campaign

This paper provides an overview of the Curiosity rover's exploration at Vera Rubin ridge (VRR) and summarizes the science results. VRR is a distinct geomorphic feature on lower Aeolis Mons (informally known as Mount Sharp) that was identified in orbital data based on its distinct texture, topographic expression, and association with a hematite spectral signature. Curiosity conducted extensive remote sensing observations, acquired data on dozens of contact science targets, and drilled three outcrop samples from the ridge, as well as one outcrop sample immediately below the ridge. Our observations indicate that strata composing VRR were deposited in a predominantly lacustrine setting and are part of the Murray formation. The rocks within the ridge are chemically in family with underlying Murray formation strata. Red hematite is dispersed throughout much of the VRR bedrock, and this is the source of the orbital spectral detection. Gray hematite is also present in isolated, gray‐colored patches concentrated toward the upper elevations of VRR, and these gray patches also contain small, dark Fe‐rich nodules. We propose that VRR formed when diagenetic event(s) preferentially hardened rocks, which were subsequently eroded into a ridge by wind. Diagenesis also led to enhanced crystallization and/or cementation that deepened the ferric‐related spectral absorptions on the ridge, which helped make them readily distinguishable from orbit. Results add to existing evidence of protracted aqueous environments at Gale crater and give new insight into how diagenesis shaped Mars' rock record.

Journal of Geophysical Research: Planets↗
USGS2015

Source mechanism of small long-period events at Mount St. Helens in July 2005 using template matching, phase-weighted stacking, and full-waveform inversion

Long-period (LP, 0.5-5 Hz) seismicity, observed at volcanoes worldwide, is a recognized signature of unrest and eruption. Cyclic LP “drumbeating” was the characteristic seismicity accompanying the sustained dome-building phase of the 2004–2008 eruption of Mount St. Helens (MSH), WA. However, together with the LP drumbeating was a near-continuous, randomly occurring series of tiny LP seismic events (LP “subevents”), which may hold important additional information on the mechanism of seismogenesis at restless volcanoes. We employ template matching, phase-weighted stacking, and full-waveform inversion to image the source mechanism of one multiplet of these LP subevents at MSH in July 2005. The signal-to-noise ratios of the individual events are too low to produce reliable waveform-inversion results, but the events are repetitive and can be stacked. We apply network-based template matching to 8 days of continuous velocity waveform data from 29 June to 7 July 2005 using a master event to detect 822 network triggers. We stack waveforms for 359 high-quality triggers at each station and component, using a combination of linear and phase-weighted stacking to produce clean stacks for use in waveform inversion. The derived source mechanism pointsto the volumetric oscillation (~10 m 3 ) of a subhorizontal crack located at shallow depth (~30 m) in an area to the south of Crater Glacier in the southern portion of the breached MSH crater. A possible excitation mechanism is the sudden condensation of metastable steam from a shallow pressurized hydrothermal system as it encounters cool meteoric water in the outer parts of the edifice, perhaps supplied from snow melt.

Washington↗
USGS2009

Development of time series stacks of Landsat images for reconstructing forest disturbance history

Forest dynamics is highly relevant to a broad range of earth science studies, many of which have geographic coverage ranging from regional to global scales. While the temporally dense Landsat acquisitions available in many regions provide a unique opportunity for understanding forest disturbance history dating back to 1972, large quantities of Landsat images will need to be analysed for studies at regional to global scales. This will not only require effective change detection algorithms, but also highly automated, high level preprocessing capabilities to produce images with subpixel geolocation accuracies and best achievable radiometric consistency, a status called imagery-ready-to-use (IRU). This paper describes a streamlined approach for producing IRU quality Landsat time series stacks (LTSS). This approach consists of an image selection protocol, high level preprocessing algorithms and IRU quality verification procedures. The high level preprocessing algorithms include updated radiometric calibration and atmospheric correction for calculating surface reflectance and precision registration and orthorectification routines for improving geolocation accuracy. These automated routines have been implemented in the Landsat Ecosystem Disturbance Adaptive System (LEDAPS) designed for processing large quantities of Landsat images. Some characteristics of the LTSS developed using this approach are discussed.

International Journal of Digital Earth↗
USGS2000

Temporal and spectral characteristics of seismicity observed at Popocatepetl volcano, central Mexico

Popocatepetl volcano entered an eruptive phase from December 21, 1994 to March 30, 1995, which was characterized by ash and fumarolic emissions. During this eruptive episode, the observed seismicity consisted of volcano-tectonic (VT) events, long-period (LP) events and sustained tremor. Before the initial eruption on December 21, VT seismicity exhibited no increase in number until a swarm of VT earthquakes was observed at 01:31 hours local time. Visual observations of the eruption occurred at dawn the next morning. LP activity increased from an average of 7 events a day in October 1994 to 22 events per day in December 1994. At the onset of the eruption, LP activity peaked at 49 events per day. LP activity declined until mid-January 1995 when no events were observed. Tremor was first observed about one day after the initial eruption and averaged 10 h per episode. By late February 1995, tremor episodes became more intermittent, lasting less than 5 min, and the number of LP events returned to pre-eruption levels (7 events per day). Using a spectral ratio technique, low-frequency oceanic microseismic noise with a predominant peak around 7 s was removed from the broadband seismic signal of tremor and LP events. Stacks of corrected tremor episodes and LP events show that both tremor and LP events contain similar frequency features with major peaks around 1.4 Hz. Frequency analyses of LP events and tremor suggest a shallow extended source with similar radiation pattern characteristics. The distribution of VT events (between 2.5 and 10 km) also points to a shallow source of the tremor and LP events located in the first 2500 m beneath the crater. Under the assumption that the frequency characteristics of the signals are representative of an oscillator we used a fluid-filled-crack model to infer the length of the resonator.

Journal of Volcanology and Geothermal Research↗