Search USGSSearch

SEARCH · Search USGS

Results for “Journal of Sedimentary Research”

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 235 records · Page 13Linked to original sources

Relationships among seismic velocity, metamorphism, and seismic and aseismic fault slip in the Salton Sea Geothermal Field region

The Salton Sea Geothermal Field is one of the most geothermally and seismically active areas in California and presents an opportunity to study the effect of high-temperature metamorphism on the properties of seismogenic faults. The area includes numerous active tectonic faults that have recently been imaged with active source seismic reflection and refraction. We utilize the active source surveys, along with the abundant microseismicity data from a dense borehole seismic network, to image the 3-D variations in seismic velocity in the upper 5 km of the crust. There are strong velocity variations, up to ~30%, that correlate spatially with the distribution of shallow heat flow patterns. The combination of hydrothermal circulation and high-temperature contact metamorphism has significantly altered the shallow sandstone sedimentary layers within the geothermal field to denser, more feldspathic, rock with higher P wave velocity, as is seen in the numerous exploration wells within the field. This alteration appears to have a first-order effect on the frictional stability of shallow faults. In 2005, a large earthquake swarm and deformation event occurred. Analysis of interferometric synthetic aperture radar data and earthquake relocations indicates that the shallow aseismic fault creep that occurred in 2005 was localized on the Kalin fault system that lies just outside the region of high-temperature metamorphism. In contrast, the earthquake swarm, which includes all of the M > 4 earthquakes to have occurred within the Salton Sea Geothermal Field in the last 15 years, ruptured the Main Central Fault (MCF) system that is localized in the heart of the geothermal anomaly. The background microseismicity induced by the geothermal operations is also concentrated in the high-temperature regions in the vicinity of operational wells. However, while this microseismicity occurs over a few kilometer scale region, much of it is clustered in earthquake swarms that last from hours to a few days and are localized near the MCF system.

Salton Sea

A review of the geologic framework of the Long Island Sound Basin, with some observations relating to postglacial sedimentation

Most of the papers in this thematic section present regional perspectives that build on more than 100 years of geologic investigation in Long Island Sound. When viewed collectively, a common theme emerges in these works. The major geologic components of the Long Island Sound basin (bedrock, buried coastal-plain strata, recessional moraines, glacial-lake deposits, and the remains of a large marine delta) interact with the water body to affect the way the modern sedimentary system functions. Previous work, along with our present knowledge of the geologic framework of the Long Island Sound basin, is comprehensively reviewed with this theme in mind. Aspects of the crystalline bedrock, and the deltaic deposits associated with glacial Lake Connecticut, are examined with respect to their influence on sedimentation along the Connecticut coast and in the northern and western Sound. We also discuss the influence of the glacial drift that mantles the coastal-plain remnant along the north shore of Long Island and in the southern Sound. A total of approximately 22.7 billion m 3 of marine sediment has accumulated in the Long Island Sound basin. A significant portion (44%) of the fine-grained marine section in the central and western basins was redistributed there from the eastern Sound, as tidal scour removed slightly over 5 billion m 3 (5.3 X 10 12 kg) of fine material from glacial lake and early-marine deposits east of the Connecticut River. The remainder of the estimated 1.2 X 10 13 kg of fine-grained marine sediment that now resides in the central and western Sound can be accounted for by riverine input over the past 13.5 ka.

Connecticut, New York

Crustal resistivity structure from magnetotelluric soundings in the Colorado Plateau-Basin and Range provinces, central and western Arizona

Resistivity structure to about 25 km depth is defined from two-dimensional modeling of 29 magnetotelluric (MT) soundings (0.002–5 Hz) that traverse 280 km of the southwestern Colorado Plateau, transition zone, and Basin and Range provinces in Arizona. From the surface to 5 km depth, the MT model suggests structural relationships between low-resistivity sedimentary and volcanic rocks (50–300 ohm m) and high-resistivity granitic and gneissic basement (500–9000 ohm m). In the Basin and Range province, the MT model and a seismic reflection section show a generally consistent distribution of supracrustal rocks that have relatively low to moderate resistivity (MT) and relatively strong, locally coherent reflectivity. The supracrustal zone defined by these physical properties is inferred to be composed of upper plate rocks above a middle Tertiary detachment fault system which has been mapped in surrounding ranges. Some low-angle fault zones inferred from seismic reflections to extend into high resistivity basement below the supracrustal rocks are not resolved by the MT model. A low-resistivity zone with a conductance of 500 S or more is modeled in the crust at an average depth of about 15 km in the Basin and Range province and transition zone and may deepen below the southwestern part of the Colorado Plateau. In the Basin and Range province, the top of the low resistivity may correspond to a reflective layer with a 6-s two-way-travel time. This deep low-resistivity zone might be caused by a small fraction of connected hydrous solutions or silicic melts.

