SOIL,SUBSOIL, AND FOREST COMPOSITION IN SOUTH-CENTRAL MICHIGAN,USA

Although soil texture is an important predictor of upland forest composition in southern Michigan, the relationship in some areas is weak and other environmental factors may exert stronger control on species composition. One of these may be concealed subsoil lithology where coarse-textured glaciofluvial sediments are buried beneath finer-textured till. Samples of forest composition, soil, and subsoil characteristics from 48 less-disturbed woodlots in the south-central Lower Peninsula reveal that upland oaks and hickories are predominant where the till veneer is thin and sand and gravel are relatively near the surface. Where till is relatively thick and outwash is deeply buried, woodlot composition is sugar maple-beech. In each of these cases, soil texture is unrelated to stand composition. The likely effect of buried near-surface glaciofluvial sediments is moisture stress such that xeric species dominate these drier sites. The influence of subsoil on forest geography is a significant new finding with regard to southern Michigan and has possible application elsewhere in glaciated eastern North America.     

Application of the Ti-in-quartz thermobarometer to rutile-free systems. Reply to: a comment on: ‘TitaniQ under pressure: the effect of pressure and temperature on the solubility of Ti in quartz’ by Thomas et al.

The premise of the Wilson et al. comment is that the Ti-in-quartz solubility calibration (Thomas et al. in Contrib Mineral Petrol 160:743–759, 2010) is fundamentally flawed. They reach this conclusion because P–T estimates using the Ti-in-quartz calibration differ from their previous interpretations for crystallization conditions of the Bishop and Oruanui rhyolites. If correct, this assertion has far-reaching implications, so a careful assessment of the Wilson et al. reasoning is warranted. Application of the Ti-in-quartz calibration as a thermobarometer in rutile-free rocks requires an estimation of TiO₂ activity in the liquid ( (liquid–rutile); referenced to rutile saturation) and an independent constraint on either P or T to obtain the crystallization temperature or pressure, respectively. The foundation of Wilson et al.’s argument is that temperature estimates obtained from Fe–Ti oxide thermometry accurately reflect crystallization conditions of quartz in the two rhyolites discussed. We maintain that our experimental approach is sound, the thermodynamic basis of the Ti-in-quartz calibration is fundamentally correct, and our experimental results are robust and reproducible. We suggest that the reason Wilson et al. obtain implausible pressure estimates is because estimates for T and they used as input values for the Ti-in-quartz calibration are demonstrably too high. Numerous studies show that Fe–Ti oxide temperature estimates of some rhyolites are substantially higher than those predicted by well-constrained phase equilibria. In this reply, we show that when reasonable input values for T and (liquid–rutile) are used, pressure estimates obtained from the Ti-in-quartz calibration are well aligned with phase equilibria and essentially identical to melt inclusion volatile saturation pressures.     

Disconnected Surface Water and Groundwater: From Theory to Practice

When describing the hydraulic relationship between rivers and aquifers, the term disconnected is frequently misunderstood or used in an incorrect way. The problem is compounded by the fact that there is no definitive literature on the topic of disconnected surface water and groundwater. We aim at closing this gap and begin the discussion with a short introduction to the historical background of the terminology. Even though a conceptual illustration of a disconnected system was published by Meinzer (1923) , it is only within the last few years that the underlying physics of the disconnection process has been described. The importance of disconnected systems, however, is not widely appreciated. Although rarely explicitly stated, many approaches for predicting the impacts of groundwater development on surface water resources assume full connection. Furthermore, management policies often suggest that surface water and groundwater should only be managed jointly if they are connected. However, although lowering the water table beneath a disconnected section of a river will not change the infiltration rate at that point, it can increase the length of stream that is disconnected. Because knowing the state of connection is of fundamental importance for sustainable water management, robust field methods that allow the identification of the state of connection are required. Currently, disconnection is identified by showing that the infiltration rate from a stream to an underlying aquifer is independent of the water table position or by identifying an unsaturated zone under the stream. More field studies are required to develop better methods for the identification of disconnection and to quantify the implications of heterogeneity and clogging processes in the streambed on disconnection.     

Overview and Assessment of Antarctic Ice-Sheet Mass Balance Estimates: 1992�C2009

