Estimation of the inland extending length of the freshwater–saltwater interface in coastal unconfined aquifers

ABSTRACT

The inland extending length of the freshwater–saltwater interface toe is useful in studies of seawater intrusion in coastal areas. The submarine fresh groundwater discharge in coastal zones is affected not only by hydraulic conductivity and hydraulic gradient of the aquifer, but also by the position of the interface. Two observation wells at different distances from the coast are required to calculate the fresh groundwater flow rate in coastal unconfined aquifers. By considering that the submarine groundwater discharge is equal to the groundwater flow rate, the length of the interface toe extending inland can be estimated when the groundwater flow is at a steady-flow state. Aquifers with horizontal and sloping confined beds and without/with unique surface vertical infiltration are considered. Examples used to illustrate the application of these methods indicate that the inland extending lengths of the interface toe in aquifers with vertical surface infiltration are much shorter than those in aquifers without vertical surface infiltration, and the length of the interface in aquifers with a horizontal confining lower bed are smaller than those in aquifers with a confining lower bed sloping towards the sea. The extent of the interface on the northwestern coast near the city of Beihai in southern Guangxi, China, on 18 January 2013 was estimated as 471–478 m.

Editor M.C. Acreman Associate editor not assigned     

EM induction in the Vancouver Island region: 3D numerical,analogue model,and field site results

Electromagnetic induction in the Vancouver Island region for a uniform inducing source field for 300 s period is investigated with the aid of three-dimensional (3-D) numerical and analogue model results and field site measurements. The thin sheet numerical model, based on the subducting Juan de Fuca plate analogue model ofDosso et al., consists of a 5km thick non-uniform thin sheet (comprising the lateral conductivity contrasts arising from the land, the varying depth ocean, and the sediment) underlain by a four-layer conductive structure. The four-layer conductive structure beneath the non-uniform thin sheet simulates the effect of the Juan de Fuca plate subducting Vancouver Island. To examine the effects of the ocean channel depth between Vancouver Island and the British Columbia (Canada) mainland, numerical results were obtained for two channel depths (0 and 600 m). The results indicate that the channel plays an important role in the geomagnetic response in the central and inner coastal regions of Vancouver Island. The general agreement of the 3-D numerical model induction arrows with the analogue model and field site induction arrows for 300 s supports the premise of a layered conductive substructure dipping at a small angle, at most, beneath Vancouver Island.Lithoprobe Publication No. 311.     

A laboratory analogue model study of ocean-wave induced magnetic fields for cases of non-uniform depths and sea-land interfaces

The laboratory analogue model described by Miles et al. (1977) for the case of a uniform depth ocean is extended to study ocean-wave induced magnetic fields for cases of non-uniform ocean depths and sea-land interfaces. The models studied are: (i) the step and shelf model, (ii) the wedge and shelf model, (iii) the dyke model, (iv) the sea mount model, (v) the sloping bay and shelf model, and (vi) the reef and shelf model. For shallow depths, the behaviour of the induced magnetic field is strongly dependent on the fluid depth over the structure. The shape of the leading edge of the submerged structure, and the contour of the coastline also significantly affect the induced horizontal magnetic field, and possible mechanisms for these effects include wave interference and current channelling. On the basis of the model results, it is apparent that an irregular ocean bottom, or a large submerged structure, should affect the induced magnetic field at the surface only for rather shallow ocean depths of less than 40 m for a fluid wavelength of 360 m.     

大地电导率横向突变处磁暴感应地电场的邻近效应

大地电性结构的横向变化会对磁暴时的感应地电流和地面电磁场产生影响.本文假设扰动地磁场变化的源为地面以上一定高度的面电流,以某一典型层状大地电导率结构为基础,构造含有电导率横向突变的地电模型.针对感应电流的方向与横向分界面平行的情形,采用伽辽金有限元法对电导率横向突变处的感应地电场进行了分析,揭示了电导率横向差异产生的趋肤效应和邻近效应的机理,针对与电性结构分界面平行的输电线路,从评估地磁感应电流的角度讨论了影响的严重程度和范围.     

Electromagnetic analogue model measurements and finite-difference numerical calculations of the response of a conducting slab to three different source fields

