Geodynamic evolution of a forearc rift in the southernmost Mariana Arc

The southernmost Mariana forearc stretched to accommodate opening of the Mariana Trough backarc basin in late Neogene time, erupting basalts at 3.7–2.7 Ma that are now exposed in the Southeast Mariana Forearc Rift (SEMFR). Today, SEMFR is a broad zone of extension that formed on hydrated, forearc lithosphere and overlies the shallow subducting slab (slab depth ≤ 30–50 km). It comprises NW–SE trending subparallel deeps, 3–16 km wide, that can be traced ≥ ∼30 km from the trench almost to the backarc spreading center, the Malaguana‐Gadao Ridge (MGR). While forearcs are usually underlain by serpentinized harzburgites too cold to melt, SEMFR crust is mostly composed of Pliocene, low‐K basaltic to basaltic andesite lavas that are compositionally similar to arc lavas and backarc basin (BAB) lavas, and thus defines a forearc region that recently witnessed abundant igneous activity in the form of seafloor spreading. SEMFR igneous rocks have low Na₈, Ti₈, and Fe₈, consistent with extensive melting, at ∼23 ± 6.6 km depth and 1239 ± 40°C, by adiabatic decompression of depleted asthenospheric mantle metasomatized by slab‐derived fluids. Stretching of pre‐existing forearc lithosphere allowed BAB‐like mantle to flow along the SEMFR and melt, forming new oceanic crust. Melts interacted with pre‐existing forearc lithosphere during ascent. The SEMFR is no longer magmatically active and post‐magmatic tectonic activity dominates the rift.     

Stratification response of soil water content during rainfall events under different rainfall patterns

Many researchers have studied the influence of rainfall patterns on soil water movement processes using rainfall simulation experiments. However, less attention has been paid to the influence under natural condition. In this paper, rainfall, soil water content (SWC), and soil temperature at 10‐, 20‐, 30‐, 40‐, and 50‐cm depths were simultaneously monitored at 1‐min intervals to measure the variation in SWC (SWCv) in response to rainfall under different rainfall patterns. First, we classified rainfall events into four patterns. During the study period, the main pattern was the advanced rainfall pattern (38% of all rainfall events), whereas the delayed, central, and uniform rainfall patterns had similar frequencies of about 20%. During natural rainfall, rainwater rapidly passed through the top soil layers (10–40 cm) and was accumulated in the bottom layer (50 cm). When a high rainfall pulse occurred, the water storage balance was disturbed, resulting in the drainage of initial soil water from the top layers into the deeper layers. Therefore, the critical function of the top layers and the bottom layers was infiltration and storage, respectively. The source of water stored in the bottom layer was not only rainfall but also the initial soil water in the upper soil layers. Changes in soil temperature at each soil depth were comonitored with SWCv to determine the movement characteristics of soil water under different rainfall patterns. Under the delayed rainfall pattern, preferential flows preferred to occur. Under the other rainfall patterns, matrix flow was the main form of soil water movement. Rainfall amount was a better indicator than rainfall intensity for SWCv in the bottom layer under the delayed rainfall pattern. These results provide insights into the responses of SWCv under different rainfall patterns in northern China.     

Model Simulations on the Variability of Particulate MSA to Non-Sea-Salt Sulfate Ratio in the Marine Environment

A box model was constructed to investigate connections between the particulate MSA to non-sea-salt sulfate ratio, R, and DMS chemistry in a clean marine boundary layer. The simulations demonstrated that R varies widely with particle size, which must be taken into account when interpreting field measurements or comparing them with each other. In addition to DMS gas-phase chemistry, R in the submicron size range was shown to be sensitive to the factors dictating sulfate production via cloud processing, to the removal of SO₂ from the boundary layer by dry deposition and sea-salt oxidation, to the entrainment of SO₂ from the free troposphere, to the relative concentration of sub- and supermicron particles, and to meteorology. Three potential explanations for the increase of R toward high-latitudes during the summer were found: larger MSA yields from DMS oxidation at high latitudes, larger DMSO yields from DMS oxidation followed by the conversion of DMSO to MSA at high latitudes, or lower ambient H₂O₂ concentrations at high latitudes leading to less efficient sulfate production in clouds. Possible reasons for the large seasonal amplitude of R at mid and high latitudes include seasonal changes in the partitioning of DMS oxidation to the OH and NO₃ initiated pathways, seasonal changes in the concentration of species participating the DMS-OH reaction pathway, or the existence of a SO₂ source other than DMS oxidation in the marine boundary layer. Even small anthropogenic perturbations were shown to have a potential to alter the MSA to non-sea-salt sulfate ratio.     

