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.     

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.     

Removal of Pyridine from Aqueous Solution by Adsorption on an Activated Carbon Cloth

Pyridine is a very toxic pollutant that has to be removed from wastewater. In this work, adsorption of pyridine on activated carbon cloth (ACC) is studied as a possible alternative for eliminating pyridine from aqueous solution. The ACC was produced from polyacrylonitrile. The adsorption equilibrium data of pyridine on ACC was obtained in a batch adsorber. The experimental data was interpreted with the isotherms of Langmuir, Freundlich, and Prausnitz‐Radke (PR), and the PR isotherm better represented the experimental data. The capacity of ACC for adsorbing pyridine was favored increasing the solution pH from 3 to 6, and this effect was due to the π–π dispersive and electrostatic interactions between the pyridine species in solution and the surface complexes of ACC. The modified Langmuir model fitted reasonably well the influence of pH on the adsorption capacity. In this model was assumed that both neutral pyridine and pyridinium were simultaneously adsorbed on ACC accordingly to the experimental results. The adsorption capacity was almost independent of temperature. The reversibility study revealed that 75% of the pyridine can be desorbed from ACC indicating that part of the pyridine was irreversibly adsorbed, and possibly chemisorbed.     

New P–T constraints on the Tamayen glaucophane‐bearing rocks,eastern Taiwan: Perple_X modelling results and geodynamic implications

New pseudosection modelling was applied to better constrain the P–T conditions and evolution of glaucophane‐bearing rocks in the Tamayen block of the Yuli belt, recognized as the world's youngest known blueschist complex. Based on the predominant clinoamphibole, textural relationships, and mineral compositions, these glaucophane‐bearing high‐P rocks can be divided into four types. We focused on the three containing garnet. The chief phase assemblages are (in decreasing mode): amphibole + quartz + epidote + garnet + chlorite + rutile/titanite (Type‐I), phengite + amphibole + quartz + garnet + chlorite + epidote + titanite + biotite + magnetite (Type‐II), and amphibole + quartz + albite + epidote + garnet + rutile + hematite + titanite (Type‐III). Amphibole exhibits compositional zoning from core to rim as follows: glaucophane → pargasitic amphibole → actinolite (Type‐I), barroisite → Mg‐katophorite/taramite → Fe‐glaucophane (Type‐II), glaucophane → winchite (Type‐III). Using petrographic data, mineral compositions and Perple_X modelling (pseudosections and superimposed isopleths), peak P–T conditions were determined as 13 ± 1 kbar and 550 ± 40 °C for Type‐I, 10.5 ± 0.5 kbar and 560 ± 30 °C for Type‐II (thermal peak) and 11 ± 1 kbar and 530 ± 30 °C for Type‐III. The calculations yield higher pressures and temperatures than previously thought; the difference is ~1–6 kbar and 50–200 °C. The three rock types record similar P–T retrograde paths with clockwise trajectories; all rocks followed trajectories with substantial pressure decrease under near‐isothermal conditions (Type‐I and Type‐III), with the probable exception of Type‐II where decompression followed colder geotherms. The P–T paths suggest a tectonic environment in which the rocks were exhumed from maximum depths of ~45 km within a subduction channel along a relative cold geothermal gradient of ~11–14 °C km^?1.     

The effect of lithology on valley width,terrace distribution,and bedload provenance in a tectonically stable catchment with flat‐lying stratigraphy