Journal of Geophysical Research

Transport slopes, sediment cover, and bedrock channel incision in the Henry Mountains, Utah

[1] Field data from channels in the Henry Mountains of Utah demonstrate that abundant coarse sediment can inhibit fluvial incision into bedrock by armoring channel beds (the cover effect). We compare several small channels that share tributary junctions and have incised into the same sedimentary bedrock unit (Navajo Sandstone) but contain differing amounts of coarse diorite clasts owing to the spatial distribution of localized sediment sources. Bedrock channels that contain abundant clasts (diorite-rich) have steeper longitudinal slopes than tributaries of these channels with smaller drainage areas and less sediment (diorite-poor). The diorite-poor tributaries have incised more deeply to lower average slopes and have more reach-scale slope variability, which may reflect bedrock properties, longitudinal sediment sorting, and incision at lower sediment supply. Diorite-rich channels have less bedrock exposed and smoother longitudinal profiles than diorite-poor channels. We find that (1) coarse sediment can mantle bedrock channel beds and reduce the efficiency of incision, validating the hypothesized cover effect in fluvial incision models; (2) the channel slope needed to transport the sediment load can be larger than that needed to erode bedrock, suggesting that the slope of incising bedrock channels can become adjusted to the sediment load; (3) when abundant sediment is available, transport capacity rather than thresholds of motion can be dominant in setting bedrock channel slope; and (4) cover effects can be important even when moderate amounts of bedrock are exposed in channel beds.

Utah

Quantifying the effect of petrogenic carbon on SOC turnover for two Rocky Mountain soils: When are petrogenic carbon corrections required?

Petrogenic organic carbon (OC petro ), derived from sedimentary rocks, is an often overlooked and poorly quantified source of soil organic carbon (SOC), which may influence measured or modeled SOC composition, age, and stability. In this study, we exploited differences in thermochemical stability between OC petro and biogenic SOC (OC bio ) using stepped elemental analysis to quantify the fractional contribution of OC petro to total SOC (f petro ), and we conducted a sensitivity analysis to estimate the effects of OC petro on modeled SOC transit times and system ages. Specifically, we compared the effects of accounting for OC petro inputs in SOC turnover modeling (using SoilR) for two montane meadow soils that are underlain by Cretaceous Mancos Shale. At these sites, we estimate that OC petro comprises 7%–9% of the total SOC stock (f petro = 0.07–0.09). However, accounting for OC petro as a mixture of inert and passive C or as completely inert C had negligible effects on SOC transit times and system ages, suggesting that there is a threshold of OC petro content under which there is minimal effect on calculated SOC turnover. Based on our sensitivity analysis, we estimate this threshold to be f petro = 0.125, further supporting that the accurate calculation of OC petro remains an important factor in estimating SOC turnover.

Colorado

On extracting sediment transport information from measurements of luminescence in river sediment

Accurately quantifying sediment transport rates in rivers remains an important goal for geomorphologists, hydraulic engineers, and environmental scientists. However, current techniques for measuring long-time scale (102–106 years) transport rates are laborious, and formulae to predict transport are notoriously inaccurate. Here we attempt to estimate sediment transport rates by using luminescence, a property of common sedimentary minerals that is used by the geoscience community for geochronology. This method is advantageous because of the ease of measurement on ubiquitous quartz and feldspar sand. We develop a model from first principles by using conservation of energy and sediment mass to explain the downstream pattern of luminescence in river channel sediment. We show that the model can accurately reproduce the luminescence observed in previously published field measurements from two rivers with very different sediment transport styles. The model demonstrates that the downstream pattern of river sand luminescence should show exponential-like decay in the headwaters which asymptotes to a constant value with further downstream distance. The parameters from the model can then be used to estimate the time-averaged virtual velocity, characteristic transport lengthscale, storage time scale, and floodplain exchange rate of fine sand-sized sediment in a fluvial system. The sediment transport values predicted from the luminescence method show a broader range than those reported in the literature, but the results are nonetheless encouraging and suggest that luminescence demonstrates potential as a sediment transport indicator. However, caution is warranted when applying the model as the complex nature of sediment transport can sometimes invalidate underlying simplifications.

