USGS Science⌕ Search

SEARCH · USGS Science

Results for “The Earth Observer”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,783 records · Page 99Linked to original sources

Silicic lunar volcanism: Testing the crustal melting model

Lunar silicic rocks were first identified by granitic fragments found in samples brought to Earth by the Apollo missions, followed by the discovery of silicic domes on the lunar surface through remote sensing. Although these silicic lithologies are thought to make up a small portion of the lunar crust, their presence indicates that lunar crustal evolution is more complex than originally thought. Models currently used to describe the formation of silicic lithologies on the Moon include in situ differentiation of a magma, magma differentiation with silicate liquid immiscibility, and partial melting of the crust. This study focuses on testing a crustal melting model through partial melting experiments on compositions representing lithologies spatially associated with the silicic domes. The experiments were guided by the results of modeling melting temperatures and residual melt compositions of possible protoliths for lunar silicic rocks using the thermodynamic modeling software, rhyolite-MELTS. Rhyolite-MELTS simulations predict liquidus temperatures of 950–1040 °C for lunar granites under anhydrous conditions, which guided the temperature range for the experiments. Monzogabbro, alkali gabbronorite, and KREEP basalt were identified as potential protoliths due to their ages, locations on the Moon (i.e., located near observed silicic domes), chemically evolved compositions, and the results from rhyolite-MELTS modeling. Partial melting experiments, using mixtures of reagent grade oxide powders representing bulk rock compositions of these rock types, were carried out at atmospheric pressure over the temperature range of 900–1100 °C. Because all lunar granite samples and remotely sensed domes have an elevated abundance of Th, some of the mixtures were doped with Th to observe its partitioning behavior. Run products show that at temperatures of 1050 and 1100 °C, melts of the three protoliths are not silicic in nature (i.e., they have <63 wt% SiO 2 ). By 1000 °C, melts of both monzogabbro and alkali gabbronorite approach the composition of granite, but are also characterized by immiscible Si-rich and Fe-rich liquids. Furthermore, Th strongly partitions into the Fe-rich, and not the Si-rich glass in all experimental runs. Our work provides important constraints on the mechanism of silicic melt formation on the Moon. The observed high-Th content of lunar granite is difficult to explain by silicate liquid immiscibility, because through this process, Th is not fractionated into the Si-rich phase. Results of our experiments and modeling suggests that silicic lunar rocks could be produced from monzogabbro and alkali gabbronorite protoliths by partial melting at T < 1000 °C. Additionally, we speculate that at higher pressures ( P ≥ 0.005 GPa), the observed immiscibility in the partial melting experiments would be suppressed.

American Mineralogist↗

The Phanerozoic impact cratering rate: Evidence from the farside of the Moon

The relatively recent (<1 b.y.) flux of asteroids and comets forming large craters on the Earth and Moon may be accurately recorded by craters with bright rays on the Moon's farside. Many previously unknown farside rayed craters are clearly distinguished in the low-phase-angle images returned by the Clementine spacecraft. Some large rayed craters on the lunar nearside are probably significantly older than 1 Ga; rays remain visible over the maria due to compositional contrasts long after soils have reached optical maturity. Most of the farside crust has a more homogeneous composition and only immature rays are visible. The size-frequency distribution of farside rayed craters is similar to that measured for Eratosthenian craters (up to 3.2 b.y.) at diameters larger than 15 km. The areal density of farside rayed craters matches that of a corrected tabulation of nearside Copernican craters. Hence the presence of bright rays due to immature soils around large craters provides a consistent time-stratigraphic basis for defining the base of the Copernican System. The density of large craters less than ∼3.2 b.y. old is ∼3.2 times higher than that of large farside rayed craters alone. This observation can be interpreted in two ways: (1) the average cratering rate has been constant over the past 3.2 b.y. and the base of the Copernican is ∼1 Ga, or (2) the cratering rate has increased in recent geologic time and the base of the Copernican is less than 1 Ga. We favor the latter interpretation because the rays of Copernicus (800–850 m.y. old) appear to be very close to optical maturity, suggesting that the average Copernican cratering rate was ∼35% higher than the average Eratosthenian rate. Other lines of evidence for an increase in the Phanerozoic (545 Ga) cratering rate are (1) the densities of small craters superimposed on Copernicus and Apollo landing sites, (2) the rates estimated from well-dated terrestrial craters (≤120 m.y.) and from present-day astronomical observations, and (3) the Proterozoic rate suggested by the crater record of Australia. The hypothesis most consistent with several key observations is that the cratering rate has increased by ∼2x during the past ∼300 m.y‥

