Stable staurolite–cordierite assemblages in K‐poor metapelitic schists in Aston and Hospitalet gneiss domes of the central Pyrenees (France,Andorra)

Staurolite–cordierite assemblages are common in mica schists of the Aston and Hospitalet gneiss domes of the central Axial Zone, Pyrenees (France, Andorra). Within a 200 m wide zone, staurolite, cordierite and andalusite porphyroblasts contain inclusion trails that preserve the same stage of development of a crenulation cleavage, strongly suggesting that all three phases are contemporaneous. Their syntectonic growth occurred during a short period at the beginning of the formation of the dominant schistosity (S₂) of the domes. Staurolite and cordierite touching each other further indicates an equilibrium relationship. Whole‐rock analyses show that some staurolite–cordierite schists are depleted in K₂O compared to post‐Archean shales (PAAS) and amphibolite facies pelites. Analysis of the st‐crd paragenesis in K‐poor schists without muscovite using KFMASH and MnNCKFMASH petrogentic grids, pseudosections and AFM compatibility diagrams predicts stable conditions at pressures of ~3.5 kbar at 575 °C. For metapelites with intermediate X_Mg values (0.7 >  X_Mg >0.48) a ‘muscovite‐out window’ exists from 550–650 °C at 3.5 kbar in the KFMASH system. Conventional thermobarometry (GB‐GASP, AvT‐AvP) and petrogenetic grids show an isobaric P–T path to peak temperatures of ~650 °C, supported by the presence of sillimanite‐K‐feldspar gneiss and migmatites. LP‐HT metamorphism in the Aston dome is related to early Carboniferous (c. 339 Ma) granitic intrusions into the dome core. As metamorphism is directly linked with the formation of the main S₂ schistosity, the temporal relations demonstrated in this study conflict with previous studies which constrained LP‐HT metamorphism and the development of flat‐lying schistosity to the late Carboniferous (315–305 Ma) – at least in the eastern Axial Zone.     

Influence of anthropogenic land‐use change on hillslope erosion in the Waipaoa River Basin,New Zealand

European settlement of the Poverty Bay Region resulted in deforestation and conversion of > 90% of the landscape to pastureland. The resulting loss of vegetation triggered a rapid increase in hillslope erosion as widespread landslide complexes and gully systems developed on weak lithologic units in the Waipaoa Basin. To quantify the rate and volume of historic hillslope degradation, we used a 1956–2010 sequence of aerial photographs for a ~16 km² catchment to map temporal changes in the spatial extent of active landslides. Then we created a ‘turf index’ based on the extent and style of pastoral ground disruption, which correlates with downslope velocity. Based on the movement of trees and other features, we assigned average velocities to the turf classes as follows: (1) minimal disrupted ground: 0.6 m/yr, (2) a mix of disrupted ground and intact blocks: 3.4 m/yr, and (3) no intact blocks or vegetation: > 6 m/yr. We then calculated the average annual sediment flux using these turf‐derived velocities, the width of the landslide‐channel intersection, and an average toe depth of 4.4 ± 1.3 m (mean ± standard deviation [SD]) from 37 field measurements. The resulting catchment averaged erosion rates are (mean ± SD): 29.9 ± 12.9 mm/yr (1956), 28.8 ± 13.7 mm/yr (1969), 13.4 ± 4.9 mm/yr (1979), 17.0 ± 6.2 mm/yr (1988), and 9.9 ± 3.6 mm/yr (2010). Compared with long‐term (post‐18 ka) erosion rates (1.6 mm/yr) and the long‐term uplift rate (~1 mm/yr) for this site, the 50‐year anthropogenically‐driven rate is an order of magnitude larger (~20 mm/yr). Previously, we measured an increase in erosion over the past 3.4 kyr (2.2 mm/yr), and here, we demonstrate this increase could be primarily due to human land‐use change – showing that a century of rapid erosion superimposed on the background geologic rate can profoundly skew the interpretation of erosion rates. Copyright © 2016 John Wiley & Sons, Ltd.     

Application of integrated back‐propagation network and self organizing map for flood forecasting

