Fast computation of Hankel Transform using orthonormal exponential approximation of complex kernel function

The computation of electromagnetic (EM) fields, for 1-D layered earth model, requires evaluation of Hankel Transform (HT) of the EM kernel function. The digital filtering is the most widely used technique to evaluate HT integrals. However, it has some obvious shortcomings. We present an alternative scheme, based on an orthonormal exponential approximation of the kernel function, for evaluating HT integrals. This approximation of the kernel function was chosen because the analytical solution of HT of an exponential function is readily available in literature. This expansion reduces the integral to a simple algebraic sum. The implementation of such a scheme requires that the weights and the exponents of the exponential function be estimated. The exponents were estimated through a guided search algorithm while the weights were obtained using Marquardt matrix inversion method. The algorithm was tested on analytical HT pairs available in literature. The results are compared with those obtained using the digital filtering technique with Anderson filters. The field curves for four types (A-, K-, H-and Q-type) of 3-layer earth models are generated using the present scheme and compared with the corresponding curves obtained using the Anderson sc heme. It is concluded that the present scheme is more accurate than the Anderson scheme     

Construction and Characterization of Transportable Calibration Pads for Calibration of Portable and Airborne Gamma Ray Spectrometers in India

Characterization of the Panandhro lignite deposits from western Indian state of Gujarat, based on the geochemical and palynological evidences, has been performed to assess the floral composition, maturity and hydrocarbon potential of the sequence. Elementally, the lignites consist of moderate carbon, low hydrogen and moderate sulfur contents. The samples are characterized by high TOC contents (lignite: av. 46.43 wt.%, resin: 62.47 wt.%). The average HI values for the lignite is 136 mg HC/g TOC, and that of the associated resin is 671 mg HC/g TOC. The highest T_max is recoded in lignite (422°C) and lowest in the resin (39°C) samples. The FTIR spectrum of lignite is characterized by highly intense OH stretching peak ~3350 cm^-1, aliphatic CH_x stretching peaks between 3000-2800 cm^-1, aromatic C=O stretching and an aromatic C=C stretching. The spectrum of resin shows strongest absorption due to aliphatic CHx stretching between 2940-2915 cm^-1 and 2870-2850 cm^-1, and deformation by the medium peak between 1450 and 1650 cm^-1. The recovered palynofloral assemblage indicates the dominance of angiosperm pollen grains with maximum abundance of Arecaceae family, and subdominant pteridophytic spores. Marine influence is indicated by the presence of abundant dinoflagellate cysts. The occurrence of flora from a variety of ecological niches suggests a luxuriant diverse vegetation pattern existed in the vicinity of depositional site under humid tropical conditions. The overall characteristics of the lignite deposits point towards their ability to generate (upon maturation) hydrocarbons as they have types III–II admixed kerogen (organic matters).     

Influences of local hydroclimatology and teleconnections on Florida's precipitation and temperature variability

Understanding the influences of local hydroclimatology and two large-scale oceanic-atmospheric oscillations (i.e., Atlantic Multidecadal Oscillation (AMO) and El Niño-Southern Oscillation (ENSO)) on seasonal precipitation (P) and temperature (T) relationships for a tropical region (i.e., Florida) is the focus of this study. The warm and cool phases of AMO and ENSO are initially identified using sea surface temperatures (SSTs). The associations of SSTs and regional minimum, maximum and average surface air temperatures (SATs) with precipitation are then evaluated. The seasonal variations in P-SATs and P-SSTs associations considering AMO and ENSO phases for sites in (1) two soil temperature regimes (i.e., thermic and hyperthermic); (2) urban and non-urban regions; and (3) regions with and without water bodies, are analysed using two monthly datasets. The analyses are carried out using trend tests, two association measures, nonparametric and parametric statistical hypothesis tests and kernel density estimates. Decreasing (increasing) trend in precipitation (SATs) is noted in the recent multi-decadal period (1985–2019) compared to the previous one (1950–1984) indicating a progression towards warmer and drier climatic conditions across Florida. Spatially and temporally non-uniform variations in the associations of precipitation with SATs and SSTs are noted. Strong positive (weak negative) P–T associations are noted during the wet (dry) season for both AMO phases and El Niño, while significant (positive) P–T associations are observed across southern Florida during La Niña in the dry season. The seasonal influences are predominant in governing the P–T relationship over the regions with and without water bodies; however, considerable variations between El Niño and La Niña are noted during the dry season. The climate variability influences on P–T correlations for hyperthermic and thermic soil zones are found to be insignificant (significant) during the wet (dry) season. Nonparametric clustering is performed to identify the spatial clusters exhibiting homogeneous P–T relationships considering seasonal and climate variability influences.     

Combined straightforward inversion of resistivity and induced polarization (time-domain) sounding data

