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.     

Sm–Nd, Rb–Sr and K–Ar geochronology of the Higo metamorphic terrane, west-central Kyushu, Japan

The Higo metamorphic terrane situated in west-central Kyushu island, southwest Japan, is composed of greenschist- to granulite-facies metamorphic rocks. The southern part of the metamorphic terrane consists mainly of garnet–biotite gneiss and garnet–cordierite–biotite gneiss, and orthopyroxene or cordierite-bearing S-type tonalite with subordinate amounts of hornblende gabbro. Rb–Sr, Sm–Nd and K–Ar isotopic ages for these rocks have been determined here. The garnet–biotite gneiss gives an Sm–Nd age of 227.1 ± 4.9 Ma and a Rb–Sr age of 101.0 ± 1.0 Ma. The hornblende gabbro has an Sm–Nd age of 257.9 ± 2.5 Ma and a K–Ar age of 103.4 ± 1.1 Ma. These age differences of the same samples are due to the difference in the closure temperature for each system and minerals. The garnet-cordierite–biotite gneiss contains coarse-grained garnet with a zonal structure conspicuously distinguished in color difference (core: dark red; rim: pink). Sm–Nd internal isochrons of the garnet core and the rim give ages of 278.8 ± 4.9 Ma (initial ¹⁴³Nd/¹⁴⁴Nd ratio = 0.512311 ± 0.000005) and 226.1 ± 28.4 Ma (0.512277 ± 0.000038), respectively. These ages are close to formation of the garnet core and the rim. It has been previously suggested that the Higo metamorphic terrane belongs to the Ryoke metamorphic belt. But Sr and Nd isotopic features of the rocks from the former are different from those of the Ryoke metamorphic rocks, and are similar to those of the granulite xenoliths contained in the Ryoke younger granite.     

Using sediment ‘fingerprints’ to assess sediment-budget errors,North Halawa Valley,Oahu, Hawaii, 1991–92

Reliable estimates of sediment-budget errors are important for interpreting sediment-budget results. Sediment-budget errors are commonly considered equal to sediment-budget imbalances, which may underestimate actual sediment-budget errors if they include compensating positive and negative errors. We modified the sediment ‘fingerprinting’ approach to qualitatively evaluate compensating errors in an annual (1991) fine (<63 μm) sediment budget for the North Halawa Valley, a mountainous, forested drainage basin on the island of Oahu, Hawaii, during construction of a major highway. We measured concentrations of aeolian quartz and ¹³⁷Cs in sediment sources and fluvial sediments, and combined concentrations of these aerosols with the sediment budget to construct aerosol budgets. Aerosol concentrations were independent of the sediment budget, hence aerosol budgets were less likely than sediment budgets to include compensating errors. Differences between sediment-budget and aerosol-budget imbalances therefore provide a measure of compensating errors in the sediment budget. The sediment-budget imbalance equalled 25 per cent of the fluvial fine-sediment load. Aerosol-budget imbalances were equal to 19 per cent of the fluvial ¹³⁷Cs load and 34 per cent of the fluvial quartz load. The reasonably close agreement between sediment- and aerosol-budget imbalances indicates that compensating errors in the sediment budget were not large and that the sediment-budget imbalance is a reliable measure of sediment-budget error. We attribute at least one-third of the 1991 fluvial fine-sediment load to highway construction. Continued monitoring indicated that highway construction produced 90 per cent of the fluvial fine-sediment load during 1992. Erosion of channel margins and attrition of coarse particles provided most of the fine sediment produced by natural processes. Hillslope processes contributed relatively minor amounts of sediment. © 1998 John Wiley & Sons, Ltd.     

A jadeite–quartz–glaucophane rock from Karangsambung, central Java, Indonesia

High-pressure metamorphic rocks are exposed in Karangsambung area of central Java, Indonesia. They form part of a Cretaceous subduction complex (Luk–Ulo Complex) with fault-bounded slices of shale, sandstone, chert, basalt, limestone, conglomerate and ultrabasic rocks. The most abundant metamorphic rock type are pelitic schists, which have yielded late Early Cretaceous K–Ar ages. Small amounts of eclogite, glaucophane rock, garnet–amphibolite and jadeite–quartz–glaucophane rock occur as tectonic blocks in sheared serpentinite. Using the jadeite–garnet–glaucophane–phengite–quartz equilibrium, peak pressure and temperature of the jadeite–quartz–glaucophane rock are P  = 22 ± 2 kbar and T  = 530 ± 40 °C. The estimated P–T conditions indicate that the rock was subducted to ca 80 km depth, and that the overall geothermal gradient was ∼ 7.0 °C/km. This rock type is interpreted to have been generated by the metamorphism of cold oceanic lithosphere subducted to upper mantle depths. The exhumation from the upper mantle to lower or middle crustal depths can be explained by buoyancy forces. The tectonic block is interpreted to be combined with the quartz–mica schists at lower or middle crustal depths.     

It's life, Jim…

Lewis Dartnell wonders how astronomers find life on this and other planets. After all, they do it on Star Trek with just a scanner…     

Dendrogeomorphic dating of rockfalls on low‐latitude,high‐elevation slopes: Rodadero,Iztaccíhuatl volcano,Mexico

Dynamics and rates of rockfalls have been repeatedly studied in mountain environments with archival records as well as lichenometric, radiocarbon or dendrogeomorphic approaches. In this study, we test the potential of conifers growing at a low‐latitude, high‐elevation site as a dendrogeomorphic tool to reconstruct to calendar dates associated rockfall activity. Analysis is based on tree‐ring records of Mexican mountain pine (Pinus hartwegii Lindl.) growing at timberline [~4000 m above sea level (a.s.l.)] and at the runout fringe of a north–northeast (NNE)‐facing slope of the dormant Iztaccíhuatl volcano (Mexico), which is subject to frequent rockfalls. The potential and limitations of tree‐ring data are demonstrated based on 67 rockfall impacts dated in the increment‐ring series of 24 trees since ad 1836. While findings of this paper are site‐specific, the study clearly shows the potential of dendrogeomorphic approaches in extra‐Alpine, low‐latitude environments and for the understanding of rockfall processes in space and time. Copyright © 2011 John Wiley & Sons, Ltd.     

The December 26, 2004, off the west coast of northern Sumatra,Indonesia, <Emphasis Type="Italic">M</Emphasis><Subscript>W</Subscript>=9.0, earthquake and the critical-point-like model of earthquake preparation

Long-term seismic activity prior to the December 26, 2004, off the west coast of northern Sumatra, Indonesia, M _W=9.0 earthquake was investigated using the Harvard CMT catalogue. It is observed that before this great earthquake, there exists an accelerating moment release (AMR) process with the temporal scale of a quarter century and the spatial scale of 1 500 km. Within this spatial range, the M _W=9.0 event falls into the piece-wise power-law-like frequency-magnitude distribution. Therefore, in the perspective of the critical-point-like model of earthquake preparation, the failure to forecast/predict the approaching and/or the size of this earthquake is not due to the physically intrinsic unpredictability of earthquakes.