This work compares experimental analogue model measurements and finite-difference numerical calculations of the electric and magnetic fields for a highly conducting slab embedded in a poorly conducting host earth for three different source field configurations. Measurements and calculations were carried out for a uniform source, a sheet current source with a y exp(?ay) current intensity distribution, and a horizontal magnetic dipole source. The results indicate reasonable agreement with some exceptions between the analogue and numerical methods. The source field is found to have an important effect on the field anomalies at the interface of the highly conducting slab and the poorly conducting host medium.     

Electromagnetic response of a buried cylindrical conductor for sheet and line current sources

A major problem in electromagnetic induction studies in regions of localized source fields, such as the auroral and equatorial electrojet regions, is the source effect. Using an analytical model, the electromagnetic response of a buried conducting cylinder to sheet current and line current excitations has been studied for the period rangeT=5 s to 24 h. The validity of the numerical results obtained from the analytical model are compared with the numerical results obtained from a finite difference model. The results show that for periods less than 30 min, there is no significant difference in the response of the cylinder to both source fields. However, significant differences are observed at longer periods. It was also observed that the equivalent height at which a uniform sheet current at 100 km above the earth's surface can be approximated by a line current varies as a function of the source period.     

利用气枪地震资料对福建及台湾海峡南部地壳三维P波速度结构研究

根据福建及台湾海峡南部海陆联测试验记录到的Pg和Pm震相走时数据,利用速度与界面联合成像方法构建地壳三维P波速度结构,揭示了该区地壳深部构造特征.结果表明：福建和台湾海峡海陆过渡带以及海峡南部地壳速度结构存在明显的不均匀性,滨海断裂两侧速度结构复杂,随深度呈现明显的分段特征,其上地壳海陆过渡带呈高速特性,台湾海峡呈低速特性;下地壳海陆过渡带呈低速特性,台湾海峡呈高速特性;研究区莫霍面的深度约为28~33 km,存在较明显差异,闽粤交接部位存在明显的地壳厚度减薄,莫霍面深度接近28 km,这与正常型华南活动地块与减薄型南海活动地块交汇致使地壳厚度减薄有关,体现了活动块体边界构造特征.历史大震主要发生在高低速异常过渡带且有深大断裂穿过的区域,现今中小震主要分布于闽粤海陆过渡带,这一特征可能与此地广泛发育的断层和华南与南海活动地块相互作用有关.     

A STUDY OF MAGNETOTELLURIC STATIC DISTORTION IN THE CONTEXT OF INTRUSIVE VOLCANISM1

MT data from a detailed 84-site grid array (28 × 12 km) in the Paraná Basin (Brazil) reveal a high degree of frequency-independent and parallel amplitude distortion. The Palaeozoic sediments across the survey area are covered by 1.2 km of flood basalts. A deep well-log provides some control regarding the emplacement of thin diabase sills but information regarding vertical feeder dikes is non-existent. The degree of parallel behaviour is identified using anisotropy ratios which quantify the extent and bandwidth of the distortion characteristics. A 2D modelling study is carried out using the concept of a horizontally layered (1D) basin with superimposed, small-scale inhomogeneities. Two likely geological distortion structures are considered. The first represents at- or near-surface (0 to 100 m) inhomogeneities. The second represents thin but vertically elongate dikes. The two distortion effects are found to produce different and characteristic behaviour in the H-polarization (TM) mode model data. The modelled E-polarization (TE) mode data are undistorted and provide control. The model developed to represent vertical dikes gives rise to a characteristic frequency dependence which is distinct from that of the undistorted response. This model also generates subparallel distortion effects which give rise to a form of variance in the response estimates as a function of location. The behaviour occurs, importantly, in both the amplitude and phase of the distorted data. A comparison of the normalized amplitude data, observed and modelled, indicates a high degree of correspondence down to a low frequency limit (< 0.1 Hz) where additional and deep contributions to the observed response become apparent. The modelling results and data characteristics indicate that the survey area is underlain by a series of closely parallel, thin intrusive dikes (a dike swarm).     

Two-dimensional scattering of SH waves in a soil layer underlain with a sloping bedrock

