Slab–mantle flow interaction: influence on subduction dynamics and duration

We investigate the influence of mantle flow relative to the lithosphere on subduction dynamics. We use 2D thermo‐mechanical models assuming incompressible non‐Newtonian fluid rheology. Different mantle flow velocities consistent with absolute plate motion models are tested, as well as both directions of flow, either sustaining or opposing slab dip. The effects of different inflow/outflow velocity profiles, slab strengths and upper–lower mantle viscosity contrasts are also evaluated. Slab dip deviations between models with opposite mantle flow directions range from 37° for relatively strong slabs (η_max = 10²⁵ Pa s) to 50° for weaker slabs (η_max = 10²⁴ Pa s), accounting for a significant amount of natural slab dip variability. For imposed mantle flow supporting the slab, the initial stage of slab steepening is followed by a stage of continuous slab dip decrease. This slab shallowing eventually leads to mantle wedge closure, subduction cessation and slab break‐off, possibly driving subduction flips.     

Effect of rainfall uncertainty on the performance of physically based rainfall–runoff models

This paper analyses the effect of rain data uncertainty on the performance of two hydrological models with different spatial structures: a semidistributed and a fully distributed model. The study is performed on a small catchment of 19.6 km² located in the north‐west of Spain, where the arrival of low pressure fronts from the Atlantic Ocean causes highly variable rainfall events. The rainfall fields in this catchment during a series of storm events are estimated using rainfall point measurements. The uncertainty of the estimated fields is quantified using a conditional simulation technique. Discharge and rain data, including the uncertainty of the estimated rainfall fields, are then used to calibrate and validate both hydrological models following the generalized likelihood uncertainty estimation (GLUE) methodology. In the storm events analysed, the two models show similar performance. In all cases, results show that the calibrated distribution of the input parameters narrows when the rain uncertainty is included in the analysis. Otherwise, when rain uncertainty is not considered, the calibration of the input parameters must account for all uncertainty in the rainfall–runoff transformation process. Also, in both models, the uncertainty of the predicted discharges increase in similar magnitude when the uncertainty of rainfall input increase.     

Modelling sediment deposition and phosphorus retention in a river floodplain

Phosphorus (P) is one of the major limiting nutrient in many freshwater ecosystems. During the last decade, attention has been focused on the fluxes of suspended sediment and particulate P through freshwater drainage systems because of severe eutrophication effects in aquatic ecosystems. Hence, the analysis and prediction of phosphorus and sediment dynamics constitute an important element for ecological conservation and restoration of freshwater ecosystems. In that sense, the development of a suitable prediction model is justified, and the present work is devoted to the validation and application of a predictive soluble reactive phosphorus (SRP) uptake and sedimentation models, to a real riparian system of the middle Ebro river floodplain. Both models are coupled to a fully distributed two‐dimensional shallow‐water flow numerical model. The SRP uptake model is validated using data from three field experiments. The model predictions show a good accuracy for SRP concentration, where the linear regressions between measured and calculated values of the three experiments were significant (r² ≥ 0.62; p ≤ 0.05), and a Nash–Sutcliffe coefficient (E) that ranged from 0.54 to 0.62. The sedimentation model is validated using field data collected during two real flooding events within the same river reach. The comparison between calculated and measured sediment depositions showed a significant linear regression (p ≤ 0.05; r² = 0.97) and an E that ranged from 0.63 to 0.78. Subsequently, the complete model that includes flow dynamics, solute transport, SRP uptake and sedimentation is used to simulate and analyse floodplain sediment deposition, river nutrient contribution and SRP uptake. According to this analysis, the main SRP uptake process appears to be the sediment sorption. The analysis also reveals the presence of a lateral gradient of hydrological connectivity that decreases with distance from the river and controls the river matter contribution to the floodplain. Copyright © 2014 John Wiley & Sons, Ltd.     

OSL signal resetting in young deposits determined with a pulsed photon‐stimulated luminescence (PPSL) unit

Optical Stimulated Luminescence (OSL) is a technique that can be used for dating geological materials deposited within the last half‐million years, including sediments transported by air, water or gravity, as well as rocks heated at high temperatures. Recently, several studies have shown that OSL can also provide information on sediment transport. The pulsed photon‐stimulated luminescence (PPSL) unit (also known as a portable OSL reader) developed by the Scottish Universities Environmental Research Centre is an instrument designed to read luminescence signals from bulk (untreated) sediment samples comprising poly‐mineral and poly‐grain fractions. In this contribution, we evaluate the potential of the PPSL unit to assess the degree of OSL signal resetting in 27 young deposits (<2 ka) transported by different geomorphic agents in volcanic, coastal and fluvial depositional settings located in Mexico. Our results are in agreement with previous findings that used the Risø TL/OSL reader, confirming that sediments transported by debrisflows contain the highest inherited luminesce signals. Infrared stimulation (IRSL) values in volcanic ash, lavas, and sand beach and dune deposits exhibit low scatter. However, with blue stimulation (BLSL) these samples reveal a large degree of scattering, attributed to charge transfer in the case of the coastal deposits and to the low sensitivity of quartz in the case of volcanic material. The luminescence signals of fluvial sediments exhibit a highly scattered distribution in both IRSL and BLSL. We conclude that the use of a PPSL unit is a simple approach to assess the degree of OSL signal resetting in deposits sourced from different geological environments. This research contributes to previous studies that have investigated new applications of the PPSL unit to assist in OSL dating of geological materials.     

