Sundararajan transform—an application to geophysical data analysis

An integral transform, called in this study as Sundararajan transform, has been used to estimate the parameters of vertical magnetic effect of a fault structure. It differs from the well-known Hilbert transform in the property of phase shift as it yields a phase shift of 270° unlike the Hilbert transform, which is a 90° phase shifter. Other properties of the Sundararajan transform remain almost the same as the Hilbert transform. The transform has been tested on synthetic data and a field example of Lachlan Foldbelts, New South Wales, Australia. The results of this technique agree with the one published in the literature. The noise analysis has been examined and showed that it still provides acceptable results. The application of this transform to geophysical interpretation illustrates its potentiality, and it may be widely applied in various disciplines mainly in the field of communication engineering, signal, and image processing. An interesting property of this transform is that two successive transforms of a function return it to its original form unlike the Hilbert transform, which returns it to the negative of the original form. The procedure discussed may be automated.     

Potential field transformation on the basis of a continuous wavelet transform

The possibility of applying the continuous wavelet analysis on the basis of a Poissonian core for the processing and interpretation of gravity and magnetic survey data is considered. Based on the wavelet spectrum, W(x,h), reconstruction of the initial signal, g(x), is possible. Additionally, based on this wavelet spectrum, the equivalent distributions of mass and magnetization that produce the initial field, g(x), can be reconstructed. These possibilities allow the wavelet transform to be used for solving such classical problems as the filtering of initial signals, continuation of the field, calculation of higher field derivatives in the upper and lower half-spaces, and reduction of magnetic anomalies to the pole.     

Interpretation of gravity anomalies due to some two-dimensional structures: A Hilbert transform technique

An interpretation procedure is formulated using the Hilbert transform for analysing the gravity effect of (a) two-dimensional horizontal circular cylinder, (b) semi-infinite thin fault block, and (c) a geologic contact. In all three cases the abscissa of the point of intersection of the gravity anomaly curve or the horizontal derivative curve and its Hilbert transform yields directly the depth of the body. The proposed method is tested on theoretical models. Also, the method is applied to the field data of a geologic contact taken over the Himalayan foothills across the Brahmaputra valley.     

Analysis of gravity gradients over a thin infinite sheet

A simple method to interpret gravity gradients over a thin infinite dipping sheet is discussed. The Hilbert transform is used to compute the vertical gradient from the horizontal gradient of the gravity field. The method is illustrated with a theoretical example.     

Geodynamics and underlying bedrock of the magnetically active crust layer of the Lut block,Eastern Iran

The Curie point depth map of Eastern Iran was constituted from spectral analysis of the aeromagnetic data. The reduction to pole (RTP) was applied to the magnetic anomaly data. The Curie point depth values from 165 overlapping blocks, 100 × 100 km in size, have been estimated. The Curie point depth method provides a relationship between the 2-D FFT power spectrum of the magnetic anomalies and the depth of magnetic sources by transforming the spatial data into the frequency domain. The centroid and top depth of the magnetic sources (respectively Z₀ and Z_t) is calculated from radially averaged log power spectrum for each block. Finally, the Curie point depth of Eastern Iran is obtained by Z_b = 2Z₀–Z_t. The highest value of 24 km is located in eastern and western boundaries of the Lut block, and the lowest value of 12 km is located at north of study area. The shallow depths in the Curie-point depth map are well correlated with the young volcanic areas and geothermal potential fields. Geothermal gradient ranging from 24 to 45°C/km. The deduced thermal structure in eastern Iran has a relationship with orogenic collapse associated with delamination of thickened lithospheric root between the Lut and Afghan continental blocks.     

规则体重力异常余弦变换谱特征及反演

提出利用余弦变换,研究规则体重力异常谱特征及反演方法。推导出球体、无限长水平圆柱体、无限和有限延伸垂直台阶的重力异常余弦变换谱及反演表达式。球体与无限长水平圆柱体的余弦变换谱相似,无限延伸垂直台阶的变化很尖锐,有限台阶的表现为上、下有规律的跳跃。在模型反演算例中,选择了适当的采样间隔和“趋势镶边”法,反演的圆柱体线密度和轴心深度误差分别为0.25%和3.8%。     

A condition for the existence of quasi-periodic nonlinear internal waves in the ocean shelf zone

Free internal gravity waves, which are typical of the shelf ocean zone are studied. A necessary condition for the existence of nonlinear wave disturbances quasi-periodic in time in a continuously stratified sea with variable depth H is found in the quasistatic and the “hard cap” approximation with respect to dissipative factors and latitudinal variation in the Coriolis parameter. The obtained assessment is equivalent to the condition obtained within the framework of linear theory for the case of the f-plane and a constant depth H.     