Journal of Geophysical Research E: Planets↗

Crustal structure and tectonics of the northern part of the Southern Granulite Terrane, India

Deep seismic reflection studies investigating the exposed Archean lower continental crust of the Southern Granulite Terrane, India, yield important constraints on the nature and evolution of the deep crust, including the formation and exhumation of granulites. Seismic reflection images along the Kuppam–Bhavani profile reveal a band of reflections that dip southward from 10.5 to 15.0 s two-way-time (TWT), across a distance of 50 km. The bottom of these reflections beneath the Dharwar craton is interpreted as the Moho. Further south, another reflection band dipping northward is observed. These bands of reflectivity constitute a divergent reflection fabric that converges at the Moho boundary observed at the Mettur shear zone. Reflection fabrics that intersect at a steep angle are interpreted as a collisional signature due to the convergence of crustal blocks, which we infer resulted in crustal thickening and the formation of granulites. Anomalous gravity and magnetic signatures are also observed across the Mettur shear zone. The gravity model derived from the Bouguer gravity data corroborates seismic results. The tectonic regime and seismic reflection profiles are combined in a 3-D representation that illustrates our evidence for paleo-subduction at a collision zone. The structural dissimilarities and geophysical anomalies suggest that the Mettur shear zone is a suture between the Dharwar craton in the north and another crustal block in the south. This study contributes significantly to our understanding of the operation of Archean plate tectonics, here inferred to involve collision and subduction. Furthermore, it provides an important link between the Gondwanaland and global granulite evolution occurring throughout the late Archean.

Earth and Planetary Science Letters↗

Locating shallow seismic sources with waves scattered by surface topography: Validation of the method at the Nevada Test Site

Accurate and robust source location is fundamental in seismology. Previously, we proposed a new full waveform location method using waves scattered by topography near the source, and we demonstrated its potential in obtaining accurate source location using synthetic data. In the work presented here, we validate this method with observed seismic data from the 1993 Non-Proliferation Experiment, a chemical explosion with a precisely known location in a region with moderate topography at the Nevada Test Site. We select the first arrivals (Pn/Pg) and their immediate codas to determine the source location and compare our solution with the known location. We use a collocated-grid finite-difference method to calculate the strain Green’s tensor in a grid-search volume containing the source, and obtain three-component synthetic waveforms at 12 broadband seismic stations at regional distances using source-receiver reciprocity. We assess the solution using a least-squares misfit between the observed and synthetic waveforms. When scattered coda waves are incorporated, the best solution is within a few hundreds of meters of the exact source location, and the estimated uncertainty of the solution is reduced compared to the waveform based solution using only the P waves. The solution is robust to the choice of the frequency content and to the addition of random velocity heterogeneity. We conclude that the full waveform source location method is effective for real seismic data, though more validation tests like this one are needed to further understand its efficacy for source-station geometry, roughness of topography, signal quality and other factors.

Nevada↗

Temporal greenness trends in stable natural land cover and relationships with climatic variability across the conterminous United States

