USGS Science⌕ Search

SEARCH · USGS Science

Results for “One Earth”

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 973 records · Page 54Linked to original sources
USGS2000

Annual bed-elevation regime in the alluvial channel of Squamish River, southwestern British Columbia Canada

The aim of this study is to examine the annual regime of channel scour and fill by monitoring bed-elevation changes in a reach of Squamish River in southwestern British Columbia, Canada. Sonar surveys of 13 river cross-sections in a sandy gravel-bed single-channel study reach were repeated biweekly over a full hydrologic year (1995/6). The survey results show that bedload movement occurs as waves or pulses forming bedwaves that appear to maintain an overall coherence with movement downstream. These bedwaves propagate downstream by a mode here termed pulse scour and pulse fill, a process distinguished from the conventional mode of scour and fill commonly associated with flood events (here termed local scour and local fill). Bedwave celerity was estimated to be about 15.5 m d-1 corresponding to a bedwave residence time in the study reach of almost one hydrologic year. The total amount of local bed-elevation change ranged between 0.22 m and 2.41 m during the period of study. Analysis of the bed-elevation and flow data reveals that, because of the bedware phenomenon, there is no simple relation between the mean bed-elevation and discharge nor any strong linear correlation among cross-sectional behaviour. The bed-elevation data also suggest that complex changes to the bed within a cross-section are masked when the bed is viewed in one dimension, although no definitive trends in bed behaviour were found in the two-dimensional analysis. Although a weak seasonal effect is evident in this study, the bed-elevation regime is dominated by sediment supply-driven fluctuations in bedload transport occurring at timescales shorter than the seasonal fluctuation in discharge. The study also indicates that bed-elevation monitoring on Squamish River, and others like it, for purposes of detecting and measuring aggradation/degradation must take into account very considerable and normal channel-bed variability operating at timescales from hours to months. Copyright (C) 2000 John Wiley and Sons, Ltd.

Earth Surface Processes and Landforms↗
USGS2025

Volcanic gases reflect magma stalling and launching depths

Many open-vent arc volcanoes display two modes in their continuous gas emissions, one with a characteristic CO 2 / S T ratio typical of periods of quiescent degassing and another punctuated by high CO 2 / S T gas emitted in the weeks before eruption, a recently recognized eruption precursor. In this study we explore the origin of the two modes of degassing revealed by time-series gas data at Turrialba volcano (Costa Rica) in the context of new melt inclusion (MI) data. To reconstruct the c[CO 2 ] of undegassed magma, we developed a rapid-quench piston-cylinder assembly to rehomogenize the vapor bubble commonly contained in MIs. We focus on olivine-hosted MIs from a mafic scoria sample erupted from Turrialba in 1864–1866. The reconstructed CO 2 contents in MIs decrease from ∼4000 to <1000 ppmw as S contents decrease from 3500 to <1000 ppmw. The highest reconstructed S and CO 2 in the MIs resulted in an initial magmatic CO 2 / S T ratio (molar) of 0.83. Informed by the MI data, we modeled the decompression degassing of Turrialba magma and vapor composition using the Sulfur_X and EVo models. Instead of being controlled by initial magmatic CO 2 /S T ratio as suggested by previous studies, we find that the quiescent gas emitted from Turrialba during 2014–2018 (CO 2 / S T = 2.3 ± 0.8, molar) appears to reflectequilibrium with magmas stored at 4–8 km (Sulfur_X) or 2 km (EVo) depth, when H 2 O is degassing extensively from the magma. A magma storage region at 4–8 km is also supported by seismic tomography. The second gas mode is noted by spikes in CO 2 / S T ∼ 7.9 ± 2 in the weeks prior to eruption. This gas reflects equilibrium with a magma at 12–18 km (Sulfur_X) or 4–8 km (EVo), where the ascending magma is saturated with a CO 2 -rich vapor. Thus, there are two important trans crustal depths beneath the volcano: one where the rate of H 2 O loss from the magma and thus magma viscosity increases, and one at greater depths where high CO 2 /S T vapor forms and may facilitate dike propagation. We interpret the shallower, H 2 O-loss region as the main site of magma stalling and storage, where quiescent gas is generated continuously. We interpret the greater depth (12–18 km) as the source of the precursory gas that precedes eruption, and where the mafic melt lastly equilibrated with a mush zone before ascending and triggering eruption weeks later. This hypothesis is ripe for testing at other volcanoes that exhibit two modes in gas geochemistry.