A Nd- and O-isotope study of the REE-rich peralkaline Strange Lake granite: implications for Mesoproterozoic A-type magmatism in the Core Zone (NE-Canada)

It is well established that A-type granites enriched in high field strength elements, such as Zr, Nb and the REE, form in anorogenic tectonic settings. The sources of these elements and the processes controlling their unusual enrichment, however, are still debated. They are addressed here using neodymium and oxygen isotope analyses of samples from the 1.24 Ga Strange Lake pluton in the Paleoproterozoic Core Zone of Québec-Labrador, an A-type granitic body characterized by hyper-enrichment in the REE, Zr, and Nb. Age-corrected εNd values for bulk rock samples and sodic amphiboles (mainly arfvedsonite) from the pluton range from ?0.6 to ?5.7, and ?0.3 to ?5.3, respectively. The εNd values for the Napeu Kainiut quartz monzonite, which hosts the pluton, range from ?4.8 to ?8.1. The ¹⁴⁷Sm/¹⁴⁴Nd ratios of the suite and the host quartz monzonite range from 0.0967 to 0.1659, large variations that can be explained by in situ fractionation of early LREE-minerals (Strange Lake), and late hydrothermal HREE remobilization. Oxygen isotope analyses of quartz of both Strange Lake and the host yielded δ¹⁸O values between +8.2 and +9.1, which are considerably higher than the mantle value of 5.7 ± 0.2‰. Bulk rock oxygen isotope analyses of biotite-gneisses in the vicinity of the Strange Lake pluton yielded δ¹⁸O values of 6.3, 8.6 and 9.6‰. The negative εNd values and positive δ¹⁸O values of the Strange Lake and Napeu Kainiut samples indicate that both magmas experienced considerable crustal contamination. The extent of this contamination was estimated, assuming that the contaminants were sedimentary-derived rocks from the underlying Archean Mistinibi (para-) gneiss complex, which is characterized by low εNd and high δ¹⁸O values. Mixing of 5–15% of a gneiss, having an εNd value of ?15 and a δ¹⁸O value of +11, with a moderately enriched mantle source (εNd = +0.9, δ¹⁸O = +6.3) would produce values similar to those obtained for the Strange Lake granites. Based on analogies between the Nain Plutonic Suite and the Gardar alkaline igneous province (SW-Greenland), we conclude that the Strange Lake pluton and associated REE-mineralized anorogenic bodies formed from a combination of subduction-induced fertilization of the sublithospheric mantle, crustal extension and in situ magma evolution.     

An analytical approach to ascertain saturation‐excess versus infiltration‐excess overland flow in urban and reference landscapes

Uncontrolled overland flow drives flooding, erosion, and contaminant transport, with the severity of these outcomes often amplified in urban areas. In pervious media such as urban soils, overland flow is initiated via either infiltration‐excess (where precipitation rate exceeds infiltration capacity) or saturation‐excess (when precipitation volume exceeds soil profile storage) mechanisms. These processes call for different management strategies, making it important for municipalities to discern between them. In this study, we derived a generalized one‐dimensional model that distinguishes between infiltration‐excess overland flow (IEOF) and saturation‐excess overland flow (SEOF) using Green–Ampt infiltration concepts. Next, we applied this model to estimate overland flow generation from pervious areas in 11 U.S. cities. We used rainfall forcing that represented low‐ and high‐intensity events and compared responses among measured urban versus predevelopment reference soil hydraulic properties. The derivation showed that the propensity for IEOF versus SEOF is related to the equivalence between two nondimensional ratios: (a) precipitation rate to depth‐weighted hydraulic conductivity and (b) depth of soil profile restrictive layer to soil capillary potential. Across all cities, reference soil profiles were associated with greater IEOF for the high‐intensity set of storms, and urbanized soil profiles tended towards production of SEOF during the lower intensity set of storms. Urban soils produced more cumulative overland flow as a fraction of cumulative precipitation than did reference soils, particularly under conditions associated with SEOF. These results will assist cities in identifying the type and extent of interventions needed to manage storm water produced from pervious areas.     

