Earth Rotation and Gravity Field Parameters from Satellite Laser Ranging

We present the results from a simultaneous estimation of the gravity field, Earth rotation parameters, and station coordinates from combined SLR solutions incorporating up to nine geodetic satellites: LAGEOS-1/2, Starlette, Stella, AJISAI, Beacon-C, Lares, Blits and LARES. These solutions cover all three pillars of satellite geodesy and ensure full consistency between the Earth rotation parameters, gravity field coefficients, and geometry-related parameters. We address benefits emerging from such an approach and discuss particular aspects and limitations of the gravity field recovery using SLR data. The current accuracy of SLR-derived polar motion, by the means of WRMS w.r.t. IERS-08-C04 series, is at a level of 118-149 μas, which corresponds to 4 to 5 mm on the Earth’s surface. The WRMS of SLR-derived Length-of-Day, when the gravity field parameters are simultaneously estimated, is 56 μs/day, corresponding to about 26 mm on the ground, and the mean bias of SLR-derived Length-of-Day is 6.3 μs/day, corresponding to 3 mm.

Variability of mesospheric water vapor above Bern in relation to the 27-day solar rotation cycle

Many studies investigated solar–terrestrial responses (thermal state, O₃ , OH, H₂O) with emphasis on the tropical upper atmosphere. In this paper the Focus is switched to water vapor in the mesosphere at a mid-latitudinal location. Eight years of water vapor profile measurements above Bern (46.88°N/7.46°E) are investigated to study oscillations with the Focus on periods between 10 and 50 days. Different spectral analyses revealed prominent features in the 27-day oscillation band, which are enhanced in the upper mesosphere (above 0.1 hPa, ∼64 km) during the rising sun spot activity of solar cycle 24. Local as well as zonal mean Aura MLS observations Support these results by showing a similar behavior. The relationship between mesospheric water and the solar Lyman-α flux is studied by comparing thesi-milarity of their temporal oscillations. The H₂O oscillation is negatively correlated to solar Lyman-α oscillation with a correlation coefficient of up to −0.3 to −0.4, and the Phase lag is 6–10 days at 0.04 hPa. The confidence level of the correlation is ≥99%. This finding supports the assumption that the 27-day oscillation in Lyman-α causes a periodical photo dissociation loss in mesospheric water. Wavelet power spectra, cross-wavelet transform and wavelet coherence analysis (WTC)complete our study. More periods of high common wavelet power of H₂O and solar Lyman-α are present when amplitudes of the Lyman-α flux increase. Since this is not a measure of physical correlation a more detailed view on WTC is necessary, where significant (two sigma level)correlations occur intermittently in the 27 and 13-day band with variable Phase lock behavior. Large Lyman-α oscillations appeared after the solar super storm in July 2012 and the H₂O oscillations show a well pronounced anticorrelation. The competition between advective transport and photo dissociation loss of mesospheric water vapor may explain the sometimes variable Phase relationship of mesospheric H₂O and solar Lyman-α oscillations. Generally, the WTC analysis indicates that solar variability causes observable photochemical and dynamical processes in the mid-latitude mesosphere.

Present-day King George Island ice cap geometry measured by ground penetrating radar in 2008/09

We present a consistent data set for the ice thickness, the bedrock topography and the ice surface topography of the King George Island ice cap (Arctowski Icefield and the adjacent central part). The data set is composed of groundbased and airborne Ground Penetrating Radar (GPR) and differential GPS (DGPS) measurements, obtained during several field campaigns. The data set incorporates groundbased measurements in the safely accessible inner parts and airborne measurements in the heavily crevassed coastal areas of the ice cap. In particular, the inclusion of airborne GPR measurements with the 30MHz BGR-P30-System developed at the Institute of Geophysics (University of Münster) completes the picture of the ice geometry substantially. The compiled digital elevation model of the bedrock shows a rough, highly variable topography with pronounced valleys, ridges, and troughs. The mean ice thickness is approx. 238m, with a maximum value of approx. 400m in the surveyed area. Noticeable are bounded areas in the bedrock topography below sea level where marine based ice exists.

Lead and strontium isotope composition, mass accumulation rates and rare earth elements of selected samples from DSDP Sites 92-597 to 92-601

