DuMu: DUNE for multi-{phase, component, scale, physics, …} flow and transport in porous media

DuMu^x is a free and open-source simulator for flow and transport processes in porous media, based on the Distributed and Unified Numerics Environment (DUNE). Its main intention is to provide a sustainable and consistent framework for the implementation and application of model concepts, constitutive relations, discretizations, and solvers. The paper provides an overview of DuMu^x with the focus on software-related aspects. Selected examples highlight the multi-scale and the parallel capabilities.     

Flood‐loss modelling for assessing impacts of flood‐frequency adjustment,Middle Mississippi River,USA

This study modelled flood losses (economic damages) along the Middle Mississippi River (MMR) (1) using current US government estimates of flow frequencies and (2) using frequencies based on the original, unaltered discharge measurements. The official flood frequencies were quantified in the Upper Mississippi River System Flow Frequency Study (UMRSFFS), but as a last step in that study, early discharges along the MMR were reduced by up to 54% to reflect a purported bias in early measurements. Subsequently, early discharge measurements were rigorously tested, and no such bias was found. Here, flood damages were quantified using a combination of one‐dimensional hydraulic modelling and flood‐loss modelling. For all recurrence intervals, damages were much less using the UMRSFFS flow frequencies compared with the frequencies based on the original discharge measurements, with differences ranging up to 79% (100‐year event) and $2.9bn (200‐year event). Annualized losses in the study area based on the UMRSFFS frequencies were just $41.6m versus $125.6m using the raw frequencies (an underestimation of 67%). These totals do not include flood losses elsewhere along the MMR, including in metropolitan St Louis. In summary, a seemingly small methodological adjustment – in this case, a single hidden adjustment, not documented anywhere within the UMRSFFS – can have dramatic societal impacts in terms of underestimation of flood probabilities and flood risk. Copyright © 2012 John Wiley & Sons, Ltd.     

Hydrologic‐response simulations for the North Fork of Caspar Creek: second‐growth,clear‐cut,new‐growth,and cumulative watershed effect scenarios

This study demonstrates that comprehensive hydrologic‐response simulation can be a useful tool for studying cumulative watershed effects. The simulations reported here were conducted with the Integrated Hydrology Model (InHM). The location of the 473 ha study site is the North Fork of the Caspar Creek Experimental Watershed, near Fort Bragg, California. Existing information from a long‐term monitoring programme and new soil‐hydraulic property measurements made for this study were used to parameterize InHM. Long‐term continuous wet‐season simulations were conducted for the North Fork catchments and main stem for second‐growth, clear‐cut and new‐growth scenarios. The simulation results show that the increases and decreases, respectively, for throughfall and potential evapotranspiration related to clear‐cutting had quantifiable impacts on the simulated hydrologic response at both the catchment and watershed scales. Model performance was best for the new‐growth simulation scenarios. To improve upon the simulations reported here would require additional soil‐hydraulic property information from across the study area. Although principally focused on the integrated hydrologic response, the effort reported here demonstrates the potential for characterizing distributed responses with physics‐based simulation. The search for a comprehensive understanding of hydrologic response will require both data‐intensive discovery and concept‐development simulation, from both integrated and distributed perspectives. Copyright © 2013 John Wiley & Sons, Ltd.     

Multi‐angle,multi‐damage fragility curves for seismic assessment of bridges

