Temporal and chemical connections between plutons and ignimbrites from the Mount Princeton magmatic center

The Mount Princeton magmatic center, located in central Colorado, consists of the epizonal Mount Princeton batholith, the nested Mount Aetna caldera, and volumetrically minor leucogranites. New CA-TIMS U/Pb zircon ages indicate the majority of the Mount Princeton batholith was emplaced during a period of regional ignimbrite quiescence. The structurally highest unit of quartz monzonite yields a ²⁰⁶Pb/²³⁸U age of 35.80 ± 0.10 Ma, and the youngest dated unit of the quartz monzonite is a porphyritic unit that yields a ²⁰⁶Pb/²³⁸U age of 35.37 ± 0.10 Ma. Using the exposed, dated volume of the quartz monzonite and new geochronology yields an estimated pluton filling rate of ~0.002 km³/a. This rate is comparable to the accumulation rates published for other plutons, and at least an order of magnitude slower than fluxes necessary to support accumulation of large eruptible magma volumes. Geochronology for the two large ignimbrites spatially associated with the batholith indicates a temporal disconnect between the vast majority of pluton building and explosive eruption of magma. The Wall Mountain Tuff erupted from a source in the same geographic area as the Mount Princeton batholith at 37.3 Ma (Ar/Ar sanidine), but no structural evidence of a caldera or temporally associated plutonic rocks is known. The Badger Creek Tuff erupted at 34.3 Ma (Ar/Ar sanidine) during the formation of the Mount Aetna caldera in the southern portion of the batholith. Our ²⁰⁶Pb/²³⁸U age for the Badger Creek Tuff is 34.47 ± 0.05. The only analyzed plutonic rocks of similar age to the Badger Creek Tuff are an extra-caldera dike with a ²⁰⁶Pb/²³⁸U age of 34.57 ± 0.08 Ma, a ring dike with a ²⁰⁶Pb/²³⁸U age of 34.48 ± 0.09 Ma, and a portion of the Mount Aetna pluton with a ²⁰⁶Pb/²³⁸U age of 34.60 ± 0.13 Ma. The small volume intrusions related to the eruption of the Badger Creek Tuff are chemically similar to the ignimbrite and show no signature of crystal–liquid separation in the shallow crust.     

Rio Tinto estuary (Spain): 5000 years of pollution

 Mining of massive sulfide deposits in southwestern Spain extending back to the Copper and Bronze Ages has resulted in the pollution of the Rio Tinto fluvial-estuarine complex, the site of Columbus' departure for the New World in 1492. Additional sources of potential pollution include the large industrial complex at Huelva near the lower portion of the estuary. Extensive analysis of surface sediment samples and cores has established that there are no geographic trends in the distribution of the pollutants, which include Cu, Fe, Pb, Zn, Ti, Ba, Cr, V and Co. These data have, however, demonstrated that tidal flux within the estuary carries phosphorus and perhaps other elements from the industrial complex at Huelva to the tidal limit of the system, several kilometers upstream from the discharge site. Radiometric analysis of short cores shows that sedimentation rates over at least the past couple of centuries have been about 0.3 cm/year. These data and that from a single deep core demonstrate that the estuary was polluted from mining activity long before the large-scale operations began in the late nineteenth century. Received: 3 May 1999 · Accepted: 18 October 1999     

Implications of spatial reservoir uncertainty for CO2 sequestration in the east Snake River Plain, Idaho (USA)

Basalt-hosted hydrogeologic systems have been proposed for geologic CO₂ sequestration based on laboratory research suggesting rapid mineralization rates. However, despite this theoretical appeal, little is known about the impacts of basalt fracture heterogeneity on CO₂ migration at commercial scales. Evaluating the suitability of basalt reservoirs is complicated by incomplete knowledge of in-situ fracture distributions at depths required for CO₂ sequestration. In this work, a numerical experiment is used to investigate the effects of spatial reservoir uncertainty for geologic CO₂ sequestration in the east Snake River Plain, Idaho (USA). Two criteria are investigated: (1) formation injectivity and (2) confinement potential. Several theoretical tools are invoked to develop a field-based approach for geostatistical reservoir characterization and their implementation is illustrated. Geologic CO₂ sequestration is simulated for 10?years of constant-rate injection at ~680,000 tons per year and modeled by Monte Carlo simulation such that model variability is a function of spatial reservoir heterogeneity. Results suggest that the spatial distribution of heterogeneous permeability structures is a controlling influence on formation injectivity. Analysis of confinement potential is less conclusive; however, in the absence of confining sedimentary interbeds within the basalt pile, rapid mineralization may be necessary to reduce the risk of escape.     

Source and supply of sediment to a shoreline salient in a fringing reef environment

