线性非平衡态下矿物-水溶液界面的迁移系数与固相表面分形特征的关系

运用分形理论讨论了线性非平衡态下矿物（岩石）—水溶液界面的迁移系数,得出了迁移系数与固相颗粒表面分形特征的关系。     

地震层析成象技术在测井和地面反射资料集中的应用研究

本文以射线的折射和反射定律为基础,采用阻尼最小二乘法（在CT技术中,常称为SIRT算法）和单井单侧多偏移距的反射资料以及地面反射资料,研究反射波层析成象技术对二维弯曲界面的成层结构进行重构的情况。数字模拟试验表明,本文提出的逐层成象方法对透镜体和尖灭地层的成象颇为有效,反演结果对误差扰动是稳定的。 本文所述的成象方法及处理技术可推广到只用地面反射资料进行地层重构的情形。     

青藏铁路多年冻土区桥墩随机地震反应分析

利用随机振动理论和动力分析的有限元方法,将场地地震动考虑为均值为零的高斯平稳随机过程,对青藏铁路多年冻土区典型桥墩进行了随机地震反应分析,计算了9度地震作用下桥墩随机地震响应的统计特性,分析了冻土层对桥墩地震反应的影响。研究结果表明,冻土层对桥墩地震反应具有显著影响。     

Wenner法探测日本海岸附近淡盐境界面的实验研究

本文在理论上对电探查法探测由数层地质构造组成的海岸带含水层的淡盐境界面的可能性进行探讨的基础上,对太平洋沿岸某地区进行了实地测量研究。在测量设置在海岸附近及分布于海岸纵深方向的观测井铅直方向电导率的同时,在观测井附近进行了电探查。结果表明,电探查法探测的淡盐境界面的深度与依据井水电导率区分的混合区域上端几乎一致,即电探查法探测的淡盐境界面的位置在混合区域的上端。因此,该方法作为一种监测海水入侵的手段是非常适宜的,且与以往的手法比较也是一种简单、快速、成本低廉的确定淡盐境界面深度的方法。     

RC pile–soil interaction analysis using a 3D‐finite element method with fibre theory‐based beam elements

To evaluate the overall response of a structural system including its foundation and surrounding soil, an equivalent finite element model with reduced degrees of freedom using fibre theory‐based beam element was proposed. The proposed model was based on investigations of the subgrade soil reaction of a single‐layer model, and was verified for the cyclic behaviour of a laterally loaded single RC pile in terms of the load–displacement relationship, pile deformation, and soil pressures on the pile surface. Also investigated was the effect of the interfacial element between pile and soil on the behaviour of the laterally loaded pile. Copyright © 2006 John Wiley & Sons, Ltd.     

Tidal effects on variations of fresh–saltwater interface and groundwater flow in a multilayered coastal aquifer on a volcanic island (Jeju Island, Korea)

We conducted various field studies at the seawater intrusion monitoring wells located in the eastern part of Jeju Island, Korea, to observe the tidal effect on groundwater–seawater flow in the coastal aquifer. Studies included monitoring the fluctuations of groundwater and tide levels, electrical and temperature logging, and 2-D heat-pulse flowmeter tests. According to time-series analysis, tidal effects on groundwater level reached up to 3 km inland from the coastline. Water-level variation was more sensitive to tidal fluctuations near the coast, and more related to rainfall toward inland areas. Temporal and spatial variations in the shape and location of the freshwater–saltwater interface were analyzed using data from nine monitoring wells. The results indicated that the interface toe is located at a distance of 6–8 km from the coastline and its location was related to geological layers present. Long-term seasonal variations revealed no major changes in the interface; minor variations were due to moving boundary conditions induced by tidal fluctuations. Using the two-dimensional heat-pulse flowmeter, groundwater flow directions and velocities at four tidal stages were measured on three monitoring wells drilled into the multilayered aquifers. This direct measurement enabled us to relate the differences of flow velocities and directions with geology and tidal fluctuations. Combining the results of EC logging and flowmeter tests, we found a zone where freshwater and saltwater moved alternately in opposite directions, as influenced by the tidal fluctuations. Integrating various physical logging and flowmeter data with water-level fluctuations improved our understanding of the behavior of fresh and seawater flow in the coastal aquifers.     

Dolomite formation in the dynamic deep biosphere: results from the Peru Margin