Turrialba volcano↗
USGS2000

Fracture process zone in granite

In uniaxial compression tests performed on Aue granite cores (diameter 50 mm, length 100 mm), a steel loading plate was used to induce the formation of a discrete shear fracture. A zone of distributed microcracks surrounds the tip of the propagating fracture. This process zone is imaged by locating acoustic emission events using 12 piezoceramic sensors attached to the samples. Propagation velocity of the process zone is varied by using the rate of acoustic emissions to control the applied axial force. The resulting velocities range from 2 mm/s in displacement-controlled tests to 2 μm/s in tests controlled by acoustic emission rate. Wave velocities and amplitudes are monitored during fault formation. P waves transmitted through the approaching process zone show a drop in amplitude of 26 dB, and ultrasonic velocities are reduced by 10%. The width of the process zone is ∼9 times the grain diameter inferred from acoustic data but is only 2 times the grain size from optical crack inspection. The process zone of fast propagating fractures is wider than for slow ones. The density of microcracks and acoustic emissions increases approaching the main fracture. Shear displacement scales linearly with fracture length. Fault plane solutions from acoustic events show similar orientation of nodal planes on both sides of the shear fracture. The ratio of the process zone width to the fault length in Aue granite ranges from 0.01 to 0.1 inferred from crack data and acoustic emissions, respectively. The fracture surface energy is estimated from microstructure analysis to be ∼2 J. A lower bound estimate for the energy dissipated by acoustic events is 0.1 J.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Regional deformation near Palmdale, California, 1973-1983 (USA)

The Tehachapi trilateration network spans the intersection of the San Andreas and Garlock faults in southern California in the “Big Bend” region of the San Andreas fault. Analysis of data from 1973–1983 shows strain differences between the northwest and southeast regions of the network and slip at depth on both faults. The Palmdale network, spanning the San Andreas fault entirely within the Tehachapi network, showed increases of about 1 μstrain in both east‐west and north‐south extension in late 1979. The Tehachapi strains also jumped at this time, but the magnitude of the increase was only about one third that of Palmdale. The principal strain rates for Tehachapi over the time interval 1973–1983 are μstrain/yr and μstrain/yr, with the 1 axis directed N76°E. Strains were also computed for two subregions. The principal strain rates for the southeast Tehachapi region are μstrain/yr and μStrain/yr, with the 1 axis directed N73°E. This result differs significantly from the principal strain rates at Palmdale ( μstrain/yr, μstrain/yr, with the 1 axis directed N71°E). For the northwest Tehachapi region, μstrain/yr and μstrain/yr, with the 1 axis directed N93°E.

California↗
USGS2021

Dating by cosmogenic nuclides

Since the 1990s, cosmogenic nuclides have revolutionized the study of Earth surface processes, particularly the understanding of rates and dates. These nuclides, including 3 He, 10 Be, 14 C, 21 Ne, 26 Al, and 36 Cl, enable dating of landforms and the measurement of erosion rates both at the scale of drainage basins and at specific locations on Earth's surface. Cosmogenic nuclides are produced at low rates (several to hundreds of atoms per gram per year) by the interaction of cosmic rays with elements both in the atmosphere and in surficial materials, including in rock and soil. Measuring nuclide concentrations requires elemental separation from source geologic material followed by counting of atoms using sensitive accelerator mass spectrometers. Because nuclide production rates have been quantified, the measured concentration of these nuclides can be interpreted as a near-surface residence time. Here, we review the systematics of commonly used cosmogenic nuclides, describe how they are extracted and measured, and then present case studies focusing on the most commonly measured cosmogenic nuclide, 10 Be. We present common applications such as dating surface features, including moraines and outcrops shaped by glaciation, the use of cosmogenic nuclides for inferring tectonic and erosion processes in drainage basins, and the use of these nuclides to trace sediment sources in drainage basins. When multiple nuclides are measured in one sample, they can be used to model burial and exposure histories in stratigraphic sections. We conclude by exploring what the future might bring in terms of measurements and applications.

Book chapter↗
USGS2013

Remote sensing survey of Chinese tallow tree in the Toledo Bend Reservoir area, Louisiana and Texas

