湖南印支期过铝质花岗岩的形成: 岩浆底侵与地壳加厚热效应的数值模拟

在合理构建华南印支期地质-物理模型的基础上, 利用FLAC软件, 模拟了该区印支期过铝质富钾花岗岩形成与基性岩浆底侵, 及陆壳变形叠置加厚两种动力学背景的可能联系. 模拟结果表明, 220 Ma±的基性岩浆底侵能导致地壳含水矿物相岩石的深熔, 但除非印支期存在大规模基性岩浆的底侵作用, 否则其热效应持续时间和热效应波及范围难以形成具大岩基规模的湖南印支期花岗岩. 在陆壳叠置加厚模型中, 地壳的叠置加厚可导致中下地壳界面温度升高到700℃以上, 引起片麻质岩石熔融, 当加厚因子达1.3, 白云母矿物脱水熔融产生的熔体达到熔体流动临界比例(≥20%), 从而形成花岗岩基; 结合印支期挤压逆冲推覆构造和同期基性火山岩极少出露的地质事实, 认为陆壳变形加厚可能是湖南印支期构造岩浆作用形成的主导因素.     

Improved NMO correction with a specific application to shallow-seismic data

The application of traditional NMO‐correction techniques in the processing of seismic data may result in severe distortion. This distortion is observed as a decrease in frequency content (NMO stretch) or even a total loss of data for steep velocity gradients. It is shown that the specific nature of the observed distortion is introduced by the particular implementation of the traditional NMO‐correction technique. It is also shown that other techniques that have been proposed in the literature suffer from similar artefacts. An alternative approach is suggested, based on the correction of tapered blocks of seismic data, followed by a coherence filter to compensate for the specific artefacts thus introduced. Correction of seismic amplitudes (e.g. geometrical spreading and attenuation) is implemented as an integral part of the NMO‐correction method thus allowing for both dynamic and kinematic reconstruction of interfering events. From the application of this method to synthetic and field data, it is concluded that the proposed technique may be particularly useful in the processing of shallow seismic (shear‐wave) data. Although not illustrated in this paper, it is emphasized that application of the proposed method is by no means restricted to shear waves; application to P‐wave data will prove equally useful, although the stretch effects involved may not be as severe as in the examples given in this paper.     

Analysis of the resolution function in seismic prestack depth imaging

We consider the problem of estimating subsurface quantities such as velocity or reflectivity from seismic measurements. Because of a limited aperture and band-limited signals, the output from a seismic prestack reconstruction method is a distorted or blurred image. This distortion can be computed using the concept of resolution function, which is a quantity readily accessible in the Fourier space of the model. The key parameter is the scattering wavenumber, which at a particular image point is defined by the incident and scattered ray directions in a given background model. Any location in any background model can be considered. In general, the resolution function will depend on the following four quantities: the background velocity model, the frequency bandwidth, the wavefield type and the acquisition geometry.
We first establish the resolution function for a general scattering model assuming local reaction. We then adapt this result for two well-known scattering models: Born and Kirchhoff. For each of these approximations the corresponding resolution function is derived and discussed. Finally, by employing a simple synthetic data example we demonstrate the ability of the resolution function to predict the image distortions.     

Wave velocities and attenuation of shaley sandstones as a function of pore pressure and partial saturation

We obtain the wave velocities and quality factors of clay‐bearing sandstones as a function of pore pressure, frequency and partial saturation. The model is based on a Biot‐type three‐phase theory that considers the coexistence of two solids (sand grains and clay particles) and a fluid mixture. Additional attenuation is described with the constant‐Q model and viscodynamic functions to model the high‐frequency behaviour. We apply a uniform gas/fluid mixing law that satisfies the Wood and Voigt averages at low and high frequencies, respectively. Pressure effects are accounted for by using an effective stress law. By fitting a permeability model of the Kozeny– Carman type to core data, the model is able to predict wave velocity and attenuation from seismic to ultrasonic frequencies, including the effects of partial saturation. Testing of the model with laboratory data shows good agreement between predictions and measurements.     

Transient foundation solution of saturated soil to impulsive concentrated loading

The transient dynamic response of saturated soil under suddenly applied normal and horizontal concentrated loading is studied in this paper. The behavior of saturated soil is governed by Biot's consolidation theory. The general solutions for Biot equations of equilibrium are derived in terms of displacements and variations of fluid volume, using Laplace–Hankel integral transforms. The solutions in the time domain can be evaluated by numerical inverse Laplace–Hankel transforms. Selected numerical results for displacements, stresses, and pore pressures are presented. Comparisons with existing closed-form solutions for the elastic half-space are made to confirm the accuracy of the present solutions. The solutions can be used to study a variety of transient wave propagation problems and dynamical interactions between saturated soil and structures.     

The El Salvador earthquakes of January and February 2001: context, characteristics and implications for seismic risk

