Earth gravity field recovery using GPS, GLONASS, and SLR satellites

The time variable Earth’s gravity field provides the information about mass transport within the system Earth, i.e., the relationship of mass transport between atmosphere, oceans, and land hydrology. We recover the low-degree parameters of the time variable gravity field using microwave observations from GPS and GLONASS satellites and from SLR data to five geodetic satellites, namely LAGEOS-1/2, Starlette, Stella, and AJISAI. GPS satellites are particularly sensitive to specific coefficients of the Earth's gravity field, because of the deep 2:1 orbital resonance with Earth rotation (two revolutions of the GPS satellites per sidereal day). The resonant coefficients cause, among other, a “secular” drift (actually periodic variations of very long periods) of the semi-major axes of up to 5.3 m/day of GPS satellites. We processed 10 years of GPS and GLONASS data using the standard orbit models from the Center of Orbit Determination in Europe (CODE) with a simultaneous estimation of the Earth gravity field coefficients and other parameters, e.g., satellite orbit parameters, station coordinates, Earth rotation parameters, troposphere delays, etc. The weekly GNSS gravity solutions up to degree and order 4/4 are compared to the weekly SLR gravity field solutions. The SLR-derived geopotential coefficients are compared to monthly GRACE and CHAMP results.

The Response of a Small Stream in the Lesni Potok Forested Catchment, Central Czech Republic, to a Short-Term In-Stream Acidification

Lesni Potok stream drains a forested headwater catchment in the central Czech Republic. It was artificially acidified with hydrochloric acid (HCl) for four hours to assess the role of stream substrate in acid-neutralisation and recovery. The pH was lowered from 4.7 to 3.2. Desorption of Ca and MP and desorption or solution of Al dominated acid-neutralisation; Al mobilisation was more important later. The stream substrate released 4.542 meq Ca, 1, 184 meq Mg, and 2,329 meq Al over a 45 in long and I in wide stream segment, smaller amounts of Be. Cd, Fe, and Mn were released. Adsorption of SO42- and desorption of F- occurred during the acidification phase of the experiment. The exchange reactions were rapidly reversible for Ca, Mg and SO42- but not symmetric as the substrate resorbed 1083, 790 and 0 meq Ca, Mg, and Al. respectively, in a 4-hour recovery period. Desorption of SO42- occurred during the resorption of Ca and Mg. These exchange and dissolution reactions delay acidification, diminish the pH depression and retard recovery from episodic acidification. The behaviour of the stream substrate-water interaction resembles that for soil-soil water interactions. A mathematical dynamic mass-balance based model, MASS (Modelling Acidification of Stream Sediments), was developed which simulates the adsorption and desorption of base cations during the experiment and was successfully calibrated to the experimental data.

Time variable Earth’s gravity field from SLR satellites

The time variable Earth’s gravity field contains information about the mass transport within the system Earth, i.e., the relationship between mass variations in the atmosphere, oceans, land hydrology, and ice sheets. For many years, satellite laser ranging (SLR) observations to geodetic satellites have provided valuable information of the low-degree coefficients of the Earth’s gravity field. Today, the Gravity Recovery and Climate Experiment (GRACE) mission is the major source of information for the time variable field of a high spatial resolution. We recover the low-degree coefficients of the time variable Earth’s gravity field using SLR observations up to nine geodetic satellites: LAGEOS-1, LAGEOS-2, Starlette, Stella, AJISAI, LARES, Larets, BLITS, and Beacon-C. We estimate monthly gravity field coefficients up to degree and order 10/10 for the time span 2003–2013 and we compare the results with the GRACE-derived gravity field coefficients. We show that not only degree-2 gravity field coefficients can be well determined from SLR, but also other coefficients up to degree 10 using the combination of short 1-day arcs for low orbiting satellites and 10-day arcs for LAGEOS-1/2. In this way, LAGEOS-1/2 allow recovering zonal terms, which are associated with long-term satellite orbit perturbations, whereas the tesseral and sectorial terms benefit most from low orbiting satellites, whose orbit modeling deficiencies are minimized due to short 1-day arcs. The amplitudes of the annual signal in the low-degree gravity field coefficients derived from SLR agree with GRACE K-band results at a level of 77 %. This implies that SLR has a great potential to fill the gap between the current GRACE and the future GRACE Follow-On mission for recovering of the seasonal variations and secular trends of the longest wavelengths in gravity field, which are associated with the large-scale mass transport in the system Earth.

