Reducing the draconitic errors in GNSS geodetic products

Systematic errors at harmonics of the GPS draconitic year have been found in diverse GPS-derived geodetic products like the geocenter $Z$ -component, station coordinates, $Y$ -pole rate and orbits (i.e. orbit overlaps). The GPS draconitic year is the repeat period of the GPS constellation w.r.t. the Sun which is about 351 days. Different error sources have been proposed which could generate these spurious signals at the draconitic harmonics. In this study, we focus on one of these error sources, namely the radiation pressure orbit modeling deficiencies. For this purpose, three GPS+GLONASS solutions of 8 years (2004–2011) were computed which differ only in the solar radiation pressure (SRP) and satellite attitude models. The models employed in the solutions are: (1) the CODE (5-parameter) radiation pressure model widely used within the International GNSS Service community, (2) the adjustable box-wing model for SRP impacting GPS (and GLONASS) satellites, and (3) the adjustable box-wing model upgraded to use non-nominal yaw attitude, specially for satellites in eclipse seasons. When comparing the first solution with the third one we achieved the following in the GNSS geodetic products. Orbits: the draconitic errors in the orbit overlaps are reduced for the GPS satellites in all the harmonics on average 46, 38 and 57 % for the radial, along-track and cross-track components, while for GLONASS satellites they are mainly reduced in the cross-track component by 39 %. Geocenter $Z$ -component: all the odd draconitic harmonics found when the CODE model is used show a very important reduction (almost disappearing with a 92 % average reduction) with the new radiation pressure models. Earth orientation parameters: the draconitic errors are reduced for the $X$ -pole rate and especially for the $Y$ -pole rate by 24 and 50 % respectively. Station coordinates: all the draconitic harmonics (except the 2nd harmonic in the North component) are reduced in the North, East and Height components, with average reductions of 41, 39 and 35 % respectively. This shows, that part of the draconitic errors currently found in GNSS geodetic products are definitely induced by the CODE radiation pressure orbit modeling deficiencies.     

Radiative forcing and climate impact resulting from SO2 injections based on a 200,000-year record of Plinian eruptions along the Central American Volcanic Arc

We present for the first time a self-consistent methodology connecting volcanological field data to global climate model estimates for a regional time series of explosive volcanic events. Using the petrologic method, we estimated SO₂ emissions from 36 detected Plinian volcanic eruptions occurring at the Central American Volcanic Arc (CAVA) during the past 200,000 years. Together with simple parametrized relationships collected from past studies, we derive estimates of global maximum volcanic aerosol optical depth (AOD) and radiative forcing (RF) describing the effect of each eruption on radiation reaching the Earth’s surface. In parallel, AOD and RF time series for selected CAVA eruptions are simulated with the global aerosol model MAECHAM5-HAM, which shows a relationship between stratospheric SO₂ injection and maximum global mean AOD that is linear for smaller volcanic eruptions (<5 Mt SO₂) and nonlinear for larger ones (≥5 Mt SO₂) and is qualitatively and quantitatively consistent with the relationship used in the simple parametrized approximation. Potential climate impacts of the selected CAVA eruptions are estimated using an earth system model of intermediate complexity by RF time series derived by (1) directly from the global aerosol model and (2) from the simple parametrized approximation assuming a 12-month exponential decay of global AOD. We find that while the maximum AOD and RF values are consistent between the two methods, their temporal evolutions are significantly different. As a result, simulated global maximum temperature anomalies and the duration of the temperature response depend on which RF time series is used, varying between 2 and 3 K and 60 and 90 years for the largest eruption of the CAVA dataset. Comparing the recurrence time of eruptions, based on the CAVA dataset, with the duration of climate impacts, based on the model results, we conclude that cumulative impacts due to successive eruptions are unlikely. The methodology and results presented here can be used to calculate approximate volcanic forcings and potential climate impacts from sulfur emissions, sulfate aerosol or AOD data for any eruption that injects sulfur into the tropical stratosphere.     

Variability over time in the sources of South Portuguese Zone turbidites: evidence of denudation of different crustal blocks during the assembly of Pangaea