We applied Hyperion sensor satellite data acquired by the National Aeronautics and Space Administration’s Earth Observing-1 (EO-1) satellite in conjunction with reconnaissance surveys to map the occurrences of the invasive Chinese tallow tree (Triadica sebifera) in the Toledo Bend Reservoir study area of northwestern Louisiana and northeastern Texas. The rationale for application of high spectral resolution EO-1 Hyperion data was based on the successful use of Hyperion data in the mapping of Chinese tallow tree in southwestern Louisiana in 2005. In contrast to the single Hyperion image used in the 2005 project, more than 20 EO-1 Hyperion and Advanced Land Imager (ALI) images of the study area were collected in 2009 and 2010 during the fall senescence when Chinese tallow tree leaves turn red. Atmospherically corrected reflectance spectra of Hyperion imagery collected at ground and aerial observation locations provided the input datasets used in the program for spectral discrimination analysis. Discrimination analysis was used to identify spectral indicator sets to best explain variance contained in the input databases. The expectation was that at least one set of Hyperion-based indicator spectra would uniquely identify occurrences of red-leaf Chinese tallow tree; however, no combination of Hyperion-based reflectance datasets produced a unique identifier. The inability to discover a unique spectral indicator resulted primarily from relatively sparse coverage by red-leaf Chinese tallow tree within the study area (percentage of coverage was less than 5 percent per 30- by 30-meter Hyperion pixel). To enhance the performance of the spectral discrimination analysis, leaf and canopy spectra of Chinese tallow tree were added to the input datasets to guide the indicator selection. In addition, input databases were segregated by land class obtained from an ALI-based landcover classification in order to reduce the input variance and to promote spectral discrimination of red-leaf Chinese tallow tree. Although no unique spectral identifier for red-leaf Chinese tallow tree was uncovered with these enhanced methods, in some cases predicted spatial patterns throughout the Hyperion images revealed alignment with vegetation associations within each land class that was often observed to contain Chinese tallow trees. These instances were associated particularly with the addition of helicopter-based spectra to the input databases. It was attempted to extend such predictions of likely occurrences of Chinese tallow tree by mapping six of the nine Hyperion swaths and four of the nine land classes, but this attempt produced uncertain results that could not be fully evaluated for accuracy. Even though the final mapping showed promise in identifying likely Chinese tallow tree occurrences, the low percentage of occurrences hindered mapping performance and validation. Results of the mapping suggested that successful detection of Chinese tallow tree in the study area would require a spectral sensor similar to the Hyperion but with a higher ground-level spatial resolution. Although the Hyperion-based spectral mapping did not provide the desired results, the associated field (ground and aerial) surveys did provide for a qualitative assessment of the overall Chinese tallow tree distribution within the study area. Ground and aerial surveys suggested that Chinese tallow tree occurrences were uncommon and were without an observed pattern in relation to proximity to the Toledo Bend Reservoir. Although uncommon and scattered, Chinese tallow trees and shrubs most commonly existed along forest edges, water edges, and fence lines, probably most in line with seed dispersal by birds. Chinese tallow trees were observed to be more densely dispersed within some scrublands and grasslands than were observed in pine, hardwood, and mixed forests.

Louisiana;Texas↗
USGS2017

Soils as relative-age dating tools

Soils develop at the earth's surface via multiple processes that act through time. Precluding burial or disturbance, soil genetic horizons form progressively and reflect the balance among formation processes, surface age, and original substrate composition. Soil morphology provides a key link between process and time (soil age), enabling soils to serve as both relative and numerical dating tools for geomorphic studies and landscape evolution. Five major factors define the contemporary state of all soils: climate, organisms, topography, parent material, and time. Soils developed on similar landforms and parent materials within a given landscape comprise what we term a soil/landform/substrate complex. Soils on such complexes that differ in development as a function of time represent a soil chronosequence. In a soil chronosequence, time constitutes the only independent formation factor; the other factors act through time. Time dictates the variations in soil development or properties (field or laboratory measured) on a soil/landform/substrate complex. Using a dataset within the chronosequence model, we can also formulate various soil development indices based upon one or a combination of soil properties, either for individual soil horizons or for an entire profile. When we evaluate soil data or soil indices mathematically, the resulting equation creates a chronofunction. Chronofunctions help quantify processes and mechanisms involved in soil development, and relate them mathematically to time. These rigorous kinds of comparisons among and within soil/landform complexes constitute an important tool for relative-age dating. After determining one or more absolute ages for a soil/landform complex, we can calculate quantitative soil formation, and or landform-development rates. Multiple dates for several complexes allow rate calculations for soil/landform-chronosequence development and soil-chronofunction calibration.