A two-dimensional analysis is applied to examine the effect that a sloping bedrock half-space has on the amplification of an anti-plane shear wave. The direct boundary integral equation method is used for the two-dimensional analysis. The particular soil–rock configuration investigated includes a homogeneous soil layer underlain by a sloping rock half-space. The rock half-space dips for a horizontal distance L and then becomes horizontal so that the overlying soil layer has a thickness H that remains constant from this point to infinity. The materials in the soil–rock configuration are considered viscoelastic except in the rock half-space below soil layer thickness H, which is considered elastic. This limitation in damping is due to the correction used for the truncation of the half-space boundary. Four cases are used to study the relationship between rock slope and surface displacement, vertical, 1:2, 1:4, 1:8. Surface displacements are determined for each of these cases for half-space incidence angles of 90, 75, and 60°. To allow for applicability to a wide range of problems, results are determined as a function of dimensionless parameters. In addition, solutions from a one-dimensional analysis are compared with the results of the two-dimensional analysis to establish limits outside of which a one-dimensional analysis suffices.     

Effect of topography and subsurface inhomogeneity on seismic rayleigh waves

The effect of topography and subsurface inhomogeneity on surface motion is investigated in the case of Rayleigh waves. In the previous paper, the same effect was investigated in the case of SV waves. Several types of topography, such as cliffs both with and without a soft layer at the foot of the slope, are considered. Computations are made using a new hybrid method combining a particle model with a finite element method. In cases of harmonic Rayleigh waves, surface motions with amplitudes as large as 1.5 to 5 times the horizontal surface displacement of the incident Rayleigh waves are produced near the slope and the sloping interface. When a Rayleigh wave propagating through a hard single-layered ground encounters a sloping interface where hard ground and soft ground make contact with each other, Rayleigh waves having two different, phase velocities are produced and they correspond to the fundamental mode, and the first mode determined by Haskell's method. In addition, the transient response when Rayleigh waves propagate through the cliff is also simulated. Assuming the vertical component of the Tokachi-oki Earthquake (1968) measured on the surface to be a Rayleigh wave, the incident Rayleigh wave can be obtained by a Fourier synthesis of eigenfunctions of Rayleigh waves.     

A laboratory electromagnetic model study of the Juan de Fuca Plate region

The behaviour of the magnetic field variations over the Juan de Fuca Plate region is studied using a scaled laboratory analogue model. The model includes a simulation of the complex Juan de Fuca Plate subducting the Vancouver Island region. The subducting plate is modelled with a profile of increasing inclination from east to west; horizontal offshore, dipping at 10° under Vancouver Island, and bending further under Georgia Strait to subduct the continent at 30° for the B.C. region and 45° for the Washington-Oregon region. The strike of the bending plate follows the general strike of the continental coastline with an abrupt change in direction (42°) in the Puget Sound area. The model substructure simulates a subducting plate, overplated by a sediment layer several kilometres thick, and underlain by a 30 km thick highly conducting upper asthenosphere. The model source frequencies used simulate periods 5–120 min in the geophysical scale. In-phase and quadrature H_x, H_y, and H_z magnetic field measurements for the modelled region are presented for an approximately uniform overhead horizontal source field for E- and H-polarizations (electric field of the source approximately parallel and perpendicular, respectively, to the west coast of Vancouver Island). The fields for three regions of the model; over Vancouver Island, over the Olympic Peninsula and over a linear portion of the U.S. coastline, are examined in detail. The general conclusion is that the effect of the dipping subducting plates is to significantly attenuate, at short periods, the maxima in the anomalies at the coastlines underlain by the 10° dipping plate, while leading to anomalous vertical and horizontal fields over ranges as large as 500 km inland over a wide period range. Anomalous fields are observed over the offshore and inland knee-bends of the subducting plates at all periods for both E- and H-polarizations. For locations on land, the in-phase induction arrows point seaward and perpendicular to the strikes of the dipping plates for all periods, while the quadrature arrows at short periods point landward and rotate to point seaward for periods greater than 20 min.     

Magnetotelluric response of a two-dimensional sloping contact by the finite element method

Summary The finite element method, with triangular elements, is used to study the effect of a two-dimensional sloping contact on the surface electromagnetic fields. It is found in the case ofH-polarization and small slopes that the electric field and the apparent resistivity near the contact, on the conductive side, are higher than their asymptotic values. In the case ofE-polarization the apparent resistivity and phase values on the conductive side fall off less rapidly to their asymptotic values with decreasing slope resulting in higher apparent resistivity and phase values on the conductive side, than those expected for a vertical contact. The peak in the amplitude and phase of the normalized vertical magnetic field shifts from the resistive side for a vertical contact to the conductive side for a sloping contact. Far from the sloping contact, on the conductive side, higher values are observed for the normalized vertical magnetic field than in the case of a vertical contact.     