This study combines geochemical and geochronological data in order to decipher the provenance of Carboniferous turbidites from the South Portuguese Zone (SW Iberia). Major and trace elements of 25 samples of graywackes and mudstones from the Mértola (Visean), Mira (Serpukhovian), and Brejeira (Moscovian) Formations were analyzed, and 363 U-Pb ages were obtained on detrital zircons from five samples of graywackes from the Mira and Brejeira Formations using LA-ICPMS. The results indicate that turbiditic sedimentation during the Carboniferous was marked by variability in the sources, involving the denudation of different crustal blocks and a break in synorogenic volcanism. The Visean is characterized by the accumulation of immature turbidites (Mértola Formation and the base of the Mira Formation) inherited from a terrane with intermediate to mafic source rocks. These source rocks were probably formed in relation to Devonian magmatic arcs poorly influenced by sedimentary recycling, as indicated by the almost total absence of pre-Devonian zircons typical of the Gondwana and/or Laurussia basements. The presence of Carboniferous grains in Visean turbidites indicates that volcanism was active at this time. Later, Serpukhovian to Moscovian turbiditic sedimentation (Mira and Brejeira Formations) included sedimentary detritus derived from felsic mature source rocks situated far from active magmatism. The abundance of Precambrian and Paleozoic zircons reveals strong recycling of the Gondwana and/or Laurussia basements. A peri-Gondwanan provenance is indicated by zircon populations with Neoproterozoic (Cadomian-Avalonian and Pan-African zircon-forming events), Paleoproterozoic, and Archean ages. The presence of late Ordovician and Silurian detrital zircons in Brejeira turbidites, which have no correspondence in the Gondwana basement of SW Iberia, indicates Laurussia as their most probable source.     

Ozone trends in the vertical structure of Upper Troposphere and Lower stratosphere over the Indian monsoon region

Ozone trends in the Upper Troposphere and Lower Stratosphere over the Indian region are investigated using three satellite data sets namely Halogen Occultation Experiment (1993–2005), Stratospheric Aerosol and Gas Experiment (1993–2005) II, and Aura Microwave Limb Sounder (MLS, 2005–2011). Estimated ozone trends using multi-variate regression analysis are compared with trends at two Indian ozonesonde stations (Delhi, 28°N, 77°E and Pune, 18°N, 73°E), and a 3-D Chemical Transport Model (CTM, SLIMCAT) for the 1993–2005 time period. Overall, all the observational data sets and model simulations indicate significant increasing trend in the upper troposphere (0–2.5 %/year). In the lower stratosphere, estimated trends are slightly positive up to 30 mb and are negative between 30 and 10 mb. Increasing trends in the upper troposphere is probably due to increasing trends in the tropospheric ozone precursor gases (e.g. CO, NO _x , NMHCs). Here, we argue that these contrasting ozone-trend profiles might be partially responsible for insignificant long-term trends in the tropical total column ozone. On seasonal scale, positive trends are observed during all the seasons in the upper troposphere while structure of trend profile varies in lower stratosphere. Seasonal variations of ozone trends and its linkages with stratospheric intrusions and increasing trends in lightning flashes in the troposphere are also discussed.     

Geochemistry and stable isotope composition of surface waters from the Ravenna plain (Italy): implications for the management of water resources in agricultural lands

This paper examines the issues related to salinization and water quality in the complex drainage system of the historical land reclamation of Lamone basin (Ravenna coastal plain, northern Italy), with the aim of guiding ongoing agricultural-related decisions. Major and trace element concentration and O–H–B–Sr isotope-ratios were measured on surface water from a network of canals and ditches. Sampling was carried out during the winter period and in summer, to assess the effects of agricultural water management on the aquatic system. Results show widespread salinization of waters over the entire crop, due to both the direct saltwater inland flushing through the canals in proximity of the mouths, and sea salt leaching from soils. During winter, salinization is partly mitigated by rainfalls, while in summer dilution processes are due to freshwater input from outside the watershed, planned to assist the agricultural water demand. In the winter season, the concentration of some elements such as Fe, Mn, Al, and Cu in waters exceeds the maximum permissible limits imposed by the Italian regulations, while during summer the concentration of these elements is significantly reduced. A seasonal cycling is established, where the bottom sediments of canals and ditches act as sinks of harmful elements, mostly through adsorption by Fe–Mn–Al oxy-hydroxide solid phases. The irrigation practice, although improving the water quality, increases water turbidity by re-mobilization of sediments which act as transport agents of contaminants, with detrimental effects that may become significant over the years.     

Geo-Characterization at selected accelerometric stations in Crete (Greece) and comparison of earthquake data recordings with EC8 elastic spectra

To estimate the seismic response according to Eurocode (EC8) and almost all other national codes, site conditions have to be properly characterized so that soil amplification and the corresponding peak ground motion can be calculated. In this work, different geophysical and geotechnical methods are combined in order to define the detailed ground conditions in selected sites of the Hellenic Accelerometric Network (HAN) in Crete. For this purpose, the geological information of the sites and shear wave velocity, calculated from surface wave measurements, is used. Additionally, ground acceleration data recorded through HAN have been utilized from intermediate depth earthquakes in the broader area of South Aegean Sea. Using the recorded ground motion data and the procedure defined in EC8, the corresponding elastic response spectrum is calculated for the selected sites. The resulting information is compared to the values defined in the corresponding EC8 spectrum for the seismic zone that includes the island of Crete. The comparison shows that accurate definition of ground type through geological, geotechnical and geophysical investigations is important. However, our current comparison focuses on the distribution of values rather than the absolute values of EC8-prescribed spectra, and the results should be considered in this context.     