Reef-associated landforms are coupled to the health of the reef ecosystem which produces the sediment that forms and maintains these landforms. However, this connection can make reef-fronted coastlines sensitive to the impacts of climate change, given that any decline in ecosystem health (e.g. decreasing sediment supply) or changes to physical processes (e.g. sea level rise, increasing wave energy) could drive the sediment budgets of these systems into a net erosive state. Therefore, knowledge of both the sediment sources and transport mechanisms is required to predict the sensitivity of reef-associated landforms to future climate change. Here, we examine the benthic habitat composition, sediment characteristics (composition, texture, and age), and transport mechanisms and pathways to understand the interconnections between coastal morphology and the reef system at Tantabiddi, Ningaloo Reef, Western Australia. Benthic surveys and sediment composition analysis revealed that although live coral accounts for less than 5% of the benthic cover, coral is the dominant sediment constituent (34% on average). Sediment ages (²³⁸U/²³⁰Th) were mostly found to be thousands of years old, suggesting that the primary sediment source is relic reef material (e.g. Holocene reef framework). Sediment transport across the lagoon was quantified through measurements of ripple migration rates, which were found to be shoreward migrating and responsible for feeding the large shoreline salient in the lee of the reef. The derived sediment fluxes were comparable with previously measured rates of sediment production by bioerosion. These results suggest that sediment budgets of systems dependent on old (>10³ years) source materials may be more resilient to climate change as present-day reef health and community composition (i.e. sources of ‘new’ carbonate production) have limited influence on sediment supply. Therefore, the vulnerability of reef-associated landforms in these systems will be dictated by future changes to mechanisms of sediment generation (e.g. bioerosion) and/or physical processes. © 2018 John Wiley & Sons, Ltd.     

Late Pleistocene to Holocene sedimentation and hydrocarbon seeps on the continental shelf of a steep,tectonically active margin,southern California,USA

Small, steep, uplifting coastal watersheds are prolific sediment producers that contribute significantly to the global marine sediment budget. This study illustrates how sedimentation evolves in one such system where the continental shelf is largely sediment-starved, with most terrestrial sediment bypassing the shelf in favor of deposition in deeper basins. The Santa Barbara–Ventura coast of southern California, USA, is considered a classic area for the study of active tectonics and of Tertiary and Quaternary climatic evolution, interpretations of which depend upon an understanding of sedimentation patterns. High-resolution seismic-reflection data over >570 km² of this shelf show that sediment production is concentrated in a few drainage basins, with the Ventura and Santa Clara River deltas containing most of the upper Pleistocene to Holocene sediment on the shelf. Away from those deltas, the major factor controlling shelf sedimentation is the interaction of wave energy with coastline geometry. Depocenters containing sediment 5–20 m thick exist opposite broad coastal embayments, whereas relict material (bedrock below a regional unconformity) is exposed at the sea floor in areas of the shelf opposite coastal headlands. Locally, natural hydrocarbon seeps interact with sediment deposition either to produce elevated tar-and-sediment mounds or as gas plumes that hinder sediment settling. As much as 80% of fluvial sediment delivered by the Ventura and Santa Clara Rivers is transported off the shelf (some into the Santa Barbara Basin and some into the Santa Monica Basin via Hueneme Canyon), leaving a shelf with relatively little recent sediment accumulation. Understanding factors that control large-scale sediment dispersal along a rapidly uplifting coast that produces substantial quantities of sediment has implications for interpreting the ancient stratigraphic record of active and transform continental margins, and for inferring the distribution of hydrocarbon resources in relict shelf deposits.     

Relationships between molecular bacterial diversity and chemistry of groundwater in the Wairarapa Valley,New Zealand

Groundwater plays an important role in New Zealand water supplies and hence monitoring activities are conducted regularly. Most monitoring programmes aim to evaluate groundwater chemistry and almost completely overlook the microbial component in this ecosystem. In our present study, the bacterial community structure of groundwater in the Wairarapa Valley was examined using the terminal restriction fragment length polymorphism (T-RFLP), and relationships between bacterial community structure and groundwater chemistry, aquifer confinement and groundwater usage were explored. In addition, the results from this study were compared with a previous T-RFLP survey of the same area in an attempt to detect changes in bacterial community structure over time. The data obtained suggested that bacterial community structure was related to groundwater chemistry, especially to redox conditions. Species composition showed minimal variation over time if groundwater chemistry remained unchanged. These findings reflect the potential of using bacterial communities as biological indicators to evaluate the health of groundwater ecosystems. We suggest that it is important to include this type of broad bacterial diversity assessment criteria into regular groundwater monitoring activities.     

Reproduction of secondary data in projection pursuit transformation

