Search USGSSearch

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Tectonics of the western Valley and Ridge foldbelt, Pendleton County, West Virginia - a summary report

A belt of high anticlines, the Nittany anticlinorium, occupies the western Valley and Ridge foldbelt in the central Appalachians. It extends southwestward from the Nittany arch of central Pennsylvania into the Virginias. An investigation of the tectonics of this anticlinorium in Pendleton County, W. Va., rules out active basement involvement in the deformation of the area. Cross-sectional models consistent with the accumulated data show that Middle Cambrian through Middle Ordovician carbonate rocks are technically stacked, shingle-fashion, from southeast to northwest below predominantly folded younger strata that have undergone less lateral shortening. Differential shortening in this area is of the proper order to balance cover deformation in the Allegheny synclinorium to the west. Field relations suggest a long period of abnormally high fluid pressures in Lower Devonian and older strata during deformation. At this time, the area was under sufficient northwest, near-horizontal compressive stress for abundant quartz deformation lamellae to form. Gravity sliding is ruled out as the deforming mechanism for this part of the Appalachian foldbelt. No significant tectonism appears to have occurred prior to Pennsylvanian time in this area.

West Virginia

Inversion of seismic refraction data in planar dipping structure

A new method is presented for the direct inversion of seismic refraction data in dipping planar structure. Three recording geometries, each consisting of two common-shot profiles, are considered: reversed, split, and roll-along profiles. Inversion is achieved via slant stacking the common-shot wavefield to obtain a delay time—slowness (tau— p ) wavefield. The tau— p curves from two shotpoints describing the critical raypath of refracted and post-critically reflected arrivals are automatically picked using coherency measurements and the two curves are jointly used to calculate velocity and dip of isovelocity lines iteratively, thereby obtaining the final two-dimensional velocity model. This procedure has been successfully applied to synthetic seismograms calculated for a dipping structure and to field data from central California. The results indicate that direct inversion of closely-spaced refraction/wide-aperture reflection data can practically be achieved in laterally inhomogeneous structures.

Geophysical Journal International

Wind shear stress measurements in a coastal marsh during Hurricane Andrew

Hurricane Andrew produced changes to the Louisiana wetlands not normally observed after lesser, more common storms. For example, the <25 m/s wind speeds generated by cold fronts and winter storms, and any accompanying storm surge, do not cause substantial, wide-spread alteration of marsh vegetation. During Hurricane Andrew, however, the wind, the wind-driven storm surge, or both produced severe, wide-spread wetland alteration, especially in areas that primarily consisted of densely vegetated floating mats. In a few hours, vegetated brackish marsh was severely torn and large areas were converted to open water, a process that takes decades when driven by geologic subsidence, human intervention, and lesser storms. During the passage of Hurricane Andrew, wind measurements were taken inside an impoundment within a brackish part of Louisiana's coastal wetlands system. At its closest point, the site lay 50 km to the right of the north-trending storm track, placing it in or near the zone of maximum wind (the eye wall). As the hurricane approached, the wind blew from the north; after it passed, the wind direction swung around to the southeast. Several hours after the eye passed, the southeasterly wind drove a 1.5-m storm surge through the area, causing the collapse of the meteorology tower. Wind shear stress calculations, based on data from two vertically stacked sensors, showed a direct correlation between wind shear stress and wind speed. The greatest increase in wind shear stress occurred when wind speed exceeded 20 m/s. Overall, wind shear stress increased more than three orders of magnitude — from approximately 0.01 N/m² at a wind speed of 6 m/s through 1 N/m² at 20 m/s to 18 N/m² at the maximum sustained speed of 43 m/s. Drag-coefficient calculations show that the open-ocean CD formulations, such as the popular WAMDI model, cannot be employed for wetland use because it overestimates CD for velocities less than approximately 20 m/s and underestimates it for higher velocities.

Louisiana

Kaguyak dome field and its Holocene caldera, Alaska Peninsula

Kaguyak Caldera lies in a remote corner of Katmai National Park, 375 km SW of Anchorage, Alaska. The 2.5-by-3-km caldera collapsed ~ 5.8 ± 0.2 ka ( 14 C age) during emplacement of a radial apron of poorly pumiceous crystal-rich dacitic pyroclastic flows (61–67% SiO 2 ). Proximal pumice-fall deposits are thin and sparsely preserved, but an oxidized coignimbrite ash is found as far as the Valley of Ten Thousand Smokes, 80 km southwest. Postcaldera events include filling the 150-m-deep caldera lake, emplacement of two intracaldera domes (61.5–64.5% SiO 2 ), and phreatic ejection of lakefloor sediments onto the caldera rim. CO 2 and H 2 S bubble up through the lake, weakly but widely. Geochemical analyses ( n = 148), including pre-and post-caldera lavas (53–74% SiO 2 ), define one of the lowest-K arc suites in Alaska. The precaldera edifice was not a stratocone but was, instead, nine contiguous but discrete clusters of lava domes, themselves stacks of rhyolite to basalt exogenous lobes and flows. Four extra caldera clusters are mid-to-late Pleistocene, but the other five are younger than 60 ka, were truncated by the collapse, and now make up the steep inner walls. The climactic ignimbrite was preceded by ~ 200 years by radial emplacement of a 100-m-thick sheet of block-rich glassy lava breccia (62–65.5% SiO 2 ). Filling the notches between the truncated dome clusters, the breccia now makes up three segments of the steep caldera wall, which beheads gullies incised into the breccia deposit prior to caldera formation. They were probably shed by a large lava dome extruding where the lake is today.

Alaska

Depth determination of the 2010 El Mayor‐Cucapah earthquake sequence (M ≥ 4.0)

The 2010 M W 7.2 El Mayor‐Cucapah earthquake ruptured a zone of ~120 km in length in northern Baja California. The geographic distribution of this earthquake sequence was well constrained by waveform relocation. The depth distribution, however, was poorly determined as it is near the edge of, or outside, the Southern California Seismic Network. Here we use two complementary methods to constrain the focal depths of moderate‐sized events ( M ≥ 4.0) in this sequence. We first determine the absolute earthquake depth by modeling the regional depth phases at high frequencies (~1 Hz). We mainly focus on Pn and its depth phases pPn and sPn , which arrive early at regional distance and are less contaminated by crustal multiples. To facilitate depth phase identification and to improve signal‐to‐noise ratio, we take advantage of the dense Southern California Seismic Network and use array analysis to align and stack Pn waveforms. For events without clear depth phases, we further determine their relative depths with respect to those with known depths using differential travel times of the Pn , direct P , and direct S phases recorded for event pairs. Focal depths of 93 out of 122 M ≥ 4.0 events are tightly constrained with absolute uncertainty of about 1 km. Aftershocks are clustered in the depth range of 3–10 km, suggesting a relatively shallow seismogenic zone, consistent with high surface heat flow in this region. Most aftershocks are located outside or near the lower terminus of coseismic high‐slip patches of the main shock, which may be governed by residual strains, local stress concentration, or postseismic slip.

Baja California