Soil loss estimation by field measurements in the badlands along Pravara river (Western India)

Nearly 4 million hectares of land in India is affected by rill and gully erosion. It is a severe form of soil erosion rendering vast tracts of lands into wastelands giving rise to the formation of badlands. Rate of land degradation through gullying is triggered recently due to the ill-advised landuse practices all over the country. The study area represents one such badland locality along the banks of Pravara river and two of its tributaries in the Deccan trap region, Maharashtra. The area is extensively reclaimed in the last two decades for agricultural practices. The rates in the cross profile changes as well as land lowering were measured with the help of a self fabricated micro-profilometer and erosion pin method. Five first order gullies were monitored using the micro-profilometer. Results of the micro-profilometer technique reveal considerable changes in the cross-sectional areas of individual gullies in this area, which shows the dynamic soil removal mechanism operating in this region. In addition to this, two well developed gully catchments were surveyed using a theodolite and erosion pins were installed in these basins in May 2007. After one year, the exposed pinheads were measured and net gain and net loss in the volume of the sediments were calculated for both the basins. The same practice continued in May 2009. After two years of monitoring the sediment yields were calculated. Average of these two years indicates an annual sediment yield of 1.79 kg/m² for sample basin 1 and 0.76 kg/m² for sample basin 2. The calculated yield was compared with the threshold values of soil tolerance limits and it was found that the area has crossed the threshold limit of soil tolerance.     

Impact of Doppler weather radar data on thunderstorm simulation during STORM pilot phase—2009

This study assesses the impact of Doppler weather radar (DWR) data (reflectivity and radial wind) assimilation on the simulation of severe thunderstorms (STS) events over the Indian monsoon region. Two different events that occurred during the Severe Thunderstorms Observations and Regional Modeling (STORM) pilot phase in 2009 were simulated. Numerical experiments—3DV (assimilation of DWR observations) and CNTL (without data assimilation)—were conducted using the three-dimensional variational data assimilation technique with the Advanced Research Weather Research and Forecasting model (WRF-ARW). The results show that consistent with prior studies the 3DV experiment, initialized by assimilation of DWR observations, performed better than the CNTL experiment over the Indian region. The enhanced performance was a result of improved representation and simulation of wind and moisture fields in the boundary layer at the initial time in the model. Assimilating DWR data caused higher moisture incursion and increased instability, which led to stronger convective activity in the simulations. Overall, the dynamic and thermodynamic features of the two thunderstorms were consistently better simulated after ingesting DWR data, as compared to the CNTL simulation. In the 3DV experiment, higher instability was observed in the analyses of thermodynamic indices and equivalent potential temperature (θ _e) fields. Maximum convergence during the mature stage was also noted, consistent with maximum vertical velocities in the assimilation experiment (3DV). In addition, simulated hydrometeor (water vapor mixing ratio, cloud water mixing ratio, and rain water mixing ratio) structures improved with the 3DV experiment, compared to that of CNTL. From the higher equitable threat scores, it is evident that the assimilation of DWR data enhanced the skill in rainfall prediction associated with the STS over the Indian monsoon region. These results add to the body of evidence now which provide consistent and notable improvements in the mesoscale model results over the Indian monsoon region after assimilating DWR fields.     

Disc precession in the M31 dipping X-ray binary Bo 158?

We present results from three XMM–Newton observations of the M31 low mass X-ray binary (LMXB) XMMU J004314.4+410726.3 (Bo 158), spaced over 3 d in 2004 July. Bo 158 was the first dipping LMXB to be discovered in M31. Periodic intensity dips were previously seen to occur on a 2.78-h period, due to absorption in material that is raised out of the plane of the accretion disc. The report of these observations stated that the dip depth was anticorrelated with source intensity. In light of the 2004 XMM–Newton observations of Bo 158, we suggest that the dip variation is due to precession of the accretion disc. This is to be expected in LMXBs with a mass ratio ≲0.3 (period ≲4 h), as the disc reaches the 3:1 resonance with the binary companion, causing elongation and precession of the disc. A smoothed particle hydrodynamics simulation of the disc in this system shows retrograde rotation of a disc warp on a period of  ∼11 P _orb  , and prograde disc precession on a period of  29 ± 1 P _orb  . This is consistent with the observed variation in the depth of the dips. We find that the dipping behaviour is most likely to be modified by the disc precession, hence we predict that the dipping behaviour repeats on an  81 ± 3 h  cycle.