Effects of the geometry of two‐dimensional fractures on their hydraulic aperture and on the validity of the local cubic law

Flow through rough fractures is investigated numerically in order to assess the validity of the local cubic law for different fracture geometries. Two‐dimensional channels with sinusoidal walls having different geometrical properties defined by the aperture, the amplitude, and the wavelength of the walls' corrugations, the corrugations asymmetry, and the phase shift between the two walls are considered to represent different fracture geometries. First, it is analytically shown that the hydraulic aperture clearly deviates from the mean aperture when the walls' roughness, the phase shift, and/or the asymmetry between the fracture walls are relatively high. The continuity and the Navier–Stokes equations are then solved by means of the finite element method and the numerical solutions compared to the theoretical predictions of the local cubic law. Reynolds numbers ranging from 0.066 to 66.66 are investigated so as to focus more particularly on the effect of flow inertial effects on the validity of the local cubic law. For low Reynolds number, typically less than 15, the local cubic law properly describes the fracture flow, especially when the fracture walls have small corrugation amplitudes. For Reynolds numbers higher than 15, the local cubic law is valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitudes, and a moderate phase lag between its walls.     

Geothermal palaeogradients and metamorphic zonation from the conodont colour alteration index (CAI)

The conodont colour alteration index (CAI) is potentially valuable in thermal history investigations, but there are problems in its use. An approach is proposed herein to overcome the problem of establishing reliable geothermal palaeogradients in areas with burial metamorphism, based on knowledge of the stratigraphic succession thickness, the boundaries of chronostratigraphic units, and the age and duration of deformation episodes. Although also problematical that a correlation between CAI and Kübler index (KI) values does not exist, there are areas where CAI and KI values are notably consistent for the anchizone boundaries; however, anomalous correlations between CAI and KI values have also been found in some areas. A CAI anchizone (or ancaizone) is defined in order to enable CAI values to be used independently to establish a metamorphic zonation.     

Lead isotope systematics of Late Cretaceous – Tertiary Andean arc magmas and associated ores between 8°N and 40°S: evidence for latitudinal mantle heterogeneity beneath the Andes

ABSTRACT Lead isotope variability of magmatic arc rocks and associated mineralization of the Central Andes is usually considered to be the result of mixing between a homogeneous mantle and heterogeneous continental crust. About 230 new lead isotope data on the Northern and Central Andes allow us to compare for the first time lead isotope systematics of the Late Cretaceous – Tertiary arc magmatism and associated mineralization along the Andean chain between 8°N and 40°S. Lead isotope compositions indicate mixing between mantle and upper crustal rocks along the whole Andean chain. Additionally, we have found that mantle end-members of the Late Cretaceous – Tertiary magmatism are heterogeneous and systematically shifted towards less radiogenic ²⁰⁶Pb/²⁰⁴Pb compositions from north to south along the Andes. This heterogeneity most likely results from mixing between a low radiogenic mantle, possibly carrying a DMM or EM I component, and a more radiogenic mantle, possibly carrying an HIMU component. Thus, our results imply that lead isotope variability of Andean magmas at the continental scale is caused not only by crustal but also by mantle heterogeneity.     

Light curve and O−C diagram analysis of RZ Cassiopeiae

The aim of the present paper is to analyze the light variation as well as the period changes of the eclipsing binary RZ Cas. New photometric elements are computed using the frequency-domain method. The possibility of the light-time effect and the apsidal motion is discussed.     

Pb, Sr and Nd isotope study of two co-existing magmas in the Nízke Tatry Mountains, Western Carpathians (Slovakia)

Summary Two co-existing plutonic rocks (diorite and granodiorite) were studied from an intrusion of Variscan age in the Raztocna Valley – Nízke Tatry Mountains, Western Carpathians. Geochemical analyses of major and trace elements constrain a volcanic arc as emplacement environment and give the first hints of a mixture of two magmatic end-members: the so-called Prasivá granodiorite and the Raztocna diorite. The ⁸⁷Sr/⁸⁶Sr₍₀₎ ratios vary between 0.7075 and 0.7118, the ε Nd₍₀₎ values range from −1.4 to −5.0. Common Pb isotopes reveal a dominant crustal source with minor influences from a mantle and a lower crustal source. Modelling based on Sr and Nd isotope data and using three component mixing calculations indicates that mixing of 2/3 of upper mantle material with 1/3 upper crustal material can produce the isotopic composition of the Raztocna diorite. Very minor amounts of lower crust were incorporated in the diorite. For the Prasivá granodiorite, the mixing ratio of upper mantle and upper crust is similar, but a lower crustal reservoir contributed about 5–10% of the source material.     

