Position and velocity perturbations for the determination of geopotential from space geodetic measurements

Although space geodetic observing systems have been advanced recently to such a revolutionary level that low Earth Orbiting (LEO) satellites can now be tracked almost continuously and at the unprecedented high accuracy, none of the three basic methods for mapping the Earth’s gravity field, namely, Kaula linear perturbation, the numerical integration method and the orbit energy-based method, could meet the demand of these challenging data. Some theoretical effort has been made in order to establish comparable mathematical modellings for these measurements, notably by Mayer-Gürr et al. (J Geod 78:462–480, 2005). Although the numerical integration method has been routinely used to produce models of the Earth’s gravity field, for example, from recent satellite gravity missions CHAMP and GRACE, the modelling error of the method increases with the increase of the length of an arc. In order to best exploit the almost continuity and unprecedented high accuracy provided by modern space observing technology for the determination of the Earth’s gravity field, we propose using measured orbits as approximate values and derive the corresponding coordinate and velocity perturbations. The perturbations derived are quasi-linear, linear and of second-order approximation. Unlike conventional perturbation techniques which are only valid in the vicinity of reference mean values, our coordinate and velocity perturbations are mathematically valid uniformly through a whole orbital arc of any length. In particular, the derived coordinate and velocity perturbations are free of singularity due to the critical inclination and resonance inherent in the solution of artificial satellite motion by using various types of orbital elements. We then transform the coordinate and velocity perturbations into those of the six Keplerian orbital elements. For completeness, we also briefly outline how to use the derived coordinate and velocity perturbations to establish observation equations of space geodetic measurements for the determination of geopotential.     

Geothermal energy potential of Tulsishyam thermal springs of Gujarat,India

Tulsishyam thermal springs are located in the Saurashtra region of Gujarat, India with discharge temperatures varying from 39 to 42 °C. The pH of these thermal springs varies from 7.1 to 7.4, indicating neutral character. Though these thermal springs propagate through the near surface layer of Deccan basalt, detailed geochemical analysis of the thermal waters using Piper diagram suggests that the water is interacting with the granitic basement rock. Silica and cation geothermometry estimates have reservoir temperature in the range of 138 to 207 °C categorizing it into a low to moderate enthalpy geothermal system. Furthermore, the area has high heat flow values of 53–90 mW/m² because of shallow Moho depth. The prevailing conditions suggest that the geothermal energy can potentially be exploited through an enhanced geothermal system (EGS). The study also indicates different mineral phases that may precipitate out of water during exploitation of geothermal energy and it should be taken into account while designing an EGS for the area.     

Assessment of H<Subscript>2</Subscript>S emission hazards into tunnels: the Nosoud tunnel case study from Iran

The presence of a gassy ground condition is an important problem in tunneling. In this study, the effects of groundwater H₂S and CH₄ emissions are investigated and characterized together with the factors that created these conditions in Nosoud tunnel in Iran. Through the geological investigations, the presence of these gasses was not detected prior to the construction of the tunnel. Groundwater sampling indicated that about 1 L of H₂S is released per 100 L of the water inflow into the Nosoud tunnel under normal conditions. However, the volume of the released gas was varying with the changes in the groundwater discharge rate. Thus, estimation of groundwater inflow into the tunnel is necessary for predicting the volume of gas emission. Based on the experience of the Nosoud tunnel excavations, there are several geological and hydrogeological factors that must be considered as the indicators of gas emissions during tunneling. Considering the importance of ground water gas emission into the tunnels located in gassy conditions, the present work was conducted to predict the H₂S seepage before the excavation using geological and hydrogeological indicators.     

Hydrogeological assessment of the Tabarteh fault zone by physicochemical analysis and multivariate statistical methods

Hydrogeologically, faults may impede, conduit, exert no influence, or may play a combination of these roles on groundwater flow. The object of this paper is to study the hydrogeological role of the Tabarteh fault, which is located on the border of Zagros and Central Iran tectonic zones in an alluvial aquifer. The recorded data of water table levels, chemical parameters, and discharge rate of wells, in addition to geological maps and geophysical results, were collected and evaluated. The outcrop of travertine in limited areas and the emergence of a few small springs within the alluvium show a barrier role of the fault in the groundwater flow. The spatial analysis of chemical components, head time series, and groundwater flow direction assessment demonstrated that the fault acts as both a barrier and a non-barrier in different sections. The multivariate statistical methods of cluster and discriminant analyses also confirm the dual role of the fault.     

Effects of historical and projected land use/cover change on runoff and sediment yield in the Netravati river basin,Western Ghats,India

In this study, the effects of changes in historical and projected land use land cover (LULC) on monthly streamflow and sediment yield for the Netravati river basin in the Western Ghats of India are explored using land use maps from six time periods (1972, 1979, 1991, 2000, 2012, and 2030) and the soil and water assessment tool (SWAT). The LULC for 2030 is projected using the land change modeller with the assumption of normal growth. The sensitivity analysis, model calibration, and validation indicated that the SWAT model could reasonably simulate streamflow and sediment yield in the river basin. The results showed that the spatial extent of the LULC classes of urban (1.80–9.96%), agriculture (31.38–55.75%), and water bodies (1.48–2.66%) increased, whereas that of forest (53.04–27.03%), grassland (11.17–4.41%), and bare land (1.09–0.16%) decreased from 1972 to 2030. The streamflow increased steadily (7.88%) with changes in LULC, whereas the average annual sediment yield decreased (0.028%) between 1972 and 1991 and increased later (0.029%) until 2012. However, it may increase by 0.43% from 2012 to 2030. The results indicate that LULC changes in urbanization and agricultural intensification have contributed to the increase in runoff, amounting to 428.65 and 58.67 mm, respectively, and sediment yield, amounting to 348 and 43 ton/km², respectively, in the catchment area from 1972 to 2030. The proposed methodology can be applied to other river basins for which temporal digital LULC maps are available for better water resource management plans.     