Journal of Geophysical Research F: Earth Surface

Coral ages and island subsidence, Hilo drill hole

A 25.8-m-thick sedimentary section containing coral fragments occurs directly below a surface lava flow (the ∼1340 year old Panaewa lava flow) at the Hilo drill hole. Ten coral samples from this section dated by accelerator mass spectrometry (AMS) radiocarbon and five by thermal infrared multispectral scanner (TIMS) 230 Th/U methods show good agreement. The calcareous unit is 9790 years old at the bottom and 1690 years old at the top and was deposited in a shallow lagoon behind an actively growing reef. This sedimentary unit is underlain by a 34-m-thick lava flow which in turn overlies a thin volcaniclastic silt with coral fragments that yield a single 14 C date of 10,340 years. The age-depth relations of the dated samples can be compared with proposed eustatic sea level curves after allowance for island subsidence is taken. Island subsidence averages 2.2 mm/yr for the last 47 years based on measurements from a tide gage near the drill hole or 2.5–2.6 mm/yr for the last 500,000 years based on the ages and depths of a series of drowned coral reefs offshore from west Hawaii. The age-depth measurements of coral fragments are more consistent with eustatic sea levels as determined by coral dating at Barbados and Albrolhos Islands than those based on oxygen isotopic data from deep sea cores. The Panaewa lava flow entered a lagoon underlain by coral debris and covered the drill site with 30.9 m of lava of which 11 m was above sea level. This surface has now subsided to 4.2 m above sea level, but it demonstrates how a modern lava flow entering Hilo Bay would not only change the coastline but could extensively modify the offshore shelf.

Journal of Geophysical Research B: Solid Earth

Diagenesis of lower Cretaceous pelagic carbonates, North Atlantic: Paleoceanographic signals obscured

The stable isotope and minor element geochemistry of Neocomian (Lower Cretaceous) pelagic carbonates of the North Atlantic Basin (Deep Sea Drilling Project Sites 105, 367, 387, 391, and 603) was examined to develop a diagenetic model for pelagic limestones. In particular, we hoped to test the fidelity of whole-rock geochemical records as paleoceanographic indicators for pelagic deposits of pre-Aptian age, in which individual microfossils are not available for analysis. Data indicate that in addition to depth of burial, rhythmic variations in primary carbonate content have strongly controlled diagenetic patterns and associated geochemical signatures in these Neocomian sequences. Samples become increasingly depleted in Sr and 18O with increasing CaCO3 content. Within individual sedimentary sections, substantial decreases in Sr/Ca ratios and ??18O values are evident over a range of 4 to 98% CaCO3. However, even over a relatively narrow range of 50 to 98% CaCO3 a 2.5%c variation in ??18O values and a change of a factor of 1.7 in Sr/Ca ratios are observed. Carbon isotope compositions do not vary as extensively with CaCO3 content, but carbonate-rich intervals tend to be relatively depleted in 13C. Petrographic analysis reveals that these geochemical patterns are related to the transfer of CaCO3 from carbonate-poor intervals (calcareous shales and marlstones) to adjacent carbonate-rich intervals (limestones) during burial compaction and pressure solution. This process results in the addition of diagenetic cement to carbonate-rich intervals to produce a bulk composition that is relatively depleted in Sr and 18O and, at the same time, enables the retention of more-or-less primary carbonate that is relatively enriched in Sr and 18O in adjacent carbonate-poor intervals. Thus, although cyclic variations in CaCO3 content are primary in the Neocomian sequences examined, measured variations in Sr/Ca ratios and ??18O values are not and, as such, do not provide reliable proxies for past variations in climate, oceanographic conditions, or global ice volume.

Journal of Foraminiferal Research

The collision zone between the North d'Entrecasteaux Ridge and the New Hebrides island arc. 2. Structure from multichannel seismic data