Assessment of temporal trends in vegetation greenness and related influences aids understanding of recent change in terrestrial ecosystems and feedbacks from weather, climate, and environment. We analyzed 1-km normalized difference vegetation index (NDVI) timeseries data (1989–2016) derived from the Advanced Very High Resolution Radiometer (AVHRR) and developed growing season time-integrated NDVI (GS-TIN) for estimating seasonal vegetation activity across stable natural land cover in the conterminous United States (CONUS). After removing areas from analysis that had experienced land cover conversion or modification, we conducted a monotonic trend analysis on the GS-TIN timeseries and found that significant positive temporal trends occurred over 35% of the area, while significant negative trends were observed over only 3.5%. Positive trends were prevalent in the forested lands of the eastern third of CONUS and far northwest, as well as in grasslands in the north central plains. We observed negative and nonsignificant trends mainly in the shrublands and grasslands across the northwest, southwest, and west central plains. To understand the relationship of climate variability with these temporal trends, we conducted partial and multiple correlation analyses on GS-TIN, growing season temperature, and water-year precipitation timeseries. The GS-TIN trends in northern forests were positively correlated with temperature. The GS-TIN trends in the central and western shrublands and grasslands were negatively correlated with temperature and positively correlated with precipitation. Our results revealed spatial patterns in vegetation greenness trends for different stable natural vegetation types across CONUS, enhancing understanding gained from prior studies based on coarser 8-km AVHRR data.

Earth Interactions↗

Recency of faulting and subsurface architecture of the San Diego Bay pull-apart basin, California, USA

In southern California, plate boundary motion between the North American and Pacific plates is distributed across several sub-parallel fault systems. The offshore faults of the California Continental Borderland (CCB) are thought to accommodate ~10-15% of the total plate boundary motion, but the exact distribution of slip and the mechanics of slip partitioning remain uncertain. The Newport-Inglewood-Rose Canyon fault is the easternmost fault within the CCB whose southern segment splays out into a complex network of faults beneath San Diego Bay. A pull-apart basin model between the Rose Canyon and the offshore Descanso fault has been used to explain prominent fault orientations and subsidence beneath San Diego Bay; however this model does not account for faults in the southern portion of the bay or faulting east of the bay. To investigate the characteristics of faulting and stratigraphic architecture beneath San Diego Bay, we combined a suite of reprocessed legacy airgun multi-channel seismic profiles and high-resolution Chirp data, with age and lithology controls from geotechnical boreholes and shallow sub-surface vibracores. This combined dataset is used to create gridded horizon surfaces, fault maps, and perform a kinematic fault analysis. The structure beneath San Diego Bay is dominated by down-to-the-east motion on normal faults that can be separated into two distinct groups. The strikes of these two fault groups can be explained with a double pull-apart basin model for San Diego Bay. In our conceptual model, the western portion of San Diego Bay is controlled by a right-step between the Rose Canyon and Descanso faults, which matches both observations and predictions from laboratory models. The eastern portion of San Diego Bay appears to be controlled by an inferred step-over between the Rose Canyon and San Miguel-Vallecitos faults and displays distinct fault strike orientations, which kinematic analysis indicates should have a significant component of strike-slip partitioning that is not detectable in the seismic data. The potential of a Rose Canyon-San Miguel-Vallecitos fault connection would effectively cut the stepover distance in half and have important implications for the seismic hazard of the San Diego-Tijuana metropolitan area (population ~3 million people).

California↗

Viking landing sites, remote-sensing observations, and physical properties of Martian surface materials

