The critical state behaviour of barodesy compared with the Matsuoka–Nakai failure criterion

Barodesy is a new approach to constitutive modelling of soil. It is based on Goldscheider's principles and maps stretching directions onto corresponding stress directions with the help of a simple exponential function. This mapping also determines a critical state surface in principal stress space. The article investigates this surface and relates it to the well‐known Matsuoka–Nakai failure criterion. It turns out that the difference between these two surfaces is negligible for practical applications. Copyright © 2011 John Wiley & Sons, Ltd.     

Chloride transfer in cement‐based materials. Part 1. Theoretical basis and modelling

In this first part of the work, we develop macroscopic models for migration and diffusion–migration of ionic species in saturated porous media, based on periodic homogenization. The prior application is chloride transport in cementitious materials. The dimensional analysis of Nernst–Planck equation lets appear to dimensionless numbers characterizing the ionic transfer in the porous medium. Using experimental data obtained from electrodiffusion tests on cement‐based materials (Part II), these dimensionless numbers are linked to the perturbation parameter ?. For a strong imposed electrical field, the asymptotic expansion of Nernst–Planck equation leads to a macroscopic model where the migration is predominant. For a weak imposed electrical field or in natural diffusion, we obtain a macroscopic model coupling diffusion and migration at the same order. In both models, the expression of the homogenized diffusion tensor is identical and only involves the geometric properties of the material microstructure. Copyright © 2012 John Wiley & Sons, Ltd.     

Phase equilibria constraints on melting of stromatic migmatites from Ronda (S. Spain): insights on the formation of peritectic garnet

Stromatic metatexites occurring structurally below the contact with the Ronda peridotite (Ojén nappe, Betic Cordillera, S Spain) are characterized by the mineral assemblage Qtz+Pl+Kfs+Bt+Sil+Grt+Ap+Gr+Ilm. Garnet occurs in low modal amount (2–5 vol.%). Very rare muscovite is present as armoured inclusions, indicating prograde exhaustion. Microstructural evidence of melting in the migmatites includes pseudomorphs after melt films and nanogranite and glassy inclusions hosted in garnet cores. The latter microstructure demonstrates that garnet crystallized in the presence of melt. Re‐melted nanogranites and preserved glassy inclusions show leucogranitic compositions. Phase equilibria modelling of the stromatic migmatite in the MnO–Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂–O₂–C (MnNCaKFMASHOC) system with graphite‐saturated fluid shows P–T conditions of equilibration of 4.5–5 kbar, 660–700 °C. These results are consistent with the complete experimental re‐melting of nanogranites at 700 °C and indicate that nanogranites represent the anatectic melt generated immediately after entering supersolidus conditions. The P–T estimate for garnet and melt development does not, however, overlap with the low‐temperature tip of the pure melt field in the phase diagram calculated for the composition of preserved glassy inclusions in garnet in the Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O (NCKFMASH) system. A comparison of measured melt compositions formed immediately beyond the solidus with results of phase equilibria modelling points to the systematic underestimation of FeO, MgO and CaO in the calculated melt. These discrepancies are present also when calculated melts are compared with low‐T natural and experimental melts from the literature. Under such conditions, the available melt model does not perform well. Given the presence of melt inclusions in garnet cores and the P–T estimates for their formation, we argue that small amounts (<5 vol.%) of peritectic garnet may grow at low temperatures (≤700 °C), as a result of continuous melting reactions consuming biotite.     

Diachronous Palaeoproterozoic deformation and metamorphism in the Committee Bay belt,Rae Province,Nunavut: insights from 40Ar–39Ar cooling ages and thermal modelling

