Upscaling permeability for three-dimensional fractured porous rocks with the multiple boundary method

Upscaling permeability of grid blocks is crucial for groundwater models. A novel upscaling method for three-dimensional fractured porous rocks is presented. The objective of the study was to compare this method with the commonly used Oda upscaling method and the volume averaging method. First, the multiple boundary method and its computational framework were defined for three-dimensional stochastic fracture networks. Then, the different upscaling methods were compared for a set of rotated fractures, for tortuous fractures, and for two discrete fracture networks. The results computed by the multiple boundary method are comparable with those of the other two methods and fit best the analytical solution for a set of rotated fractures. The errors in flow rate of the equivalent fracture model decrease when using the multiple boundary method. Furthermore, the errors of the equivalent fracture models increase from well-connected fracture networks to poorly connected ones. Finally, the diagonal components of the equivalent permeability tensors tend to follow a normal or log-normal distribution for the well-connected fracture network model with infinite fracture size. By contrast, they exhibit a power-law distribution for the poorly connected fracture network with multiple scale fractures. The study demonstrates the accuracy and the flexibility of the multiple boundary upscaling concept. This makes it attractive for being incorporated into any existing flow-based upscaling procedures, which helps in reducing the uncertainty of groundwater models.     

Combined effects of pumped-storage operation and climate change on thermal structure and water quality

Climatic Change - The assessment of ecological impacts of pumped-storage (PS) hydropower plants on the two connected water bodies is usually based on present climatic conditions. However,...     

Microbial photosynthesis in coral reef sediments (Heron Reef, Australia)

We investigated microphytobenthic photosynthesis at four stations in the coral reef sediments at Heron Reef, Australia. The microphytobenthos was dominated by diatoms, dinoflagellates and cyanobacteria, as indicated by biomarker pigment analysis. Conspicuous algae firmly attached to the sand grains (ca. 100 μm in diameter, surrounded by a hard transparent wall) were rich in peridinin, a marker pigment for dinoflagellates, but also showed a high diversity based on cyanobacterial 16S rDNA gene sequence analysis. Specimens of these algae that were buried below the photic zone exhibited an unexpected stimulation of respiration by light, resulting in an increase of local oxygen concentrations upon darkening. Net photosynthesis of the sediments varied between 1.9 and 8.5 mmol O₂ m⁻² h⁻¹ and was strongly correlated with Chl a content, which lay between 31 and 84 mg m⁻². An estimate based on our spatially limited dataset indicates that the microphytobenthic production for the entire reef is in the order of magnitude of the production estimated for corals. Photosynthesis stimulated calcification at all investigated sites (0.2–1.0 mmol Ca²⁺ m⁻² h⁻¹). The sediments of at least three stations were net calcifying. Sedimentary N₂-fixation rates (measured by acetylene reduction assays at two sites) ranged between 0.9 to 3.9 mmol N₂ m⁻² h⁻¹ and were highest in the light, indicating the importance of heterocystous cyanobacteria. In coral fingers no N₂-fixation was measurable, which stresses the importance of the sediment compartment for reef nitrogen cycling.     

Capture into first-order resonances and long-term stability of pairs of equal-mass planets

Massive planets form within the lifetime of protoplanetary disks, and therefore, they are subject to orbital migration due to planet–disk interactions. When the first planet reaches the inner edge of the disk, its migration stops and consequently the second planet ends up locked in resonance with the first one. We detail how the resonant trapping works comparing semi-analytical formulae and numerical simulations. We restrict to the case of two equal-mass coplanar planets trapped in first-order resonances, but the method can be easily generalized. We first describe the family of resonant stable equilibrium points (zero-amplitude libration orbits) using series expansions up to different orders in eccentricity as well as a non-expanded Hamiltonian. Then we show that during convergent migration the planets evolve along these families of equilibrium points. Eccentricity damping from the disk leads to a final equilibrium configuration that we predict precisely analytically. The fact that observed multi-exoplanetary systems are rarely seen in resonances suggests that in most cases the resonant configurations achieved by migration become unstable after the removal of the protoplanetary disk. Here we probe the stability of the resonances as a function of planetary mass. For this purpose, we fictitiously increase the masses of resonant planets, adiabatically maintaining the low-amplitude libration regime until instability occurs. We discuss two hypotheses for the instability, that of a low-order secondary resonance of the libration frequency with a fast synodic frequency of the system, and that of minimal approach distance between planets. We show that secondary resonances do not seem to impact resonant systems at low amplitude of libration. Resonant systems are more stable than non-resonant ones for a given minimal distance at close encounters, but we show that the latter nevertheless play the decisive role in the destabilization of resonant pairs. We show evidence that as the planetary mass increases and the minimal distance between planets gets smaller in terms of mutual Hill radius, the region of stability around the resonance center shrinks, until the equilibrium point itself becomes unstable.     

