Evaluation of two satellite-based products against ground-based observation for drought analysis in the southern part of Iran

Natural Hazards - Water stress or more specifically drought assessment plays a key role in water management, especially in extreme climate conditions. Basically, globally gridded satellite-based...     

Oxygen isotope speedometry in granulite facies garnet recording fluid/melt–rock interaction (Sør Rondane Mountains,East Antarctica)

In situ analysis of a garnet porphyroblast from a granulite facies gneiss from Sør Rondane Mountains, East Antarctica, reveals discontinuous step‐wise zoning in phosphorus and large δ¹⁸O variations from the phosphorus‐rich core to the phosphorus‐poor rim. The gradually decreasing profile of oxygen isotope from the core (δ¹⁸O = ~15‰) to the rim (δ¹⁸O = ~11‰) suggests that the ¹⁸O/¹⁶O zoning was originally step‐wise, and modified by diffusion after the garnet rim formation at ~800°C and 0.8 GPa. Fitting of the ¹⁸O/¹⁶O data to the diffusion equation constrains a duration of the high‐T event (~800°C) to c. 0.5–40 Ma after the garnet rim formation. The low δ¹⁸O value of the garnet rim, together with the previously reported low δ¹⁸O values in metacarbonates, indicates regional infiltration, probably along a detachment fault, of low δ¹⁸O fluid/melt possibly derived from meta‐mafic to ultramafic rocks.     

Simulating urban growth in a megalopolitan area using a patch‐based cellular automata

Although traditional cellular automata (CA)‐based models can effectively simulate urban land‐use changes, they typically ignore the spatial evolution of urban patches, due to their use of cell‐based simulation strategies. This research proposes a new patch‐based CA model to incorporate a spatial constraint based on the growth patterns of urban patches into the conventional CA model for reducing the uncertainty of the distribution of simulated new urban patches. In this model, the growth pattern of urban patches is first estimated using a developed indicator that is based on the local variations in existing urban patches. The urban growth is then simulated by integrating the estimated growth pattern and land suitability using a pattern‐calibrated method. In this method, the pattern of new urban patches is gradually calibrated toward the dominant growth pattern through the steps of the CA model. The proposed model is applied to simulate urban growth in the Tehran megalopolitan area during 2000–2006–2012. The results from this model were compared with two common models: cell‐based CA and logistic‐patch CA. The proposed model yields a degree of patch‐level agreement that is 23.4 and 7.5% higher than those of these pre‐existing models, respectively. This reveals that the patch‐based CA model simulates actual development patterns much better than the two other models.     

Basaltic and Solution Reference Materials for Iron,Copper and Zinc Isotope Measurements

Iron, Cu and Zn stable isotope systems are applied in constraining a variety of geochemical and environmental processes. Secondary reference materials have been developed by the Institute of Geology, Chinese Academy of Geological Sciences (CAGS), in collaboration with other participating laboratories, comprising three solutions (CAGS‐Fe, CAGS‐Cu and CAGS‐Zn) and one basalt (CAGS‐Basalt). These materials exhibit sufficient homogeneity and stability for application in Fe, Cu and Zn isotopic ratio determinations. Reference values were determined by inter‐laboratory analytical comparisons involving up to eight participating laboratories employing MC‐ICP‐MS techniques, based on the unweighted means of submitted results. Isotopic compositions are reported in per mil notation, based on reference materials IRMM‐014 for Fe, NIST SRM 976 for Cu and IRMM‐3702 for Zn. Respective reference values of CAGS‐Fe, CAGS‐Cu and CAGS‐Zn solutions are as follows: δ⁵⁶Fe = 0.83 ± 0.07 and δ⁵⁷Fe = 1.20 ± 0.13, δ⁶⁵Cu = 0.57 ± 0.06, and δ⁶⁶Zn = ?0.79 ± 0.12 and δ⁶⁸Zn = ?1.65 ± 0.24, respectively. Those of CAGS‐Basalt are δ⁵⁶Fe = 0.15 ± 0.07, δ⁵⁷Fe = 0.22 ± 0.10, δ⁶⁵Cu = 0.12 ± 0.08, δ⁶⁶Zn = 0.17 ± 0.13, and δ⁶⁸Zn = 0.34 ± 0.26 (2s).     

Clinoform geometry,geomorphology, facies character and stratigraphic architecture of a sand‐rich subaqueous delta: Jurassic Sognefjord Formation,offshore Norway