Important problems that confront future scientific exploration of Mars include the physical properties of Martian surface materials and the geologic processes that formed the materials. The design of landing spacecraft, roving vehicles, and sampling devices and the selection of landing sites, vehicle traverses, and sample sites will be, in part, guided by the physical properties of the materials. Four materials occur in the sample fields of the Viking landers: (1) drift, (2) crusty to cloddy, (3) blocky, and (4) rock. The first three are soillike. Drift materials is weak, loose, and porous. We estimate that it has a dielectric constant near 2.4 and a thermal inertia near 1 &#xD7; 10 &#x2212;3 to 3 &#xD7; 10 &#x2212;3 ( cal cm &#x2212;2 sec 1 2 K &#x2212;1 ) "> 1 × 10 −3 to 3 × 10 −3 (cal cm −2 sec 12 K −1 ) because of its low bulk density, fine grain size, and small cohesion. Crusty to cloddy material is expected to have a dielectric constant near 2.8 and a thermal inertia near 4 × 10 −3 to 7 × 10 −3 because of its moderate bulk density and cementation of grains. Blocky material should have a dielectric constant near 3.3 and a thermal inertia near 7 × 10 −3 to 9 × 10 −3 because of its moderate bulk density and cementation. Common basaltic rocks have dielectric constants near 8 and thermal inertias near 30 × 10 −3 to 60 × 10 −3 . Comparisons of estimated dielectric constants and thermal inertias of the materials at the landing sites with those obtained remotely by Earth-based radars and Viking Orbiter thermal sensors suggest that the materials at the landing sites are good analogs for materials elsewhere on Mars. Correlation of remotely estimated dielectric constant and thermal inertias indicates two modal values for paired values of dielectric constants and thermal inertias near (A) 2 and 2 × 10 −3 and (B) 3 and 6 × 10 −3 , respectively. These two modes are comparable to the dielectric constants and thermal inertias for drift and crusty to cloddy material, respectively. Dielectric constants and thermal inertias for blocky material are larger but consistent with values in the northern plains. Our interpretations are compatible with an aeolian origin for drift and similar materials elsewhere on Mars. The postulate that moderate dielectric constants and thermal inertias larger than 3 or 4 × 10 −3 are produced by cementation of soillike materials is partly consistent with the data. The average dielectric constant and thermal inertia and their correlation with one another suggest that most of the surface of Mars should present few difficulties to future surface exploration, but some surfaces may present difficulties for spacecraft that are not suitably designed.

Icarus↗

Very long period conduit oscillations induced by rockfalls at Kilauea Volcano, Hawaii

Eruptive activity at the summit of Kilauea Volcano, Hawaii, beginning in 2010 and continuing to the present time is characterized by transient outgassing bursts accompanied by very long period (VLP) seismic signals triggered by rockfalls from the vent walls impacting a lava lake in a pit within the Halemaumau pit crater. We use raw data recorded with an 11-station broadband network to model the source mechanism of signals accompanying two large rockfalls on 29 August 2012 and two smaller average rockfalls obtained by stacking over all events with similar waveforms to improve the signal-to-noise ratio. To determine the source centroid location and source mechanism, we minimize the residual error between data and synthetics calculated by the finite difference method for a point source embedded in a homogeneous medium that takes topography into account. We apply a new waveform inversion method that accounts for the contributions from both translation and tilt in horizontal seismograms through the use of Green's functions representing the seismometer response to translation and tilt ground motions. This method enables a robust description of the source mechanism over the period range 1–1000 s. The VLP signals associated with the rockfalls originate in a source region ∼1 km below the eastern perimeter of the Halemaumau pit crater. The observed waveforms are well explained by a simple volumetric source with geometry composed of two intersecting cracks including an east striking crack (dike) dipping 80° to the north, intersecting a north striking crack (another dike) dipping 65° to the east. Each rockfall is marked by a similar step-like inflation trailed by decaying oscillations of the volumetric source, attributed to the efficient coupling at the source centroid location of the pressure and momentum changes induced by the rock mass impacting the top of the lava column. Assuming a simple lumped parameter representation of the shallow magmatic system, the observed pressure and volume variations can be modeled with the following attributes: rockfall volume (200–4500 m 3 ), length of magma column (120–210 m), diameter of pipe connecting the Halemaumau pit crater to the subjacent dike system (6 m), average thickness of the two underlying dikes (3–6 m), and effective magma viscosity (30–210 Pa s). Most rockfalls occur during episodes of sustained deflation of the Kilauea summit. The mass loss rate in the shallow magmatic system is estimated to be 1400–15,000 kg s −1 based on measurements of the temporal variation of VLP period in the two large rockfalls that occurred on 29 August 2012.

Hawaii↗

Silicate melt inclusions and glasses in lunar soil fragments from the Luna 16 core sample

More than 2000 fragments were studied microscopically, and electron microprobe analyses were made of 39 selected areas, from a few square mm of polished surface, through 75- to 425-μm fragments of lunar soil from two samples of the Luna 16 core. The silicate melt inclusions and glasses differ in important details from those observed earlier in the Apollo samples. Melt inclusions in olivine contain epitaxially oriented daughter crystals, but also show a similar epitaxy around the outside of the crystals not observed in previous lunar samples. Melt inclusions in ilmenite suggest trapping at successive stages in a differentiation sequence. There is abundant evidence for late-stage silicate liquid immiscibility, with melt compositions similar but not identical to those from Apollo 11 and 12. A comparison of the alkali ratio of any given bulk rock analysis with that of its late-stage, high-silica melt shows gross differences for different rocks. This is pertinent to understanding late-stage differentiation processes. Glass fragments and spherules exhibit a wide range of crystallization textures, reflecting their wide range of compositions and cooling histories. No significant differences were found between the two portions of core examined (Zones A and D).