Book chapter↗
USGS1969

Some triple-filament lead isotope ratio measurements and an absolute growth curve for single-stage leads

Triple-filament analyses of three standard lead samples are used to calibrate a mass spectrometer in an absolute sense. The bias we measure is 0.0155 percent per mass unit, and the precision (for 95% confidence limits) is ±0.13% or less for all ratios relative to 204 Pb. Although its precision is not quite so good as that of the lead-tetramethyl method in the analysis of large samples, the triple-filament method is less complex and is an attractive alternative for smaller sample sizes down to 500 μg. Triple-filament data are presented for six possibly single-stage lead ores and one feldspar. These new data for ores are combined with corrected tetramethyl data for stratiform lead deposits to compute absolute parameters for a universal single-stage lead isotope growth curve. Absolute isotopic ratios for primeval lead have been determined by Oversby and because all the previous data for both meteorites and lead ores were similarly fractionated, the absolute value of 238 U 204 Pb= 9.09 ± 0.06 for stratiform leads is little different from the value 8.99 ± 0.05 originally computed by Ostic, Russell and Stanton. Absolute values for lead isotope ratios for all interlaboratory standard samples presently available from the literature are tabulated.

Earth and Planetary Science Letters↗
USGS2014

Post-earthquake relaxation using a spectral element method: 2.5-D case

The computation of quasi-static deformation for axisymmetric viscoelastic structures on a gravitating spherical earth is addressed using the spectral element method (SEM). A 2-D spectral element domain is defined with respect to spherical coordinates of radius and angular distance from a pole of symmetry, and 3-D viscoelastic structure is assumed to be azimuthally symmetric with respect to this pole. A point dislocation source that is periodic in azimuth is implemented with a truncated sequence of azimuthal order numbers. Viscoelasticity is limited to linear rheologies and is implemented with the correspondence principle in the Laplace transform domain. This leads to a series of decoupled 2-D problems which are solved with the SEM. Inverse Laplace transform of the independent 2-D solutions leads to the time-domain solution of the 3-D equations of quasi-static equilibrium imposed on a 2-D structure. The numerical procedure is verified through comparison with analytic solutions for finite faults embedded in a laterally homogeneous viscoelastic structure. This methodology is applicable to situations where the predominant structure varies in one horizontal direction, such as a structural contrast across (or parallel to) a long strike-slip fault.

Geophysical Journal International↗
USGS2000

Mesoproterozoic graphite deposits, New Jersey Highlands: Geologic and stable isotopic evidence for possible algal origins

Graphite deposits of Mesoproterozoic age are locally abundant in the eastern New Jersey Highlands, where they are hosted by sulphidic biotite–quartz–feldspar gneiss, metaquartzite, and anatectic pegmatite. Gneiss and metaquartzite represent a shallow marine shelf sequence of locally organic-rich sand and mud. Graphite from massive deposits within metaquartzite yielded δ13C values of –26 ± 2‰ (1σ), and graphite from massive deposits within biotite–quartz–feldspar gneiss yielded δ13C values of –23 ± 4‰. Disseminated graphite from biotite–quartz–feldspar gneiss country rock was –22 ± 3‰, indistinguishable from the massive deposits hosted by the same lithology. Anatectic pegmatite is graphitic only where generated from graphite-bearing host rocks; one sample gave a δ13C value of –15‰. The δ34S values of trace pyrrhotite are uniform within individual deposits, but vary from 0 to 9‰ from one deposit to another. Apart from pegmatitic occurrences, evidence is lacking for long-range mobilization of carbon during Grenvillian orogenesis or post-Grenvillian tectonism. The field, petrographic, and isotope data suggest that massive graphite was formed by granulite-facies metamorphism of Proterozoic accumulations of sedimentary organic matter, possibly algal mats. Preservation of these accumulations in the sedimentary environment requires anoxic basin waters or rapid burial. Anoxia would also favour the accumulation of dissolved ferrous iron in basin waters, which may explain some of the metasediment-hosted massive magnetite deposits in the New Jersey Highlands.

New Jersey↗
USGS2002

Rare earth element variations resulting from inversion of pigeonite and subsolidus reequilibration in lunar ferroan anorthosites

