Runoff and sediment yield from land under various uses in a Mediterranean mountain area: long‐term results from an experimental station

In this study we analyzed runoff and sediment yield from land under various traditional and current land uses in Mediterranean mountain areas, using long‐term data from an experimental station in the Aísa Valley, Central Spanish Pyrenees. Monitoring at this station has provided 20 years of data that can help explain the hydrological and geomorphological changes that have been observed at larger spatial scales, and also the changes that have occurred to some of the most characteristic landscapes of the Mediterranean middle mountains. In spite of the problems associated with the use of small experimental plots, the results obtained are consistent with other studies in the Mediterranean region, and confirm the strong influence of land use changes on runoff generation and sediment yield. The results indicate that: (i) cereal cultivation on steep slopes (both alternating cereal cultivation and fallow on sloping fields and shifting agriculture on the steepest slopes) represents a major problem for soil conservation. This explains the occurrence throughout the Mediterranean mountains of many degraded hillslopes, which show evidence of sheet wash erosion, rilling, gullying and shallow landsliding; (ii) farmland abandonment has led to a marked reduction in runoff and sediment yield as a consequence of rapid plant recolonization, particularly by dense shrubs; (iii) the natural transformation of abandoned fields into grazing meadows has reduced runoff and sediment yield. Land use trends in the Mediterranean mountains are mainly characterized by generalized farmland abandonment and a decrease in livestock pressure. From a hydrological and geomorphological point of view the main consequences have been a reduction in overland flow from the hillslopes, and a reduction in sediment sources, with differences up to one order of magnitude in sediment yield from dense shrub cover and grazing meadow areas compared with areas under shifting agriculture. Copyright © 2012 John Wiley & Sons, Ltd.     

Strong ground motion record selection for the reliable prediction of the mean seismic collapse capacity of a structure group

How to select a limited number of strong ground motion records (SGMRs) is an important challenge for the seismic collapse capacity assessment of structures. The collapse capacity is considered as the ground motion intensity measure corresponding to the drift‐related dynamic instability in the structural system. The goal of this paper is to select, from a general set of SGMRs, a small number of subsets such that each can be used for the reliable prediction of the mean collapse capacity of a particular group of structures, i.e. of single degree‐of‐freedom systems with a typical behaviour range. In order to achieve this goal, multivariate statistical analysis is first applied, to determine what degree of similarity exists between each selected small subset and the general set of SGMRs. Principal Component analysis is applied to identify the best way to group structures, resulting in a minimum number of SGMRs in a proposed subset. The structures were classified into six groups, and for each group a subset of eight SGMRs has been proposed. The methodology has been validated by analysing a first‐mode‐dominated three‐storey‐reinforced concrete structure by means of the proposed subsets, as well as the general set of SGMRs. The results of this analysis show that the mean seismic collapse capacity can be predicted by the proposed subsets with less dispersion than by the recently developed improved approach, which is based on scaling the response spectra of the records to match the conditional mean spectrum. Copyright © 2010 John Wiley & Sons, Ltd.     

A simple model accounting for the uptake,transport, and deposition of wind‐eroded mineral particles in the hyperarid coastal Atacama Desert of northern Chile

As previously observed in marine sediments collected downwind of African or South American continental sources, recent studies of sediment cores collected at the bottom of Mejillones Bay in north Chile (23°S) show a laminated structure in which the amount of particles of aeolian origin and their size create significant differences between the layers. This suggests inter‐annual to inter‐decadal variations in the strength of the local southerly winds responsible for (1) the erosion of the adjacent hyperarid surface of the Mejillones Pampa, and (2) the subsequent transport of the eroded particles towards the bay. A simple model accounting for the vertical uptake, transport, and deposition of the particles initially set into motion by wind at the surface of the pampa is proposed. This model, which could be adapted to other locations, assumes that the initial rate of (vertical) uptake is proportional to the (horizontal) saltation flux quantified by means of White's equation, that particles are lifted to a height (H), increasing with the magnitude of turbulence, and that sedimentation progressively removes the coarsest particles from the air column as it moves towards the bay. In this model, the proportionality constant (A) linking the vertical flux of particles with the horizontal flux, and the injection height (H) control the magnitude and size distribution of the deposition flux in the bay. Their values are determined using the wind speed measured over the pampa and the size distribution of particles collected in sediment traps deployed in the bay as constraints. After calibration, the model is used to assess the sensitivity of the deposition flux to the wind intensity variations. The possibility of performing such quantitative studies is necessary for interpreting precisely the variability of the aeolian material in the sediment cores collected at the bottom of Mejillones Bay. Copyright © 2010 John Wiley & Sons, Ltd.     

Zircon M127 – A Homogeneous Reference Material for SIMS U–Pb Geochronology Combined with Hafnium,Oxygen and,Potentially, Lithium Isotope Analysis