The role of hydrothermal fluids in sedimentation in saline alkaline lakes: Evidence from Nasikie Engida,Kenya Rift Valley

Saline alkaline lakes that precipitate sodium carbonate evaporites are most common in volcanic terrains in semi‐arid environments. Processes that lead to trona precipitation are poorly understood compared to those in sulphate‐dominated and chloride‐dominated lake brines. Nasikie Engida (Little Magadi) in the southern Kenya Rift shows the initial stages of soda evaporite formation. This small shallow (<2 m deep; 7 km long) lake is recharged by alkaline hot springs and seasonal runoff but unlike neighbouring Lake Magadi is perennial. This study aims to understand modern sedimentary and geochemical processes in Nasikie Engida and to assess the importance of geothermal fluids in evaporite formation. Perennial hot‐spring inflow waters along the northern shoreline evaporate and become saturated with respect to nahcolite and trona, which precipitate in the southern part of the lake, up to 6 km from the hot springs. Nahcolite (NaHCO₃) forms bladed crystals that nucleate on the lake floor. Trona (Na₂CO₃·NaHCO₃·2H₂O) precipitates from more concentrated brines as rafts and as bottom‐nucleated shrubs of acicular crystals that coalesce laterally to form bedded trona. Many processes modify the fluid composition as it evolves. Silica is removed as gels and by early diagenetic reactions and diatoms. Sulphate is depleted by bacterial reduction. Potassium and chloride, of moderate concentration, remain conservative in the brine. Clastic sedimentation is relatively minor because of the predominant hydrothermal inflow. Nahcolite precipitates when and where pCO₂ is high, notably near sublacustrine spring discharge. Results from Nasikie Engida show that hot spring discharge has maintained the lake for at least 2 kyr, and that the evaporite formation is strongly influenced by local discharge of carbon dioxide. Brine evolution and evaporite deposition at Nasikie Engida help to explain conditions under which ancient sodium carbonate evaporites formed, including those in other East African rift basins, the Eocene Green River Formation (western USA), and elsewhere.     

Bistable plant–soil dynamics and biogenic controls on the soil production function

Soil formation results from opposite processes of bedrock weathering and erosion, whose balance may be altered by natural events and human activities, resulting in reduced soil depth and function. The impacts of vegetation on soil production and erosion and the feedbacks between soil formation and vegetation growth are only beginning to be explored quantitatively. Since plants require suitable soil environments, disturbed soil states may support less vegetation, leading to a downward spiral of increased erosion and decline in ecosystem function. We explore these feedbacks with a minimal model of the soil–plant system described by two coupled nonlinear differential equations, which include key feedbacks, such as plant‐driven soil production and erosion inhibition. We show that sufficiently strong positive plant–soil feedback can lead to a ‘humped’ soil production function, a necessary condition for soil depth bistability when erosion is assumed to vary monotonically with vegetation biomass. In bistable plant–soil systems, the sustainable soil condition engineered by plants is only accessible above a threshold vegetation biomass and occurs in environments where the high potential rate of erosion exerts a strong control on soil production and erosion. Vegetation removal for agriculture reduces the stabilizing effect of vegetation and lowers the system resilience, thereby increasing the likelihood of transition to a degraded soil state. Copyright © 2016 John Wiley & Sons, Ltd.     