We present results of a secondary ion mass spectrometry study of the rare earth elements (REEs) in the minerals of two samples of lunar ferroan anorthosite, and the results are applicable to studies of REEs in all igneous rocks, no matter what their planet of origin. Our pyroxene analyses are used to determine solid-solid REE distribution coefficients (D = C REE in low-Ca pyroxene/C REE in augite) in orthopyroxene-augite pairs derived by inversion of pigeonite. Our data and predictions from crystal-chemical considerations indicate that as primary pigeonite inverts to orthopyroxene plus augite and subsolidus reequilibration proceeds, the solid-solid Ds for orthopyroxene-augite pairs progressively decrease for all REEs; the decrease is greatest for the LREEs. The REE pattern of solid-solid Ds for inversion-derived pyroxene pairs is close to a straight line for Sm-Lu and turns upward for REEs lighter than Sm; the shape of this pattern is predicted by the shapes of the REE patterns for the individual minerals. Equilibrium liquids calculated for one sample from the compositions of primary phases, using measured or experimentally determined solid-liquid Ds, have chondrite-normalized REE patterns that are very slightly enriched in LREEs. The plagioclase equilibrium liquid is overall less rich in REEs than pyroxene equilibrium liquids, and the discrepancy probably arises because the calculated plagioclase equilibrium liquid represents a liquid earlier in the fractionationsequence than the pyroxene equilibrium liquids. “Equilibrium” liquids calculated from the compositions of inversion-derived pyroxenes or orthopyroxene derived by reaction of olivine are LREE depleted (in some cases substantially) in comparison with equilibrium liquids calculated from the compositions of primary phases. These discrepancies arise because the inversion-derived and reaction-derived pyroxenes did not crystallize directly from liquid, and the use of solid-liquid Ds is inappropriate. The LREE depletion of the calculated liquids is a relic of formation of these phases from primary LREE-depleted minerals. Thus, if one attempts to calculate the compositions of equilibrium liquids from pyroxene compositions, it is important to establish that the pyroxenes are primary. In addition, our data suggest that experimental studies have underestimated solid-liquid Ds for REEs in pigeonite and that REE contents of liquids calculated using these Ds are overestimates. Our results have implications for Sm-Nd age studies. Our work shows that if pigeonite inversion and/or subsolidus reequilibration between augite and orthopyroxene occurred significantly after crystallization, and if pyroxene separates isolated for Sm-Nd studies do not have the bulk composition of the primary pyroxenes, then the Sm-Nd isochron age and ε Nd will be in error.

Geochimica et Cosmochimica Acta↗
USGS2004

Quantifying probabilities of volcanic events: The example of volcanic hazard at Mount Vesuvius

We describe an event tree scheme to quantitatively estimate both long- and short-term volcanic hazard. The procedure is based on a Bayesian approach that produces a probability estimation of any possible event in which we are interested and can make use of all available information including theoretical models, historical and geological data, and monitoring observations. The main steps in the procedure are (1) to estimate an a priori probability distribution based upon theoretical knowledge, (2) to modify that using past data, and (3) to modify it further using current monitoring data. The scheme allows epistemic and aleatoric uncertainties to be dealt with in a formal way, through estimation of probability distributions at each node of the event tree. We then describe an application of the method to the case of Mount Vesuvius. Although the primary intent of the example is to illustrate the methodology, one result of this application merits special mention. The present emergency response plan for Mount Vesuvius is referenced to a maximum expected event (MEE), the largest out of all the possible eruptions within the next few decades. Our calculation suggest that there is a nonnegligible (1-20%) chance that the next eruption could be larger than that stipulated in the present MEE. The methodology allows all assumptions and thresholds to be clearly identified and provides a rational means for their revision if new data or information are obtained. Copyright 2004 by the American Geophysical Union.

Mount Vesuvius↗
USGS2022

Quantifying large-scale surface change using SAR amplitude images: Crater morphology changes during the 2019-2020 Shishaldin Volcano eruption