Much of the nonlinearity and uncertainty regarding the flood process is because hydrologic data required for estimation are often tremendously difficult to obtain. This study employed a back‐propagation network (BPN) as the main structure in flood forecasting to learn and to demonstrate the sophisticated nonlinear mapping relationship. However, a deterministic BPN model implies high uncertainty and poor consistency for verification work even when the learning performance is satisfactory for flood forecasting. Therefore, a novel procedure was proposed in this investigation which integrates linear transfer function (LTF) and self‐organizing map (SOM) to efficiently determine the intervals of weights and biases of a flood forecasting neural network to avoid the above problems. A SOM network with classification ability was applied to the solutions and parameters of the BPN model in the learning stage, to classify the network parameter rules and to obtain the winning parameters. The outcomes from the previous stage were then used as the ranges of the parameters in the recall stage. Finally, a case study was carried out in Wu‐Shi basin to demonstrate the effectiveness of the proposal. Copyright © 2009 John Wiley & Sons, Ltd.     

Granular element method for three‐dimensional discrete element calculations

This paper endows the recently‐proposed granular element method (GEM) with the ability to perform 3D discrete element calculations. By using non‐uniform rational B‐Splines to accurately represent complex grain geometries, we proposed an alternative approach to clustering‐based and polyhedra‐based discrete element methods whereby the need for complicated and ad hoc approaches to construct 3D grain geometries is entirely bypassed. We demonstrate the ability of GEM in capturing arbitrary‐shaped 3D grains with great ease, flexibility, and without excessive geometric information. Furthermore, the applicability of GEM is enhanced by its tight integration with existing non‐uniform rational B‐Splines modeling tools and ability to provide a seamless transition from binary images of real grain shapes (e.g., from 3D X‐ray CT) to modeling and discrete mechanics computations.Copyright © 2013 John Wiley & Sons, Ltd.     

Main features of the deep structure by local earthquake tomography and active tectonics: case of Rif Mountain (morocco) and Betic Cordillera (Spain)

Within the Spain and Moroccan networks, a large volume of seismic data has been collected and used for investigating the lithosphere in the Betic–Rif Cordillera. The present study has two main goals: (1) Use the most actual seismological data from recent earthquakes in the Betic–Rif arc for investigating the lithosphere through the application of seismic local tomography techniques. (2) Define the possible structural blocks and explain the GPS velocities perturbation in this region. The resolution tests results indicate that the calculated images gave a close true structure for the studied regions from 5- to 60-km depth. The resulting tomographic image shows that the presence of two upper crust body (velocity 6.5 km/s) at 3- to 13-km depth between Iberian Betic and Moroccan Rif in the western and in the middle of Alboran Sea also shows the low velocity favoring the presence of melt in the base of these two bodies. The crustal bodies forms tectonic blocks in the Central Rif and in the Central Betic Cordillera.     

Growth fold controls on carbonate distribution in mixed foreland basins: insights from the Amiran foreland basin (NW Zagros,Iran) and stratigraphic numerical modelling

The evolution from Late Cretaceous to early Eocene of the well dated Amiran foreland basin in the NW Iranian Zagros Mountains is studied based on the reconstruction of successive thickness, palaeobathymetry and subsidence maps. These maps show the progressive forelandwards migration of the mixed carbonate‐siliciclastic system associated with a decrease in creation of accommodation. Carbonate facies variations across the basin suggest a structural control on the carbonate distribution in the Amiran foreland basin, which has been used as initial constraint to study the control exerted by syndepositional folding in basin architecture and evolution by means of stratigraphic numerical modelling. Modelled results show that shallow bathymetries on top of growing folds enhance carbonate production and basin compartmentalization. As a consequence, coarse clastics become restricted to the internal parts of the basin and only the fine sediments can by‐pass the bathymetric highs generated by folding. Additionally, the development of extensive carbonate platforms on top of the anticlines favours the basinwards migration of the depositional system, which progrades farther with higher fold uplift rates. In this scenario, build‐ups on top of anticlines record its growth and can be used as a dating method. Extrapolation of presented modelling results into the Amiran foreland basin is in agreement with an early folding stage in the SE Lurestan area, between the Khorramabad and Kabir Kuh anticlines. This folding stage would enhance the development of carbonate platforms on top of the anticlines, the south‐westward migration of the system and eventually, the complete filling of the basin north of the Chenareh anticline at the end of the Cuisian. Incremental thickness maps are consistent with a thin (0.4–2 km) ophiolite complex in the source area of the Amiran basin.     

A Compilation of Silicon,Rare Earth Element and Twenty‐One other Trace Element Concentrations in the Natural River Water Reference Material SLRS‐5 (NRC‐CNRC)