Mass balance estimates for the Antarctic Ice Sheet (AIS) in the 2007 report by the Intergovernmental Panel on Climate Change and in more recent reports lie between approximately +50 to ?250 Gt/year for 1992 to 2009. The 300 Gt/year range is approximately 15% of the annual mass input and 0.8 mm/year Sea Level Equivalent (SLE). Two estimates from radar altimeter measurements of elevation change by European Remote-sensing Satellites (ERS) (+28 and ?31 Gt/year) lie in the upper part, whereas estimates from the Input-minus-Output Method (IOM) and the Gravity Recovery and Climate Experiment (GRACE) lie in the lower part (?40 to ?246 Gt/year). We compare the various estimates, discuss the methodology used, and critically assess the results. We also modify the IOM estimate using (1) an alternate extrapolation to estimate the discharge from the non-observed 15% of the periphery, and (2) substitution of input from a field data compilation for input from an atmospheric model in 6% of area. The modified IOM estimate reduces the loss from 136 Gt/year to 13 Gt/year. Two ERS-based estimates, the modified IOM, and a GRACE-based estimate for observations within 1992?C2005 lie in a narrowed range of +27 to ?40 Gt/year, which is about 3% of the annual mass input and only 0.2 mm/year SLE. Our preferred estimate for 1992?C2001 is ?47 Gt/year for West Antarctica, +16 Gt/year for East Antarctica, and ?31 Gt/year overall (+0.1 mm/year SLE), not including part of the Antarctic Peninsula (1.07% of the AIS area). Although recent reports of large and increasing rates of mass loss with time from GRACE-based studies cite agreement with IOM results, our evaluation does not support that conclusion.     

Evolution of a Carboniferous carbonate-hosted sphalerite breccia deposit,Isle of Man

A newly discovered, extensive sphalerite-bearing breccia (~7.5 wt.% Zn) is hosted in dolomitised Carboniferous limestones overlying Ordovician–Silurian metasedimentary rocks on the Isle of Man. Although base metal sulphide deposits have been mined historically on the island, they are nearly all quartz vein deposits in the metamorphic basement. This study investigates the origin of the unusual sphalerite breccia and its relationship to basement-hosted deposits, through a combination of petrographic, cathodoluminescence, fluid inclusion, stable isotope and hydrogeologic modelling techniques. Breccia mineralisation comprises four stages, marked by episodes of structural deformation and abrupt changes in fluid temperature and chemistry. In stage I, high-temperature (T _h > 300°C), high-salinity (20–45 wt.% equiv. NaCl) fluid of likely basement origin deposited a discontinuous quartz vein. This vein was subsequently dismembered during a major brecciation event. Stages II–IV are dominated by open-space filling sphalerite, quartz and dolomite, respectively. Fluid inclusions in these minerals record temperatures of ~105–180°C and salinities of ~15–20 wt.% equiv. NaCl. The δ³⁴S values of sphalerite (6.5–6.9‰ Vienna-Canyon Diablo troilite) are nearly identical to those of ore sulphides from mines in the Lower Palaeozoic metamorphic rocks. The δ¹⁸O values for quartz and dolomite indicate two main fluid sources in the breccia’s hydrothermal system, local Carboniferous-hosted brines (~0.5–6.0‰ Vienna standard mean ocean water) and basement-involved fluids (~5.5–11.5‰). Ore sulphide deposition in the breccia is compatible with the introduction and cooling of a hot, basement-derived fluid that interacted with local sedimentary brines.     

Chronostratigraphy and lake-level changes of Laguna Cari-Laufquén, Río Negro, Argentina

Evidence from shoreline and deep-lake sediments show Laguna Cari-Laufquén, located at 41°S in central Argentina, rose and fell repeatedly during the late Quaternary. Our results show that a deep (> 38 m above modern lake level) lake persisted from no later than 28 ka to 19 ka, with the deepest lake phase from 27 to 22 ka. No evidence of highstands is found after 19 ka until the lake rose briefly in the last millennia to 12 m above the modern lake, before regressing to present levels. Laguna Cari-Laufquén broadly matches other regional records in showing last glacial maximum (LGM) highstands, but contrasts with sub-tropical lake records in South America where the hydrologic maximum occurred during deglaciation (17–10 ka). Our lake record from Cari-Laufquén mimics that of high-latitude records from the Northern Hemisphere. This points to a common cause for lake expansions, likely involving some combination of temperature depression and intensification of storminess in the westerlies belt of both hemispheres during the LGM.     

Isotopic tracers of paleohydrologic change in large lakes of the Bolivian Altiplano