Sediment mobilization deposits from episodic subsurface fluid flow—A new tool to reveal long-term earthquake records?

Subsurface fluid flow can be affected by earthquakes; increased spring activity, mud vol- cano eruptions, groundwater fluctuations, changes in geyser frequency, and other forms of altered subsurface fluid flow have been documented during, after, or even prior to seismic shaking. Recently discovered giant pockmarks on the bottom of Lake Neuchâtel, Switzerland, are the lake-floor expression of subsurface fluid flow. They discharge groundwater from the Jura Mountains karstic aquifers and experience episodically increased subsurface fluid flow documented by subsurface sediment mobilization deposits at the levees of the pockmarks. In this study, we present the spatio-temporal distribution of event deposits from these phases of sediment expulsion and of multiple time-correlative mass-transport deposits. We report five striking instances of concurrent multiple subsurface sediment deposits and multiple mass- transport deposits since late glacial times, for which we propose past earthquakes as a trigger. Comparison of this new event catalogue with historic earthquakes and other independent paleoseismic records suggests that initiation of sediment expulsion requires a minimum mac- roseismic intensity of VII. Thus, our study presents for the first time sedimentary deposits resulting from increased subsurface fluid flow as a paleoseismic proxy.

Gravity potential spherical harmonic series

Time variable gravity fields, reflecting variations of mass distribution in the system Earth is one of the key parameters to understand the changing Earth. Mass variations are caused either by redistribution of mass in, on or above the Earth's surface or by geophysical processes in the Earth's interior. The first set of observations of monthly variations of the Earth gravity field was provided by the US/German GRACE satellite mission beginning in 2002. This mission is still providing valuable information to the science community. However, as GRACE has outlived its expected lifetime, the geoscience community is currently seeking successor missions in order to maintain the long time series of climate change that was begun by GRACE. Several studies on science requirements and technical feasibility have been conducted in the recent years. These studies required a realistic model of the time variable gravity field in order to perform simulation studies on sensitivity of satellites and their instrumentation. This was the primary reason for the European Space Agency (ESA) to initiate a study on ''Monitoring and Modelling individual Sources of Mass Distribution and Transport in the Earth System by Means of Satellites''. The goal of this interdisciplinary study was to create as realistic as possible simulated time variable gravity fields based on coupled geophysical models, which could be used in the simulation processes in a controlled environment. For this purpose global atmosphere, ocean, continental hydrology and ice models were used. The coupling was performed by using consistent forcing throughout the models and by including water flow between the different domains of the Earth system. In addition gravity field changes due to solid Earth processes like continuous glacial isostatic adjustment (GIA) and a sudden earthquake with co-seismic and post-seismic signals were modelled. All individual model results were combined and converted to gravity field spherical harmonic series, which is the quantity commonly used to describe the Earth's global gravity field. The result of this study is a twelve-year time-series of 6-hourly time variable gravity field spherical harmonics up to degree and order 180 corresponding to a global spatial resolution of 1 degree in latitude and longitude. In this paper, we outline the input data sets and the process of combining these data sets into a coherent model of temporal gravity field changes. The resulting time series was used in some follow-on studies and is available to anybody interested.

137Cs concentrations in the southeastern Baltic Sea in 1997-2000

Variations of 137Cs concentration in the southeastern Baltic Sea were investigated over the period 1997-2000, i.e. in 11-14 years after the Chernobyl Nuclear Power Plant accident. Rate of "self-cleaning" proved to be very slow. Some results obtained in 1999 were almost the same as those measured after the accident, in 1986. Calculated results showed that "Chernobyl" caesium-137 would be "cleaned" in the Baltic Sea by 2020-2022. In 2000 average concentration had to be about 50-60 Bq/m**3. Sometimes mentioned concentrations were observed. In some cases higher concentrations averaging from 67 to 80 Bq/m**3 were registered in the southeastern Baltic Sea in 1999; and in some samples 137Cs concentrations were very high. They varied from 110 to 212 Bq/m**3. No steady correlation was observed between 137Cs concentration, salinity and temperature in surface water of the area. Distribution of radionuclide concentration sometimes depends on direction of water mass transport. Abnormally high concentrations of 137Cs in the southeastern Baltic Sea may result from additional radioactive waste discharge.