The scope of this study is to investigate the effect of the direction of seismic excitation on the fragility of an already constructed, 99‐m‐long, three‐span highway overpass. First, the investigation is performed at a component level, quantifying the sensitivity of local damage modes of individual bridge components (namely, piers, bearings, abutments, and footings) to the direction of earthquake excitation. The global vulnerability at the system level is then assessed for a given angle of incidence of the earthquake ground motion to provide a single‐angle, multi‐damage probabilistic estimate of the bridge overall performance. A multi‐angle, multi‐damage, vulnerability assessment methodology is then followed, assuming uniform distribution for the angle of incidence of seismic waves with respect to the bridge axis. The above three levels of investigation highlight that the directivity of ground motion excitation may have a significant impact on the fragility of the individual bridge components, which shall not be a priori neglected. Most importantly, depending on the assumptions made for the component to the system level transition, this local sensitivity is often suppressed. It may be therefore necessary, based on the ultimate purpose of the vulnerability or the life cycle analysis, to obtain a comprehensive insight on the multiple damage potential of all individual structural and foundation components under multi‐angle excitation, to quantify the statistical correlation among the distinct damage modes and to identify the components that are both most critical and sensitive to the direction of ground motion and carefully define their limit states which control the predicted bridge fragility. Copyright © 2015 John Wiley & Sons, Ltd.     

N<Subscript>2</Subscript>, Ar,and He as a tool for discriminating sources of volcanic fluids with application to Vulcano,Italy

A detailed analysis of published data on the N₂, Ar, and He content and Ar and He isotopic composition of fumarolic fluids from Vulcano crater (south Italy) supports a model with two endmembers comprising magmatic and hydrothermal fluids with correspondingly low and high H₂O content. The magmatic component with the highest ³He/⁴He and highest absolute concentrations of N₂, Ar, and He also has the lowest N₂/Ar and N₂/He ratios (∼300 and ∼500, respectively). In contrast, the hydrothermal endmember, with the lower ³He/⁴He and lower absolute N₂, Ar, and He abundances, has high N₂/Ar (∼1,000) and high N₂/He (>3,000) ratios. The hydrothermal component is also characterized by the highest ⁴⁰Ar/³⁶Ar ratios (>1,000) and is proposed to be the main carrier of metamorphic gases from the arc crust.     

Chemical geothermometry and fluid–mineral equilibria for the Ömer–Gecek thermal waters, Afyon area, Turkey

Thermal waters of the Ömer–Gecek geothermal field, Turkey, with temperatures ranging from 32 to 92°C vary in chemical composition and TDS contents. They are generally enriched in Na–Cl–HCO₃ and suggest deep water circulation. Silica and cation geothermometers applied to the Ömer–Gecek thermal waters yield reservoir temperatures of 75–155°C. The enthalpy–chloride mixing model, which approximates a reservoir temperature of 125°C for the Ömer–Gecek field, accounts for the diversity in the chemical composition and temperature of the waters by a combination of processes including boiling and conductive cooling of deep thermal water and mixing of the deep thermal water with cold water. It is also determined that the solubility of silica in most of the waters is controlled by the chalcedony phase. Equilibrium states of the Ömer–Gecek thermal waters studied by means of the Na–K–Mg triangular diagram, Na–K–Mg–Ca diagram, K–Mg–Ca geoindicator diagram, activity diagrams in the systems composed of Na₂O–CaO–K₂O–Al₂O₃–SiO₂–CO₂–H₂O phases, log SI diagrams, and finally the alteration mineralogy indicate that most of the spring and low-temperature well waters in the area can be classified as shallow or mixed waters which are likely to be equilibrated with calcite, chalcedony and kaolinite at predicted temperature ranges similar to those calculated from the chemical geothermometers. It was also observed that mineral equilibrium in the Ömer–Gecek waters is largely controlled by CO₂ concentrations.     

Network‐scale dynamics of grain‐size sorting: implications for downstream fining,stream‐profile concavity,and drainage basin morphology

