宁南深层岩溶地下水系统三维数值模拟

宁南“南北古脊梁”岩溶裂隙水系统的储水空间以岩溶裂隙为主,岩溶水的空间分布明显受南北向大型断裂构造控制。根据水动力场和水化学场特征,划分了3个相对独立的岩溶水子系统。通过三维数值模拟计算,确定了马渠—洪河岩溶水子系统天然资源量和可开采资源量,为本区的地下水合理开发利用提供了科学依据。     

Crystal morphologies in pillow basalts: implications for mid-ocean ridge processes

We present the results of a detailed petrological study of a sparsely phyric basalt (MAPCO CH98-DR11) dredged along the Mid-Atlantic Ridge (30°41′N). The sample contains microphenocrysts of olivine that display four different rapid-growth morphologies. Comparison of these morphologies with those obtained in dynamic crystallization experiments allows us to constrain the thermal history of the sample. The dendritic morphology (swallowtail, chain and lattice olivine) is directly related to the final quenching during magma–seawater interaction. In contrast, the three other morphologies, namely the complex polyhedral crystal, the closed hopper and the complex swallowtail morphology result from several cycles of cooling–heating (corresponding to a maximum degree of undercooling of 20–25°C) during crystal growth. These thermal variations occurred before eruption and are interpreted to be the result of turbulent convection in a small magmatic body beneath the ridge. The results suggest that the Mid-Atlantic Ridge is underlain by a mush zone that releases batches of liquid during tectonic segregation. Aphyric basalts are emitted during eruptions controlled by the tectonic activity, whereas phyric basalts correspond to small fractions of magma from the mush zone mobilized by reinjections of primitive magmas.     

Magmatic effects of the lunar late heavy bombardment

Large volumes of mare basalts are present on the surface of the moon, located preferentially in large impact basins. Mechanisms relating impact basins and mare basalt eruptions have previously been suggested: lunar impacts removed low-density material that may have inhibited eruption, and created cracks for fluid flow [Icarus 139 (1999) 246], and lunar basins have long been described as catchments for magma (e.g., [Rev. Geophys. Space Phys. 18 (1980) 107] and references therein). We present a new model for melt creation under near side lunar basins that is triggered by the impacts themselves. Magma can be produced in two stages. First, crater excavation depressurizes underlying material such that it may melt in-situ. Second, the cratered lithosphere rises isostatically, warping isotherms at the lithosphere-asthenosphere boundary which may initiate convection, in which adiabatic melting can occur. The first stage produces by far the largest volume of melt, but convective melting can continue for up to 350 Ma. We propose that giant impacts account for a large portion of the volume and longevity of mare basalt volcanism, as well as for several compositional groups, including high alumina, high titanium, KREEP-rich, and picritic magmas.     

The influence of thermochemical convection on the fixity of mantle plumes

A general feature of both isochemical and thermochemical studies of mantle convection is that horizontal plume velocities tend to be smaller than typical convective velocities, however, it is not clear which system leads to a greater fixity of mantle plumes. We perform two- and three-dimensional numerical calculations and compare both thermochemical and isochemical cases with similar convective vigor to determine whether presence of a dense component in the mantle can lead to smaller ratios of horizontal plume velocity to surface velocity. We investigate different viscosity and density contrasts between chemical components in the thermochemical calculations, and we perform isochemical calculations with both free-slip and no-slip bottom boundary conditions. We then compare both visually and quantitatively the results of the thermochemical and isochemical calculations to determine which leads to greater plume fixity. We find that horizontal plume velocities for thermochemical calculations are similar to those from isochemical calculations with no-slip bottom boundary conditions. In addition, we find that plumes tend to be more fixed for isochemical cases with free-slip bottom boundary conditions for two-dimensional calculations, however, in three dimensions, we find that plume fixity is similar to that observed in thermochemical calculations.     

Formation and shortening deformation of a back-arc rift basin revealed by deep seismic profiling, central Japan

The northern Fossa Magna (NFM) basin is a Miocene rift system produced in the final stages of the opening of the Sea of Japan. It divides the major structure of Japan into two regions, with north-trending geological structures to the NE of the basin and EW trending structures to the west of the basin. The Itoigawa-Shizuoka Tectonic Line (ISTL) bounds the western part of the northern Fossa Magna and forms an active fault system that displays one of the largest slip rates (4–9 mm/year) in the Japanese islands. Deep seismic reflection and refraction/wide-angle reflection profiling were undertaken in 2002 across the northern part of ISTL in order to delineate structures in the crust, and the deep geometry of the active fault systems. The seismic images are interpreted based on the pattern of reflectors, the surface geology and velocities derived from refraction analysis. The 68-km-long seismic section suggests that the Miocene NFM basin was formed by an east dipping normal fault with a shallow flat segment to 6 km depth and a deeper ramp penetrating to 15 km depth. This low-angle normal fault originated as a comparatively shallow brittle/ductile detachment in a high thermal regime present in the Miocene. The NFM basin was filled by a thick (>6 km) accumulation of sediments. Shortening since the late Neogene is accommodated along NS to NE–SE trending thrust faults that previously accommodated extension and produce fault-related folds on their hanging wall. Based on our balanced geologic cross-section, the total amount of Miocene extension is ca. 42 km and the total amount of late Neogene to Quaternary shortening is ca. 23 km.     

Nonlinear free thermal convection in a spherical shell:A variable viscosity model

Introduction The velocity field of surface plate motion can be split into a poloidal and a toroidal parts.At the Earth′s surface,the toroidal component is manifested by the existence of transform faults,and the poloidal component by the presence of convergence and divergence,i.e.spreading and subduc-tion zones.They have coupled each other and completely depicted the characteristics of plate tec-tonic motions.The mechanism of poloidal field has been studied fairly clearly which is related to …     

Numerical model of convection in sills

The numerical model of convection in magma sills is developed. The model is based on a full system of equations of fluid dynamics and includes heat transfer, buoyancy effects and diffusion of some minor component (marker). Solidification is treated as a phase transition. The results indicate that there are some qualitative differences between very thin sills with Rayleigh number Ra = 10⁵ and thin sills with Ra = 10⁶. For a basaltic magma the first case corresponds to the thickness of the sills of approximately 30 cm and the second case corresponds to the thickness of 60 cm. In the first case mixing is inefficient and conduction is the dominant form of heat transfer. In the second case mixing is efficient and convection is the dominant form of heat transfer. Some of the results can be scaled for the more viscous magmas in thicker sills.