Regional‐scale ⁴⁰Ar–³⁹Ar data presented in this paper reveal significant across‐strike and along‐strike age differences in the Committee Bay belt (CBb), Rae Province, Nunavut, Canada, that complement variations in observed monazite ages. ⁴⁰Ar–³⁹Ar hornblende ages are c. 1795, 1775, and 1750 Ma in the western, eastern and central parts of the Prince Albert Group (PAG) domain respectively. The migmatite domain and Walker Lake intrusive complex are characterized by c. 1750–1730 ⁴⁰Ar–³⁹Ar hornblende ages without significant along‐strike variation. The ⁴⁰Ar–³⁹Ar data provide important constraints on the cooling history and on thermal modelling that elucidates the controls on diachroneity and metamorphic patterns within the belt. In the western CBb, prograde monazite growth occurred 26 ± 10 Myr earlier in the migmatite domain (1864 ± 9 Ma; peak P–T = 5 kbar?700 °C) than in the PAG domain (1838 ± 5 Ma; peak P–T = 5 kbar?580 °C). Calculations indicate that this earlier monazite growth results from tectonic thickening of higher heat productivity Archean lithologies in the migmatite domain, which undergoes more rapid prograde heating than the less radiogenetic and lower grade rocks of the PAG domain. Granite generation via biotite dehydration melting at 800 °C and 20 km depth is predicted to occur c. 1835 Ma, in agreement with geochronological constraints. The tectonic burial of crustal domains with contrasting radiogenic properties also explains the general congruence of lower to upper amphibolite facies metamorphic zones generated during the two main orogenic cycles (i.e. M₂–D₁ and M₃–D₂). The modelled timing of prograde monazite growth in the migmatite domain suggests that D₂ tectonic thickening began at 1872 ± 9 Ma, some 8 ± 3 Myr before monzazite growth, coeval with the inferred time of collision of the Meta Incognita terrane with the southern Rae Province. Along‐strike diachroneity, reflected in 25 Myr younger monazite and ⁴⁰Ar–³⁹Ar hornblende ages in the eastern relative to the western PAG domain, cannot be accounted for by heat productivity contrasts along the belt. Instead the younger deformation and metamorphism in the eastern CBb was driven by its proximity to the eastern promontory of the Superior Province which collided with the Rae Province at c. 1820 Ma. The ⁴⁰Ar–³⁹Ar data presented here support the interpretation that the youngest monazite in the CBb crystallized at c. 1790 Ma in the central CBb when this part of the belt was downfolded into a gentle synformal structure while the western part of the belt cooled through ⁴⁰Ar–³⁹Ar hornblende closure. The results of this study illustrate the important influence of contrasting rock properties on the thermal evolution of orogenic belts and on the temporal record of this evolution.     

Parameterization and multi‐criteria calibration of a distributed storm flow model applied to a Mediterranean agricultural catchment

The principal challenge in the parameterization of storm flow models for agricultural catchments with an artificial drainage network and fields with different degrees of tillage lies in the parsimonious definition of distributed model parameters in a way that reduces the number of calibration parameters to a justifiable minimum. This paper presents a comprehensive case study for the parameter estimation of a distributed storm flow model applied to an agricultural catchment (0.91 km²) in the Mediterranean region. Model parameterization was combined with procedures for multi‐criteria, multi‐storm calibration, where we automatically calibrated three parameters related to flow velocity and infiltration, and compared single and multi‐storm criteria that are based on discharge volume, peak flow, and the Nash–Sutcliffe coefficient. Multi‐storm calibration yielded a set of parameter values for the simulation batch with best multi‐storm overall performance, which are close to the median values in the pre‐calibration of individual storms. Our results suggest that flow velocities and proportionality of the channel infiltration rate do not vary significantly over the course of 11 years. Copyright © 2012 John Wiley & Sons, Ltd.     

A hydro‐economic modelling framework for flood damage estimation and the role of riparian vegetation

A modelling framework for the quick estimate of flood inundation and the resultant damages is developed in this paper. The model, called the flood economic impact analysis system (FEIAS), can be applied to a river reach of any hydrogeological river basin. For the development of the integrated modelling framework, three models were employed: (1) a modelling scheme based on the Hydrological Simulation Program FORTRAN model that was developed for any geomorphological river basin, (2) a river flow/floodplain model, and (3) a flood loss estimation model. The first sub‐model of the flood economic impact analysis system simulates the hydrological processes for extended periods of time, and its output is used as input to a second component, the river/floodplain model. The hydraulic model MIKE 11 (quasi‐2D) is the river/floodplain model employed in this study. The simulated flood parameters from the hydraulic model MIKE 11 (quasi‐2D) are passed, at the end of each time step, to a third component, the flood loss model for the estimation of flood damage. In the present work, emphasis was given to the seasonal variation of Manning's coefficient (n), which is an important parameter for the determination of the flood inundation in hydraulic modelling. High values of Manning's coefficient for a channel indicate high flow resistance. The riparian vegetation can have a large impact on channel resistance. The modelling framework developed in this paper was used to investigate the role of riparian vegetation in reducing flood damage. Moreover, it was used to investigate the influence of cutting riparian vegetation scenarios on the flow characteristics. The proposed framework was applied to the downstream part of the Koiliaris River basin in Crete, Greece, and was tested and validated with historical data. Copyright © 2012 John Wiley & Sons, Ltd.     