In this article, we document a detailed analytical characterisation of zircon M127, a homogeneous 12.7 carat gemstone from Ratnapura, Sri Lanka. Zircon M127 has TIMS‐determined mean U–Pb radiogenic isotopic ratios of 0.084743 ± 0.000027 for ²⁰⁶Pb/²³⁸U and 0.67676 ± 0.00023 for ²⁰⁷Pb/²³⁵U (weighted means, 2s uncertainties). Its ²⁰⁶Pb/²³⁸U age of 524.36 ± 0.16 Ma (95% confidence uncertainty) is concordant within the uncertainties of decay constants. The δ¹⁸O value (determined by laser fluorination) is 8.26 ± 0.06‰ VSMOW (2s), and the mean ¹⁷⁶Hf/¹⁷⁷Hf ratio (determined by solution ICP‐MS) is 0.282396 ± 0.000004 (2s). The SIMS‐determined δ⁷Li value is ?0.6 ± 0.9‰ (2s), with a mean mass fraction of 1.0 ± 0.1 μg g^?1 Li (2s). Zircon M127 contains ~ 923 μg g^?1 U. The moderate degree of radiation damage corresponds well with the time‐integrated self‐irradiation dose of 1.82 × 10¹⁸ alpha events per gram. This observation, and the (U–Th)/He age of 426 ± 7 Ma (2s), which is typical of unheated Sri Lankan zircon, enable us to exclude any thermal treatment. Zircon M127 is proposed as a reference material for the determination of zircon U–Pb ages by means of SIMS in combination with hafnium and stable isotope (oxygen and potentially also lithium) determination.     

Top-down and bottom-up modelling to support low-carbon scenarios: climate policy implications

Bottom-up and top-down models are used to support climate policies, to identify the options required to meet GHG abatement targets and to evaluate their economic impact. Some studies have shown that the GHG mitigation options provided by economic top-down and technological bottom-up models tend to vary. One reason for this is that these models tend to use different baseline scenarios. The bottom-up TIMES_PT and the top-down computable general equilibrium GEM-E3_PT models are examined using a common baseline scenario to calibrate them, and the extend of their different mitigation options and its relevant to domestic policy making are assessed. Three low-carbon scenarios for Portugal until 2050 are generated, each with different GHG reduction targets. Both models suggest close mitigation options and locate the largest mitigation potential to energy supply. However, the models suggest different mitigation options for the end-use sectors: GEM-E3_PT focuses more on energy efficiency, while TIMES_PT relies on decrease carbon intensity due to a shift to electricity. Although a common baseline scenario cannot be ignored, the models’ inherent characteristics are the main factor for the different outcomes, thereby highlighting different mitigation options.

Policy relevance

The relevance of modelling tools used to support the design of domestic climate policies is assessed by evaluating the mitigation options suggested by a bottom-up and a top-down model. The different outcomes of each model are significant for climate policy design since each suggest different mitigation options like end-use energy efficiency and the promotion of low-carbon technologies. Policy makers should carefully select the modelling tool used to support their policies. The specific modelling structures of each model make them more appropriate to address certain policy questions than others. Using both modelling approaches for policy support can therefore bring added value and result in more robust climate policy design. Although the results are specific for Portugal, the insights provided by the analysis of both models can be extended to, and used in the climate policy decisions of, other countries.     

Effects of mesoscale sea-surface temperature fronts on the marine atmospheric boundary layer

A numerical modelling study is presented focusing on the effects of mesoscale sea-surface temperature (SST) variability on surface fluxes and the marine atmospheric boundary-layer structure. A basic scenario is examined having two regions of SST anomaly with alternating warm/cold or cold/warm water regions. Conditions upstream from the anomaly region have SST values equal to the ambient atmosphere temperature, creating an upstream neutrally stratified boundary layer. Downstream from the anomaly region the SST is also set to the ambient atmosphere value. When the warm anomaly is upstream from the cold anomaly, the downstream boundary layer exhibits a more complex structure because of convective forcing and mixed layer deepening upstream from the cold anomaly. An internal boundary layer forms over the cold anomaly in this case, generating two distinct layers over the downstream region. When the cold anomaly is upstream from the warm anomaly, mixing over the warm anomaly quickly destroys the shallow cold layer, yielding a more uniform downstream boundary-layer vertical structure compared with the warm-to- cold case. Analysis of the momentum budget indicates that turbulent momentum flux divergence dominates the velocity field tendency, with pressure forcing accounting for only about 20% of the changes in momentum. Parameterization of surface fluxes and boundary-layer structure at these scales would be very difficult because of their dependence on subgrid-scale SST spatial order. Simulations of similar flow over smaller scale fronts (<5 km) suggest that small-scale SST variability might be parameterized in mesoscale models by relating the effective heat flux to the strength of the SST variance.     

Effects of forestry on summertime low flows and physical fish habitat in snowmelt‐dominant headwater catchments of the Pacific Northwest