Early diagenetic dolomite beds were sampled during the Ocean Drilling Programme (ODP) Leg 201 at four reoccupied ODP Leg 112 sites on the Peru continental margin (Sites 1227/684, 1228/680, 1229/681 and 1230/685) and analysed for petrography, mineralogy, δ¹³C, δ¹⁸O and ⁸⁷Sr/⁸⁶Sr values. The results are compared with the chemistry, and δ¹³C and ⁸⁷Sr/⁸⁶Sr values of the associated porewater. Petrographic relationships indicate that dolomite forms as a primary precipitate in porous diatom ooze and siliciclastic sediment and is not replacing the small amounts of precursor carbonate. Dolomite precipitation often pre‐dates the formation of framboidal pyrite. Most dolomite layers show ⁸⁷Sr/⁸⁶Sr‐ratios similar to the composition of Quaternary seawater and do not indicate a contribution from the hypersaline brine, which is present at a greater burial depth. Also, the δ¹³C values of the dolomite are not in equilibrium with the δ¹³C values of the dissolved inorganic carbon in the associated modern porewater. Both petrography and ⁸⁷Sr/⁸⁶Sr ratios suggest a shallow depth of dolomite formation in the uppermost sediment (<30 m below the seafloor). A significant depletion in the dissolved Mg and Ca in the porewater constrains the present site of dolomite precipitation, which co‐occurs with a sharp increase in alkalinity and microbial cell concentration at the sulphate–methane interface. It has been hypothesized that microbial ‘hot‐spots’, such as the sulphate–methane interface, may act as focused sites of dolomite precipitation. Varying δ¹³C values from −15‰ to +15‰ for the dolomite are consistent with precipitation at a dynamic sulphate–methane interface, where δ¹³C of the dissolved inorganic carbon would likewise be variable. A dynamic deep biosphere with upward and downward migration of the sulphate–methane interface can be simulated using a simple numerical diffusion model for sulphate concentration in a sedimentary sequence with variable input of organic matter. Thus, the study of dolomite layers in ancient organic carbon‐rich sedimentary sequences can provide a useful window into the palaeo‐dynamics of the deep biosphere.     

Modelling of short-interval silent slip events in deeper subduction interfaces considering the frictional properties at the unstable–stable transition regime

Recent high-resolution observations of crustal movements have revealed silent slip events (SSEs) with propagation velocities of around 10–15 km d⁻¹ and with intervals of 3–14 months along the deeper parts of the Cascadia and Nankai subduction zones. This study develops 2-D and 3-D models of these short-interval SSEs considering the frictional behaviour that was confirmed experimentally by Shimamoto for the unstable–stable transition regime. To represent this frictional behaviour, a small cut-off velocity to an evolution effect is introduced in a rate- and state-dependent friction law. When the cut-off velocity to the evolution effect is significantly smaller than that to a direct effect, steady-state friction exhibits velocity weakening at low slip velocities and velocity strengthening at high slip velocities. At the deeper Cascadia and Nankai subduction interfaces, the pore pressure is inferred to be high because of the dehydration of materials in the descending plate. Under conditions where the pore-fluid pressure is nearly equal to the lithostatic pressure and the critical weakening displacement is very small, short-interval SSEs with propagation velocities and slip velocities of 4–8 km d⁻¹ and  2 − 4 × 10⁻⁷  m s⁻¹, respectively, can be reproduced. The propagation velocity of short-interval SSEs is in proportion to the slip velocity. The results also show that during the nucleation process of large earthquakes, the occurrence of short-interval SSEs becomes irregular because of the accelerated slips that occur at the bottom of the seismogenic zone. Our results suggest that monitoring of short-interval SSEs might be useful for forecasting the main earthquakes.     

Models of Earth's main magnetic field incorporating flux and radial vorticity constraints

We describe a new technique for implementing the constraints on magnetic fields arising from two hypotheses about the fluid core of the Earth, namely the frozen-flux hypothesis and the hypothesis that the core is in magnetostrophic force balance with negligible leakage of current into the mantle. These hypotheses lead to time-independence of the integrated flux through certain 'null-flux patches' on the core surface, and to time-independence of their radial vorticity. Although the frozen-flux hypothesis has received attention before, constraining the radial vorticity has not previously been attempted. We describe a parametrization and an algorithm for preserving topology of radial magnetic fields at the core surface while allowing morphological changes. The parametrization is a spherical triangle tesselation of the core surface. Topology with respect to a reference model (based on data from the Oersted satellite) is preserved as models at different epochs are perturbed to optimize the fit to the data; the topology preservation is achieved by the imposition of inequality constraints on the model, and the optimization at each iteration is cast as a bounded value least-squares problem. For epochs 2000, 1980, 1945, 1915 and 1882 we are able to produce models of the core field which are consistent with flux and radial vorticity conservation, thus providing no observational evidence for the failure of the underlying assumptions. These models are a step towards the production of models which are optimal for the retrieval of frozen-flux velocity fields at the core surface.