The need for operational reasoning in data‐driven rating curve prediction of suspended sediment

The use of data‐driven modelling techniques to deliver improved suspended sediment rating curves has received considerable interest in recent years. Studies indicate an increased level of performance over traditional approaches when such techniques are adopted. However, closer scrutiny reveals that, unlike their traditional counterparts, data‐driven solutions commonly include lagged sediment data as model inputs, and this seriously limits their operational application. In this paper, we argue the need for a greater degree of operational reasoning underpinning data‐driven rating curve solutions and demonstrate how incorrect conclusions about the performance of a data‐driven modelling technique can be reached when the model solution is based upon operationally invalid input combinations. We exemplify the problem through the re‐analysis and augmentation of a recent and typical published study, which uses gene expression programming to model the rating curve. We compare and contrast the previously published solutions, whose inputs negate their operational application, with a range of newly developed and directly comparable traditional and data‐driven solutions, which do have operational value. Results clearly demonstrate that the performance benefits of the published gene expression programming solutions are dependent on the inclusion of operationally limiting, lagged data inputs. Indeed, when operationally inapplicable input combinations are discounted from the models and the analysis is repeated, gene expression programming fails to perform as well as many simpler, more standard multiple linear regression, piecewise linear regression and neural network counterparts. The potential for overstatement of the benefits of the data‐driven paradigm in rating curve studies is thus highlighted. Copyright © 2011 John Wiley & Sons, Ltd.     

Ice‐jam flood risk assessment and mapping

In northern regions, river ice‐ jam flooding can be more severe than open‐water flooding causing property and infrastructure damages, loss of human life and adverse impacts on aquatic ecosystems. Very little has been performed to assess the risk induced by ice‐related floods because most risk assessments are limited to open‐water floods. The specific objective of this study is to incorporate ice‐jam numerical modelling tools (e.g. RIVICE, Monte‐Carlo simulation) into flood hazard and risk assessment along the Peace River at the Town of Peace River (TPR) in Alberta, Canada. Adequate historical data for different ice‐jam and open‐water flooding events were available for this study site and were useful in developing ice‐affected stage‐frequency curves. These curves were then applied to calibrate a numerical hydraulic model, which simulated different ice jams and flood scenarios along the Peace River at the TPR. A Monte‐Carlo analysis was then carried out to acquire an ensemble of water level profiles to determine the 1 : 100‐year and 1 : 200‐year annual exceedance probability flood stages for the TPR. These flood stages were then used to map flood hazard and vulnerability of the TPR. Finally, the flood risk for a 200‐year return period was calculated to be an average of $32/m²/a ($/m²/a corresponds to a unit of annual expected damages or risk). Copyright © 2016 John Wiley & Sons, Ltd.     

Quantifying the P–T–t conditions of north–south Lhasa terrane accretion: new insight into the pre‐Himalayan architecture of the Tibetan plateau

An integrated field, petrological and geochronological study of the Basong Tso region of south‐eastern Tibet has constrained the timing and P–T conditions of north–south Lhasa terrane accretion and provides new insight into the tectonothermal evolution of the Tibetan plateau. Two distinct high‐grade metamorphic belts are recognized in the region: a southern belt (the Basong Tso complex) that consists of sheared schist and orthogneiss; and a northern belt (the Zhala complex) that comprises paragneiss and granite. Combined pseudosection modelling and U–Pb geochronology of monazite and zircon indicates that the Basong Tso complex records peak metamorphic conditions of 9 ± 0.5 kbar and 690 ± 25 °C at c. 204–201 Ma, whereas the Zhala complex experienced peak metamorphic conditions of 5.0 ± 1.0 kbar and 740 ± 40 °C at c. 198–192 Ma. Microstructural analysis suggests that the two belts share a common early prograde history, after which the Basong Tso complex attained peak conditions following rapid burial, and the Zhala complex approached peak conditions along an isobaric path. Overall it is inferred that the Basong Tso and Zhala complexes represent the lower and upper structural levels of an evolving orogen that underwent Barrovian‐type metamorphism following collision (M1), followed by Buchan‐style overprinting at higher structural levels due to heat advection by syn‐tectonic granites (M2). Mylonitization (sensu lato) of the Basong Tso complex and juxtaposition of the two units occurred after attainment of peak conditions. The dominance of Mesozoic regional metamorphism across most of the Tibetan plateau indicates that Cenozoic crustal thickening processes, where present, are only manifested at depth.