Caldera formation by magma withdrawal from a reservoir beneath a volcanic edifice

Caldera formation has been explained by magma withdrawal from a crustal reservoir, but little is known about the conditions that lead to the critical reservoir pressure for collapse. During an eruption, the reservoir pressure is constrained to lie within a finite range: it cannot exceed the threshold value for eruption, and cannot decrease below another threshold value such that feeder dykes get shut by the confining pressure, which stops the eruption. For caldera collapse to occur, the critical reservoir pressure for roof failure must therefore be within this operating range. We use an analytical elastic model to evaluate the changes of reservoir pressure that are required for failure of roof rocks above the reservoir with and without a volcanic edifice at Earth's surface. With no edifice at Earth's surface, faulting in the roof region can only occur in the initial phase of reservoir inflation and affects a very small part of the focal area. Such conditions do not allow caldera collapse. With a volcanic edifice, large tensile stresses develop in the roof region, whose magnitude increase as the reservoir deflates during an eruption. The edifice size must exceed a threshold value for failure of the roof region before the end of eruption. The largest tensile stresses are reached at Earth's surface, indicating that faulting starts there. Failure affects an area whose horizontal dimensions depend on edifice and chamber dimensions. For small and deep reservoirs, failure conditions cannot be achieved even if the edifice is very large. Quantitative predictions are consistent with observations on a number of volcanoes.     

Amphibole Ar/Ar Geochronology from the Okcheon Metamorphic Belt, South Korea and its Tectonic Implications

An understanding of the Okcheon Metamorphic Belt (OMB) in South Korea is central to unraveling the tectono-metamorphic evolution of East Asia. Amphibole-bearing rocks in the OMB occur as calcsilicate layers and lenses in psammitic rocks, in the psammitic rocks themselves, and in the mafic volcanic layers and intrusives. Most amphiboles fail to show ⁴⁰Ar/³⁹Ar plateau ages; those that do have ages ranging from 132 to 975 Ma. The disturbed age pattern and wide variation in ⁴⁰Ar/³⁹Ar ages can be related to metamorphic grade, retrograde chemical reactions, excess Ar and amphibole composition. The oldest age (975 Ma) can be interpreted either as an old igneous or metamorphic age predating sedimentation or a false age caused by excess Ar. The youngest age of 132 Ma and the disturbed age pattern found in amphiboles from rocks located close to Jurassic granitoids are the result of retrograde thermal metamorphic effects accompanying intrusion of the granitoids. Some medium- or coarse-grained amphiboles in the calcsilicates are aggregates of fine-grained crystals. As a result, they are heterogeneous and prove to be readily affected by excess Ar. A disturbed age pattern in amphiboles from the calcsilicates occurring in the high-grade metamorphic zone may also be the product of excess Ar. On the other hand, the disturbed pattern of amphiboles present in the calcsilicates from the low-grade metamorphic zone could arise from both excess Ar and mixed ages. However, amphiboles from psammitic rocks and some calcsilicates in the high-grade metamorphic zone and in intrusive metabasites display real plateau ages of 237 to 261 Ma. The temperature conditions in the high-grade metamorphic zone were higher than the argon closing temperature for amphibole, and the amphiboles in this zone give plateau ages only when they are homogeneous in composition, lack excess Ar, and have not been thermally affected by intrusion of the granitoids. The unmodified ⁴⁰Ar/³⁹Ar ages prove rather younger than the age of the Late Paleozoic metamorphic event of 280 to 300 Ma, but they are close to muscovite K-Ar ages of 263 to 277 Ma. These ⁴⁰Ar/³⁹Ar amphibole ages are interpreted as the time of cooling that followed the main regional, intermediate-P/T metamorphic climax. The results demonstrate that interpretation of ⁴⁰Ar/³⁹Ar amphibole ages in an area subjected to several metamorphic events can be accomplished only by undertaking a thorough tectono-metamorphic study, accompanied by detailed chemical analysis of the amphiboles.     

Development of a soil moisture‐based distributed hydrologic model for determining hydrologically based critical source areas

A simple grid cell‐based distributed hydrologic model was developed to provide spatial information on hydrologic components for determining hydrologically based critical source areas. The model represents the critical process (soil moisture variation) to run‐off generation accounting for both local and global water balance. In this way, it simulates both infiltration excess run‐off and saturation excess run‐off. The model was tested by multisite and multivariable evaluation on the 50‐km² Little River Experimental Watershed I in Georgia, U.S. and 2 smaller nested subwatersheds. Water balance, hydrograph, and soil moisture were simulated and compared to observed data. For streamflow calibration, the daily Nash‐Sutcliffe coefficient was 0.78 at the watershed outlet and 0.56 and 0.75 at the 2 nested subwatersheds. For the validation period, the Nash‐Sutcliffe coefficients were 0.79 at the watershed outlet and 0.85 and 0.83 at the 2 subwatersheds. The per cent bias was less than 15% for all sites. For soil moisture, the model also predicted the rising and declining trends at 4 of the 5 measurement sites. The spatial distribution of surface run‐off simulated by the model was mainly controlled by local characteristics (precipitation, soil properties, and land cover) on dry days and by global watershed characteristics (relative position within the watershed and hydrologic connectivity) on wet days when saturation excess run‐off was simulated. The spatial details of run‐off generation and travel time along flow paths provided by the model are helpful for watershed managers to further identify critical source areas of non‐point source pollution and develop best management practices.     