Most of the Pb isotope data for the Leg 92 metalliferous sediments (carbonate-free fraction) form approximately linear arrays in the conventional isotopic plots, extending from the middle of the field for mid-ocean ridge basalts (MORB) toward the field for Mn nodules. These arrays are directed closely to the average values of Mn nodules, the composition of which reflects the Pb isotope composition of seawater (Reynolds and Dasch, 1971). Since the Leg 92 samples are almost devoid of continentally derived detritus, it can be inferred that the more radiogenic end-member is seawater. The less radiogenic end-member lies in the very middle of the MORB field, and hence can be considered to reflect the Pb isotope composition of typical ocean-ridge basalt. The array of data lying between these two end-members is most readily interpreted in terms of simple linear mixing of Pb from the two different end-member sources. According to this model, eight samples from Sites 599 to 601 contain 50 to 100% basaltic Pb. Five of these samples have compositions that are identical within the uncertainty of the analyses. We use the average of these five values to define our unradiogenic end-member in the linear mixing model. The ratios used for this average are 206Pb/204Pb = 18.425 ± 0.010; 207Pb/204Pb = 15.495 ± 0.018; 208Pb/204Pb = 37.879 ± 0.068. These values should approximate the average Pb isotope composition of discharging hydrothermal solutions, and therefore also that of the basaltic crust, over the period of time represented by these samples ( 4 m.y., from 4 to 8 Ma). Sr isotope ratios show a significant range of values, from 0.7082 to 0.7091. The lower ratios are well outside the value of 0.70910 ± 6 for modern-day seawater (Burke et al., 1982). However, most values correspond very closely to the curve of 87Sr/86Sr versus age for seawater, with older samples having progressively lower 87Sr/86Sr ratios. The simplest explanation for this progressive reduction is that recrystallization of the abundant biogenic carbonate in the sediments released older seawater Sr which was incorporated into ferromanganiferous phases during diagenesis. Leg 92 metalliferous sediments have total rare earth element (REE) contents that range on a carbonate-free basis from 131 to 301 ppm, with a clustering between 167 and 222 ppm. The patterns have strong negative Ce anomalies. Samples from Sites 599 to 601 display a slight but distinct enrichment in the heavy REE relative to the light REE, whereas those from Sites 597 to 598 show almost no heavy REE enrichment. The former patterns (those for Sites 599 to 601) are interpreted as indicating moderate diagenetic alteration of metalliferous sediments originating at the EPR axis; the latter reflect more complete diagenetic modification.

Geochemistry of rare earth elements in basalts from Walvis Ridge

Selected basalts from a suite of dredged and drilled samples (IPOD sites 525, 527, 528 and 530) from the Walvis Ridge have been analysed to determine their rare earth element (REE) contents in order to investigate the origin and evolution of this major structural feature in the South Atlantic Ocean. All of the samples show a high degree of light rare earth element (LREE) enrichment, quite unlike the flat or depleted patterns normally observed for normal mid-ocean ridge basalts (MORBs). Basalts from Sites 527, 528 and 530 show REE patterns characterised by an arcuate shape and relatively low (Ce/Yb)N ratios (1.46-5.22), and the ratios show a positive linear relationship to Nb content. A different trend is exhibited by the dredged basalts and the basalts from Site 525, and their REE patterns have a fairly constant slope, and higher (Ce/Yb)N ratios (4.31-8.50). These differences are further reflected in the ratios of incompatible trace elements, which also indicate considerable variations within the groups. Mixing hyperbolae for these ratios suggest that simple magma mixing between a 'hot spot' type of magma, similar to present-day volcanics of Tristan da Cunha, and a depleted source, possibly similar to that for magmas being erupted at the Mid-Atlantic Ridge, was an important process in the origin of parts of the Walvis Ridge, as exemplified by Sites 527, 528 and 530. Site 525 and dredged basalts cannot be explained by this mixing process, and their incompatible element ratios suggest either a mantle source of a different composition or some complexity to the mixing process. In addition, the occurrence of different types of basalt at the same location suggests there is vertical zonation within the volcanic pile, with the later erupted basalts becoming more alkaline arid more enriched in incompatible elements. The model proposed for the origin and evolution of the Walvis Ridge involves an initial stage of eruption in which the magma was essentially a mixture of enriched and depleted end-member sources, with the N-MORB component being small. The dredged basalts and Site 525, which represent either later-stage eruptives or those close to the hot spot plume, probably result from mixing of the enriched mantle source with variable amounts and variable low degrees of partial melting of the depleted mantle source. As the volcano leaves the hot spot, these late-stage eruptives continue for some time. The change from tholeiitic to alkalic volcanism is probably related either to evolution in the plumbing system and magma chamber of the individual volcano, or to changes in the depth of origin of the enriched mantle source melt, similar to processes in Hawaiian volcanoes.