Evaluation of the glacial isostatic adjustment (GIA) models for Antarctica based on GPS vertical velocities

Due to the scarcity of data, modeling the glacial isostatic adjustment(GIA) for Antarctica is more difficult than it is for the ancient ice sheet area in North America and Northern Europe. Large uncertainties are observed in existing GIA models for Antarctica. Modern space-based geodetic measurements provide checks and constraints for GIA models. The present-day uplift velocities of global positioning system(GPS) stations at 73 stations in Antarctica and adjacent regions from 1996 to 2014 have been estimated using GAMIT/GLOBK version 10.5 with a colored noise model. To easily analyze the effect of difference sources on the vertical velocities, and for easy comparison with both GIA model predictions and GPS results from Argus et al.(2014) and Thomas et al.(2011), seven sub-regions are divided. They are the northern Antarctic Peninsula, the Filchner-Ronne Ice Shelf, the Amundsen Sea coast, the Ross Ice Shelf, Mount Erebus, inland Southwest Antarctica and the East Antarctic coast,respectively. The results show that the fast uplift in the north Antarctic Peninsula and Pine Island Bay regions may be caused by the elastic response to snow and ice mass loss. The fast subsidence near Mount Erebus may be related to the activity of a magma body. The uplift or subsidence near the East Antarctic coast is very slow while the uplift for the rest regions is mainly caused by GIA. By analyzing the correlation and the associated weighted root mean square(WRMS) between the GIA predictions and the GPS velocities, we found that the ICE-6G_C(VM5a) model and the Geruo 13 model show the most consistency with our GPS results, while the W12a and IJ05_R2 series models show poor consistency with our GPS results. Although improved greatly in recent years, the GIA modeling in Antarctica still lags behind the modeling of the North American. Some GPS stations, for example the Bennett Nunatak station(BENN), have observed large discrepancies between GIA predictions and GPS velocities.Because of the large uncertainties in calculating elastic responses due to the significant variations of ice and snow loads, the GPS velocities still cannot be used as a precise constraint on GIA models.     

Fast Imaging of TDEM data based on S-inversion

Fast S-inversion is a method of interpretation of time-domain electromagnetic (TDEM) sounding data using the thin sheet model approach. Within the framework of this method, the electromagnetic (EM) response measured at the surface of the earth at every time delay is matched with that of a thin sheet model. The conductivity change with depth is obtained using the conductance, S, and depth, d, of the equivalent thin sheet. We analyze two different numerical techniques, the differential S-transformation and the regularized S-inversion, to determine the parameters of the thin sheet. The first technique is a direct differential transformation of the observed data into conductance and depth values. It is fast and requires no iterations or starting model. The second technique uses a regularized inversion scheme to fit the measured response with that of a thin sheet. In both techniques, the retrieved conductance values are differentiated with respect to depth to obtain the conductivity change with depth. We apply S-inversion to three-dimensional synthetic data and we successfully locate the local conductors. We also demonstrate a case history by interpreting TDEM data obtained at the Nojima fault zone in Japan. The results clearly indicate the location of the fault zone.     

Vertical and horizontal components of the lake Baikal electrotelluric field and their relation to the electric conductivity

Synchronous observations of the electrotelluric field were performed on the lake ice and the coast. The ice-based vertical/horizontal field ratio at periods of a few seconds to a few minutes ranges from a few hundredths to a few thousandths, which is comparable to the noise level. The results of these observations favor a plane wave model. At the coastal station, the vertical and horizontal fields are similar, which is related to geoelectric heterogeneity of the medium. An electric tipper is used to control the behavior of heterogeneity in the spectral-time domain. Anomalous changes associated with earthquakes of K ≥ 12 are revealed from data of the tipper monitoring at epicentral distances of up to 300 km. A longitudinal MTS curve free from the effect of local geoelectric heterogeneities is obtained from ice-based data. Deep conducting layers are recognized as minimums in this curve.     

Physiographic and geological characteristics of shelves in north and west of Taiwan