Earth and Planetary Science Letters↗

Rapid fluid disruption: A source for self-potential anomalies on volcanoes

Self-potential (SP) anomalies observed above suspected magma reservoirs, dikes, etc., on various volcanoes (Kilauea, Hawaii; Mount Unzen, Japan; Piton de la Fournaise, Reunion Island, Miyake Jima, Japan) result from transient surface electric fields of tens of millivolts per kilometer and generally have a positive polarity. These SP anomalies are usually attributed to electrokinetic effects where properties controlling this process are poorly constrained. We propose an alternate explanation that contributions to electric fields of correct polarity should be expected from charge generation by fluid vaporization/disruption. As liquids are vaporized or removed as droplets by gas transport away from hot dike intrusions, both charge generation and local increase in electrical resistivity by removal of fluids should occur. We report laboratory observations of electric fields in hot rock samples generated by pulses of fluid (water) through the rock at atmospheric pressure. These indicate the relative amplitudes of rapid fluid disruption (RFD) potentials and electrokinetic potentials to be dramatically different and the signals are opposite in sign. Above vaporization temperatures, RFD effects of positive sign in the direction of gas flow dominate, whereas below these temperatures, effects of negative sign dominate. This suggests that the primary contribution to observed self-potential anomalies arises from gas-related charge transport processes at temperatures high enough to produce vigorous boiling and vapor transport. At lower temperatures, the primary contribution is from electrokinetic effects modulated perhaps by changing electrical resistivity and RFD effects from high-pressure but low-temperature CO2 and SO2 gas flow ripping water molecules from saturated crustal rocks. If charge generation is continuous, as could well occur above a newly emplaced dike, positive static potentials will be set up that could be sustained for many years, and the simplest method for identifying these hot, active regions would be to identify the SP anomalies they generate.

Hawai'i↗

Frequency and initiation of debris flows in Grand Canyon, Arizona

Debris flows from 740 tributaries transport sediment into the Colorado River in Grand Canyon, Arizona, creating rapids that control its longitudinal profile. Debris flows mostly occur when runoff triggers failures in colluvium by a process termed “the fire hose effect.” Debris flows originate from a limited number of geologic strata, almost exclusively shales or other clay‐rich, fine‐grained formations. Observations from 1984 through 2003 provide a 20 year record of all debris flows that reached the Colorado River in Grand Canyon, and repeat photography provides a 100 year record of debris flows from 147 tributaries. Observed frequencies are 5.1 events/year from 1984 to 2003, and historic frequencies are 5.0 events/year from 1890 to 1983. Logistic regression is used to model historic frequencies based on drainage basin parameters observed to control debris flow initiation and transport. From 5 to 7 of the 16 parameters evaluated are statistically significant, including drainage area, basin relief, and the height of and gradient below debris flow source areas, variables which reflect transport distance and potential energy. The aspect of the river channel, which at least partially reflects storm movement within the canyon, is also significant. Model results are used to calculate the probability of debris flow occurrence at the river over a century for all 740 tributaries. Owing to the variability of underlying geomorphic controls, the distribution of this probability is not uniform among tributaries of the Colorado River in Grand Canyon.

Arizona↗

Numerical simulation of sound-side barrier-island inundation and breaching during Hurricane Dorian (2019)