The d'Entrecasteaux zone (DEZ) collides with the central New Hebrides island arc and consists of two subparallel ridges that strike east-west, stand 1–2 km above the surrounding oceanic plate, and subduct obliquely (15°) northward beneath the arc. Rocks dredged from the north ridge as well as reflections evident in multichannel seismic reflection data indicate that this ridge has a volcanic origin. Crystalline volcanic rocks are common along the lower flank of the ridge, but sedimentary, probably volcaniclastic, rock caps the ridge. Seismic reflection data collected over the lower arc slope reveal that mass wasting deposits locally make up most of the accretionary wedge. These deposits appear to form discrete bodies, suggesting that mass wasting occurred episodically. Large anticlines and thrust faults having large vertical separation are not readily evident where the colliding ridge intersects the arc slope; apparently, slope rocks have low strength so that mass wasting deposits formed instead of large-relief structures. Mass wasting is thought to occur as the accretionary wedge is uplifted in response to the northward oblique subduction of the north ridge. The toe of the north ridge flank marks an abrupt transition in the lithologies that make up the footwall of the interplate decollement. Footwall lithologies change from ocean basin to volcaniclastic ridge material, and this transition probably marks a discontinuity in friction along the decollement or in rock mechanical properties because north of the transition, thrust faults deform the accretionary wedge whereas south of the transition, steep reverse faults crosscut the wedge and pierce the north flank of the ridge. This piercement means that the decollement at least locally lies within the ridge and that ridge material exotic to the New Hebrides arc may be incorporated into the accretionary wedge.

Journal of Geophysical Research

Episodic fluid flow in the Nankai accretionary complex: Timescale, geochemistry, flow rates, and fluid budget

Down-hole geochemical anomalies encountered in active accretionary systems can be used to constrain the timing, rates, and localization of fluid flow. Here we combine a coupled flow and solute transport model with a kinetic model for smectite dehydration to better understand and quantify fluid flow in the Nankai accretionary complex offshore of Japan. Compaction of sediments and clay dehydration provide fluid sources which drive the model flow system. We explicitly include the consolidation rate of underthrust sediments in our calculations to evaluate the impact that variations in this unknown quantity have on pressure and chloride distribution. Sensitivity analysis of steady state pressure solutions constrains bulk and flow conduit permeabilities. Steady state simulations with 30% smectite in the incoming sedimentary sequence result in minimum chloride concentrations at site 808 of 550 m M , but measured chlorinity is as low as 447 m M . We simulate the transient effects of hydrofracture or a strain event by assuming an instantaneous permeability increase of 3–4 orders of magnitude along a flow conduit (in this case the décollement), using steady state results as initial conditions. Transient results with an increase in décollement permeability from 10 −16 m 2 to 10 −13 m 2 and 20% smectite reproduce the observed chloride profile at site 808 after 80–160 kyr. Modeled chloride concentrations are highly sensitive to the consolidation rate of underthrust sediments, such that rapid compaction of underthrust material leads to increased freshening. Pressures within the décollement during transient simulations rise rapidly to a significant fraction of lithostatic and remain high for at least 160 kyr, providing a mechanism for maintaining high permeability. Flow rates at the deformation front for transient simulations are in good agreement with direct measurements, but steady state flow rates are 2–3 orders of magnitude smaller than observed. Fluid budget calculations indicate that nearly 71% of the incoming water in the sediments leaves the accretionary wedge via diffuse flow out the seafloor, 0–5% escapes by focused flow along the décollement, and roughly 1% is subducted.

Journal of Geophysical Research B: Solid Earth

Carbon and sediment accumulation in the Everglades (USA) during the past 4000 years: rates, drivers, and sources of error

Tropical and sub-tropical wetlands are considered to be globally important sources for greenhouse gases but their capacity to store carbon is presumably limited by warm soil temperatures and high rates of decomposition. Unfortunately, these assumptions can be difficult to test across long timescales because the chronology, cumulative mass, and completeness of a sedimentary profile are often difficult to establish. We therefore made a detailed analysis of a core from the principal drainage outlet of the Everglades of South Florida, to assess these problems and determine the factors that could govern carbon accumulation in this large sub-tropical wetland. Accelerator mass spectroscopy dating provided direct evidence for both hard-water and open-system sources of dating errors, whereas cumulative mass varied depending upon the type of method used. Radiocarbon dates of gastropod shells, nevertheless, seemed to provide a reliable chronology for this core once the hard-water error was quantified and subtracted. Long-term accumulation rates were then calculated to be 12.1 g m -2 yr -1 for carbon, which is less than half the average rate reported for northern and tropical peatlands. Moreover, accumulation rates remained slow and relatively steady for both organic and inorganic strata, and the slow rate of sediment accretion ( 0.2 mm yr -1 ) tracked the correspondingly slow rise in sea level (0.35 mm yr -1 ) reported for South Florida over the past 4000 years. These results suggest that sea level and the local geologic setting may impose long-term constraints on rates of sediment and carbon accumulation in the Everglades and other wetlands.