Estimation of Pullout Capacity of Vertical Plate Anchors in Cohesionless Soil Using MARS

Vertical plate anchors provide an economical solution to safely resist the large horizontal forces experienced by the foundation of different structures such as bulkheads, sheet piles, retaining walls and so forth. This paper develops a multivariate adaptive regression spline (MARS) model-based approach for the determination of horizontal pullout capacity (P _u ) of vertical plate anchors buried in cohesionless soil by utilizing experimental results reported by different researchers. Based on the collection of forty different pullout experimental test results reported in the literature for anchors buried in loose to dense cohesionless soil with an embedment ratio ranges from 1 to 5, a predictive approach for P _u of vertical plate anchors has been developed in terms of non-dimensional pullout coefficient (M _γq ). The capability of the proposed MARS model for estimating the values of M _γq is examined by comparing the results obtained in the present study with those methods available in the literature. Using different statistical error measure criteria, this study indicates that the present approach is efficient in estimating the horizontal pullout capacity of vertical plate anchors as compared to other methods. The sensitivity analysis indicates that the embedment ratio (H/h, where H = embedment depth of anchor, and h = height of anchor) and internal friction angle (?) of soil mass are the two most important parameters for the evaluation of non-dimensional pullout coefficient (M _γq ) using the proposed MARS model.     

Common features of the tectonosphere deep structure in the Western Pacific Margins (Northeast Asia Region and Australia)

This paper analyses the spatial distributions of gravity sources and density contrast of geological media between the centers of mass of density heterogeneities and surfaces of some layers, which is reflected by the values of parameter μ_z, in the Northeast Asia Region and Australia. Statistical images of complex distributions of density heterogeneities are generally consistent with the shallowest tectonic structures, wave velocity, and electric anomalies. The correlation of μ_z minimums with high heat flow suggests melted magmas beneath the crust base and in the asthenosphere. In the lower crust layer, large μ_z maximums correlate with major continental megaelements, which are divided by extended linear μ_z minimums, high conductivity anomalies, and Vp/Vs maximums marking zones of lower viscosity. Deep boundaries and a thickness of lithospheric plates and asthenosphere lenses have been defined in the Northeast Asia Region. In Australia, two roughly NS-striking low-viscosity zones crossing the continent reflect the deep setting of the Lasseter (in the west) and the Tasman (in the east) lines. These zones are displayed in a wide depth range (15–120 km) and are confined to the Archean–Proterozoic and Proterozoic–Phanerozoic boundaries in the Australian lithosphere. Statistical gravity models derived without external (relatively to gravity) geological and geophysical information prove the multilayered tectonosphere structure of both continents and the similar deep structure features of their margins.     

Two-dimensional interpretation of gravity anomalies over sedimentary basins with an exponential decrease of density contrast with depth

The decrease of density contrast with depth in sedimentary basins is approximated by an exponential function. The anomaly equation, in frequency domain, of a prismatic model with an exponential density function is derived. The method has been extended to derive the Fourier transforms of the gravity anomalies of the sedimentary basin, wherein the basin is viewed as vertical prisms placed in juxtaposition. The gravity anomalies of the sedimentary basin are obtained by taking the inverse Fourier transforms. Filon’s method has been extended for calculating accurate inverse Fourier transforms. The accuracy of the method has been tested using a synthetic example. A combination of space and frequency domain methods have been developed for inversion of gravity anomalies over the sedimentary basin. The method has been applied to interpret one synthetic profile and one field profile over the Godavari basin. The method developed in this paper to calculate the inverse Fourier transforms yields continuous spectrum with accurate values. The maximum depth deduced from the gravity anomalies is of the same order as the depth encountered to the basement at the Aswaraopeta borewell.     

Inversion of residual gravity anomalies using neural network

A new approach is presented in order to interpret residual gravity anomalies from simple geometrically shaped bodies such as horizontal cylinder, vertical cylinder, and sphere. This approach is mainly based on using modular neural network (MNN) inversion for estimating the shape factor, the depth, and the amplitude coefficient. The sigmoid function has been used as an activation function in the MNN inversion. The new approach has been tested first on synthetic data from different models using only one well-trained network. The results of this approach show that the parameter values estimated by the modular inversion are almost identical to the true parameters. Furthermore, the noise analysis has been examined where the results of the inversion produce satisfactory results up to 10% of white Gaussian noise. The reliability of this approach is demonstrated through two published real gravity field anomalies taken over a chromite deposit in Camaguey province, Cuba and over sulfide ore body, Nornada, Quebec, Canada. A comparable and acceptable agreement is obtained between the results derived by the MNN inversion method and those deduced by other interpretation methods. Furthermore, the depth obtained by the proposed technique is found to be very close to that obtained by drilling information.     