How rock resistance or erodibility affects fluvial landforms and processes is an outstanding question in geomorphology that has recently garnered attention owing to the recognition that the erosion rates of bedrock channels largely set the pace of landscape evolution. In this work, we evaluate valley width, terrace distribution, and bedload provenance in terms of reach scale variation in lithology in the study reach and discuss the implications for landscape evolution in a catchment with relatively flat‐lying stratigraphy and very little uplift. A reach of the Buffalo National River in Arkansas was partitioned into lithologic reaches and the mechanical and chemical resistance of the main lithologies making up the catchment was measured. Valley width and the spatial distribution of terraces were compared among the different lithologic reaches. The surface grain size and provenance of coarse (2–90 mm) sediment of both modern gravel bars and older terrace deposits that make up the former bedload were measured and defined. The results demonstrate a strong impact of lithology upon valley width, terrace distribution, and bedload provenance and therefore, upon landscape evolution processes. Channel down‐cutting through different lithologies creates variable patterns of resistance across catchments and continents. Particularly in post‐tectonic and non‐tectonic landscapes, the variation in resistance that arises from the exhumation of different rocks in channel longitudinal profiles can impact local base levels, initiating responses that can be propagated through channel networks. The rate at which that response is transmitted through channels is potentially amplified and/or mitigated by differences between the resistance of channel beds and bedload sediment loads. In the study reach, variation in lithologic resistance influences the prevalence of lateral and vertical processes, thus producing a spatial pattern of terraces that reflects rock type rather than climate, regional base level change, or hydrologic variability. Copyright © 2017 John Wiley & Sons, Ltd.     

A novel approach for estimating karst groundwater recharge in mountainous regions and its application in Switzerland

A pragmatic and simple approach for estimating the groundwater recharge of karst aquifers in mountainous regions by extrapolation of the hydrological regimes of gauged and well‐documented systems is presented. Specific discharge rates are derived using annual precipitation and spring measurements by taking into account catchment size and elevation, which are assumed to be the dominant factors. Reference sites with high data reliability are used for calibration and regional extrapolation. This is performed with normalized values employing spatial precipitation deviations and correlation with the elevation of the catchment areas. A tiered step procedure provides minimum and maximum normalized gradients for the relationship between recharge quantity and elevation for karst regions. The normalized recharge can therefore be obtained and extrapolated for any location using the spatial precipitation variability to provide an estimate of annual groundwater recharge. The approach was applied to Switzerland (approximately 7500 km² of karst terrain situated between 200 and over 4000 m a.s.l.) using annual precipitation data from meteorological stations for the years 2000 to 2011. Results show that the average recharge rates of different Swiss karst domains range from 20 to 46 L/km²s, which corresponds to an infiltration ratio between 0.6 and 0.9 of total precipitation. Despite uncertainties inherent in the approach, these results provide a benchmark for renewable karst groundwater resources in Switzerland of about 8.4 km³/year. The approach can be applied to any other mountainous karst region, that is, where a clear relationship between elevation, precipitation and recharge can be assumed. Copyright © 2015 John Wiley & Sons, Ltd.     

Vegetation in western Central Mexico during the last 50 000 years: modern analogs and climate in the Zacapu Basin

The late glacial and the transition towards the Holocene marked a period of significant environmental change at a global scale. In western Central Mexico, few records span beyond the Holocene, and little is known about Pleistocene climatic and environmental variability. Here we report on the pollen record of a composite sequence made of three cores that cover the last 50 ka (thousands of calibrated radiocarbon years before present) in the Zacapu Basin, western Central Mexico (～1970 m asl). The conjunction of modern pollen rain from the area with the composite fossil pollen sequence from cores Cantabria, Cantabria 1 and CEMCA Point‐1 provided a detailed history of the vegetation of the area. Recognition of modern environmental patterns through individual modern pollen taxa proved difficult, but multivariate analysis separated different vegetation types. No‐modern‐analog communities relative to the modern vegetation cover of the area were evident, especially during periods of drought. Apparently, summer precipitation was significant before 40 ka and after the deglaciation. However, this moisture source declined from ～40 to 10 ka because tradewinds from the Gulf of Mexico were blocked by the eastern mountains of Central Mexico as a result of lower temperatures. By contrast, westerly and northerly winds were enhanced by maximum ice cover from ～19 to 17 ka, providing abundant winter precipitation and therefore preventing extreme drought in the area. Copyright © 2012 John Wiley & Sons, Ltd.