Analytical modeling on consolidation of stiffened deep mixed column-reinforced soft soil under embankment

Stiffened deep mixed (SDM) column is a new ground improvement technique to improve soft soil, which can be used to increase bearing capacity, reduce deformation, and enhance stability of soft soil. This technique has been successfully adopted to support the highway and railway embankments over soft soils in China and other countries. However, there have been limited investigations on its consolidation under embankment loading. This paper developed an analytical solution for the consolidation of embankment over soft soil with SDM column in which core pile is equal to or shorter than outer DM column. The consolidation problem was simplified as a consolidation of composite soil considering the load shear effect of core pile. The developed solution was verified by a comparison with the results computed by three-dimensional (3-D) finite element analysis. A parametric study based on the derived solution was conducted to investigate influence factors—length of core pile, diameter of core pile, diameter of SDM column, modulus of DM column, and permeability coefficient of DM column—on the consolidation behavior of SDM column-supported embankment over soft soil. The developed solution was applied to a case history of SDM column-supported embankment, and a good agreement was found between the predictions and the field measurements.     

A cavity expansion–based solution for interpretation of CPTu data in soils under partially drained conditions

A cavity expansion–based solution is proposed in this paper for the interpretation of CPTu data under a partially drained condition. Variations of the normalized cone tip resistance, cone factor, and undrained-drained resistance ratio are examined with different initial specific volume and overconsolidation ratio, based on the exact solutions of both undrained and drained cavity expansion in CASM, which is a unified state parameter model for clay and sand. A drainage index is proposed to represent the partially drained condition, and the critical state after expansion and stress paths of cavity expansion are therefore predicted by estimating a virtual plastic region and assuming a drainage-index–based mapping technique. The stress paths and distributions of stresses and specific volume are investigated for different values of drainage index, which are also related to the penetration velocity with comparisons of experimental data and numerical results. The subsequent consolidation after penetration is thus predicted with the assumption of constant deviatoric stress during dissipation of the excess pore pressure. Both spherical and cylindrical consolidations are compared for dissipation around the cone tip and the probe shaft, respectively. The effects of overconsolidation ratio on the stress paths and the distributions of excess pore pressure and specific volume are then thoroughly investigated. The proposed solution and the findings would contribute to the interpretation of CPTu tests under a random drained condition, as well as the analysis of pile installation and the subsequent consolidation.     

Spatial and temporal variability of stable isotopes (δ18O and δ2H) in surface waters of arid,mountainous Central Asia

Characterization of spatial and temporal variability of stable isotopes (δ¹⁸O and δ²H) of surface waters is essential to interpret hydrological processes and establish modern isotope–elevation gradients across mountainous terrains. Here, we present stable isotope data for river waters across Kyrgyzstan. River water isotopes exhibit substantial spatial heterogeneity among different watersheds in Kyrgyzstan. Higher river water isotope values were found mainly in the Issyk‐Kul Lake watershed, whereas waters in the Son‐Kul Lake watershed display lower values. Results show a close δ¹⁸O–δ²H relation between river water and the local meteoric water line, implying that river water experiences little evaporative enrichment. River water from the high‐elevation regions (e.g., Naryn and Son‐Kul Lake watershed) had the most negative isotope values, implying that river water is dominated by snowmelt. Higher deuterium excess (average d = 13.9‰) in river water probably represents the isotopic signature of combined contributions from direct precipitation and glacier melt in stream discharge across Kyrgyzstan. A significant relationship between river water δ¹⁸O and elevation was observed with a vertical lapse rate of 0.13‰/100 m. These findings provide crucial information about hydrological processes across Kyrgyzstan and contribute to a better understanding of the paleoclimate/elevation reconstruction of this region.     

Alteration of hydrological processes and streamflow with juniper (Juniperus virginiana) encroachment in a mesic grassland catchment

To understand the effect of woody plant encroachment on hydrological processes of mesic grasslands, we quantified infiltration capacity in situ, the temporal changes in soil water storage, and streamflow of a grassland catchment and a catchment heavily encroached by juniper (Juniperus virginiana, eastern redcedar) in previously cultivated, non‐karst substrate grasslands in north‐central Oklahoma for 3 years. The initial and steady‐state infiltration rates under the juniper canopy were nearly triple to that of the grassland catchment and were intermediate in the intercanopy spaces within the encroached catchment. Soil water content and soil water storage on the encroached catchment were generally lower than on the grassland catchment, especially when preceding the seasons of peak rainfall in spring and fall. Frequency and magnitude of streamflow events were reduced in the encroached catchment. Annual runoff coefficients for the encroached catchment averaged 2.1%, in contrast to 10.6% for the grassland catchment. Annual streamflow duration ranged from 80 to 250 h for the encroached catchment compared with 600 to 800 h for the grassland catchment. Our results showed that the encroachment of juniper into previously cultivated mesic grasslands fundamentally alters catchment hydrological function. Rapid transformation of mesic grassland to a woodland state with juniper encroachment, if not confined, has the potential to drastically reduce soil water, streamflow and flow duration of ephemeral streams in the Southern Great Plains. Copyright © 2013 John Wiley & Sons, Ltd.     

