Soil property variation mapping through data mining of soil category maps

Spatial information on soil properties is an important input to hydrological models. In current hydrological modelling practices, soil property information is often derived from soil category maps by the linking method in which a representative soil property value is linked to each soil polygon. Limited by the area‐class nature of soil category maps, the derived soil property variation is discontinuous and less detailed than high resolution digital terrain or remote sensing data. This research proposed dmSoil, a data‐mining‐based approach to derive continuous and spatially detailed soil property information from soil category maps. First, the soil–environment relationships are extracted through data mining of a soil map. The similarity of the soil at each location to different soil types in the soil map is then estimated using the mined relationships. Prediction of soil property values at each location is made by combining the similarities of the soil at that location to different soil types and the representative soil property values of these soil types. The new approach was applied in the Raffelson Watershed and Pleasant Valley in the Driftless Area of Wisconsin, United States to map soil A horizon texture (in both areas) and depth to soil C horizon (in Pleasant Valley). The property maps from the dmSoil approach capture the spatial gradation and details of soil properties better than those from the linking method. The new approach also shows consistent accuracy improvement at validation points. In addition to the improved performances, the inputs for the dmSoil approach are easy to prepare, and the approach itself is simple to deploy. It provides an effective way to derive better soil property information from soil category maps for hydrological modelling. Copyright © 2014 John Wiley & Sons, Ltd.     

Pre‐Alpine discordant granitic dikes in the metamorphic core of the Betic Cordillera: tectonic implications

Pegmatite dikes bearing andalusite crosscut foliation S₂ in Alpujarride gneisses and schists. Post‐S₂ andalusite is transposed by a foliation S₃, defined by fibrolite, which affects the dikes. The dikes represent highly differentiated granitic magmas with low REE and Zr contents and a positive Eu anomaly. U‐Pb SHRIMP dating of magmatic zircons provided Pan‐African ages (cores) and late Variscan ages (rims). However, U‐rich rims also provided metamorphic Alpine ages, supporting a polyorogenic tectonometamorphic history for pre‐Mesozoic Alpujarride rocks.     

The multi-component Hekla Ö Tephra,Iceland: a complex widespread mid-Holocene tephra layer

Large Plinian eruptions from Hekla volcano, Iceland, produce compositionally zoned tephra used as key markers in tephrochronology. However, spatial variations in chemical composition of a tephra layer may complicate its identification. An example is the 5950–6180 cal a bp Hekla Ö tephra layer, which shows compositional spread from rhyolite, dacite and andesite to basalt. In soil sections north of Hekla, the SiO₂ content of the tephra glass reaches 76 wt% in the lowest unit of the Hekla Ö deposit and decreases to 62–63 wt% in the uppermost unit. Intermingled within the whole deposit are basalt tephra grains having 46–47 wt% SiO₂. The composition of the basalt glass includes primitive basalt and a more evolved basalt (MgO >6 and <6 wt%, respectively). Together with literature data, the Hekla Ö tephra and the so-called T-Tephra/Hekla-T are most likely from contemporaneous eruptions of different vents on the Hekla volcanic system, forming a single important marker tephra (Hekla ÖT) deposited over 80% of Iceland. Identification is complicated by its spatial compositional heterogeneity, such as systematic decrease in SiO₂ content from the east to the west of Hekla volcano. Consequently, an individual tephra layer from a large explosive eruption can have different composition at different locations. © 2020 John Wiley & Sons, Ltd.     

Does stress transmission in forelands depend on structural style? Distinctive stress magnitudes during Sevier thin‐skinned and Laramide thick‐skinned layer‐parallel shortening in the Bighorn Basin (USA) revealed by stylolite and calcite twinning paleopiezometry

The Sheep Mountain‐Little Sheep Mountain Anticlines, Bighorn Basin (USA) formed as basement‐cored Laramide structures in the formerly undeformed foreland of the thin‐skinned Sevier orogen. We take advantage of the well‐constrained microstructural network there to reconstruct differential stress magnitudes that prevailed during both Sevier and Laramide layer‐parallel shortening (LPS), before the onset of large‐scale folding. Differential stress magnitudes determined from tectonic stylolites are compared and combined to previous stress estimates from calcite twinning paleopiezometry in the same formations. During stress loading related to LPS, differential stress magnitudes transmitted from the distant Sevier thin‐skinned orogen into the sedimentary cover of the Bighorn basin (11–43 MPa) are higher than the differential stress magnitudes accompanying the early stage of LPS related to the thick‐skinned Laramide deformation (2–19 MPa). This study illustrates that the tectonic style of an orogen affects the transmission of early orogenic stress into the stable continental interior.     

Crustal segments in the North Patagonian Massif,Patagonia: An integrated perspective based on Sm–Nd isotope systematics