Importance of mixotrophic nanoplankton in Aysén Fjord (Southern Chile) during austral winter

Mixotrophy, the combination of autotrophic and heterotrophic nutrition in the same organism, is widespread in planktonic algae. Several reports from temperate and high-latitude fjords in Scandinavia suggest the occurrence of a niche in late summer and autumn during post-bloom conditions in which mixotrophic algae can become important grazers in pelagic ecosystems, accessing the nutrients bound in their prey to overcome nutrient limitation. Here, we experimentally determined the trophic modes and bacterivory rates for the nanoplankton community (2–20 μm) in Aysén Fjord located in the Chilean Northern Patagonia during two contrasting seasons: winter and spring. While mixotrophic nanoplankton was virtually absent from the system in spring, in winter at occasions it even constituted the dominant trophic group of the nanoplankton with abundances of >900 cells mL^?1. This indicates a second niche for mixotrophs in winter, when mixotrophy allows overcoming light limitation.     

Calc-alkaline rear-arc magmatism in the Fuegian Andes: Implications for the mid-cretaceous tectonomagmatic evolution of southernmost South America

The magmatic arc of the Fuegian Andes is composed mostly of Upper Mesozoic to Cenozoic calc-alkaline plutons and subordinated lavas. To the rear arc, however, isolated mid-Cretaceous monzonitic plutons and small calc-alkaline dykes and sills crop out. This calc-alkaline unit (the Ushuaia Peninsula Andesites, UPA) includes hornblende-rich, porphyritic quartz meladiorites, granodiorites, andesites, dacites and lamprophyres. Radiometric dating and cross-cutting relationships indicate that UPA is younger than the monzonitic suite. The geochemistry of UPA is medium to high K, with high LILE (Ba 500–2000 ppm, Sr 800–1400 ppm), HFSE (Th 7–23 ppm, Nb 7–13 ppm, Ta 0.5–1.1 ppm) and LREE (La 16–51 ppm) contents, along with relatively low HREE (Yb 1.7–1.3 ppm) and Y (9–19 ppm). The similar mineralogy and geochemistry of all UPA rocks suggest they evolved from a common parental magma, by low pressure crystal fractionation, without significant crustal assimilation. A pure Rayleigh fractionation model indicates that 60–65% of crystal fractionation of 60% hornblende + 34% plagioclase + 4% clinopyroxene + 1% Fe-Ti oxide, apatite and sphene (a paragenesis similar of UPA mafic rocks) can explain evolution from lamprophyres to dacites. The UPA has higher LILE, HFSE and LREE, and lower HREE and Y than the calc-alkaline plutons and lavas of the volcanic front. The HREE and Y are lower than in the potassic plutons as well. High concentrations of Th, Nb, Ta, Zr, Hf, LREE and Ce/Pb, and low U/Th, Ba/Th ratios in UPA, even in the least differentiated samples, suggest contributions from subducted sediments to the mantle source. On the other hand, relatively low HREE and Y, high LREE/HREE (La/Yb 11–38) ratios and Nb-Ta contents can be interpreted as mantle metasomatism by partial melts of either subducted garnetiferous oceanic sediment or basalt as well. Additionally, high LILE content in UPA, similar to the potassic plutons, suggests also a mantle wedge previously metasomatized by potassic parental magmas in their route to crustal levels. Therefore, UPA represents a unique suite in the Fuegian arc generated in a multiple hybridized source. UPA generation is related to a transition from normal to flat subduction which additionally caused the widening and landward migration of the magmatic arc, as well as crustal deformation. Rear-arc magmatism endured ca. 22 m.y.; afterwards, calc-alkaline magmatism remained at the volcanic front.     

Latitudinal patterns of export production recorded in surface sediments of the Chilean Patagonian fjords (41–55°S) as a response to water column productivity