Hurricane-induced morphological changes and associated community hazards along sandy, barrier-island coastlines have been studied primarily from the perspective of ocean-side attack by storm-driven ocean surge and large waves. Thus, our understanding of long-term barrier island morphological change focuses on beach erosion, overwash, and inlet formation. In contrast, outwash events with inundation from the sound side, such as one that occurred in Cape Lookout National Seashore, North Carolina, USA during Hurricane Dorian (September 2019), are understudied. Studying such events can improve understanding of barrier island response and stability for a broader range of conditions. Here, we model the hydrodynamics and morphological evolution of a barrier island using a coupled wave-current-sediment transport modeling system. Wind-driven surge in Pamlico Sound led to overtopping from the sound side, which eroded outwash channels and transported sediment seaward into the nearshore. Simulations reproduce the channel features observed with aerial imagery and provide information not available from the remote-sensing observations, including channel depths (>2 m) and the fate of the eroded sand. We found that >99% of the eroded sand was deposited in the nearshore, within 1,000 m of the shoreline in depths <10 m, suggesting that the deposited sediment remains available for littoral transport and beach recovery. Simulations with combinations of coarse or fine sediment and vegetated or unvegetated landcover indicate that channel position did not vary with grain size or vegetation, while volume of erosion and channel morphology were more responsive to variations in grain size and less responsive to presence of vegetation.

North Carolina↗

Tide gage measurements of uplift along the south coast of Alaska

Annual mean sea levels along the south coast of Alaska are used to measure uplift along the Alaska‐Aleutian subduction zone. Oceanographic effects are removed from the observed annual mean sea levels by subtracting a correction that is proportional to the sea level fluctuations observed in southeast Alaska. That correction is effective in reducing fluctuations in the observed, annual mean sea level as far west as the tip of Alaska peninsula. Additional corrections to remove the eustatic rise in sea level and the apparent fall in sea level due to postglacial isostatic rebound of the land are introduced. This corrected sea level record should provide a measure of tectonic subsidence. In the area affected by the 1964 Alaska earthquake, postseismic uplift occurs where coseismic subsidence was observed, and postseismic subsidence occurs where coseismic uplift was observed. The immediate postseismic response is damped out within the first decade, and the subsequent uplift rates appear to be steady over the 1974–1989 interval. However, some of those rates seem to be too high to be sustained over the ∼1000 year earthquake recurrence interval appropriate to this area if the interseismic deformation is only to recover the coseismic displacement. Thus a long‐term ( ∼100 years) relaxation in uplift rates is postulated. The immediate (time constant ∼5 years) postseismic relaxation is attributed to postseismic slip on the plate interface directly downdip from the coseismic rupture. The long‐term (time constant ∼100 years) relaxation is attributed to flow in the asthenosphere.

Alaska↗

The group separation of the rare-earth elements and yttrium from geologic materials by cation-exchange chromatography

Demand is increasing for the determination of the rare-earth elements (REE) and yttrium in geologic materials. Due to their low natural abundance in many materials and the interferences that occur in many methods of determination, a separation procedure utilizing gradient strong-acid cation-exchange chromatography is often used to preconcentrate and isolate these elements from the host-rock matrix. Two separate gradient strong-acid cation-exchange procedures were characterized and the major elements as well as those elements thought to provide the greatest interference for the determination of the REE in geologic materials were tested for separation from the REE. Simultaneous inductively coupled argon plasma-atomic emission spectroscopy (ICAP-AES) measurements were used to construct the chromatograms for the elution studies, allowing the elution patterns of all the elements of interest to be determined in a single fraction of eluent. As a rock matrix, U.S. Geological Survey standard reference BCR-1 basalt was digested using both an acid decomposition procedure and a lithium metaborate fusion. Hydrochloric and nitric acids were tested as eluents and chromatograms were plotted using the ICAP-AES data; and we observed substantial differences in the elution patterns of the REE and as well as in the solution patterns of Ba, Ca, Fe and Sr. The nitric acid elution required substantially less eluent to elute the REE and Y as a group when compared to the hydrochloric acid elution, and provided a clearer separation of the REE from interfering and matrix elements. ?? 1984.

Chemical Geology↗

Tracing chlorine sources of thermal and mineral springs along and across the Cascade Range using halogen and chlorine isotope compositions