Florida

High-resolution gamma-ray spectrometry in uranium exploration

Sedimentary-type uranium deposits accumulate at favorable sites along a migration path which may be kilometers in length. Their source is a large volume of rock from which the uranium has been leached. The geochemical mobilities and half lives of uranium and its daughter products vary widely so that they are transported from the source rocks, at different rates, along the migration path to their ultimate site. The radioactive disequilibrium resulting from this process has been well documented in the immediate vicinity of ore deposits, and disequilibrium is commonly recorded on gamma-ray logs up the hydraulic gradient from uranium ore. Little is known about the state of secular equilibrium in the leached host rocks, which often represent the only part of the migration path that is at or near the surface and is thus most accessible to the exploration geophysicist. High-resolution gamma-ray spectrometry provides a means of investigating the disequilibrium associated with uranium leaching and migration. Direct measurement of uranium can be made by this method, and the equivalent weight percents can be determined for six of the seven daughter-product decay groups that characterize the state of radioactive equilibrium. The technique has been used quantitatively in laboratory studies, where the results compare favorably with radiochemical analyses; field experiments suggest that semi-quantitative data may be obtained at the outcrop.

Journal of Research of the U.S. Geological Survey

CRUST 5.1: A global crustal model at 5° x 5°

We present a new global model for the Earth's crust based on seismic refraction data published in the period 1948-1995 and a detailed compilation of ice and sediment thickness. An extensive compilation of seismic refraction measurements has been used to determine the crustal structure on continents and their margins. Oceanic crust is modeled with both a standard model for normal oceanic crust, and variants for nonstandard regions, such as oceanic plateaus. Our model (CRUST 5.1) consists of 2592 5° × 5° tiles in which the crust and uppermost mantle are described by eight layers: (1) ice, (2) water, (3) soft sediments, (4) hard sediments, (5) crystalline upper, (6) middle, (7) lower crust, and (8) uppermost mantle. Topography and bathymetry are adopted from a Mandard database (ETOPO-5). Compressional wave velocity in each layer is based on field measurements, and shear wave velocity and density are estimated using recently published empirical V p -V s and V p -density relationships. The crustal model differs from previous models in that (1) the thickness and seismic/density structure of sedimentary basins is accounted for more completely, (2) the velocity structure of unmeasured regions is estimated using statistical averages that are based on a significantly larger database of crustal structure, (3) the compressional wave, shear wave, and density structure have been explicitly specified using newly available constraints from field and laboratory studies. Thus - this global crustal model is based on substantially more data than previous models and differs from them in many important respects. A new map of the thickness of the Earth's crust is presented, and we illustrate the application of this model by using it to provide the crustal correction for surface wave phase velocity maps. Love waves at 40 s are dominantly sensitive to crustal structure, and there is a very close correspondence between observed phase velocities at this period and those predicted by CRUST 5.1. We find that the application of crustal corrections to long-period (167 s) Rayleigh waves significantly increases the variance in the phase velocity maps and strengthens the upper mantle , velocity anomalies beneath stable continental regions. A simple calculation of crustal isostacy indicates significant lateral variations in upper mantle density. The model CRUST 5.1 provides a complete description of the physical properties of the Earth's crust at a scale of 5° × 5° and can be used for a wide range of seismological and nonseismological problems.

Journal of Geophysical Research B: Solid Earth

The storage time, age, and erosion hazard of laterally accreted sediment on the floodplain of a simulated meandering river

A sediment particle traversing the fluvial system may spend the majority of the total transit time at rest, stored in various sedimentary deposits. Floodplains are among the most important of these deposits, with the potential to store large amounts of sediment for long periods of time. The virtual velocity of a sediment grain depends strongly on the amount of time spent in storage, but little is known about sediment storage times. Measurements of floodplain vegetation age have suggested that storage times are exponentially distributed, a case that arises when all the sediment on a floodplain is equally vulnerable to erosion in a given interval. This assumption has been incorporated into sediment routing models, despite some evidence that younger sediment is more likely to be eroded from floodplains than older sediment. We investigate the relationship between sediment age and erosion, which we term the “erosion hazard,” with a model of a meandering river that constructs its floodplain by lateral accretion. We find that the erosion hazard decreases with sediment age, leading to a storage time distribution that is not exponential. We propose an alternate model that requires that channel motion is approximately diffusive and results in a heavy tailed distribution of storage time. The model applies to timescales over which the direction of channel motion is uncorrelated. We speculate that the lower end of this range of time is set by the meander cutoff timescale and the upper end is set by processes that limit the width of the meander belt.

