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).     

Beyond dual-porosity modeling for the simulation of complex flow mechanisms in shale reservoirs

The state of the art of modeling fluid flow in shale reservoirs is dominated by dual-porosity models which divide the reservoirs into matrix blocks that significantly contribute to fluid storage and fracture networks which principally control flow capacity. However, recent extensive microscopic studies reveal that there exist massive micro- and nano-pore systems in shale matrices. Because of this, the actual flow mechanisms in shale reservoirs are considerably more complex than can be simulated by the conventional dual-porosity models and Darcy’s law. Therefore, a model capturing multiple pore scales and flow can provide a better understanding of the complex flow mechanisms occurring in these reservoirs. This paper presents a micro-scale multiple-porosity model for fluid flow in shale reservoirs by capturing the dynamics occurring in three porosity systems: inorganic matter, organic matter (mainly kerogen), and natural fractures. Inorganic and organic portions of shale matrix are treated as sub-blocks with different attributes, such as wettability and pore structures. In kerogen, gas desorption and diffusion are the dominant physics. Since the flow regimes are sensitive to pore size, the effects of nano-pores and micro-pores in kerogen are incorporated into the simulator. The multiple-porosity model is built upon a unique tool for simulating general multiple-porosity systems in which several porosity systems may be tied to each other through arbitrary connectivities. This new model allows us to better understand complex flow mechanisms and eventually is extended into the reservoir scale through upscaling techniques. Sensitivity studies on the contributions of the different flow mechanisms and kerogen properties give some insight as to their importance. Results also include a comparison of the conventional dual-porosity treatment and show that significant differences in fluid distributions and dynamics are obtained with the improved multiple-porosity simulation.     

Adaptive particle swarm optimization for optimal orbital elements of binary stars

The paper presents an adaptive particle swarm optimization (APSO) as an alternative method to determine the optimal orbital elements of the star η Bootis of MK type G0 IV. The proposed algorithm transforms the problem of finding periodic orbits into the problem of detecting global minimizers as a function, to get a best fit of Keplerian and Phase curves. The experimental results demonstrate that the proposed approach of APSO generally more accurate than the standard particle swarm optimization (PSO) and other published optimization algorithms, in terms of solution accuracy, convergence speed and algorithm reliability.     

Micrometeorological Observations of a Microburst in Southern Finland

On the afternoon of 3 July 2004 in Hyytiälä (Juupajoki, Finland), convective cells produced a strong downburst causing forest damage. The SMEAR II field station, situated near the damage site, enabled a unique micrometeorological analysis of a microburst with differences above and inside the canopy. At the time of the event, a squall line associated with a cold front was crossing Hyytiälä with a reflectivity maximum in the middle of the squall line. A bow echo, rear-inflow notch, and probable mesovortex were observed in radar data. The bow echo moved west-north-west, and its apex travelled just north of Hyytiälä. The turbulence data were analysed at two locations above the forest canopy and at one location at sub-canopy. At 1412 EET (Eastern European Time, UTC+2), the horizontal and vertical wind speed increased and the wind veered, reflecting the arrival of a gust front. At the same time, the carbon dioxide concentration increased due to turbulent mixing, the temperature decreased due to cold air flow from aloft and aerosol particle concentration decreased due to rain scavenging. An increase in the number concentration of ultra-fine particles (< 10 nm) was detected, supporting the new particle formation either from cloud outflow or due to rain. Five minutes after the gust front (1417 EET), strong horizontal and downward vertical wind speed gusts occurred with maxima of 22 and 15 m s^?1, respectively, reflecting the microburst. The turbulence spectra before, during and after the event were consistent with traditional turbulence spectral theory.     

Phytoplankton productivity in the western subarctic gyre of the North Pacific in early summer 2006

We deployed a profiling buoy system incorporating a fast repetition rate fluorometer in the western subarctic Pacific and carried out time-series observations of phytoplankton productivity from 9 June to 15 July 2006. The chlorophyll a (Chl a) biomass integrated over the euphotic layer was as high as 45–50 mg Chl a m⁻² in the middle of June and remained in the 30–40 mg Chl a m⁻² range during the rest of observation period; day-to-day variation in Chl a biomass was relatively small. The daily net primary productivity integrated over the euphotic layer ranged from 144 to 919 mg C m⁻² day⁻¹ and varied greatly, depending more on insolation rather than Chl a biomass. In addition, we found that part of primary production was exported to a 150-m depth within 2 days, indicating that the variations in primary productivity quickly influenced the organic carbon flux from the upper ocean. Our results suggest that the short-term variability in primary productivity is one of the key factors controlling the carbon cycle in the surface ocean in the western subarctic Pacific.     

Hydrothermal plume mapping as a prospecting tool for seafloor sulfide deposits: a case study at the Zouyu-1 and Zouyu-2 hydrothermal fields in the southern Mid-Atlantic Ridge

Seafloor hydrothermal polymetallic sulfide deposits are a new type of resource, with great potential economic value and good prospect development. This paper discusses turbidity, oxidation–reduction potential, and temperature anomalies of hydrothermal plumes from the Zouyu-1 and Zouyu-2 hydrothermal fields on the southern Mid-Atlantic Ridge. We use the known location of these vent fields and plume data collected in multiple years (2009, 2011, 2013) to demonstrate how real-time plume exploration can be used to locate active vent fields, and thus associated sulfide deposits. Turbidity anomalies can be detected 10 s of km from an active source, but the location precision is no better than a few kilometers because fine-grained particles are quasi-conservative over periods of many days. Temperature and oxidation–reduction potential anomalies provide location precision of a few hundred meters. Temperature anomalies are generally weak and difficult to reliably detect, except by chance encounters of a buoyant plume. Oxidation–reduction potential is highly sensitive (nmol concentrations of reduced hydrothermal chemicals) to discharges of all temperatures and responds immediately to a plume encounter. Real-time surveys using continuous tows of turbidity and oxidation–reduction potential sensors offer the most efficient and precise surface ship exploration presently possible.     

Timing of the initial collision between the Indian and Asian continents

There exist three mainstream opinions regarding the timing of the initial collision between the Indian and Eurasian continents,namely,65±5,45±5,and 30±5 Ma.Five criteria are proposed for determining which tectonic event was related to the initial collision between India and Asia:the rapid decrease in the rate of plate motion,the cessation of magmatic activity originating from the subduction of oceanic crust,the end of sedimentation of oceanic facies,the occurrence of intracontinental deformation,and the exchange of sediments sourced from two continents.These criteria are used to constrain the nature of these tectonic events.It is proposed that the 65±5 Ma tectonic event is consistent with some of the criteria,but the upshot of this model is that the magmatic activity originating from the Tethyan subduction since the Mesozoic restarted along the southern margin of the Asian continent in this time after a brief calm,implying that the subduction of the Neotethys slab was still taking place.The magmatic activity that occurred along the southern margin of the Asian continent had a 7-Myr break during 72-65 Ma,which in this study is interpreted as having resulted from tectonic transformation from subduction to transform faulting,indicating that the convergence between the Indian and Asian continents was once dominated by strike-slip motion.The 30±5 Ma tectonic event resulted in the uplift of the Tibetan Plateau,which was related to the late stage of the convergence between these two continents,namely,a hard collision.The 45±5 Ma tectonic event is in accordance with most of the criteria,corresponding to the initial collision between these two continents.