A challenge when working with multivariate data in a geostatistical context is that the data are rarely Gaussian. Multivariate distributions may include nonlinear features, clustering, long tails, functional boundaries, spikes, and heteroskedasticity. Multivariate transformations account for such features so that they are reproduced in geostatistical models. Projection pursuit as developed for high dimensional data exploration can also be used to transform a multivariate distribution into a multivariate Gaussian distribution with an identity covariance matrix. Its application within a geostatistical modeling context is called the projection pursuit multivariate transform (PPMT). An approach to incorporate exhaustive secondary variables in the PPMT is introduced. With this approach the PPMT can incorporate any number of secondary variables with any number of primary variables. A necessary alteration to the approach to make this numerically practical was the implementation of a continuous probability estimator that relies on Bernstein polynomials for the transformation that takes place in the projections. Stopping criteria were updated to incorporate a bootstrap t test that compares data sampled from a multivariate Gaussian distribution with the data undergoing transformation.     

Using Diurnal Temperature Signals to Infer Vertical Groundwater‐Surface Water Exchange

Heat is a powerful tracer to quantify fluid exchange between surface water and groundwater. Temperature time series can be used to estimate pore water fluid flux, and techniques can be employed to extend these estimates to produce detailed plan‐view flux maps. Key advantages of heat tracing include cost‐effective sensors and ease of data collection and interpretation, without the need for expensive and time‐consuming laboratory analyses or induced tracers. While the collection of temperature data in saturated sediments is relatively straightforward, several factors influence the reliability of flux estimates that are based on time series analysis (diurnal signals) of recorded temperatures. Sensor resolution and deployment are particularly important in obtaining robust flux estimates in upwelling conditions. Also, processing temperature time series data involves a sequence of complex steps, including filtering temperature signals, selection of appropriate thermal parameters, and selection of the optimal analytical solution for modeling. This review provides a synthesis of heat tracing using diurnal temperature oscillations, including details on optimal sensor selection and deployment, data processing, model parameterization, and an overview of computing tools available. Recent advances in diurnal temperature methods also provide the opportunity to determine local saturated thermal diffusivity, which can improve the accuracy of fluid flux modeling and sensor spacing, which is related to streambed scour and deposition. These parameters can also be used to determine the reliability of flux estimates from the use of heat as a tracer.     

Predicting the displacement of triple pendulum™ bearings in a full‐scale shaking experiment using a three‐dimensional element

The accuracy of a series spring model to predict the peak displacement and displacement history of Triple Pendulum? (TP) bearings in a strongly shaken, full‐scale building is evaluated in this paper. The series spring model was implemented as a self‐contained three‐dimensional TP bearing element in OpenSees and is now available for general use. The TP bearing element contains the option for constant friction or a generalized friction model that accounts for the effect of instantaneous velocity and compression load on the friction coefficient. Comparison between numerical simulation and experimental data of a five‐story steel moment frame building shows that the peak displacement of isolation system can generally be predicted with confidence using a constant friction coefficient model. The friction coefficient model accounting for the effect of axial load and velocity leads to minor improvement over the constant friction coefficient models in some cases. Copyright © 2013 John Wiley & Sons, Ltd.     

Modeling Variably Saturated Subsurface Solute Transport with MODFLOW‐UZF and MT3DMS

The MT3DMS groundwater solute transport model was modified to simulate solute transport in the unsaturated zone by incorporating the unsaturated‐zone flow (UZF1) package developed for MODFLOW. The modified MT3DMS code uses a volume‐averaged approach in which Lagrangian‐based UZF1 fluid fluxes and storage changes are mapped onto a fixed grid. Referred to as UZF‐MT3DMS, the linked model was tested against published benchmarks solved analytically as well as against other published codes, most frequently the U.S. Geological Survey's Variably‐Saturated Two‐Dimensional Flow and Transport Model. Results from a suite of test cases demonstrate that the modified code accurately simulates solute advection, dispersion, and reaction in the unsaturated zone. Two‐ and three‐dimensional simulations also were investigated to ensure unsaturated‐saturated zone interaction was simulated correctly. Because the UZF1 solution is analytical, large‐scale flow and transport investigations can be performed free from the computational and data burdens required by numerical solutions to Richards' equation. Results demonstrate that significant simulation runtime savings can be achieved with UZF‐MT3DMS, an important development when hundreds or thousands of model runs are required during parameter estimation and uncertainty analysis. Three‐dimensional variably saturated flow and transport simulations revealed UZF‐MT3DMS to have runtimes that are less than one tenth of the time required by models that rely on Richards' equation. Given its accuracy and efficiency, and the wide‐spread use of both MODFLOW and MT3DMS, the added capability of unsaturated‐zone transport in this familiar modeling framework stands to benefit a broad user‐ship.