VSP wavefield separation using the gray-scale Hough transform: synthetic data

The principal objective in Vertical Seismic Profile (VSP) processing is the separation of the downgoing and upgoing wavefields. Several methods have been suggested in this area. This paper presents a new approach based on the gray-scale Hough transform (GSHT) which is an extension of the conventional Hough transform used to detect straight lines and other curves. The technique, we suggest here, directly maps the gray-scale VSP image, including the downgoing and upgoing linear events, in image coordinate space (x,t,g) to the gray Hough parameter counting space (θ,ρ,g), where θ and ρ are the polar parameters and g is the gray-scale value. In this new space, the downgoing events appear in the negative angles θ quadrant and the upgoing in the positive quadrant, owning to their opposite apparent velocities. The inverse GSHT algorithm, we developed in this study, is performed for extracting separately these two wavefields by considering the straight lines that satisfy the corresponding filtering conditions. The experimental results on synthetic VSP datasets are convincing. The wave separation is well performed, even in the presence of loud noise levels, with signal to noise ratio improvement and amplitude preservation, in contrast to median filtering.     

Identification of pre-seismic anomalies of soil radon-222 signal using Hilbert–Huang transform

Concentration of Rn-222 in soil has been monitored continuously at Ravangla in the Sikkim Himalayan Region of eastern India for about 7 months from October 2015 to May 2016 to detect earthquake-induced anomalies. The recorded data clearly show that various physical and meteorological parameters influence the soil radon concentration, leading to very complex soil Rn-222 time series. The components due to such external influences have been removed from the present time series, and Hilbert–Huang transform (HHT) applied for analysis of the data. Two radon anomalies caused due to earthquakes of magnitude M _b = 5.0 that occurred on 19 November 2015 and 5 April 2016 within an epicentral distance of 500 km from the monitoring station have been identified on the soil Rn-222 time series. These two precursory anomalies occurred 9 and 10 days, respectively, before the occurrence of the earthquakes. The absence of spurious signals or missing anomalies demonstrates that HHT is advantageous for analysis of nonlinear non-stationary data, and hence, it is a promising technique to analyse soil radon behaviour for predicting the possibility of occurrence of earthquakes.     

计算重力异常垂向二阶导数的DCT法——模型实验

基于离散余弦变换（DCT）的重力异常垂询二阶导数的计算方法是笔者提出的新方法。以无限长水平圆柱体为例研究了DCT法的计算规律：模型实验证实,该法计算的无限长水平圆柱体重力异常垂向二阶导数的精度与采样闻隔、剩余密度及圆柱体半径大小无关。与圆柱体轴心埋藏深度成正比。通过分析给出了基于DCT的合适的线性滤波方式,滤波后的重力异常垂向二阶导数与理论二阶导数拟合效果好,具有较高的计算精度。     

小波变换在位场资料去噪和位场分离中的应用

根据小波变换的原理,研究了小波阙值去噪方法和位场分离方法在某实测布格重力异常的应用及效果。多矩形棱柱体模型试验表明:选取合适的小波基和小波参数,可以更好地提取出重力异常中的有用信息,为下一步在实测重力数据中的应用提供理论基础。实测布格重力资料处理结果表明:基于小波变换的去噪方法和位场分离方法,能够较好地去除随机噪声,更好地区别出局部异常和区域异常,提高数据处理与解释中有用信号识别的精度。     

Peculiarities and variations in the optical spectrum of the post-AGB star V448 Lac = IRAS 22223+4327

Multi-epoch observations with high spectral resolution acquired in 1998–2008 are used to study the time behavior of the spectral-line profiles and velocity fields in the atmosphere and circumstellar shell of the post-AGB star V448 Lac. Asymmetry of the profiles of the strongest absorption lines with lower-level excitation potentials χ _low < 1 eV and time variations of these profiles have been detected, most prominently the profiles of the resonance lines of BaII, YII, LaII, SiII. The peculiarities of these profiles can be explained using a superposition of stellar absorption line and shell emission lines. Emission in the (0; 1) 5635 Å Swan system band of the C₂ molecule has been detected in the spectrum of V448 Lac for the first time. The core of the Hα line displays radial-velocity variations with an amplitude of ΔV _r ≈ 8 km/s. Radial-velocity variations displayed by weak metallic lines with lower amplitudes, ΔV _r ≈ 1–2 km/s, may be due to atmospheric pulsations. Differential line shifts, ΔV _r = 0–8 km/s have been detected on various dates. The position of the molecular spectrum is stationary in time, indicating a constant expansion velocity of the circumstellar shell, V _exp = 15.2 km/s, as derived from the C₂ and NaI lines.     