We have developed an ⁸⁷Sr/⁸⁶Sr, ²³⁴U/²³⁸U, and δ¹⁸O data set from carbonates associated with late Quaternary paleolake cycles on the southern Bolivian Altiplano as a tool for tracking and understanding the causes of lake-level fluctuations. Distinctive groupings of ⁸⁷Sr/⁸⁶Sr ratios are observed. Ratios are highest for the Ouki lake cycle (120-95 ka) at 0.70932, lowest for Coipasa lake cycle (12.8-11.4 ka) at 0.70853, and intermediate at 0.70881 to 0.70884 for the Salinas (95-80 ka), Inca Huasi (~ 45 ka), Sajsi (24-20.5 ka), and Tauca (18.1-14.1 ka) lake cycles. These Sr ratios reflect variable contributions from the eastern and western Cordilleras. The Laca hydrologic divide exerts a primary influence on modern and paleolake ⁸⁷Sr/⁸⁶Sr ratios; waters show higher ⁸⁷Sr/⁸⁶Sr ratios north of this divide. Most lake cycles were sustained by slightly more rainfall north of this divide but with minimal input from Lake Titicaca. The Coipasa lake cycle appears to have been sustained mainly by rainfall south of this divide. In contrast, the Ouki lake cycle was an expansive lake, deepest in the northern (Poópo) basin, and spilling southward. These results indicate that regional variability in central Andean wet events can be reconstructed using geochemical patterns from this lake system.     

The Gradenbach Observatory—monitoring deep-seated gravitational slope deformation by geodetic, hydrological, and seismological methods

The Gradenbach mass movement (GMM) is an example of DGSD (deep-seated gravitational slope deformation) in crystalline rocks of the Eastern Alps (12.85°E, 47.00°N). The main body of the GMM covers an area of 1.7 km² and its volume is about 120?×?10⁶ m³. A reconstruction of the deformation history yields a mean displacement of?～?22 m from 1962 to 2011. In 1965/66, 1975, 2001, and 2009 high sliding velocities, exceeding several meters per year, interrupt the quasi-stationary periods of slow movement (≤0.3 m/year). Since 1999 the displacement of the main body of the GMM has been observed by GPS. Time series of extensometer readings, precipitation, snow cover water equivalent, water discharge, and hydrostatic water level observed in boreholes were re-processed and are presented in this paper. Continuous recording of seismic activity by a seismic monitoring network at the GMM began in the summer of 2006. Deformation has been monitored since 2007 by an embedded strain rosette based on fiber optics technology and a local conventional geodetic deformation network. The velocity of the GMM could be modeled to a large extent by a quantitative relation to hydro-meteorological data. During the phase of high sliding velocity in spring 2009, the seismic activity in the area increased significantly. Several types of seismic events were identified with some of them preceding the acceleration of the main body by about 6 weeks. The potential inherent in the Gradenbach Observatory data to supply early warning and hazard estimation is discussed.     

Geophysical Finite-Element Simulation Tool (GeoFEST): Algorithms and Validation for Quasistatic Regional Faulted Crust Problems

GeoFEST (Geophysical Finite Element Simulation Tool) is a two- and three-dimensional finite element software package for the modeling of solid stress and strain in geophysical and other continuum domain applications. It is one of the featured high-performance applications of the NASA QuakeSim project. The program is targeted to be compiled and run on UNIX systems, and is running on diverse systems including sequential and message-passing parallel systems. Solution to the elliptical partial differential equations is obtained by finite element basis sampling, resulting in a sparse linear system primarily solved by conjugate gradient iteration to a tolerance level; on sequential systems a Crout factorization for the direct inversion of the linear system is also supported. The physics models supported include isotropic linear elasticity and both Newtonian and power-law viscoelasticity, via implicit quasi-static time stepping. In addition to triangular, quadrilateral, tetrahedral and hexahedral continuum elements, GeoFEST supports split-node faulting, body forces, and surface tractions. This software and related mesh refinement strategies have been validated on a variety of test cases with rigorous comparison to analytical solutions. These include a box-shaped domain with imposed motion on one surface, a pair of strike slip faults in stepover arrangement, and two community-agreed benchmark cases: a strike slip fault in an enclosing box, and a quarter-domain circular fault problem. Scientific applications of the code include the modeling of static and transient co- and post-seismic earth deformation, Earth response to glacial, atmospheric and hydrological loading, and other scenarios involving the bulk deformation of geologic media.