Geophysical constraints on deep weathering and water storage potential in the Southern Sierra Critical Zone Observatory

The conversion of bedrock to regolith marks the inception of critical zone processes, but the factors that regulate it remain poorly understood. Although the thickness and degree of weathering of regolith are widely thought to be important regulators of the development of regolith and its water‐storage potential, the functional relationships between regolith properties and the processes that generate it remain poorly documented. This is due in part to the fact that regolith is difficult to characterize by direct observations over the broad scales needed for process‐based understanding of the critical zone. Here we use seismic refraction and resistivity imaging techniques to estimate variations in regolith thickness and porosity across a forested slope and swampy meadow in the Southern Sierra Critical Zone Observatory (SSCZO). Inferred seismic velocities and electrical resistivities image a weathering zone ranging in thickness from 10 to 35 m (average = 23 m) along one intensively studied transect. The inferred weathering zone consists of roughly equal thicknesses of saprolite (P‐velocity < 2 km s^?1) and moderately weathered bedrock (P‐velocity = 2–4 km s^?1). A minimum‐porosity model assuming dry pore space shows porosities as high as 50% near the surface, decreasing to near zero at the base of weathered rock. Physical properties of saprolite samples from hand augering and push cores are consistent with our rock physics model when variations in pore saturation are taken into account. Our results indicate that saprolite is a crucial reservoir of water, potentially storing an average of 3 m³ m^?2 of water along a forested slope in the headwaters of the SSCZO. When coupled with published erosion rates from cosmogenic nuclides, our geophysical estimates of weathering zone thickness imply regolith residence times on the order of 10⁵ years. Thus, soils at the surface today may integrate weathering over glacial–interglacial fluctuations in climate. Copyright © 2013 John Wiley & Sons, Ltd.     

On the accuracy of the GPS L2 observable for ionospheric monitoring

The introduction of the unencrypted global positioning system (GPS) L2 civil (L2C) signal has the potential to improve measurements made with the L2 frequency, an important observable in GPS-based ionospheric research and monitoring. Recent work has shown significant differences between the legacy L2P(Y) and L2C-derived total electron content rate of change index (ROTI). This difference is observed between L2P(Y) and L2C-derived ROTI with certain receiver models and between zero-baseline receiver pairs. We discuss the likely cause for these differences: L1-aided tracking used to track both the L2P(Y) and L2C signals. We also present L2C data that are confirmed to be from tracking independent of L1. Using the ionospheric-free linear combination, we show that the independently tracked carrier phase dynamics are significantly more accurate than the L1-aided observables. This result is confirmed by comparing the behavior of the L2C and L2P(Y) carrier phase observables upon a sudden antenna rotation.     

Earthquake relocation in the Central Alborz region of Iran using a non-linear probabilistic method

In this study, we calculate accurate absolute locations for nearly 3,000 shallow earthquakes (≤20 km depth) that occurred from 1996 to 2010 in the Central Alborz region of northern Iran using a non-linear probabilistic relocation algorithm on a local scale. We aim to produce a consistent dataset with a realistic assessment of location errors using probabilistic hypocenter probability density functions. Our results indicate significant improvement in hypocenter locations and far less scattering than in the routine earthquake catalog. According to our results, 816 earthquakes have horizontal uncertainties in the 0.5–3.0 km range, and 981 earthquakes are relocated with focal-depth errors less than 3.0 km, even with a suboptimal network geometry. Earthquake relocated are tightly clustered in the eastern Tehran region and are mainly associated with active faults in the study area (the Mosha and Garmsar faults). Strong historical earthquakes have occurred along the Mosha and Garmsar faults, and the relocated earthquakes along these faults show clear north-dipping structures and align along east–west lineations, consistent with the predominant trend of faults within the study region. After event relocation, all seismicity lies in the upper 20 km of the crust, and no deep seismicity (>20 km depth) has been observed. In many circumstances, the seismicity at depth does not correlate with surface faulting, suggesting that the faulting at depth does not directly offset overlying sediments.