Physical properties and geochemical pore-water analysis of sediment core GeoB16423-1

Several studies indicate that the 2011 Tohoku-Oki earthquake (Mw 9.0) off the Pacific coast of Japan has induced slip to the trench and triggered landslides in the Japan Trench. In order to better understand these processes, detailed mapping and shallow-coring landslides at the trench as well as Integrated Ocean Drilling Program (IODP) deep drilling to recover the plate boundary décollement (Japan Trench Fast Earthquake Drilling Project, JFAST) have been conducted. In this study we report sediment core data from the rapid response R/V SONNE cruise (SO219A) to the Japan Trench, evidencing a Mass Transport Deposit (MTD) in the uppermost section later drilled at this JFAST-site during IODP Expedition 343. A 8.7 m long gravity core (GeoB16423-1) recovered from ~7,000 m water depth reveals a 8 m sequence of semi-consolidated mud clast breccias embedded in a distorted chaotic sediment matrix. The MTD is covered by a thin veneer of 50 cm hemipelagic, bioturbated diatomaceous mud. This stratigraphic boundary can be clearly distinguished by using physical properties data from Multi Sensor Core Logging and from fall-cone penetrometer shear strength measurements. The geochemical analysis of the pore-water shows undisturbed linear profiles measured from the seafloor downcore across the stratigraphic contact between overlying younger background-sediment and MTD below. This indicates that the investigated section has not been affected by a recent sediment destabilization in the course of the giant Tohoku-Oki earthquake event. Instead, we report an older landslide which occurred between 700 and 10,000 years ago, implying that submarine mass movements are dominant processes along the Japan Trench. However, they occur on local sites and not during each megathrust earthquake.

Geochemistry of frost flowers and snow obtained near Ahlmann summit, Vestfonna ice cap

The chemistry of snow and ice cores from Svalbard is influenced by variations in local sea ice margin and distance to open water. Snow pits sampled at two summits of Vestfonna ice cap (Nordaustlandet, Svalbard), exhibit spatially heterogeneous soluble ions concentrations despite similar accumulation rates, reflecting the importance of small-scale weather patterns on this island ice cap. The snow pack on the western summit shows higher average values of marine ions and a winter snow layer that is relatively depleted in sulphate. One part of the winter snow pack exhibits a [SO4-/Na+] ratio reduced by two thirds compared with its ratio in sea water. This low sulphate content in winter snow is interpreted as the signature of frost flowers, which are formed on young sea ice when offshore winds predominate. Frost flowers have been described as the dominant source of sea salt to aerosol and precipitation in ice cores in coastal Antarctica but this is the first time their chemical signal has been described in the Arctic. The eastern summit does not show any frost flower signature and we interpret the unusually dynamic ice transport and rapid formation of thin ice on the Hinlopen Strait as the source of the frost flowers.

A numerical discrete element study on the influence of an embedded viscoelastic-plastic layer at different viscosity values

The major aim of this study was to examine the influence of an embedded viscoelastic-plastic layer at different viscosity values on accretionary wedges at subduction zones. To quantify the effects of the layer viscosity, we analysed the wedge geometry, accretion mode, thrust systems and mass transport pattern. Therefore, we developed a numerical 2D 'sandbox' model utilising the Discrete Element Method. Starting with a simple pure Mohr Coulomb sequence, we added an embedded viscoelastic-plastic layer within the brittle, undeformed 'sediment' package. This layer followed Burger's rheology, which simulates the creep behaviour of natural rocks, such as evaporites. This layer got thrusted and folded during the subduction process. The testing of different bulk viscosity values, from 1 × 10**13 to 1 × 10**14 (Pa s), revealed a certain range where an active detachment evolved within the viscoelastic-plastic layer that decoupled the over- and the underlying brittle strata. This mid-level detachment caused the evolution of a frontally accreted wedge above it and a long underthrusted and subsequently basally accreted sequence beneath it. Both sequences were characterised by specific mass transport patterns depending on the used viscosity value. With decreasing bulk viscosities, thrust systems above this weak mid-level detachment became increasingly symmetrical and the particle uplift was reduced, as would be expected for a salt controlled forearc in nature. Simultaneously, antiformal stacking was favoured over hinterland dipping in the lower brittle layer and overturning of the uplifted material increased. Hence, we validated that the viscosity of an embedded detachment strongly influences the whole wedge mechanics, both the respective lower slope and the upper slope duplex, shown by e.g. the mass transport pattern.