We explore the link between channel‐bed texture and river basin concavity in equilibrium catchments using a numerical landscape evolution model. Theory from homogeneous sediment transport predicts that river basin concavity directly increases with bed sediment size. If the effective grain size on a river bed governs its concavity, then natural phenomena such as grain‐size sorting and channel armouring should be linked to concavity. We examine this hypothesis by allowing the bed sediment texture to evolve in a transport‐limited regime using a two grain‐size mixture of sand and gravel. Downstream ?ning through selective particle erosion is produced in equilibrium. As the channel‐bed texture adjusts downstream so does the local slope. Our model predicts that it is not the texture of the original sediment mixture that governs basin concavity. Rather, concavity is linked to the texture of the sorted surface layer. Two different textural regimes are produced in the experiments: a transitional regime where the mobility of sand and gravel changes with channel‐bed texture, and a sand‐dominated region where the mobility of sand and gravel is constant. The concavity of these regions varies depending on the median gravel‐ or sand‐grain size, erosion rate, and precipitation rate. The results highlight the importance of adjustments in both surface texture and slope in natural rivers in response to changes in ?uvial and sediment inputs throughout a drainage network. This adjustment can only be captured numerically using multiple grain sizes or empirical downstream ?ning rules. Copyright © 2004 John Wiley & Sons, Ltd.     

Conditions for generation of fire‐related debris flows,Capulin Canyon,New Mexico

Comparison of the responses of three drainage basins burned by the Dome fire of 1996 in New Mexico is used to identify the hillslope, channel and fire characteristics that indicate a susceptibility specifically to wildfire‐related debris flow. Summer thunderstorms generated three distinct erosive responses from each of three basins. The Capulin Canyon basin showed widespread erosive sheetwash and rilling from hillslopes, and severe flooding occurred in the channel; the North Tributary basin exhibited extensive erosion of the mineral soil to a depth of 5 cm and downslope movement of up to boulder‐sized material, and at least one debris flow occurred in the channel; negligible surface runoff was observed in the South Tributary basin. The negligible surface runoff observed in the South Tributary basin is attributed to the limited extent and severity of the fire in that basin. The factors that best distinguish between debris‐flow producing and flood‐producing drainages are drainage basin morphology and lithology. A rugged drainage basin morphology, an average 12 per cent channel gradient, and steep, rough hillslopes coupled with colluvium and soil weathered from volcaniclastic and volcanic rocks promoted the generation of debris flows. A less rugged basin morphology, an average gradient of 5 per cent, and long, smooth slopes mantled with pumice promoted flooding. Flood and debris‐flow responses were produced without the presence of water‐repellent soils. The continuity and severity of the burn mosaic, the condition of the riparian vegetation, the condition of the fibrous root mat, accumulations of dry ravel and colluvial material in the channel and on hillslopes, and past debris‐flow activity, appeared to have little bearing on the distinctive responses of the basins. Published in 2000 by John Wiley & Sons, Ltd.     

A 13 ± 3 ka age determination of a tholeiite,Pinacate volcanic field,Mexico, and improved methods for Ar/Ar dating of young basaltic rocks

Among the youngest lava flows of the Pinacate volcanic field, Sonora, Mexico, is a large outpouring of tholeiite, the Ives flow. This tube-fed pahoehoe flow contrasts sharply with other Pinacate lavas in its great volume, alkali-poor composition and morphologic features, which include novel small structures named “spatter tubes.” Despite its K-poor character, young age, and the presence of excess ⁴⁰Ar, we determined a ⁴⁰Ar/³⁹Ar age on samples of this flow at 13 ± 3 ka. Such an age determination is made possible via careful monitoring of the mass discrimination of the mass spectrometer and by stacking results from multiple incremental-heating experiments into a single, composite isochron. This age is among the youngest ever to be determined with such precision by the ⁴⁰Ar/³⁹Ar method on a K-poor tholeiite.     

Proposed GSSP for the base of Cambrian Stage 9, coinciding with the first appearance of <Emphasis Type="Italic">Agnostotes orientalis</Emphasis>, at Duibian,Zhejiang, China