The integration of core sedimentology, seismic stratigraphy and seismic geomorphology has enabled interpretation of delta‐scale (i.e. tens of metres high) subaqueous clinoforms in the upper Jurassic Sognefjord Formation of the Troll Field. Mud‐prone subaqueous deltas characterized by a compound clinoform morphology and sandy delta‐scale subaqueous clinoforms are common in recent tide‐influenced, wave‐influenced and current‐influenced settings, but ancient examples are virtually unknown. The data presented help to fully comprehend the criteria for the recognition of other ancient delta‐scale subaqueous clinoforms, as well as refining the depositional model of the reservoir in the super‐giant Troll hydrocarbon field. Two 10 to 60 m thick, overall coarsening‐upward packages are distinguished in the lower Sognefjord Formation. Progressively higher energy, wave‐dominated or current‐dominated facies occur from the base to the top of each package. Each package corresponds to a set of seismically resolved, westerly dipping clinoforms, the bounding surfaces of which form the seismic ‘envelope’ of a clinoform set and the major marine flooding surfaces recognized in cores. The packages thicken westwards, until they reach a maximum where the clinoform ‘envelope’ rolls over to define a topset–foreset–toeset geometry. All clinoforms are consistently oriented sub‐parallel to the edge of the Horda Platform (N005–N030). In the eastern half of the field, individual foresets are relatively gently dipping (1° to 6°) and bound thin (10 to 30 m) clinothems. Core data indicate that these proximal clinothems are dominated by fine‐grained, hummocky cross‐stratified sandstones. Towards the west, clinoforms gradually become steeper (5° to 14°) and bound thicker (15 to 60 m) clinothems that comprise medium‐grained, cross‐bedded sandstones. Topsets are consistently well‐developed, except in the westernmost area. No seismic or sedimentological evidence of subaerial exposure is observed. Deposition created fully subaqueous, near‐linear clinoforms that prograded westwards across the Horda Platform. Subaqueous clinoforms were probably fed by a river outlet in the north‐east and sculpted by the action of currents sub‐parallel to the clinoform strike.     

Revised age of proximal deposits in the Zagros foreland basin and implications for Cenozoic evolution of the High Zagros

The regionally extensive, coarse-grained Bakhtiyari Formation represents the youngest synorogenic fill in the Zagros foreland basin of Iran. The Bakhtiyari is present throughout the Zagros fold-thrust belt and consists of conglomerate with subordinate sandstone and marl. The formation is up to 3000 m thick and was deposited in foredeep and wedge-top depocenters flanked by fold-thrust structures. Although the Bakhtiyari concordantly overlies Miocene deposits in foreland regions, an angular unconformity above tilted Paleozoic to Miocene rocks is expressed in the hinterland (High Zagros).

The Bakhtiyari Formation has been widely considered to be a regional sheet of Pliocene–Pleistocene conglomerate deposited during and after major late Miocene–Pliocene shortening. It is further believed that rapid fold growth and Bakhtiyari deposition commenced simultaneously across the fold-thrust belt, with limited migration from hinterland (NE) to foreland (SW). Thus, the Bakhtiyari is generally interpreted as an unmistakable time indicator for shortening and surface uplift across the Zagros. However, new structural and stratigraphic data show that the most-proximal Bakhtiyari exposures, in the High Zagros south of Shahr-kord, were deposited during the early Miocene and probably Oligocene. In this locality, a coarse-grained Bakhtiyari succession several hundred meters thick contains gray marl, limestone, and sandstone with diagnostic marine pelecypod, gastropod, coral, and coralline algae fossils. Foraminiferal and palynological species indicate deposition during early Miocene time. However, the lower Miocene marine interval lies in angular unconformity above ~ 150 m of Bakhtiyari conglomerate that, in turn, unconformably caps an Oligocene marine sequence. These relationships attest to syndepositional deformation and suggest that the oldest Bakhtiyari conglomerate could be Oligocene in age.

The new age information constrains the timing of initial foreland-basin development and proximal Bakhtiyari deposition in the Zagros hinterland. These findings reveal that structural evolution of the High Zagros was underway by early Miocene and probably Oligocene time, earlier than commonly envisioned. The age of the Bakhtiyari Formation in the High Zagros contrasts significantly with the Pliocene–Quaternary Bakhtiyari deposits near the modern deformation front, suggesting a long-term (> 20 Myr) advance of deformation toward the foreland.     

Towards more systematic perceptual model development: a case study using 3 Luxembourgish catchments

The synthesis of experimental understanding of catchment behaviour and its translation into qualitative perceptual models is an important objective of hydrological sciences. We explore this challenge by examining the cumulative understanding of the hydrology of three experimental catchments and how it evolves through the application of different investigation techniques. The case study considers the Huewelerbach, Weierbach and Wollefsbach headwater catchments of the Attert basin in Luxembourg. Subsurface investigations including bore holes and pits, analysis of soil samples and Electrical Resistivity Tomography measurements are presented and discussed. Streamflow and tracer data are used to gain further insights into the streamflow dynamics of the catchments, using end‐member mixing analysis and hydrograph separation based on dissolved silica and electrical conductivity. We show that the streamflow generating processes in all three catchments are controlled primarily by the subsolum and underlying bedrock. In the Huewelerbach, the permeable sandstone formation supports a stable groundwater component with little seasonality, which reaches the stream through a series of sources at the contact zone with the impermeable marls formation. In the Weierbach, the schist formation is relatively impermeable and supports a ‘fill and spill’‐type of flow mechanism; during wet conditions, it produces a delayed response dominated by pre‐event water. In the Wollefsbach, the impermeable marls formation is responsible for a saturation‐excess runoff generating process, producing a fast and highly seasonal response dominated by event water. The distinct streamflow generating processes of the three catchments are represented qualitatively using perceptual models. The perceptual models are in turn translated into quantitative conceptual models, which simulate the hydrological processes using networks of connected reservoirs and transfer functions. More generally, the paper illustrates the evolution of perceptual models based on experimental fieldwork data, the translation of perceptual models into conceptual models and the value of different types of data for processes understanding and model representation. Copyright © 2014 John Wiley & Sons, Ltd.     