Morphological processes often induce meter-scale elevation changes. When a volcano erupts, tracking such processes provides insights into the style and evolution of eruptive activity and related hazards. Compared to optical remote-sensing products, synthetic aperture radar (SAR) observes surface change during inclement weather and at night. Differential SAR interferometry estimates phase change between SAR acquisitions and is commonly applied to quantify deformation. However, large deformation or other coherence loss can limit its use. We develop a new approach applicable when repeated digital elevation models (DEMs) cannot be otherwise retrieved. Assuming an isotropic radar cross-section, we estimate meter-scale vertical morphological change directly from SAR amplitude images via an optimization method that utilizes a high-quality DEM. We verify our implementation through simulation of a collapse feature that we modulate onto topography. We simulate radar effects and recover the simulated collapse. To validate our method, we estimate elevation changes from TerraSAR-X stripmap images for the 2011–2012 eruption of Mount Cleveland. Our results reproduce those from two previous studies; one that used the same dataset, and another based on thermal satellite data. By applying this method to the 2019–2020 eruption of Shishaldin Volcano, Alaska, we generate elevation change time series from dozens of co-registered TerraSAR-X high-resolution spotlight images. Our results quantify previously unresolved cone growth in November 2019, collapses associated with explosions in December–January, and further changes in crater elevations into spring 2020. This method can be used to track meter-scale morphology changes for ongoing eruptions with low latency as SAR imagery becomes available.

Alaska↗
USGS2001

U-Pb age of the Diana Complex and Adirondack granulite petrogenesis

U-Pb isotopic analyses of eight single and multi-grain zircon fractions separated from a syenite of the Diana Complex of the Adirondack Mountains do not define a single linear array, but a scatter along a chord that intersects the Concordia curve at 1145 ± 29 and 285 ± 204 Ma. For the most concordant analyses, the 207 Pb/ 206 Pb ages range between 1115 and 1150 Ma. Detailed petrographic studies revealed that most grains contained at least two phases of zircon growth, either primary magmatic cores enclosed by variable thickness of metamorphic overgrowths or magmatic portions enclosing presumably older xenocrystic zircon cores. The magmatic portions are characterized by typical dipyramidal prismatic zoning and numerous black inclusions that make them quite distinct from adjacent overgrowths or cores when observed in polarizing light microscopy and in backscattered electron micrographs. Careful handpicking and analysis of the “best” magmatic grains, devoid of visible overgrowth of core material, produced two nearly concordant points that along with two of the multi-grain analyses yielded an upper-intercept age of 1118 ± 2.8 Ma and a lowerintercept age of 251 ± 13 Ma. The older age is interpreted as the crystallization age of the syenite and the younger one is consistent with late stage uplift of the Appalachian region. The 1118 Ma age for the Diana Complex, some 35 Ma younger than previously believed, is now approximately synchronous with the main Adirondack anorthosite intrusion, implying a cogenetic relationship among the various meta-igneous rocks of the Adirondacks. The retention of a high-temperature contact metamorphic aureole around Diana convincingly places the timing of Adirondack regional metamorphism as early as 1118 Ma. This result also implies that the sources of anomalous hightemperature during granulite metamorphism are the syn-metamorphic intrusions, such as the Diana Complex.

New York↗
USGS2011

Arsenic associations in sediments from shallow aquifers of northwestern Hetao Basin, Inner Mongolia

Understanding the mechanism of arsenic mobilization from sediments to groundwater is important for water quality management in areas of endemic arsenic poisoning, such as the Hetao Basin in Inner Mongolia, northern China. Aquifer geochemistry was characterized at three field sites (SH, HF, TYS) in Hangjinhouqi County of northwestern Hetao Basin. The results of bulk geochemistry analysis of sediment samples indicated that total As concentrations have a range of 6. 8-58. 5 mg/kg, with a median of 14. 4 mg/kg. The highest As concentrations were found at 15-25 m depth. In the meanwhile, the range of As concentration in the sediments from background borehole is 3-21. 8 mg/kg, with a median value of 9 mg/kg. The As sediments concentrations with depth from the SH borehole were correlated with the contents of Fe, Sb, B, V, total C and total S. Generally, the abundance of elements varied with grain size, with higher concentrations in finer fractions of the sediments. Distinct lithology profile and different geochemical characteristics of aquifer sediments indicate the sediments are associated with different sources and diverse sedimentary environments. Up to one third of arsenic in the sediments could be extracted by ammonium oxalate, suggesting that Fe oxyhydroxides may be the major sink of As in the aquifer. Sequential extraction results indicate that arsenic occurs as strongly adsorbed on and/or co-precipitated with amorphous Fe oxyhydroxides in sediments accounting for 35 and 20%, respectively, of the total contents of arsenic. The release of As into groundwater may occur by desorption from the mineral surface driven by reductive dissolution of the Fe oxide minerals. Furthermore, small proportions of As associated with iron sulfides occur in the reductive sediments. ?? 2011 Springer-Verlag.