It is proposed that the Straightforward Inversion Scheme (SIS) developed by the authors for 1D inversion of resistivity sounding and magneto-telluric sounding data can also be used in similar fashion for time-domain induced polarization sounding data. The necessary formulations based on dynamic dipole theory are presented. It is shown that by using induced polarization potential, measured at the instant when steady state current is switched off, an equation can be developed for apparent ‘chargeability–resistivity’ which is similar to the one for apparent resistivity. The two data sets of apparent resistivity and apparent chargeability–resistivity can be inverted in a combined manner, using SIS for a common uniform thickness layer earth model to estimate the respective subsurface distributions of resistivity and chargeability–resistivity. The quotient of the two profiles will give the sought after chargeability profile. A brief outline of SIS is provided for completeness. Three theoretical models are included to confirm the efficacy of SIS software by inverting only the synthetic resistivity sounding data. Then one synthetic data set based on a geological model and three field data sets (combination of resistivity and IP soundings) from diverse geological and geographical regions are included as validation of the proposal. It is hoped that the proposed scheme would complement the resistivity interpretation with special reference to shaly sand formations.     

Evidence of topographic disequilibrium in the Subarnarekha River Basin,India: A digital elevation model based analysis

Cratonic areas experience complex process-response changes due to their operative endogenic and exogenic forces varying in intensity and spatiality over long timescales. Unlike zones of active deformation, the surface expression of the transient signals in relatively tectonically stable areas are usually scant. The Subarnarekha River Basin, in eastern India, is a prime example of a Precambrian cratonic landscape, overlain in places by Tertiary and Quaternary deposits. A coupled quantitative-qualitative approach is employed towards deciphering tectonic and geological influences across linear and areal aspects, at the basin and sub-basin scale. Within this landscape, the transient erosional signatures are explored, as recorded in the disequilibrium conditions of the longitudinal profiles of the major streams, which are marked by a number of waterfalls at structural and lithological boundaries. Mathematical expressions derived from the normalized longitudinal profiles of these streams are used to ascertain their stage of development. Cluster analysis and chi plots provide significant interpretations of the role of vertical displacements or litho-structural variations within the basin. These analyses suggest that a heterogeneous, piece-meal response to the ongoing deformation exists in the area, albeit, determining the actual rate of this deformation or its temporal variation is difficult without correlated chronological datasets.     

Potential Usage of Remote Sensing Data in Studying the Behaviour of Shore Drift Along Kerala Coast, India

Landsat images were analysed to examine direction, amount and behaviour of long-shore drift and its contribution towards deciding areas vulnerable to coastal erosion and accretion along the coast of Kerala. Two approaches are adopted while utilizing the remote sensing data. Net shore-drift direction is determined firstly by studying various coastal landform indicators and secondly by studying offshore turbidity distribution patterns. Landform indicators study suggest that during the south-west monsoon period, strong southerly currents are eroding protruded sectors and depositing eroded material, though partially, along varying sectors. Whereas during the rest of the year, under the influence of a northerly current, accretion is taking place along retreating sectors which has been reflected in the development of beach ridges. Turbidity pattern distribution exhibited by satellite image suggests that for January and February, long-shore drift is from south to north. The results of the two studies were found to be comparable. It may be concluded that shore-drift direction, which can be determined effectively by remote sensing, is season-dependent and plays a significant role in deciding the areas of coastal erosion and accretion.     

Prediction of sea ice edge in the Antarctic using GVF Snake model

Antarctic sea ice cover plays an important role in shaping the earth’s climate, primarily by insulating the ocean from the atmosphere and increasing the surface albedo. The convective processes accompanied with the sea ice formation result bottom water formation. The cold and dense bottom water moves towards the equator along the ocean basins and takes part in the global thermohaline circulation. Sea ice edge is a potential indicator of climate change. Additionally, fishing and commercial shipping activities as well as military submarine operations in the polar seas need reliable ice edge information. However, as the sea ice edge is unstable in time, the temporal validity of the estimated ice edge is often shorter than the time required to transfer the information to the operational user. Hence, an accurate sea ice edge prediction as well as determination is crucial for fine-scale geophysical modeling and for near-real-time operations. In this study, active contour modelling (known as Snake model) and non-rigid motion estimation techniques have been used for predicting the sea ice edge (SIE) in the Antarctic. For this purpose the SIE has been detected from sea ice concentration derived using special sensor microwave imager (SSM/I) observations. The 15% sea ice concentration pixels are being taken as the edge pixel between ice and water. The external force, gradient vector flow (GVF), of SIE for total the Antarctic region is parameterised for daily as well as weekly data set. The SIE is predicted at certain points using a statistical technique. These predicted points have been used to constitute a SIE using artificial intelligence technique, the gradient vector flow (GVF). The predicted edge has been validated with that of SSM/I. It is found that all the major curvatures have been captured by the predicated edge and it is in good agreement with that of the SSM/I observation.