Quick-look assessments to identify optimal CO2 EOR storage sites

A newly developed, multistage quick-look methodology allows for the efficient screening of an unmanageably large number of reservoirs to generate a workable set of sites that closely match the requirements for optimal CO₂ enhanced oil recovery (EOR) storage. The objective of the study is to quickly identify miscible CO₂ EOR candidates in areas that contain thousands of reservoirs and to estimate additional oil recovery and sequestration capacities of selected top options through dimensionless modeling and reservoir characterization. Quick-look assessments indicate that the CO₂ EOR resource potential along the US Gulf Coast is 4.7 billion barrels, and CO₂ sequestration capacity is 2.6 billion metric tons. In the first stage, oil reservoirs are screened and ranked in terms of technical and practical feasibility for miscible CO₂ EOR. The second stage provides quick estimates of CO₂ EOR potential and sequestration capacities. In the third stage, a dimensionless group model is applied to a selected set of sites to improve the estimates of oil recovery and storage potential using appropriate inputs for rock and fluid properties, disregarding reservoir architecture and sweep design. The fourth stage validates and refines the results by simulating flow in a model that describes the internal architecture and fluid distribution in the reservoir. The stated approach both saves time and allows more resources to be applied to the best candidate sites.     

Quantitative sinkhole hazard assessment. A case study from the Ebro Valley evaporite alluvial karst (NE Spain)

Quantitative sinkhole hazard assessments in karst areas allow calculation of the potential sinkhole risk and the performance of cost-benefit analyses. These estimations are of practical interest for planning, engineering, and insurance purposes. The sinkhole hazard assessments should include two components: the probability of occurrence of sinkholes (sinkholes/km² year) and the severity of the sinkholes, which mainly refers to the subsidence mechanisms (progressive passive bending or catastrophic collapse) and the size of the sinkholes at the time of formation; a critical engineering design parameter. This requires the compilation of an exhaustive database on recent sinkholes, including information on the: (1) location, (2) chronology (precise date or age range), (3) size, and (4) subsidence mechanisms and rate. This work presents a hazard assessment from an alluvial evaporite karst area (0.81 km²) located in the periphery of the city of Zaragoza (Ebro River valley, NE Spain). Five sinkholes and four locations with features attributable to karstic subsidence where identified in an initial investigation phase providing a preliminary probability of occurrence of 0.14 sinkholes/km² year (11.34% in annual probability). A trenching program conducted in a subsequent investigation phase allowed us to rule out the four probable sinkholes, reducing the probability of occurrence to 0.079 sinkholes/km² year (6.4% in annual probability). The information on the severity indicates that collapse sinkholes 10–15 m in diameter may occur in the area. A detailed study of the deposits and deformational structures exposed by trenching in one of the sinkholes allowed us to infer a modern collapse sinkhole approximately 12 m in diameter and with a vertical throw of 8 m. This collapse structure is superimposed on a subsidence sinkhole around 80 m across that records at least 1.7 m of synsedimentary subsidence. Trenching, in combination with dating techniques, is proposed as a useful methodology to elucidate the origin of depressions with uncertain diagnosis and to gather practical information with predictive utility about particular sinkholes in alluvial karst settings: precise location, subsidence mechanisms and magnitude, and timing and rate of the subsidence episodes.     

Rain-induced debris and mudflow triggering factors assessment in the Santiago cordilleran foothills,Central Chile

Debris and mudflows are some of the main geological hazards in the mountain foothills of the Chilean capital city of Santiago. There, the risk of flows triggered in the basins of ravines that drain the range into the city increases with time due to the city growth. A multivariate statistical study based on the logistic regression method is presented. The model provides equations that allow the computation of combined meteorological triggering factors associated with a probability of rain-induced flow occurrence. Daily rainfall, accumulated rainfall and the snowfall level, traditionally considered as the relevant factors, are analysed for a 25-year period. The results show a strong relevance of the rainfall on the day of the flow event over the other factors. However, the relatively low probabilities returned for some real flow events suggest that the model does not capture all the significant variables and the problem is more complex than as it has been traditionally assumed, and further investigations are needed to develop predictive models of flow triggering.