Exposed in natural outcrops near the Duibian Village, Jiangshan County, Zhejiang Province, China, the Duibian B section is proposed as the boundary stratotype for the base of an unnamed stage termed provisionally Cambrian Stage 9. The proposed position of the GSSP is 108.12 m above the base of the Huayansi Formation, at a horizon coinciding with the first appearance of the cosmopolitan agnostoid trilobite Agnostotes orientalis. This horizon coincides also with the first appearance of the cosmopolitan polymerid trilobite Irvingella angustilimbata. The section fulfills all the requirements for a GSSP, and the horizon can be constrained not only by the primary stratigraphic marker (A. orientalis) but also with secondary biostratigraphic, chemostratigraphic, and sequence-stratigraphic correlation tools. The first appearance of A. orientalis is one of the most readily recognizable levels in Cambrian, and can be correlated with precision to all paleocontinents. Supported by Chinese Academy of Sciences (Grant No. KZCX2-YW-122), the Ministry of Science and Technology of China (Grant Nos. 2006FY120300, 2006CB806400), National Natural Science Foundation of China (Grant Nos. 40072003, 40023002, 40332018, 40672023, 40602002), and the U.S. National Science Foundation (Grant No. EAR-0106883)     

Synthesis,Mechanism, and Performance Assessment of Zero‐Valent Iron for Metal‐Contaminated Water Remediation: A Review

Recently, increased industrial and agriculture activities have resulted in toxic metal ions, which has increased public concern about the quality of surface and groundwater. Various types of physical, biological, and chemical approaches have been developed to remove surface and groundwater metal ions contaminants. Among these practices, zero‐valent iron (ZVI) is the most studied reactive material for environmental clean‐up over the last two decade and so. Although ZVI can remove the contaminants even more efficiently than any other reactive materials. However, low reactivity due to its intrinsic passive layer, narrow working pH, and the loss of hydraulic conductivity due to the precipitation of metal hydroxides and metal carbonates limits its wide‐scale application. The aim of this work is to document properties, synthesis, and reaction mechanism of ZVI for the treatment of metal ions from the surface and groundwater in recent 10 years (2008–2018). So far, different modified techniques such as conjugation with support, bimetal alloying, weak magnetic field, and ZVI/oxidant coupling system have been developed to facilitate the use of ZVI in various environmental remediation scenarios. However, some challenges still remain to be addressed. Therefore, development and research in this field are needed to overcome or mitigate these limitations.     

Joint inversion of GNSS and teleseismic data for the rupture process of the 2017 <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript>6.5 Jiuzhaigou,China, earthquake

The spatio-temporal slip distribution of the earthquake that occurred on 8 August 2017 in Jiuzhaigou, China, was estimated from the teleseismic body wave and near-field Global Navigation Satellite System (GNSS) data (coseismic displacements and high-rate GPS data) based on a finite fault model. Compared with the inversion results from the teleseismic body waves, the near-field GNSS data can better restrain the rupture area, the maximum slip, the source time function, and the surface rupture. The results show that the maximum slip of the earthquake approaches 1.4 m, the scalar seismic moment is ~ 8.0 × 10¹⁸ N·m (M_w?≈?6.5), and the centroid depth is ~ 15 km. The slip is mainly driven by the left-lateral strike-slip and it is initially inferred that the seismogenic fault occurs in the south branch of the Tazang fault or an undetectable fault, a NW-trending left-lateral strike-slip fault, and belongs to one of the tail structures at the easternmost end of the eastern Kunlun fault zone. The earthquake rupture is mainly concentrated at depths of 5–15 km, which results in the complete rupture of the seismic gap left by the previous four earthquakes with magnitudes >?6.0 in 1973 and 1976. Therefore, the possibility of a strong aftershock on the Huya fault is low. The source duration is ~ 30 s and there are two major ruptures. The main rupture occurs in the first 10 s, 4 s after the earthquake; the second rupture peak arrives in ~ 17 s. In addition, the Coulomb stress study shows that the epicenter of the earthquake is located in the area where the static Coulomb stress change increased because of the 12 May 2017 M_w7.9 Wenchuan, China, earthquake. Therefore, the Wenchuan earthquake promoted the occurrence of the 8 August 2017 Jiuzhaigou earthquake.