Assessment of alternative adsorption models and global sensitivity analysis to characterize hexavalent chromium loss from soil to surface runoff

We investigate our ability to assess transfer of hexavalent chromium, Cr(VI), from the soil to surface runoff by considering the effect of coupling diverse adsorption models with a two‐layer solute transfer model. Our analyses are grounded on a set of two experiments associated with soils characterized by diverse particle size distributions. Our study is motivated by the observation that Cr(VI) is receiving much attention for the assessment of environmental risks due to its high solubility, mobility, and toxicological significance. Adsorption of Cr(VI) is considered to be at equilibrium in the mixing layer under our experimental conditions. Four adsorption models, that is, the Langmuir, Freundlich, Temkin, and linear models, constitute our set of alternative (competing) mathematical formulations. Experimental results reveal that the soil samples characterized by the finest grain sizes are associated with the highest release of Cr(VI) to runoff. We compare the relative abilities of the four models to interpret experimental results through maximum likelihood model calibration and four model identification criteria (i.e., the Akaike information criteria [AIC and AIC_C] and the Bayesian and Kashyap information criteria). Our study results enable us to rank the tested models on the basis of a set of posterior weights assigned to each of them. A classical variance‐based global sensitivity analysis is then performed to assess the relative importance of the uncertain parameters associated with each of the models considered, within subregions of the parameter space. In this context, the modelling strategy resulting from coupling the Langmuir isotherm with a two‐layer solute transfer model is then evaluated as the most skilful for the overall interpretation of both sets of experiments. Our results document that (a) the depth of the mixing layer is the most influential factor for all models tested, with the exception of the Freundlich isotherm, and (b) the total sensitivity of the adsorption parameters varies in time, with a trend to increase as time progresses for all of the models. These results suggest that adsorption has a significant effect on the uncertainty associated with the release of Cr(VI) from the soil to the surface runoff component.     

Effect of hysteretic smoothness on inelastic response spectra with constant‐ductility

This paper presents an efficient methodology for computing constant‐ductility inelastic response spectra. The computation of constant‐ductility spectra involves numerical root‐finding algorithms to find the strongest structure providing a desired ductility response. Smooth inelastic structural behavior is modeled using a first‐order nonlinear differential equation and the transient structural response is solved using an implicit algorithm requiring Newton iterations at each time step. For structural models with smooth hysteretic behavior (not piece‐wise linear), a simple root‐finding method involving a combination of hyperbolic fits, linear interpolation, and Newton's method converges upon the highest strength (conservative) solution with a small number of iterations. The effect of the hysteretic smoothness on the occurrence of multiple roots is examined for two near‐fault and two far‐fault earthquake records, and for two measures of ductility and for normalized hysteretic energy. The results indicate how the smoothness of the hysteretic behavior affects ductility demand and constant‐ductility response spectra. Copyright © 2010 John Wiley & Sons, Ltd.     

Stream channel erosion in a rapidly urbanizing region of the US–Mexico border: documenting the importance of channel hardpoints with Structure‐from‐Motion photogrammetry

Urbanization can lead to accelerated stream channel erosion, especially in areas experiencing rapid population growth, unregulated urban development on erodible soils, and variable enforcement of environmental regulations. A combination of field surveys and Structure‐from‐Motion (SfM) photogrammetry techniques was used to document spatial patterns in stream channel geometry in a rapidly urbanizing watershed, Los Laureles Canyon (LLCW), in Tijuana, Mexico. Ground‐based SfM photogrammetry was used to map channel dimensions with 1 to 2 cm vertical mean error for four stream reaches (100–300 m long) that were highly variable and difficult to survey with a differential GPS. Regional channel geometry curves for LLCW had statistically larger slopes and intercepts compared with regional curves developed for comparable, undisturbed reference channels. Cross‐sectional areas of channels downstream of hardpoints, such as concrete reaches or culverts, were up to 64 times greater than reference channels, with enlargement persisting, in some cases, up to 230 m downstream. Percentage impervious cover was not a good predictor of channel enlargement. Proximity to upstream hardpoint, and lack of riparian and bank vegetation paired with highly erodible bed and bank materials may account for the instability of the highly enlarged and unstable cross‐sections. Channel erosion due to urbanization accounts for approximately 25–40% of the total sediment budget for the watershed, and channel erosion downstream of hardpoints accounts for one‐third of all channel erosion. Channels downstream of hardpoints should be stabilized to prevent increased inputs of sediment to the Tijuana Estuary and local hazards near the structures, especially in areas with urban settlements near the stream channel. Copyright © 2017 John Wiley & Sons, Ltd.