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.     

Alluvial architecture of the Holocene Rhine–Meuse delta (the Netherlands)

Ancient fluvial successions often act as hydrocarbon reservoirs. Sub‐surface data on the alluvial architecture of fluvial successions are often incomplete and modelling is performed to reconstruct the stratigraphy. However, all alluvial architecture models suffer from the scarcity of field data to test and calibrate them. The purposes of this study were to quantify the alluvial architecture of the Holocene Rhine–Meuse delta (the Netherlands) and to determine spatio‐temporal trends in the architecture. Five north–south orientated cross‐sections, perpendicular to the general flow direction, were compiled for the fluvial‐dominated part of the delta. These sections were used to calculate the width/thickness ratios of fluvial sandbodies (SBW/SBT) and the proportions of channel‐belt deposits (CDP), clastic overbank deposits (ODP) and organic material (OP) in the succession. Furthermore, the connectedness ratio (CR) between channel belts was calculated for each cross‐section. Distinct spatial and temporal trends in the alluvial architecture were found. SBW/SBT ratios decrease by a factor of ca 4 in a downstream direction. CDP decreases from ca 0·7 (upstream) to ca 0·3 (downstream). OP increases from less than 0·05 in the upstream part of the delta to more than 0·25 in the downstream delta. ODP is approximately constant (0·4). CR is ca 0·25 upstream, which is approximately two times larger than in the downstream part of the delta. Furthermore, CDP in the downstream Rhine–Meuse delta increases after 3000 cal yr BP. These trends are attributed to variations in available accommodation space, floodplain geometry and channel‐belt size. For instance, channel belts tend to narrow in a downstream direction, which reduces SBW/SBT, CDP and CR. Tectonics cause local deviations in the general architectural trends. In addition, the positive correlation between avulsion frequency and the ratio of local to regional aggradation rate probably influenced alluvial architecture in the Rhine–Meuse delta. The Rhine–Meuse data set can be a great resource when developing more sophisticated models for alluvial architecture simulation, which eventually could lead to better characterizations of hydrocarbon reservoirs. To aid such usage of the Rhine–Meuse data set, constraints for relevant parameters are provided at the end of the paper.     

First seismic record for intra-arc strike-slip tectonics along the Liquiñe-Ofqui fault zone at the obliquely convergent plate margin of the southern Andes

A temporal seismic network recorded local seismicity along a 130 km long segment of the transpressional dextral strike-slip Liquiñe-Ofqui fault zone (LOFZ) in southern Chile. Seventy five shallow crustal events with magnitudes up to M_w 3.8 and depths shallower than 25 km were observed in an 11-month period mainly occurring in different clusters. Those clusters are spatially related to the LOFZ, to the volcanoes Chaitén, Michinmahuida and Corcovado, and to active faulting on secondary faults. Further activity along the LOFZ is indicated by individual events located in direct vicinity of the surface expression of the LOFZ. Focal mechanisms were calculated using deviatoric moment tensor inversion of body wave amplitude spectra which mostly yield strike-slip mechanisms indicating a NE–SW direction of the P-axis for the LOFZ at this latitude. The seismic activity reveals the present-day activity of the fault zone. The recent M_w 6.2 event near Puerto Aysén, Southern Chile at 45.4°S on April 21, 2007 shows that the LOFZ is also capable of producing large magnitude earthquakes and therefore imposing significant seismic hazard to this region.     

Age estimates of coastal terraces in the Andaman and Nicobar Islands and their tectonic implications

The great Indian Ocean earthquake of December 26, 2004 caused significant vertical changes in its rupture zone. About 800 km of the rupture is along the Andaman and Nicobar Islands, which forms the outer arc ridge of the subduction zone. Coseismic deformation along the exposed land could be observed as uplift/subsidence. Here we analyze the morphological features along the coast of the Andaman and Nicobar Islands, in an effort to reconstruct the past tectonics, taking cues from the coseismic effects. We obtained radiocarbon dates from coastal terraces of the island belt and used them to compute uplift rates, which vary from 1.33 mm yr^− 1 in the Little Andaman to 2.80 mm yr^− 1 in South Andaman and 2.45 mm yr^− 1 in the North Andaman. Our radiocarbon dates converge on  600 yr and  1000 yr old coastal uplifts, which we attribute to the level changes due to two major previous subduction earthquakes in the region.     