Lithology, stratigraphy, organic carbon content, extract and liquid chromatography yields, and Rock-Eval pyrolysis at DSDP Holes 75-530A and 75-532

Seventeen sediment samples of Albian-Cenomanian to early Pliocene age from DSDP Hole 530A in the Angola Basin and six sediment samples of early Pliocene to late Pleistocene age from the Walvis Ridge were investigated by organic geochemical methods, including organic carbon determination, Rock-Eval pyrolysis, gas chromatography and combined gas chromatography/mass spectrometry of extractable hydrocarbons, and kerogen microscopy. The organic matter in all samples is strongly influenced by a terrigenous component from the nearby continent. The amount of marine organic matter present usually increases with the total organic carbon content, which reaches an extreme value of more than 10% in a Cenomanian black shale from Hole 530A. At Site 530 the extent of preservation of organic matter in the deep sea sediments is related to mass transport down the continental slope, whereas the high organic carbon contents in the sediments from Site 532 reflect both high bioproductivity in the Benguela upwelling regime and considerable supply of terrigenous organic matter. The maturation level of the organic matter is low in all samples.

Basal dolomitic sediments ofthe Tyrrhenian Sea

Basal dolomitic sediments were recovered at three drill sites in the Tyrrhenian Sea during Ocean Drilling Program (ODP) Leg 107 (Sites 650, 651, and 655). These sediments overlie the basaltic basement complex and are enriched in iron, and in some instances, also in manganese. The manganese enrichments, together with a very slight enrichment in trace transition elements, strongly suggest that the basal sediments have an affinity to deep-sea metalliferous deposits of hydrothermal origin. At Sites 651 and 655, the dolostones contain variable amounts of authigenic palygorskite, a Mgrich clay mineral. At Site 651, the basal sediments are 40 m thick and contain nonstoichiometric dolomite, sometimes Ca rich, but primarily Mg rich. The occurrence of Mg-rich dolomite with excess Mg up to 4% is unusual for the deep-sea environment; it may be associated with a hydrothermally driven flux of altered sea water through the directly underlying basement complex, which comprises basalt, dolerite, and serpentinized peridotite. Low-temperature alteration of the basement complex could produce solutions enriched in Mg. Oxygen-isotope equilibrium temperatures indicate that all of the studied dolomites formed under low-temperature conditions (i.e., < 70?C). The carbon-isotope compositions, together with the strong isotopic covariance, suggest that the Mg-rich dolomite precipitated more rapidly than the Carich dolomite. We suggest that the low-temperature, hydrothermal convection of Mg-rich solutions through the basal sediments in this back-arc basin environment (1) overcame kinetic problems related to the formation of massive dolostones, and (2) provided a mass-transport mechanism for dolomitization.

Using ground-penetrating radar, topography and classification of vegetation to model the sediment and active layer thickness in a periglacial lake catchment, Western Greenland, link to shapefiles

The geometries of a catchment constitute the basis for distributed physically based numerical modeling of different geoscientific disciplines. In this paper results from ground-penetrating radar (GPR) measurements, in terms of a 3D model of total sediment thickness and active layer thickness in a periglacial catchment in western Greenland, is presented. Using the topography, thickness and distribution of sediments is calculated. Vegetation classification and GPR measurements are used to scale active layer thickness from local measurements to catchment scale models. Annual maximum active layer thickness varies from 0.3 m in wetlands to 2.0 m in barren areas and areas of exposed bedrock. Maximum sediment thickness is estimated to be 12.3 m in the major valleys of the catchment. A method to correlate surface vegetation with active layer thickness is also presented. By using relatively simple methods, such as probing and vegetation classification, it is possible to upscale local point measurements to catchment scale models, in areas where the upper subsurface is relatively homogenous. The resulting spatial model of active layer thickness can be used in combination with the sediment model as a geometrical input to further studies of subsurface mass-transport and hydrological flow paths in the periglacial catchment through numerical modelling.