In north of Taiwan the East China Sea Shelf extends from the coastline of mainland China to the southern Okinawa Trough, showing a shelf-slope-basin bathymetric profile. Sea floor of the East China Sea Shelf is generally smooth and flat. It is wide (230 km) and relatively shallow with an average shelf-break depth of 120 m, reflecting the effect of Quaternary glaciation about 15000 years ago. The Pliocene-Quaternary East China Sea Shelf is underlain by thin sequences of about 1000 m thick shallow marine sediments mainly derived from China, forming an eastward dipping sedimentary wedge due to regional tilting and subsidence and representing the youngest part of the Cenozoic rift margin in southeastern China. The Taiwan Strait Shelf has two phases in development: one is the early phase of the Paleocene-Miocene rift margin and the other is the late stage of the Pliocene-Quaternary foreland sedimentation. It is a foreland shelf and its present morphology results from the combined effects of tectonic subsidence and sedimentation overprinting that of the Late Pleistocene glaciation about 15000 years ago.     

台湾海峡大容量气枪震源海陆联测初探

本文利用在我国台湾海峡采用大容量气枪震源开展海陆联测获得的广角地震测线HX9, 采用二维射线追踪法反演得到了HX9剖面的地壳二维速度结构和地壳界面形态, 初步探明了福建—台湾海峡海陆过渡带的深部构造． 结果表明: HX9剖面的地壳内存在两个速度间断面, 即C界面和莫霍面, 其中: C界面为上、 下地壳的分界面, 是一个小的速度不连续面, 速度变化值达0.08—0.16 km/s; 而地壳底部的莫霍面则有较大的速度反差, 变化值达1.02—1.29 km/s, 莫霍面上、 下的速度分别为6.75—6.97 km/s和8.00—8.07 km/s． 沿剖面的地壳界面形态总体起伏不大, 陆域上、 下地壳的厚度和界面变化趋势均相似, 从陆域到海域呈微倾斜变化趋势, 表现为减薄陆壳的特征． 莫霍面陆域埋深约为31.6 km, 向福建东南沿海逐渐减薄至27.4 km左右．      

A FIRST-ORDER CORRECTION TO THE THEORY FOR THE ELECTROMAGNETIC RESPONSE OF A THIN CONDUCTING SHEET1

The eddy currents induced in a thin conducting sheet by a sinusoidally varying primary magnetic field are investigated in the low-frequency limit when the depth of penetration of the primary field is much greater than the thickness of the sheet. A perturbation scheme in terms of the small parameter, δ, which is the ratio of the thickness of the sheet to the length scale of the primary field is initiated. The zero-order terms in δ give the familiar results that were first obtained by Maxwell. The first-order terms in δ depend on the tangential as well as the normal component of the primary field and give eddy-current distributions that vary quadratically across the sheet. The boundary conditions that determine the accompanying secondary magnetic field are derived. Detailed results are given when the primary field is that due to a dipole.     

Barotropic instability of a boundary jet on a sloping bottom

Abstract

The stability of a shear flow on a sloping bottom in a homogeneous, rotating system was investigated by means of a laboratory experiment.

The basic flow was driven near a vertical wall of a circular container by a ring-shaped plate that contacted with a free surface of the working fluid and rotated relative to the fluid container. The velocity profile was asymmetric in the radial direction and had only one inflection point. The velocity profile was well expressed by a linear theory for the vertical shear layer.

The effect of the circular geometry was checked by comparing experimental results obtained in two fluid systems in which only the sign of the curvature was opposite and it was confirmed that circular geometry was not essential for the shear flow on the sloping bottom in this experiment.

It was found that the sloping bottom stabilizes the basic flow only when the drift direction of the topographic Rossby wave is opposite to that of the basic flow. The viscous dissipation in both the Ekman layer and the interior region was also important in determining the critical Rossby number.

The eddy fields caused by the instability can be classified into two types: One is the stationary eddy field in which a row of eddies moves along the basic flow without changing form. The other is the flow pattern in which eddies have finite life times and their configuration is not well organized. When the sloping bottom does not stabilize the basic flow, the former flow pattern is realized, otherwise the latter flow pattern appears.

The wave numbers of the eddies in the regular flow pattern were observed as a function of the Rossby number. The relation did not fit to linear preferred modes predicted by an eigenvalue problem.     

Analogue model measurements for a horizontal magnetic dipole embedded within a conducting medium

The validity of an analogue model method employing a horizontal magnetic dipole source situated within a conducting layer for the cases of a poorly conducting model earth and a highly conducting model ocean is studied by comparing model magnetic field measurements with theoretical calculations. The model is then used to study one example for each case; the response of a conducting cylinder simulating an ore body embedded in the earth, and the response of a conducting wedge simulating a shelving ocean.