Checking vulnerability to pure loss of stability in long crested following waves: A probabilistic approach

This paper presents a probabilistic methodology for the analysis of the vulnerability of a ship to the risk of inception of pure loss of stability events. A pure loss of stability failure is modelled as the persistence of the metacentric height below a critical level for a too long time. The metacentric height is modelled as a stationary Gaussian process with a spectrum obtained from the sea elevation spectrum. The time dependent failure index is obtained under the assumption of filtered Poisson process for the occurrence of critical events. The analysis separates cases where the fluctuation of the metacentric height is narrow-band from those where the bandwidth of the spectrum is wide, with an intermediate blending. In case of narrow-band processes appropriate approximate solutions to the problem are provided, while in the wide-band cases an exponential distribution for the persistence time below the critical level is employed. A rational development for the critical persistence time is also provided considering an approximation of the roll dynamics during periods of time where the metacentric height is negative. Monte Carlo simulations are performed to check the developed approximate distributions for the persistence time, and examples of application are provided for a sample ship.     

Hydrodynamic and isotopic investigations for evaluating the mechanisms and amount of groundwater seepage through a rockslide dam

This study addresses the influence of landslide dams on surface water drainage and groundwater flow. In the study area of Scanno Lake and Sagittario River (Central Italy), a limestone rockslide‐avalanche formed a lake, which has an outlet that is occasionally active, showing infiltration into the rockslide dam. Several springs are present at the lake's base and are partly fed by seepage through the rockslide debris. Piezometric surveys, discharge measurements, pumping tests and chemical analyses are tools used to build a conceptual model of the groundwater flow and to evaluate the flow through the rockslide debris. Seasonal water isotopic signatures validate the assumed model, showing a mixing of infiltration recharge and groundwater seepage throughout the rockslide debris. Various recharge areas have been found for springs, pointing out those directly fed by the rockslide debris aquifer. Hypotheses about seasonal groundwater mixing between the regional carbonate aquifer and the rockslide debris aquifer are supported by isotope results. Seasonal changes in groundwater table level due to recharge and surface losses from seasonal outlet have been correlated with isotopic groundwater composition from the rockslide debris aquifer and the downstream springs; this relationship highlights the role of the rockslide dam body on the hydrodynamics of the studied area. Relationships between surface waters and groundwater in the area have been completely understood on the basis of water isotopic fingerprinting, finally obtaining a complete evaluation of groundwater renewable resources and its regimen. Copyright © 2010 John Wiley & Sons, Ltd.     

Paläomagnetische Untersuchungen im Schilfsandstein (Trias,Km2) Westeuropas

Zusammenfassung Von insgesamt 46 beprobten Profilen aus dem Schilfsandstein (Trias, Km2) Westeuropas wurden bislang 41 Profile mit 343 Einheiten und 870 Handstücken mit 1897 Einzelkernen auf ihren paläomagnetischen Inhalt untersucht.Thermische Entmagnetisierungsversuche an 277 Proben zeigen, daß die natürliche remanente Magnetisierung (NRM) von mehreren, an verschiedene Magnetisierungsträger gebundenen Komponenten getragen wird.Als charakteristische remanente Magnetisierung (CARM), von der angenommen wird, daß sie das zur Ablagerungszeit herrschende Paläofeld abbildet, kristallisiert sich für die Sedimentationsräume Weserbergland und Franken eine Deklination (DR=Deklination rotiert) von 60° und eine Inklination (IR=Inklination rotiert) von 30° (DR/IR 60/30) und für Schwaben und die Nordschweiz DR/IR 240/–30 heraus.Neben der stabilen Paläofeldrichtung werden auch sedimentäre und tektonische Gefüge als weitere stabile Magnetisierungsrichtungen abgebildet. Diese Richtungen liegen der Achse maximaler Suszeptibilität (K1 max) parallel.Unter Annahme einer normalen Polarität des Erdmagnetfeldes zur Zeit der Schilfsandsteinablagerung im Weserbergland und Franken und einer reversen Polarität zur Zeit der Ablagerung in Schwaben und der Nordschweiz errechnet sich ein Paläonordpol bei 111 E 28 N; er fällt damit in die Nähe des Triaspoles der Europäischen Polwanderungskurve (McElhinny, 1973).Für den Ablagerungszeitraum des Raibler Sandsteins in den Südalpen wird ein Paläonordpol bei 128 W 48 N errechnet.Die magnetostratigraphische Zeitskala der Trias (Burek, 1970) kann um die Einstufung des Schilfsandsteins in das höhere obere Karn mit beginnender normaler Polarität und dann reverser Polarität ergänzt werden.