Present day and paleo-geographic coordinates, thicknesses and isopach area of the Early Eocene +19 ash layer (Table 1)

23 layers of altered volcanic ash (bentonites) originating from the North Atlantic Igneous Province have been recorded in early Eocene deposits of the Austrian Alps, about 1,900 km away from the source area. The Austrian bentonites are distal equivalents of the ''main ash-phase'' in Denmark and the North Sea basin. We have calculated the total eruption volume of this series as 21,000 km**3, which occurred in 600,000 years. The most powerful single eruption of this series took place 54.0 million years ago (Ma) and ejected ca. 1,200 km**3 of ash material, which makes it one of the largest basaltic pyroclastic eruptions in geological history. The clustering of eruptions must have significantly affected the incoming solar radiation in the early Eocene by the continuous production of stratospheric dust and aerosol clouds. This hypothesis is corroborated by oxygen isotope values, which indicate a global decrease of sea surface temperatures between 1 and 2 C during this major phase of explosive volcanism.

Characterising the three-dimensional ozone distribution of a tidally locked Earth-like planet

We simulate the 3D ozone distribution of a tidally locked Earth-like exoplanet using the high-resolution, 3D chemistry climate model CESM1(WACCM) and study how the ozone layer of a tidally locked Earth (TLE) (ΩTLE = 1/365 days) differs from that of our present-day Earth (PDE) (ΩPDE = 1/1 day). The middle atmosphere reaches a steady state a symptotically within the first 80 days of the simulation. An upwelling, centred on the subsolar point, is present on the day side while a downwelling, centred on the antisolar point, is present on the night side. In the mesosphere, we find similar global ozone distributions for the TLE and the PDE, with decreased ozone on the day side and enhanced ozone on the night side. In the lower mesosphere, a jet stream transitions into a large-scale vortex around a low-pressure system, located at low latitudes of the TLE night side. In the middle stratosphere, the concentration of odd oxygen is approximately equal to that of the ozone [(Ox) ≈ (O3)]. At these altitudes, the lifetime of odd oxygen is ~16 h and the transport processes significantly contribute to the global distribution of stratospheric ozone. Compared to the PDE, where the strong Coriolis force acts as a mixing barrier between low and high latitudes, the transport processes of the TLE are governed by jet streams variable in the zonal and meridional directions. In the middle stratosphere of the TLE, we find high ozone values on the day side, due to the increased production of atomic oxygen on the day side, where it immediately recombines with molecular oxygen to form ozone. In contrast, the ozone is depleted on the night side, due to changes in the solar radiation distribution and the presence of a downwelling. As a result of the reduced Coriolis force, the tropical and extratropical air masses are well mixed and the global temperature distribution of the TLE stratosphere has smaller horizontal gradients than the PDE. Compared to the PDE, the total ozone column global mean is reduced by ~19.3 %. The day side and the night side total ozone column means are reduced by 23.21 and 15.52 %, respectively. Finally, we present the total ozone column (TOC) maps as viewed by a remote observer for four phases of the TLE during its revolution around the star. The mean TOC values of the four phases of the TLE vary by up to 23 %.

Global non-dynamic refinement of radial orbit error for altimetric earth satellites

The accuracy of altimetrically derived oceanographic and geophysical information is limited by the precision of the radial component of the satellite ephemeris. A non-dynamic technique is proposed as a method of reducing the global radial orbit error of altimetric satellites. This involves the recovery of each coefficient of an analytically derived radial error correction through a refinement of crossover difference residuals. The crossover data is supplemented by absolute height measurements to permit the retrieval of otherwise unobservable geographically correlated and linearly combined parameters. The feasibility of the radial reduction procedure is established upon application to the three day repeat orbit of SEASAT. The concept of arc aggregates is devised as a means of extending the method to incorporate longer durations, such as the 35 day repeat period of ERS-1. A continuous orbit is effectively created by including the radial misclosure between consecutive long arcs as an infallible observation. The arc aggregate procedure is validated using a combination of three successive SEASAT ephemerides. A complete simulation of the 501 revolution per 35 day repeat orbit of ERS-1 is derived and the recovery of the global radial orbit error over the full repeat period is successfully accomplished. The radial reduction is dependent upon the geographical locations of the supplementary direct height data. Investigations into the respective influences of various sites proposed for the tracking of ERS-1 by ground-based transponders are carried out. The potential effectiveness on the radial orbital accuracy of locating future tracking sites in regions of high latitudinal magnitude is demonstrated.