Morphological analysis of the drainage system in the Eastern Alps

We study the morphology of the major rivers draining the Eastern Alps to test whether the active tectonics of this part of the orogen is reflected in the shape of channel profiles of the river network. In our approach we compare channel profiles measured from digital elevation models with numerically modelled channel profiles using a stream power approach. It is shown that regions of high stream power coincide largely with regions of highest topography and largest uplift rates, while the forelands and the Pannonian Basin are characterised by a significantly lower stream power. From stream power modelling we conclude that there is young uplift at the very east of the Eastern Alps, in the Bohemian Massif and in the Pohorje Range. The impact of the Pleistocene glaciations is explored by comparing properties of rivers that drain in proximal and distal positions relative to the ice sheet during the last glacial maximum. Our analysis shows that most knick points, wind gaps and other non-equilibrium features of catchments covered by ice during the last glaciations (Salzach, Enns) can be correlated with glacial processes. In contrast the ice free catchments of the Mur and Drava are characterized by channels in morphological equilibrium at the first approximation and are showing only weak evidence of the strong tectonic activity within these catchments. Finally, the channel profiles of the Adige and the divide between the upper Rhine and Danube catchments differ significantly from the other catchments. We relate this to the fact that the Adige and the Rhine respond to different base levels from the remainder of the Eastern Alps: The Adige may preserve a record from the Messininan base level change and the Rhine is subject to the base level lowering in the Rhine Graben.     

Migration of the Ganga river and its implication on hydro-geological potential of Varanasi area, U.P., India

Borehole data reveals that during Late Quaternary, the Ganga river was non-existent in its present location near Varanasi. Instead, it was flowing further south towards peripheral craton. Himalayan derived grey micaceous sands were being carried by southward flowing rivers beyond the present day water divide of Ganga and mixed with pink arkosic sand brought by northward flowing peninsular rivers. Subsequently, the Ganga shifted to its present position and got incised. Near Varanasi, the Ganga river is flowing along a NW-SE tectonic lineament. The migration of Ganga river is believed to have been in response to basin expansion caused due to Himalayan tectonics during Middle Pleistocene times. Multi-storied sand bodies generated as a result of channel migration provide excellent aquifers confined by a thick zone of muddy sediments near the surface. Good quality potable water is available at various levels below about 70 m depth in sandy aquifers. Craton derived gravelly coarse-to-medium grained sand forms the main aquifer zones of tens of meter thickness with enormous yield. In contrast, the shallow aquifers made up of recycled interfluve silt and sandy silt occur under unconfined conditions and show water-level fluctuation of a few meters during pre-and post-monsoon periods.     

Triassic Nb-enriched basalts, magnesian andesites, and adakites of the Qiangtang terrane (Central Tibet): evidence for metasomatism by slab-derived melts in the mantle wedge

New chronological, geochemical, and isotopic data are reported for Triassic (219–236 Ma) adakite-magnesian andesite-Nb-enriched basaltic rock associations from the Tuotuohe area, central Qiangtang terrane. The adakites and magnesian andesites are characterized by high Sr/Y (25–45), La/Yb (14–42) and Na₂O/K₂O (12–49) ratios, high Al₂O₃ (15.34–18.28 wt%) and moderate to high Sr concentrations (220–498 ppm) and ε_ND (t) (+0.86 to +1.21) values. Low enrichments of Th, Rb relative to Nb, and subequal normalized Nb and La contents, and enrichments of light rare earth elements combine to distinguish a group of Nb-enriched basaltic rocks (NEBs). They have positive ε_ND (t) (+2.57 to +5.16) values. Positive correlations between Th, La and Nb and an absence of negative Nb anomalies on mantle normalized plots indicate the NEBs are products of a mantle source metasomatized by a slab melt rather than by hydrous fluids. A continuous compositional variation between adakites and magnesian andesites confirms slab melt interaction with mantle peridotite. The spatial association of the NEBs with adakites and magnesian andesites define an “adakitic metasomatic volcanic series” recognized in many demonstrably subduction-related environments (e.g., Mindanao arc, Philippines; Kamchatka arc, Russia; and southern Baja California arc, Mexico). The age of the Touhuohe suite, and its correlation with Triassic NEB to the north indicates that volcanism derived from subduction-modified mantle was abundant prior to 220 Ma in the central Qiangtang terrane.     