In order to provide constraints on the sources of chlorine in spring waters associated with arc volcanism, the major/minor element concentrations and stable isotope compositions of chlorine, oxygen, and hydrogen were measured in 28 thermal and mineral springs along the Cascade Range in northwestern USA. Chloride concentrations in the springs range from 64 to 19,000 mg/L and View the MathML source values range from +0.2&permil; to +1.9&permil; (average=+1.0&plusmn;0.4&permil;), with no systematic variation along or across the arc, nor correlations with their presumed underlying basement lithologies. Additionally, nine geochemically well-characterized lavas from across the Mt. St. Helens/Mt. Adams region of the Cascade Range (Leeman et al., 2004 and Leeman et al., 2005) were analyzed for their halogen concentrations and Cl isotope compositions. In the arc lavas, Cl and Br concentrations from the volcanic front are higher than in lavas from the forearc and backarc. F and I concentrations progressively decrease from forearc to backarc, similar to the trend documented for B in most arcs. View the MathML source values of the lavas range from &minus;0.1 to +0.8&permil; (average = +0.4&plusmn;0.3&permil;). Our results suggest that the predominantly positive View the MathML source values observed in the springs are consistent with water interaction with underlying 37Cl-enriched basalt and/or altered oceanic crust, thereby making thermal spring waters a reasonable proxy for the Cl isotope compositions of associated volcanic rocks in the Cascades. However, waters with View the MathML source values >+1.0&permil; also suggest additional contributions of chlorine degassed from cooling magmas due to subsurface vapor&ndash;liquid HCl fractionation in which Cl is lost to the aqueous fluid phase and 37Cl is concentrated in the ascending magmatic HCl vapor. Future work is necessary to better constrain Cl isotope behavior during volcanic degassing and fluid&ndash;rock interaction in order to improve volatile flux estimates through subduction zones.

Cascade Range↗

Trend analysis of ground-water levels and spring discharge in the Yucca Mountain Region, Nevada and California, 1960-2000

Ground-water level and discharge data from 1960 to 2000 were analyzed for the Yucca Mountain region of southern Nevada and eastern California. Included were water-level data from 37 wells and a fissure (Devils Hole) and discharge data from five springs and from a flowing well. Data were evaluated for variability and for upward, downward, or cyclic trends with an emphasis on the period 1992-2000. Potential factors causing trends in water levels and discharge include ground-water withdrawal, infiltration of precipitation, earthquakes, evapotranspiration, barometric pressure, and earth tides. Statistically significant trends in ground-water levels or spring discharge from 1992 to 2000 were upward at 12 water-level sites and downward at 14 water-level sites and 1 spring-discharge site. In general, the magnitude of the change in water level from 1992 to 2000 was small (less than 2 feet), except where influenced by pumping or local effects such as possible equilibration from well construction or diversion of nearby surface water. Seasonal trends are superimposed on some of the long-term (1992-2000) trends in water levels and discharge. Factors causing seasonal trends include barometric pressure, evapotranspiration, and pumping. The magnitude of seasonal change in water level can vary from as little as 0.05 foot in regional aquifers to greater than 5 feet in monitoring wells near large supply wells in the Amargosa Farms area. Three major episodes of earthquake activity affected water levels in wells in the Yucca Mountain region between 1992 and 2000: the Landers/Little Skull Mountain, Northridge, and Hector Mine earthquakes. The Landers/Little Skull Mountain earthquakes, in June 1992, had the largest observed effect on water levels and on discharge during the study period. Monthly measurements of wells in the study network show that earthquakes affected water levels from a few tenths of a foot to 3.5 feet. In the Ash Meadows area, water levels remained relatively stable from 1992 to 2000, with some water levels showing small rising trends and some declining slightly. Possible reasons for water-level fluctuations at sites AD-6 (Tracer Well 3), AM-5 (Devils Hole Well), and AM-4 (Devils Hole) from 1960 to 2000 include climate change, local and regional ground-water withdrawals, and tectonic activity. In Jackass Flats, water levels from 1992 to 2000 in six wells adjacent to Fortymile Wash displayed either small upward trends or no upward or downward trend. Comparison of trends in water levels from 1983 to 2000 for these six wells shows good correlations between all wells and suggests a common mechanism controlling water levels in the area. Of the likely controls on the system--precipitation or pumping in Jackass Flats--precipitation appears to be the predominant factor controlling water levels near Fortymile Wash. Water levels in the heavily pumped Amargosa Farms area declined from about 10 to 30 feet from 1964 to 2000. Water-level declines accelerated beginning in the early 1990's as pumping rates increased substantially. Pumping in the Amargosa Farms area may affect water levels in some wells as far away as 5-14 miles. The water level at site DV-3 (Travertine Point 1 Well) and discharge at site DV-2 (Navel Spring), both in the Death Valley hydrographic area, had downward trends from 1992 to 2000. The cause of these downward trends may be linked to earthquakes, pumping in the Amargosa Farms area, or both.