Environmental Earth Sciences↗
USGS2002

Prehistoric ruptures of the Gurvan Bulag fault, Gobi Altay, Mongolia

The 1957 Gobi Altay M8.3 earthquake in southern Mongolia was associated with the simultaneous rupture of several faults, including the Gurvan Bulag reverse fault, which is located about 25 km south of the main strike-slip Bogd fault. Our study of paleoseismic excavations across the Gurvan Bulag fault suggests that the penultimate surface rupture occurred after 6.0 ka, most likely between 2.6 and 4.4 ka, and a possible earlier rupture occurred after 7.3 ka. Our interpretation of the stratigraphic relations in one of the exposures suggests that at least five earthquakes have generated surface rupture of the Gurvan Bulag fault since the abandonment of an ancient alluvial fan surface. Luminescence dating of sediment associated with this surface indicates that it formed either 26.6 ± 2.1 ka or 16.1 ± 2.0 ka. These data imply that the recurrence intervals for surface faulting on the Gurvan Bulag and Bogd faults are similar, on the order of several thousands of years, but that the penultimate surface ruptures of the two faults did not occur during the same earthquake.

Gobi Altay↗
USGS2012

The Galapagos Spreading Centre at 86° W: A detailed geothermal field study

We report here measurements of the heat flow field of the Galapagos Spreading Center on crust of age less than 1.0 m.y. The 443 measurements in an area of about 570 km 2 reveal the general planform of the geothermal flux and permit the first truly areal estimate of the near-axis conductive heat flux. The intrusion process and associated hydrothermal circulation dominate the surface heat flow pattern, with circulation apparently continuing beyond the limits of our survey. The areal average of the conductive heat flux is 7.1 ± 0.8 HFU (295 ± 33 mW/m 2 ), about one-third the heat flux predicted by plate models. The remaining heat is apparently removed by venting of hydrothermal waters at the spreading axis and through basalt outcrops and hydrothermal mounds off axis. The pattern of surface heat flux is lineated parallel to the axis and the strongly lineated topography. Sharp lateral gradients in the heat flow, greater than 10 HFU/km near escarpments and commonly expressed as high heat flow at the tops of the scarps and lower heat flow in the valleys, may indicate a local concentration of the circulation by surface fault systems and/or variable sediment thickness.

Journal of Geophysical Research Solid Earth↗
USGS2001

Resolution analysis of finite fault source inversion using one- and three-dimensional Green's functions 1. Strong motions

We develop a methodology to perform finite fault source inversions from strong motion data using Green's functions (GFs) calculated for a three-dimensional (3-D) velocity structure. The 3-D GFs are calculated numerically by inserting body forces at each of the strong motion sites and then recording the resulting strains along the target fault surface. Using reciprocity, these GFs can be recombined to represent the ground motion at each site for any (heterogeneous) slip distribution on the fault. The reciprocal formulation significantly reduces the required number of 3-D finite difference computations to at most 3NS, where NS is the number of strong motion sites used in the inversion. Using controlled numerical resolution tests, we have examined the relative importance of accurate GFs for finite fault source inversions which rely on near-source ground motions. These experiments use both 1-D and 3-D GFs in inversions for hypothetical rupture models in order (1) to analyze the ability of the 3-D methodology to resolve trade-offs between complex source phenomena and 3-D path effects, (2) to address the sensitivity of the inversion results to uncertainties in the 3-D velocity structure, and (3) to test the adequacy of the 1-D GF method when propagation effects are known to be three-dimensional. We find that given "data" from a prescribed 3-D Earth structure, the use of well-calibrated 3-D GFs in the inversion provides very good resolution of the assumed slip distribution, thus adequately separating source and 3-D propagation effects. In contrast, using a set of inexact 3-D GFs or a set of hybrid 1-D GFs allows only partial recovery of the slip distribution. These findings suggest that in regions of complex geology the use of well-calibrated 3-D GFs has the potential for increased resolution of the rupture process relative to 1-D GFs. However, realizing this full potential requires that the 3-D velocity model and associated GFs should be carefully validated against the true 3-D Earth structure before performing the inverse problem with actual data.

Journal of Geophysical Research B: Solid Earth↗