The Chilean Patagonian fjords region (41–56°S) is characterized by highly complex geomorphology and hydrographic conditions, and strong seasonal and latitudinal patterns in precipitation, freshwater discharge, glacier coverage, and light regime; all of these directly affect biological production in the water column. In this study, we compiled published and new information on water column properties (primary production, nutrients) and surface sediment characteristics (biogenic opal, organic carbon, molar C/N, bulk sedimentary δ¹³C_org) from the Chilean Patagonian fjords between 41°S and 55°S, describing herein the latitudinal pattern of water column productivity and its imprint in the underlying sediments. Based on information collected at 188 water column and 118 sediment sampling sites, we grouped the Chilean fjords into four main zones: Inner Sea of Chiloé (41° to ～44°S), Northern Patagonia (44° to ～47°S), Central Patagonia (48–51°S), and Southern Patagonia (Magellan Strait region between 52° and 55°S). Primary production in the Chilean Patagonian fjords was the highest in spring–summer, reflecting the seasonal pattern of water column productivity. A clear north–south latitudinal pattern in primary production was observed, with the highest average spring and summer estimates in the Inner Sea of Chiloé (2427 and 5860 mg C m^?2 d^?1) and Northern Patagonia (1667 and 2616 mg C m^?2 d^?1). This pattern was closely related to the higher availability of nutrients, greater solar radiation, and extended photoperiod during the productive season in these two zones. The lowest spring value was found in Caleta Tortel, Central Patagonia (91 mg C m^?2 d^?1), a site heavily influenced by glacier meltwater and river discharge loaded with glacial sediments. Biogenic opal, an important constituent of the Chilean fjord surface sediments (Si_OPAL ～1–13%), reproduced the general north–south pattern of primary production and was directly related to water column silicic acid concentrations. Surface sediments were also rich in organic carbon content and the highest values corresponded to locations far away from glacier influence, sites within fjords, and/or semi-enclosed and protected basins, reflecting both autochthonous (water column productivity) and allochthonous sources (contribution of terrestrial organic matter from fluvial input to the fjords). A gradient was observed from the more oceanic sites to the fjord heads (west–east) in terms of bulk sedimentary δ¹³C_org and C/N ratios; the more depleted (δ¹³C_org ?26‰) and higher C/N (23) values corresponded to areas close to rivers and glaciers. A comparison of the Chilean Patagonian fjords with other fjord systems in the world revealed high variability in primary production for all fjord systems as well as similar surface sediment geochemistry due to the mixing of marine and terrestrial organic carbon.     

A multi-tracer approach to delineate groundwater dynamics in the Rio Actopan Basin,Veracruz State,Mexico

Geochemistry and environmental tracers were used to understand groundwater resources, recharge processes, and potential sources of contamination in the Rio Actopan Basin, Veracruz State, Mexico. Total dissolved solids are lower in wells and springs located in the basin uplands compared with those closer to the coast, likely associated with rock/water interaction. Geochemical results also indicate some saltwater intrusion near the coast and increased nitrate near urban centers. Stable isotopes show that precipitation is the source of recharge to the groundwater system. Interestingly, some high-elevation springs are more isotopically enriched than average annual precipitation at higher elevations, indicating preferential recharge during the drier but cooler winter months when evapotranspiration is reduced. In contrast, groundwater below 1,200 m elevation is more isotopically depleted than average precipitation, indicating recharge occurring at much higher elevation than the sampling site. Relatively cool recharge temperatures, derived from noble gas measurements at four sites (11–20 °C), also suggest higher elevation recharge. Environmental tracers indicate that groundwater residence time in the basin ranges from 12,000 years to modern. While this large range shows varying groundwater flowpaths and travel times, ages using different tracer methods (¹⁴C, ³H/³He, CFCs) were generally consistent. Comparing multiple tracers such as CFC-12 with CFC-113 indicates piston-flow to some discharge points, yet binary mixing of young and older groundwater at other points. In summary, groundwater within the Rio Actopan Basin watershed is relatively young (Holocene) and the majority of recharge occurs in the basin uplands and moves towards the coast.     

A robust debris-flow and GLOF risk management strategy for a data-scarce catchment in Santa Teresa,Peru

The town of Santa Teresa (Cusco Region, Peru) has been affected by several large debris-flow events in the recent past, which destroyed parts of the town and resulted in a resettlement of the municipality. Here, we present a risk analysis and a risk management strategy for debris-flows and glacier lake outbursts in the Sacsara catchment. Data scarcity and limited understanding of both physical and social processes impede a full quantitative risk assessment. Therefore, a bottom-up approach is chosen in order to establish an integrated risk management strategy that is robust against uncertainties in the risk analysis. With the Rapid Mass Movement Simulation (RAMMS) model, a reconstruction of a major event from 1998 in the Sacsara catchment is calculated, including a sensitivity analysis for various model parameters. Based on the simulation results, potential future debris-flows scenarios of different magnitudes, including outbursts of two glacier lakes, are modeled for assessing the hazard. For the local communities in the catchment, the hazard assessment is complemented by the analysis of high-resolution satellite imagery and fieldwork. Physical, social, economic, and institutional vulnerability are considered for the vulnerability assessment, and risk is eventually evaluated by crossing the local hazard maps with the vulnerability. Based on this risk analysis, a risk management strategy is developed, consisting of three complementing elements: (i) standardized risk sheets for the communities; (ii) activities with the local population and authorities to increase social and institutional preparedness; and (iii) a simple Early Warning System. By combining scientific, technical, and social aspects, this work is an example of a framework for an integrated risk management strategy in a data scarce, remote mountain catchment in a developing country.