Dissolved inorganic carbon evolution in neutral discharge from mine tailings piles

We measured the concentrations of dissolved inorganic carbon (DIC) and major ions and the stable carbon isotope ratios of DIC (δ¹³C_DIC) in two creeks discharging from carbonate‐rich sulphide‐containing mine tailings piles. Our aim was to assess downstream carbon evolution of the tailings discharge as it interacted with the atmosphere. The discharge had pH of 6.5–8.1 and was saturated with respect to carbonates. Over the reach of one creek, the DIC concentrations decreased by 1.1 mmol C/l and δ¹³C_DIC increased by ~4.0‰ 200 m from the seep source. The decrease in the DIC concentrations was concomitant with decreases in the partial pressure of CO_2(aq) because of the loss of excess CO_2(aq) from the discharge. The corresponding enrichment in the δ¹³C_DIC is because of kinetic isotope fractionation accompanying the loss of CO_2(g). Over the reach of the other creek, there was no significant decrease in the DIC concentrations or notable changes in the δ¹³C_DIC. The insignificant change in the DIC concentrations and the δ¹³C_DIC is because the first water sample was collected 160 m away from the discharge seep, not accessible during this research. In this case, most of the excess CO_2(aq) was lost before our first sampling station. Our results indicate that neutral discharges from tailings piles quickly lose excess CO_2(aq) to the atmosphere and the DIC becomes enrich in ¹³C. We suggest that a significant amount of carbon cycling in neutral discharges from tailings piles occur close to the locations where the discharge seeps to the surface. Copyright © 2015 John Wiley & Sons, Ltd.     

Constitutive modelling of fine-grained gassy soil: A composite approach

Fine-grained marine sediments containing large undissolved gas bubbles are widely distributed around the world. Presence of the bubbles could degrade the undrained shear strength (s_u ) of the soil, when the gas pressure u_g is relatively high as compared with the effective stress in the saturated soil matrix. Meanwhile, the addition of bubbles may also increase s_u when the difference between u_g and pore water pressure u_w becomes smaller than the water entry value, causing partial water drainage from the saturated matrix into the bubbles (bubble flooding) during globally undrained shearing. A new constitutive model for describing the two competing effects on the stress-strain relationship of fine-grained gassy soil is proposed within the framework of critical state soil mechanics. The gassy soil is considered as a three-phase composite material with compressible cavities, which allows water entry from the saturated matrix. Bubble flooding is modelled by introducing an additional positive volumetric strain increment of the saturated clay matrix, which is dependent on the difference between pore gas and pore water pressure based on experimental observations. A modified hardening law based on that of the modified Cam clay model is employed, which in conjunction with the expression for bubble flooding, can describe both the detrimental and beneficial effects of gas bubbles on soil strength and plastic hardening in shear. Only two extra parameters in addition to those in the modified Cam clay model are used. It is shown that the key features of the stress-strain relationship of three fine-grained gassy soils can be reproduced satisfactorily.     

Tropical precipitation anomalies and d‐excess evolution during El Niño 2014‐16

The last 2014‐16 El Niño event was among the three strongest episodes on record. El Niño considerably changes annual and seasonal precipitation across the tropics. Here, we present a unique stable isotope data set of daily precipitation collected in Costa Rica prior to, during, and after El Niño 2014‐16, in combination with Lagrangian moisture source and precipitation anomaly diagnostics. δ²H composition ranged from ‐129.4 to +18.1 (‰) while δ¹⁸O ranged from ‐17.3 to +1.0 (‰). No significant difference was observed among δ¹⁸O (P=0.186) and δ²H (P=0.664) mean annual compositions. However, mean annual d‐excess showed a significant decreasing trend (from +13.3 to +8.7 ‰) (P<0.001) with values ranging from +26.6 to ‐13.9 ‰ prior to and during the El Niño evolution. The latter decrease in d‐excess can be partly explained by an enhanced moisture flux convergence across the southeastern Caribbean Sea coupled with moisture transport from northern South America by means of an increased Caribbean Low Level Jet regime. During 2014‐15, precipitation deficit across the Pacific domain averaged 46% resulting in a very severe drought; while a 94% precipitation surplus was observed in the Caribbean domain. Understanding these regional moisture transport mechanisms during a strong El Niño event may contribute to a) better understanding of precipitation anomalies in the tropics and b) re‐evaluate past stable isotope interpretations of ENSO events in paleoclimatic archives within the Central America region.