Water-Resources Investigations Report↗

Thermal and petrologic constraints on lower crustal melt accumulation under the Salton Sea Geothermal Field

In the Salton Sea region of southern California (USA), concurrent magmatism, extension, subsidence, and sedimentation over the past 0.5 to 1.0 Ma have led to the creation of the Salton Sea Geothermal Field (SSGF)—the second largest and hottest geothermal system in the continental United States—and the small-volume rhyolite eruptions that created the Salton Buttes. In this study, we determine the flux of mantle-derived basaltic magma that would be required to produce the elevated average heat flow and sustain the magmatic roots of rhyolite volcanism observed at the surface of the Salton Sea region. We use a 2D thermal model to show that a lower-crustal, partially molten mush containing &lt; 20 &#x2013; 40 % "> <20–40% interstitial melt develops over a &#x223C; 10 5 "> ∼105 -yr timescale for basalt fluxes of 0.008 to 0.010 m 3 / m 2 / yr "> 0.010m3/m2/yr (∼0.0008 to ∼0.001 km / 3 yr "> /3yr injection rate) given extension rates at or below the current value of ∼0.01 m/yr ( Brothers et al., 2009 ). These regions of partial melt are a natural consequence of a thermal regime that scales with average surface heat flow in the Salton Trough, and are consistent with seismic observations. Our results indicate limited melting and assimilation of pre-existing rocks in the lower crust. Instead, we find that basalt fractionation in the lower crust produces derivative melts of andesitic to dacitic composition. Such melts are then expected to ascend and accumulate in the upper crust, where they further evolve to give rise to small-volume rhyolite eruptions (Salton Buttes) and fuel local spikes in surface heat flux as currently seen in the SSGF. Such upper crustal magma evolution, with limited assimilation of hydrothermally altered material, is required to explain the slight decrease in &#x3B4; 18 O "> δ18O values of zircons (and melts) that have been measured in these rhyolites.

California↗

Hydrate morphology: Physical properties of sands with patchy hydrate saturation

The physical properties of gas hydrate-bearing sediments depend on the volume fraction and spatial distribution of the hydrate phase. The host sediment grain size and the state of effective stress determine the hydrate morphology in sediments; this information can be used to significantly constrain estimates of the physical properties of hydrate-bearing sediments, including the coarse-grained sands subjected to high effective stress that are of interest as potential energy resources. Reported data and physical analyses suggest hydrate-bearing sands contain a heterogeneous, patchy hydrate distribution, whereby zones with 100% pore-space hydrate saturation are embedded in hydrate-free sand. Accounting for patchy rather than homogeneous hydrate distribution yields more tightly constrained estimates of physical properties in hydrate-bearing sands and captures observed physical-property dependencies on hydrate saturation. For example, numerical modeling results of sands with patchy saturation agree with experimental observation, showing a transition in stiffness starting near the series bound at low hydrate saturations but moving toward the parallel bound at high hydrate saturations. The hydrate-patch size itself impacts the physical properties of hydrate-bearing sediments; for example, at constant hydrate saturation, we find that conductivity (electrical, hydraulic and thermal) increases as the number of hydrate-saturated patches increases. This increase reflects the larger number of conductive flow paths that exist in specimens with many small hydrate-saturated patches in comparison to specimens in which a few large hydrate saturated patches can block flow over a significant cross-section of the specimen.

Journal of Geophysical Research B: Solid Earth↗