---

In the Schilfsandstein (Triassic, Km2) of Western Europe 41 sections with 343 units and 1897 cores were sampled for palaeomagnetic investigations. By thermal demagnetization experiments different magnetic minerals were determined as carriers of the natural rémanent magnetization (NRM).The characteristic remanent magnetization (CARM) reflecting the palaeofield direction is described in the Weserbergland and Franken by DR/IR 60/30 (DR = declination after tectonic correction/IR=inclination after tectonic correction) and in Schwaben and Northern Switzerland by DR/IR 240/–30.In addition stable directions are reported which are related to sedimentary and tectonic fabric. These directions are parallel to the axis of maximum susceptibility (Klmax).Assuming normal polarity for the magnetic field during deposition of the Schilfsandstein in the Weserbergland and Franken and assuming reversed polarity in Schwaben and Northern Switzerland a palaeonorthpole is estimated at 111 E 28 N. This poleposition fits near the Triassic poleposition of the European polar wandering curve (McElhinny, 1973).The poleposition of the Bergamasker Alps (Raibler Sandstein, Southern Alps) is figured out at 128 W 48 N.The triassic magnetostratigraphic time scale ofBurek (1970) can be completed as far as the Schilfsandstein is to arrange into the highest Upper Carnian starting with normal polarity going on with reversed polarity.

---

Résumé Dans le Grès à roseaux (Trias, Km 2) de l'Europe occidentale, 343 prélèvements provenant de 41 profils et 1897 carrottes extraites de 870 échantillons ont fait l'objet de recherches paléomagnétiques. — Des essais de démagnétisation thermique sur 277 échantillons ont révélé différents minéraux magnétiques comme étant les porteurs de l'aimantation rémanente naturelle (NRM). — L'aimentation rémanente caractéristique (CARM) reflétant la direction du paléochamp magnétique au moment de la sédimentation donne des valeurs de 60/30 pour DR/IR (DR=déclinaison après correction tectonique, — IR=inclinaison après correction tectonique) dans le Weserbergland et en Franconie et 240/– 30 pour DR/IR en Souabe et dans la Suisse septentrionale. — En plus de la direction stable du paléo champ, on a trouvé d'autres directions d'aimentation stables reliées aux structures sédimentaires et tectoniques. Ces directions sont parallèles à l'axe de la susceptibilité magnétique maximale (K1 max). — En supposant une polarité normale du champ magnétique terrestre au moment de la sédimentation du Grès à roseaux dans le Weserbergland et la Franconie, et une polarité inversée pendant sa formation en Souabe et en Suisse septentrionale, il en résulte une position géographique de 111 E et 28° N pour le paléopôle Nord, ce qui coïncide à peu près avec la position du pôle triasique sur la courbe polaire de l'Europe (McElhinny, 1973). Pendant la sédimentation du Grès de Raibler dans les Alpes bergamasques, le paléopôle Nord se trouvait à 128° O 48° N. — L'echelle magnétostratigraphique pour le Trias (Burek, 1970) peut maintenant être complétée en datant le Grès à roseaux dans la partie supérieure du Carnien supérieur, débutant par une aimantation normale et continuant par une aimantation inversée.

---

343 41 ( ) 1897 870 . 277 , — NRM — . — CARM —, , , , , (DR = Deklination rotiert) 60°, (IR = Inklination rotiert) 30° (DR/IR 60/30), DR/IR 240/-30. , . /_1makc. , 111 E 28 N, . ., . 128 W 48 N. (Burek, 1970) . . , .