The small Central American republic of El Salvador has experienced, on average, one destructive earthquake per decade during the last hundred years. The latest events occurred on 13 January and 13 February 2001, with magnitudes M_w 7.7 and 6.6, respectively. The two events, which were of different tectonic origin, follow the patterns of the seismicity of the region although neither event has a known precedent in the earthquake catalogue in terms of size and location. The earthquakes caused damage to thousands of traditionally built houses and triggered hundreds of landslides, which were the main causes of fatalities. The earthquakes have clearly demonstrated trends of increasing seismic risk in El Salvador due to rapid population expansion in areas of high shaking and landslide hazard, exacerbated by deforestation and uncontrolled urbanisation. The institutional mechanisms required for the control of land use and building practice are very weak and present a major obstacle to risk mitigation.     

Fractal properties of medium and seismoelectric phenomena

Electrokinetic phenomena in a water-porous medium with a fractal structure above percolation threshold are theoretically investigated. Fracture zone with space-variable porosity is considered as a model of an earthquake hypocenter zone in which the electrokinetic current results from fluid filtration in a fractal pore network. A critical exponent of the streaming potential coefficient is found to depend on both the transport critical exponent and correlation length critical exponent. In this model, logarithmic dependence of electric field amplitude E on the earthquake magnitude M is derived which is compatible with the one observed by the VAN group. Without fractal properties, this form of dependence contradicts the empirical data. The electromagnetic field far from the hypocenter is calculated, which leads to the prediction of weak magnetic field variations. To explain the observed amplitude of VAN's Seismic Electric Signals (SES), the electric source must be at a distance of about 10 km from the registration point if the medium is homogeneous. Therefore, some conductive channel(s) are needed to explain the long distance selective SES transmission.     

Paleolatitudes of the Kerguelen hotspot: new paleomagnetic results and dynamic modeling

The Kerguelen Plateau, a Large Igneous Province in the southern Indian Ocean, was formed as a product of the Kerguelen hotspot in several eruptive phases during the last 120 Myr. We obtained new paleolatitudes for the central and northern Kerguelen Plateau from paleomagnetic investigations on basalts, which were drilled during ODP Leg 183 to the Kerguelen Plateau-Broken Ridge. The paleolatitudes coincide with paleolatitudes from previous investigations at the Kerguelen Plateau and Ninetyeast Ridge (the track of the Kerguelen hotspot) and indicate a difference between paleolatitudes and present position at 49°S of the Kerguelen hotspot. We show that true polar wander, the global motion between the mantle and the rotation axis, cannot explain this difference in latitudes. We present numerical model results of plume conduit motion in a large-scale mantle flow and the resulting surface hotspot motion. A large number of models all predict southward motion between 3° and 10° for the Kerguelen hotspot during the last 100 Myr, which is consistent with our paleomagnetic results.     

Development and testing of numerical methods for two-way nested air pollution modelling

A new 3-D model REGINA (REGIonal high resolutioN Air pollution model) is under development at the National Environmental Research Institute (NERI). The model is based on pieces from several models developed over the last decades at NERI. The aim of the work is to develop a nested model which can operate with very high resolution in both space and time. To fulfill this aim the choice and implementation of accurate numerical methods is crucial. The model will be applied for studying air pollution phenomena (both monitoring, forecasting and scenarios) over Denmark. The present paper is focussed on the modification, implementation and testing of a numerical method for treating the horizontal advection in the model as well as the implementation of two-way nesting techniques. The horizontal transport in the model is solved using an accurate space derivative algorithm. This method is traditionally implemented with periodic boundary conditions, however, this is not an option for nested modelling. A new method for calculating non-periodic boundary conditions has been developed in order to overcome this problem. Extensive testing of the numerical solution of the advection and the coupling of the solution of advection and chemistry in the model using Molenkamp–Crowley rotation tests have been carried out. The results show that the model with the current implementation of numerical methods is suitable for calculating air pollution levels with high resolution.     

Spatial and temporal variability of pore‐water pressures in residual soil slopes in a tropical climate

It is critical to understand and quantify the temporal and spatial variability in hillslope hydrological data in order to advance hillslope hydrological studies, evaluate distributed parameter hydrological models, analyse variability in hydrological response of slopes and design efficient field data sampling networks. The spatial and temporal variability of field‐measured pore‐water pressures in three residual soil slopes in Singapore was investigated using geostatistical methods. Parameters of the semivariograms, namely the range, sill and nugget effect, revealed interesting insights into the spatial structure of the temporal situation of pore‐water pressures in the slopes. While informative, mean estimates have been shown to be inadequate for modelling purposes, indicator semivariograms together with mean prediction by kriging provide a better form of model input. Results also indicate that significant temporal and spatial variability in pore‐water pressures exists in the slope profile and thereby induces variability in hydrological response of the slope. Spatial and temporal variability in pore‐water pressure decreases with increasing soil depth. The variability decreases during wet conditions as the slope approaches near saturation and the variability increases with high matric suction development following rainfall periods. Variability in pore‐water pressures is greatest at shallow depths and near the slope crest and is strongly influenced by the combined action of microclimate, vegetation and soil properties. Copyright © 2002 John Wiley & Sons, Ltd.