Journal of Geophysical Research F: Earth Surface

Geometry and subsidence history of the Dead Sea basin: A case for fluid-induced mid-crustal shear zone?

Pull‐apart basins are narrow zones of crustal extension bounded by strike‐slip faults that can serve as analogs to the early stages of crustal rifting. We use seismic tomography, 2‐D ray tracing, gravity modeling, and subsidence analysis to study crustal extension of the Dead Sea basin (DSB), a large and long‐lived pull‐apart basin along the Dead Sea transform (DST). The basin gradually shallows southward for 50 km from the only significant transverse normal fault. Stratigraphic relationships there indicate basin elongation with time. The basin is deepest (8–8.5 km) and widest (∼15 km) under the Lisan about 40 km north of the transverse fault. Farther north, basin depth is ambiguous, but is 3 km deep immediately north of the lake. The underlying pre‐basin sedimentary layer thickens gradually from 2 to 3 km under the southern edge of the DSB to 3–4 km under the northern end of the lake and 5–6 km farther north. Crystalline basement is ∼11 km deep under the deepest part of the basin. The upper crust under the basin has lower P wave velocity than in the surrounding regions, which is interpreted to reflect elevated pore fluids there. Within data resolution, the lower crust below ∼18 km and the Moho are not affected by basin development. The subsidence rate was several hundreds of m/m.y. since the development of the DST ∼17 Ma, similar to other basins along the DST, but subsidence rate has accelerated by an order of magnitude during the Pleistocene, which allowed the accumulation of 4 km of sediment. We propose that the rapid subsidence and perhaps elongation of the DSB are due to the development of inter‐connected mid‐crustal ductile shear zones caused by alteration of feldspar to muscovite in the presence of pore fluids. This alteration resulted in a significant strength decrease and viscous creep. We propose a similar cause to the enigmatic rapid subsidence of the North Sea at the onset the North Atlantic mantle plume. Thus, we propose that aqueous fluid flux into a slowly extending continental crust can cause rapid basin subsidence that may be erroneously interpreted as an increased rate of tectonic activity.

Dead Sea shoreline

Predicted pH at the domestic and public supply drinking water depths, Central Valley, California