New insights on the Paleozoic evolution of the continental crust in the North Patagonian Massif are presented based on the analysis of Sm–Nd systematics. New evidence is presented to constrain tectonic models for the origin of Patagonia and its relations with the South American crustal blocks. Geologic, isotopic and tectonic characterization of the North Patagonian Massif and comparison of the Nd parameters lead us to conclude that: (1) The North Patagonian Massif is a crustal block with bulk crustal average ages between 2.1 and 1.6 Ga T_DM (Nd) and (2) At least three metamorphic episodes could be identified in the Paleozoic rocks of the North Patagonian Massif. In the northeastern corner, Famatinian metamorphism is widely identified. However field and petrographic evidence indicate a Middle to Late Cambrian metamorphism pre-dating the emplacement of the ca. 475 Ma granitoids. In the southwestern area, are apparent 425–420 Ma (?) and 380–360 Ma metamorphic peaks. The latter episode might have resulted from the collision of the Antonia terrane; and (3) Early Paleozoic magmatism in the northeastern area is coeval with the Famatinian arc. Nd isotopic compositions reveal that Ordovician magmatism was associated with attenuated crust. On the southwestern border, the first magmatic recycling record is Devonian. Nd data shows a step by step melting of different levels of the continental crust in the Late Palaeozoic. Between 330 and 295 Ma magmatism was likely the product of a crustal source with an average 1.5 Ga T_DM (Nd). Widespread magmatism represented by the 295–260 Ma granitoids involved a lower crustal mafic source, and continued with massive shallower-acid plutono volcanic complexes which might have recycled an upper crustal segment of the Proterozoic continental basement, resulting in a more felsic crust until the Triassic. (4) Sm–Nd parameters and detrital zircon age patterns of Early Paleozoic (meta)-sedimentary rocks from the North Patagonian Massif and those from the neighboring blocks, suggest crustal continuity between Eastern Sierras Pampeanas, southern Arequipa-Antofalla and the northeastern sector of the North Patagonian Massif by the Early Paleozoic. This evidence suggests that, at least, this corner of the North Patagonian Massif is not allochthonous to Gondwana. A Late Paleozoic frontal collision with the southwestern margin of Gondwana can be reconcilied in a para-autochthonous model including a rifting event from a similar or neighbouring position to its post-collision location. Possible Proterozoic or Early Paleozoic connections of the NPM with the Kalahari craton or the western Antartic blocks should be investigated.     

Polygonal grain‐based distinct element modeling for mechanical behavior of brittle rock

The microstructure of rock was numerically reproduced by a polygonal grain‐based model, and its mechanical behavior was examined by performing the uniaxial compression test and Brazilian tests via the Universal Distinct Element Code. The numerical results of the model demonstrated good agreement with the experimental results obtained with rock specimens in terms of the stress–strain behavior, strength characteristics, and brittle fracture phenomenon. An encouraging result is that the grain‐based model‐Universal Distinct Element Code model can reproduce a low ratio of tensile to compressive strength of 1/20 to 1/10 without the need for an additional process. This finding is ascribed to the fact that the geometrical features of polygons can effectively capture the effects of angularity, finite rotation, and interlocking of grains that exist in reality. A numerical methodology to monitor the evolution of micro‐cracks was developed, which enabled us to examine the progressive process of the failure and distinguish the contribution of tensile cracking to the process from that of shear cracking. From the observations of the micro‐cracking process in reference to the stress–strain relation, crack initiation stress, and crack damage stress, it can be concluded that the failure process of the model closely resembles the microscopic observations of rock. We also carried out a parametric study to examine the relationships between the microscopic properties and the macroscopic behavior of the model. Depending on the micro‐properties, the model exhibited a variety of responses to the external load in terms of the strength and deformation characteristics, the evolution of micro‐cracks, and the post‐peak behavior. Copyright © 2016 John Wiley & Sons, Ltd.     

The influence of coupled physical swelling and chemical reactions on deformable geomaterials

Coupled thermo‐hydro‐mechanical‐chemical modelling has attracted attention in past decades due to many contemporary geotechnical engineering applications (e.g., waste disposal, carbon capture and storage). However, molecular‐scale interactions within geomaterials (e.g., swelling and dissolution/precipitation) have a significant influence on the mechanical behaviour, yet are rarely incorporated into existing Thermal‐Hydro‐Mechanical‐Chemical (THMC) frameworks. This paper presents a new coupled hydro‐mechanical‐chemical constitutive model to bridge molecular‐scale interactions with macro‐physical deformation by combining the swelling and dissolution/precipitation through an extension of the new mixture‐coupling theory. Entropy analysis of the geomaterial system provides dissipation energy, and Helmholtz free energy gives the relationship between solids and fluids. Numerical simulation is used to compare with the selected recognized models, which demonstrates that the swelling and dissolution/precipitation processes may have a significant influence on the mechanical deformation of the geomaterials.     

Including spatial distribution in a data‐driven rainfall‐runoff model to improve reservoir inflow forecasting in Taiwan

Multi‐step ahead inflow forecasting has a critical role to play in reservoir operation and management in Taiwan during typhoons as statutory legislation requires a minimum of 3‐h warning to be issued before any reservoir releases are made. However, the complex spatial and temporal heterogeneity of typhoon rainfall, coupled with a remote and mountainous physiographic context, makes the development of real‐time rainfall‐runoff models that can accurately predict reservoir inflow several hours ahead of time challenging. Consequently, there is an urgent, operational requirement for models that can enhance reservoir inflow prediction at forecast horizons of more than 3 h. In this paper, we develop a novel semi‐distributed, data‐driven, rainfall‐runoff model for the Shihmen catchment, north Taiwan. A suite of Adaptive Network‐based Fuzzy Inference System solutions is created using various combinations of autoregressive, spatially lumped radar and point‐based rain gauge predictors. Different levels of spatially aggregated radar‐derived rainfall data are used to generate 4, 8 and 12 sub‐catchment input drivers. In general, the semi‐distributed radar rainfall models outperform their less complex counterparts in predictions of reservoir inflow at lead times greater than 3 h. Performance is found to be optimal when spatial aggregation is restricted to four sub‐catchments, with up to 30% improvements in the performance over lumped and point‐based models being evident at 5‐h lead times. The potential benefits of applying semi‐distributed, data‐driven models in reservoir inflow modelling specifically, and hydrological modelling more generally, are thus demonstrated. Copyright © 2012 John Wiley & Sons, Ltd.