The natural river water certified reference material SLRS‐5 (NRC‐CNRC) was routinely analysed in this study for major and trace elements by ten French laboratories. Most of the measurements were made using ICP‐MS. Because no certified values are assigned by NRC‐CNRC for silicon and 35 trace element concentrations (rare earth elements, Ag, B, Bi, Cs, Ga, Ge, Li, Nb, P, Rb, Rh, Re, S, Sc, Sn, Th, Ti, Tl, W, Y and Zr), or for isotopic ratios, we provide a compilation of the concentrations and related uncertainties obtained by the participating laboratories. Strontium isotopic ratios are also given.     

Sorption and Chemical Reactions of PAHs with Dissolved Humic Substances and Related Model Polymers

Sorption coefficients measured for PAHs on dissolved humic substances by SPME and FQT were found to be inevitably different and method‐dependent – SPME provides activity‐based and FQT concentration‐based sorption coefficients. Poly(acrylic acid) esters as well‐defined model polymers were used in sorption experiments, leading to the conclusion that short aliphatic chains are more effective in binding PAHs than aromatic moieties. FQT was inappropriate to measure sorption coefficients for the interaction of pyrene with poly(acrylic acid) esters but the experiments revealed a characteristic shift in the fluorescence spectrum. Using pyrene as a probe for the molecular environment in the sorbed state, the observed spectral shift indicated a highly hydrophobic microenvironment. The empirical relationships between lg K_DOC and lg K_OW were generalized on the basis of a modified Flory‐Huggins concept. Introducing only one sorbent‐specific parameter, the solubility parameter δ_DOM, the calculation of sorption coefficients became possible for a wide range of HOCs using fundamental data readily available from the literature. Long‐term experiments showed that reactive PAHs (such as acenaphthylene and 9‐methylanthracene) are able to react with HAs under strictly abiotic and anoxic conditions, whereas less reactive PAHs (such as naphthalene and dihydroanthracene) do not form bound residues. The HA reveals two functions in the interaction, behaving as a reaction partner and as a protecting ligand.     

Frequency‐domain acoustic‐elastic coupled waveform inversion using the Gauss‐Newton conjugate gradient method

We developed a frequency‐domain acoustic‐elastic coupled waveform inversion based on the Gauss‐Newton conjugate gradient method. Despite the use of a high‐performance computer system and a state‐of‐the‐art parallel computation algorithm, it remained computationally prohibitive to calculate the approximate Hessian explicitly for a large‐scale inverse problem. Therefore, we adopted the conjugate gradient least‐squares algorithm, which is frequently used for geophysical inverse problems, to implement the Gauss‐Newton method so that the approximate Hessian is calculated implicitly. Thus, there was no need to store the Hessian matrix. By simultaneously back‐propagating multi‐components consisting of the pressure and displacements, we could efficiently extract information on the subsurface structures. To verify our algorithm, we applied it to synthetic data sets generated from the Marmousi‐2 model and the modified SEG/EAGE salt model. We also extended our algorithm to the ocean‐bottom cable environment and verified it using ocean‐bottom cable data generated from the Marmousi‐2 model. With the assumption of a hard seafloor, we recovered both the P‐wave velocity of complicated subsurface structures as well as the S‐wave velocity. Although the inversion of the S‐wave velocity is not feasible for the high Poisson's ratios used to simulate a soft seafloor, several strategies exist to treat this problem. Our example using multi‐component data showed some promise in mitigating the soft seafloor effect. However, this issue still remains open.     

Application of singular spectrum analysis to identify the degrading structure using deteriorating distributed element model

Reinforced concrete structure may exhibit significant inelastic hysteretic behavior when subject to strong earthquake excitation. To investigate such an inelastic behavior, in this study, a new system identification technique is applied by using the deteriorating distributed element (DDE) model to simulate the hysteretic behavior of a degrading structure. Through the advanced signal processing technique, the multiple singular spectrum analysis (SSA) and the nonlinear SSA, the recorded inelastic restoring force of a deteriorating structure can be decomposed into several independent additive components in its sequentially degrading order and with decreasing weight. With each decomposed hysteresis loop, the model parameters of the DDE model are identified. The evolutionary properties of the progressive stiffness degradation behavior of reinforced concrete structure can be observed from the identified model parameters. Finally, comparison on the physical properties of the identified DDE model with respect to the seismic response data of the deteriorating structure is also discussed. The result shows that the proposed identification technique can have a good estimation on the seismic behavior of the degrading structure. Copyright © 2012 John Wiley & Sons, Ltd.