Periods of summertime low flows are often critical for fish. This study quantified the impacts of forest clear‐cutting on summertime low flows and fish habitat and how they evolved through time in two snowmelt‐dominant headwater catchments in the southern interior of British Columbia, Canada. A paired‐catchment analysis was applied to July–September water yield, the number of days each year with flow less than 10% of mean annual discharge, and daily streamflow for each calendar day. The postharvest time series were divided into treatment periods of approximately 6–10 years, which were analysed independently to evaluate how the effects of forestry changed through time. An instream flow assessment using a physical habitat simulation‐style approach was used to relate streamflow to the availability of physical habitat for resident rainbow trout. About two decades after the onset of logging and as the extent of logging increased to approximately 50% of the catchments, reductions in daily summertime low flows became more significant for the July–September yield (43%) and for the analysis by calendar day (11–68%). Reductions in summertime low flows were most pronounced in the catchment with the longest postharvest time series. On the basis of the temporal patterns of response, we hypothesize that the delayed reductions in late‐summer flow represent the combined effects of a persistent advance in snowmelt timing in combination with at least a partial recovery of transpiration and interception loss from the regenerating forests. These results indicate that asymptotic hydrological recovery as time progresses following logging is not suitable for understanding the impacts of forest harvesting on summertime low flows. Additionally, these reductions in streamflow corresponded to persistent decreases in modelled fish habitat availability that typically ranged from 20% to 50% during the summer low‐flow period in one of the catchments, suggesting that forest harvest may have substantial delayed effects on rearing salmonids in headwater streams.     

Phase equilibria modelling of retrograde amphibole and clinozoisite in mafic eclogite from the Tso Morari massif,northwest India: constraining the P–T–M(H2O) conditions of exhumation

Phase equilibria modelling of post‐peak metamorphic mineral assemblages in (ultra)high‐P mafic eclogite from the Tso Morari massif, Ladakh Himalaya, northwest India, has provided new insights into the potential behaviour and source of metamorphic fluid during exhumation, and constrained the P–T conditions of hydration. A series of P–M(H₂O) pseudosections constructed in the Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O (NCKFMASHTO) system show that a number of petrographically distinct hydration episodes occurred during exhumation from peak P–T conditions (~640 °C, 27–28 kbar), resulting in the formation of abundant compositionally zoned amphibole and minor clinozoisite poikiloblasts at the expense of a peak assemblage dominated by garnet and omphacite. Initial hydration is interpreted to have occurred as a result of the destabilization of talc following isothermal decompression to ~23 kbar, which led to the formation of barroisite–winchite amphibole core domains. An episode of fluid infiltration from an external source at ~19 kbar, with or without syn‐decompressional cooling to ~560 °C, resulted in further barroisitic–winchitic amphibole growth, followed by the formation of clinozoisite poikiloblasts. Continued buoyancy‐driven exhumation to the base of the lower crust is constrained to have taken place with no additional fluid input. A final hydration event is characterized by the formation of magnesiohornblende rims on the barroisite–winchite cores, with the former interpreted to have formed during later prograde overprinting in the middle crust associated with the final stages of exhumation. Notably, the vast majority of externally sourced H₂O, comprising just over half of the current bulk rock fluid content, was added during this later hydration event. In a middle crustal setting, this is interpreted as the result of devolatilization reactions occurring in migmatitic host orthogneiss and/or metasedimentary units, or following the crystallization of partial melt.     

Simulation of less‐mobile porosity dynamics in contrasting sediment water interface porous media

Considering heterogeneity in porous media pore size and connectivity is essential to predicting reactive solute transport across interfaces. However, exchange with less‐mobile porosity is rarely considered in surface water/groundwater recharge studies. Previous research indicates that a combination of pore‐fluid sampling and geoelectrical measurements can be used to quantify less‐mobile porosity exchange dynamics using the time‐varying relation between fluid and bulk electrical conductivity. For this study, we use macro‐scale (10 s of cm) advection–dispersion solute transport models linked with electrical conduction in COMSOL Multiphysics to explore less‐mobile porosity dynamics in two different types of observed sediment water interface porous media. Modeled sediment textures contrast from strongly layered streambed deposits to poorly sorted lakebed sands and cobbles. During simulated ionic tracer perturbations, a lag between fluid and bulk electrical conductivity, and the resultant hysteresis, is observed for all simulations indicating differential loading of pore spaces with tracer. Less‐mobile exchange parameters are determined graphically from these tracer time series data without the need for inverse numerical model simulation. In both sediment types, effective less‐mobile porosity exchange parameters are variable in response to changes in flow direction and fluid flux. These observed flow‐dependent effects directly impact local less‐mobile residence times and associated contact time for biogeochemical reaction. The simulations indicate that for the sediment textures explored here, less‐mobile porosity exchange is dominated by variable rates of advection through the domain, rather than diffusion of solute, for typical low‐to‐moderate rate (approximately 3–40 cm/day) hyporheic fluid fluxes. Overall, our model‐based results show that less‐mobile porosity may be expected in a range of natural hyporheic sediments and that changes in flowpath orientation and magnitude will impact less‐mobile exchange parameters. These temporal dynamics can be assessed with the geoelectrical experimental tracer method applied at laboratory and field scales.