This scientific investigations map is a product of the U.S. Geological Survey (USGS) National Water-Quality Assessment (NAWQA) project modeling and mapping team. The prediction grids depicted in this map are of continuous pH and are intended to provide an understanding of groundwater-quality conditions at the domestic and public supply drinking water zones in the groundwater of the Central Valley of California. The chemical quality of groundwater and the fate of many contaminants is often influenced by pH in all aquifers. These grids are of interest to water-resource managers, water-quality researchers, and groundwater modelers concerned with the occurrence of natural and anthropogenic contaminants related to pH. In this work, the median well depth categorized as domestic supply was 30 meters below land surface, and the median well depth categorized as public supply is 100 meters below land surface. Prediction grids were created using prediction modeling methods, specifically boosted regression trees (BRT) with a Gaussian error distribution within a statistical learning framework within the computing framework of R ( http://www.r-project.org/ ). The statistical learning framework seeks to maximize the predictive performance of machine learning methods through model tuning by cross validation. The response variable was measured pH from 1,337 wells and was compiled from two sources: USGS National Water Information System (NWIS) database (all data are publicly available from the USGS: http://waterdata.usgs.gov/ca/nwis/nwis ) and the California State Water Resources Control Board Division of Drinking Water (SWRCB-DDW) database (water quality data are publicly available from the SWRCB: http://www.waterboards.ca.gov/gama/geotracker_gama.shtml ). Only wells with measured pH and well depth data were selected, and for wells with multiple records, only the most recent sample in the period 1993–2014 was used. A total of 1,003 wells (training dataset) were used to train the BRT model, and 334 wells (hold-out dataset) were used to validate the prediction model. The training r-squared was 0.70, and the root-mean-square error (RMSE) in standard pH units was 0.26. The hold-out r-squared was 0.43, and RMSE in standard pH units was 0.37. Predictor variables consisting of more than 60 variables from 7 sources were assembled to develop a model that incorporates regional-scale soil properties, soil chemistry, land use, aquifer textures, and aquifer hydrology. Previously developed Central Valley model outputs of textures (Central Valley Textural Model, CVTM; Faunt and others, 2010) and MODFLOW-simulated vertical water fluxes and predicted depth to water table (Central Valley Hydrologic Model, CVHM; Faunt, 2009) were used to represent aquifer textures and groundwater hydraulics, respectively. In this work, wells were attributed to predictor variable values in ArcGIS using a 500-meter buffer. Faunt, C.C., ed., 2009, Groundwater availability in the Central Valley aquifer, California: U.S. Geological Survey Professional Paper 1776, 225 p., accessed at https://pubs.usgs.gov/pp/1766/ . Faunt, C.C., Belitz, K., and Hanson, R.T., 2010, Development of a three-dimensional model of sedimentary texture in valley-fill deposits of Central Valley, California, USA: Hydrogeology Journal, v. 18, no. 3, p. 625–649, https://doi.org/10.1007/s10040-009-0539-7 .

California

Veneers, rinds, and fracture fills: Relatively late alteration of sedimentary rocks at Meridiani Planum, Mars

Veneers and thicker rinds that coat outcrop surfaces and partially cemented fracture fills formed perpendicular to bedding document relatively late stage alteration of ancient sedimentary rocks at Meridiani Planum, Mars. The chemistry of submillimeter thick, buff-colored veneers reflects multiple processes at work since the establishment of the current plains surface. Veneer composition is dominated by the mixing of silicate-rich dust and sulfate-rich outcrop surface, but it has also been influenced by mineral precipitation, including NaCl, and possibly by limited physical or chemical weathering of sulfate minerals. Competing processes of chemical alteration (perhaps mediated by thin films of water or water vapor beneath blanketing soils) and sandblasting of exposed outcrop surfaces determine the current distribution of veneers. Dark-toned rinds several millimeters thick reflect more extensive surface alteration but also indicate combined dust admixture, halite precipitation, and possible minor sulfate removal. Cemented fracture fills that are differentially resistant to erosion occur along the margins of linear fracture systems possibly related to impact. These appear to reflect limited groundwater activity along the margins of fractures, cementing mechanically introduced fill derived principally from outcrop rocks. The limited thickness and spatial distribution of these three features suggest that aqueous activity has been rare and transient or has operated at exceedingly low rates during the protracted interval since outcropping Meridiani strata were exposed on the plains surface.

Journal of Geophysical Research E: Planets

Petrogenesis of garnet two-mica granites in the Ruby Mountains, Nevada

Garnet-two-mica granites of Jurassic (160 ± 3 m.y.) and Cretaceous (83 ± 1.3 m.y.) ages intrude amphibolite facies metasedimentary rocks of Precambrian age and lower Paleozoic sedimentary rocks in the northern part of the Ruby Mountains, Nevada. High initial 87 Sr/ 86 Sr and high values of δ 18 O for minerals suggest a dominantly sedimentary source for the magmas. Total rare-earth contents are low and less fractionated (Ce N /Yb N = 10) than granitic rocks of similar composition in the Sierra Nevada and show a negative Eu anomaly (Eu/Eu* = 0.57). These data are consistent with melting in the stability field of plagioclase with little contribution to the melt from refractory minerals such as zircon, sphene, and apatite. Minor euhedral garnets are manganese-rich and magnesium-poor (al 60 sp 34 py 3.6 gr 2.8 ) compared to garnet (al 79 sp 6 py 13 gr 2 ) in the intruded metamorphic rocks and show an increase in manganese and a decrease in calcium in the outer few microns. Fractionation of magnesium-iron between garnet and biotite, and of sodium-potassium between K-feldspar and plagioclase in the granites suggest submagmatic equilibration temperatures, in the range of 365° to 505°C, whereas oxygen isotope equilibration temperatures for two granites are higher at 480° and 570°C. The compositions of the granites plot near minimum melting compositions in the water-saturated quartz-albite-orthoclase system. Comparison of muscovite plus quartz stability, the water-saturated granite solidus, and experimental garnet-melt equilibria suggests pressures of crystallization no lower than about 3.5 kbar. This pressure is compatible with that estimated from garnet-plagioclase-sillimanite-quartz equilibria for the intruded metasedimentary rocks.

Nevada