Estimating satellite phase fractional cycle biases based on Kalman filter

Phase fractional cycle biases (FCBs) originating from satellites and receivers destroy the integer nature of PPP carrier phase ambiguities. To achieve integer ambiguity resolution of PPP, FCBs of satellites are required. In former work, least squares methods are commonly adopted to isolate FCBs from a network of reference stations. However, it can be extremely time consuming concerning the large number of observations from hundreds of stations and thousands of epochs. In addition, iterations are required to deal with the one-cycle inconsistency among FCB measurements. We propose to estimate the FCB based on a Kalman filter. The large number of observations are handled epoch by epoch, which significantly reduces the dimension of the involved matrix and accelerates the computation. In addition, it is also suitable for real-time applications. As for the one-cycle inconsistency, a pre-elimination method is developed to avoid iterations and posterior adjustments. A globally distributed network consisting of about 200 IGS stations is selected to determine the GPS satellite FCBs. Observations recorded from DoY 52 to 61 in 2016 are processed to verify the proposed approach. The RMS of wide lane (WL) posterior residuals is 0.09 cycles while that of the narrow lane (NL) is about 0.05 cycles, which indicates a good internal accuracy. The estimated WL FCBs also have a good consistency with existing WL FCB products (e.g., CNES-GRG, WHU-SGG). The RMS of differences with respect to GRG and SGG products are 0.03 and 0.05 cycles. For satellite NL FCB estimates, 97.9% of the differences with respect to SGG products are within ±?0.1 cycles. The RMS of the difference is 0.05 cycles. These results prove the efficiency of the proposed approach.     

Estimation of Above Ground Biomass for Central Indian Deciduous Forests Using ALOS PALSAR L-Band Data

Reliable and accurate estimates of tropical forest above ground biomass (AGB) are important to reduce uncertainties in carbon budgeting. In the present study we estimated AGB of central Indian deciduous forests of Madhya Pradesh (M.P.) state, India, using Advanced Land Observing Satellite – Phased Array type L-band Synthetic Aperture Radar (ALOS-PALSAR) L-band data of year 2010 in conjunction with field based AGB estimates using empirical models. Digital numbers of gridded 1?×?1° dual polarization (HH & HV) PALSAR mosaics for the study area were converted to normalized radar cross section (sigma naught - σ⁰). A total of 415 sampling plots (0.1 ha) data collected over the study area during 2009–10 was used in the present study. Plot-level AGB estimates using volume equations representative to the study area were computed using field inventory data. The plot-level AGB estimates were empirically modeled with the PALSAR backscatter information in HH, HV and their ratios from different forest types of the study area. The HV backscatter information showed better relation with field based AGB estimates with a coefficient of determination (R²) of 0.509 which was used to estimate spatial AGB of the study area. Results suggested a total AGB of 367.4 Mt for forests of M.P. state. Further, validation of the model was carried out using observed vs. predicted AGB estimates, which suggested a root mean square error (RMSE) of ±19.32 t/ha. The model reported robust and defensible relation for observed vs. predicted AGB values of the study area.     

Heat tracing to examine seasonal groundwater flow beneath a low-gradient stream in rural central Illinois,USA

The thermal profile of a streambed is affected by a number of factors including: temperatures of stream water and groundwater, hydraulic conductivity, thermal conductivity, heat capacity of the streambed, and the geometry of hyporheic flow paths. Changes in these parameters over time cause changes in thermal profiles. In this study, temperature data were collected at depths of 30, 60, 90 and 150 cm at six streambed wells 5 m apart along the thalweg of Little Kickapoo Creek, in rural central Illinois, USA. This is a third-order low-gradient baseflow-fed stream. A positive temperature gradient with inflection at 90-cm depth was observed during the summer period. A negative temperature gradient with inflection at 30 cm was observed during the winter period, which suggests greater influence of stream-water temperatures in the substrate during the summer. Thermal models of the streambed were built using VS2DHI to simulate the thermal profiles observed in the field. Comparison of the parameters along with analysis of temperature envelopes and Peclet numbers suggested greater upwelling and stability in temperatures during the winter than during the summer. Upwelling was more pronounced in the downstream reach of the pool in the riffle and pool sequence.     

Using radon-222 to study coastal groundwater/surface-water interaction in the Crau coastal aquifer (southeastern France)

Radon has been used to determine groundwater velocity and groundwater discharge into wetlands at the southern downstream boundary of the Crau aquifer, southeastern France. This aquifer constitutes an important high-quality freshwater resource exploited for agriculture, industry and human consumption. An increase in salinity occurs close to the sea, highlighting the need to investigate the water balance and groundwater behavior. Darcy velocity was estimated using radon activities in well waters according to the Hamada “single-well method” (involving comparison with radon in groundwater in the aquifer itself). Measurements done at three depths (7, 15 and 21 m) provided velocity ranging from a few mm/day to more than 20 cm/day, with highest velocities observed at the 15-m depth. Resulting hydraulic conductivities agree with the known geology. Waters showing high radon activity and high salinity were found near the presumed shoreline at 3,000 years BP, highlighting the presence of ancient saltwater. Radon activity has also been measured in canals, rivers and ponds, to trace groundwater discharges and evaluate water balance. A model of the radon spatial evolution explains the observed radon activities. Groundwater discharge to surface water is low in pond waters (4 % of total inputs) but significant in canals (55 l/m²/day).