Geoblocks Recognition and Delineation of the Earthquake Prone Areas in the Tuan Giao Area (Northwest Vietnam)

The goal of the paper is to reveal the hierarchical block structure in the Tuan Giao area (Northwest Vietnam) and, on that base, to identify areas prone to earthquakes with M ≥ 5. Four large geoblocks of the second rank have been delineated on the basis of the joint analysis of geological, geophysical, geomorphic, and remote sensing data. The second rank geoblocks have been divided into smaller sub-blocks of the third rank. The recent geodynamics of the geoblocks have been characterized using geomorphic, seismological, gravity, and GPS data.The system of the delineated geoblocks is viewed as the Geodynamic Blocks model. The areas prone to earthquakes with M ≥ 5 have been identified using the pattern recognition algorithm CORA-3. The objects of the recognition were defined as circular areas, where boundaries of the geoblocks intersect each other. The recognition results confirm high seismic potential of the study region and provide information on potential earthquake sources for seismic hazard assessment: a number of boundary intersections have been identified as areas prone to earthquakes with M ≥ 5, where events of such magnitude have not been recorded up to date. Mordern geothermal activity and anomalies of dissolved methane (up to 10000 nL/L), helium and hydrogen in the hot mineral water springs prove high permeability of the study area and its active seismo-tectonic state.     

Time constant of hydraulic-head response in aquifers subjected to sudden recharge change: application to large basins

Analytical formulae are proposed to describe the first-order temporal evolution of the head in large groundwater systems (such as those found in North Africa or eastern Australia) that are subjected to drastic modifications of their recharge conditions (such as those in Pleistocene and Holocene times). The mathematical model is based on the hydrodynamics of a mixed-aquifer system composed of a confined aquifer connected to an unconfined one with a large storage capacity. The transient behaviour of the head following a sudden change of recharge conditions is computed with Laplace transforms for linear one-dimensional and cylindrical geometries. This transient evolution closely follows an exponential trend exp(?t/τ). The time constant τ is expressed analytically as a function of the various parameters characterizing the system. In many commonly occurring situations, τ depends on only four parameters: the width a _c of the main confined aquifer, its transmissivity T _c, the integrated storage situated upstream in the unconfined aquifer M?=?S _u a _u, and a curvature parameter accounting for convergence/divergence effects. This model is applied to the natural decay of large aquifer basins of the Sahara and Australia following the end of the mid-Holocene humid period. The observed persistence of the resource is discussed on the basis of the time constant estimated with the system parameters. This comparison confirms the role of the upstream water reserve, which is modelled as an unconfined aquifer, and highlights the significant increase of the time constant in case of converging flow.     

梯度法解释复杂二维断裂重力异常

对以往的二维断裂重力异常的梯度解释方法作了几点改进:①采用最小二乘拟合方法求取台阶倾角α;②直接利用台阶倾角与台阶的水平边界位置来计算台阶的埋深;③将重力二维梯度解释方法推广到复杂断裂模型的解释。文中详细论述了这种方法的原理, 给出了其具体实现步骤。该方法比起以往重力梯度解释方法, 应用更加简便而且解释结果更加稳定。用该法成功地解释了某油气盆地的实测重力资料。     

Emerging grassroots resilience and flood responses in informal settlements in Accra,Ghana

In the last five decades, many informal communities in Accra, Ghana have suffered from annual flood hazards. Residents of these communities appear to have successfully resisted evictions by city authorities; survived flood hazards and poor environmental health conditions. These flood affected households continue to survive with increasing housing and population densities in the face of these annual floods. Are they becoming resilient? Have residents built adaptive capacities through learning experiences from previous flood occurrences and evictions attempts? What has produced and continued to shape their responses to flooding? What can be learned from this supposed grassroots resilience to inform flood management in urban Africa? Using case studies of three informal communities of Glefe, Agbogbloshie and Old Fadama, this paper explores the gradual and evolving adaptive capacities and social resilience to flood hazards among poor urban dwellers. The paper reveals the depth of understanding and embodied nature of flood experiences among affected slum dwellers and how these are gradually being transformed into adaptive capacities and shaping their responses. In the absence of efficient state flood interventions, there are emerging and enduring flood responses and adaptation practices shaped by residents’ social networks, political alliances and sense of place. These responses translates into continuous re-structuring of housing units, construction of communal drains, creation of local evacuation teams and safe havens. Urban policy contributions that can be learned from these emerging grassroots capacities for flood vulnerability management have been proposed.