Structural elements of the Makran region, Oman sea and their potential relevance to tsunamigenisis

The character of convergence along the Arabian–Iranian plate boundary changes radically eastward from the Zagros ranges to the Makran region. This appears to be due to collision of continental crust in the west, in contrast to subduction of oceanic crust in the east. The Makran subduction zone with a length of about 900 km display progressively older and highly deformed sedimentary units northward from the coast, together with an increase in elevation of the ranges. North of the Makran ranges are large subsiding basins, flanked to the north by active volcanoes. Based on 2D seismic reflection data obtained in this study, the main structural provinces and elements in the Gulf of Oman include: (i) the structural elements on the northeastern part of the Arabian Plate and, (ii) the Offshore Makran Accretionary Complex. Based on detailed analysis of these data on the northeastern part of the Arabian Plate five structural provinces and elements—the Musendam High, the Musendam Peneplain, the Musendam Slope, the Dibba Zone, and the Abyssal Plain have been identified. Further, the Offshore Makran Accretionary Complex shown is to consist Accretionary Prism and the For-Arc Basin, while the Accretionary Prism has been subdivided into the Accretionary Wedge and the Accreted/Colored Mélange. Lastly, it is important to note that the Makran subduction zone lacks the trench. The identification of these structural elements should help in better understanding the seismicity of the Makran region in general and the subduction zone in particular. The 1945 magnitude 8.1 tsunamigenic earthquake of the Makran and some other historical events are illustrative of the coastal region’s vulnerability to future tsunami in the area, and such data should be of value to the developing Indian Ocean Tsunami Warning System.     

深水沉积层序特点及构成要素

本文在回顾当前国际上深水沉积研究热点的基础上,结合在墨西哥湾深水研究的成果系统描述了深水沉积的定义、形成机理、深水沉积层序及深水沉积构成要素的特点.深水沉积主要是在重力流作用下深水环境的沉积,主要形成于相对水平面下降和早期上升的时期,主要分布在低位体系域中.深水层序以凝缩段为边界,块状搬运沉积最早形成并直接位于层序界面上,其上被河道-天然堤沉积所覆盖.典型深水沉积的要素主要由河道、天然堤及越岸沉积、板状砂、块状搬运沉积等构成,这些沉积要素时空上有序地分布.深水河道是物源的主要通道和沉积的重要场所,从上游至下游河道弯曲度增加,能量逐渐减弱.侧向迁移明显,垂向上由富砂的顺直河道演化为相对富泥的弯曲河道.天然堤及越岸沉积以泥质为主,天然堤沿河道呈楔状分布,其近端砂岩含量高,地层厚且倾角较陡;远端砂岩含量低,地层薄且平缓,侧向连续性好但垂向连续性差.板状砂主要为深水扇前缘非限制性沉积,可分为块型和层型.块型侧向连续性好,同时垂向连通性高.层型侧向连续性好,垂向连通性差.块状搬运沉积主要是低水位期坡上沉积物失稳形成的各类滑塌体及碎屑流,其对下伏地层侵蚀明显,分布广泛,变形构造常见,可作为油气良好的封盖层.     

伶仃洋沉积动力特点的研究

位于珠江三角洲东侧的伶仃洋，因径流下泄与潮流进退的流向不一，使它各分流口门的出口水道都有主槽和支槽之分，即都有主干水道和分汊水道。陆架高盐海水入侵又使伶仃洋内沉积动力过程在空间分布上发生差异，如沉积物分布有粗－细－稍粗之分；而水体中的密度、速度差异，常常产生锋带，对水下地形的发展有不可忽视的影响。因而全面认识发生在伶仃洋内的沉积动力作用，对深水航道的选线极为重要。