A plume tectonics model for the Tharsis province, Mars

Morphological and structural data from the whole Tharsis province suggest that a number of shallow grabens radially oriented about the Tharsis bulge on Mars are underlain by dykes, which define giant radiating swarms similar to, e.g. the Mackenzie dyke swarm of the Canadian shield. Mechanisms for graben formation are proposed, and the depth, width, and height of the associated dykes are estimated. Structural mapping leads to define successive stages of dyke emplacement, and provide stress-trajectory maps that indicate a steady source of the regional stress during the whole history of the Tharsis province. A new tectonic model of Tharsis is presented, based on an analogy with dyke swarms on the Earth that form inside hot spots. This model successfully matches the following features: (1) the geometry of the South Tharsis Ridge Belt, which may have been a consequence of the compressional stress field at the boundary between the uplifted and non-uplifted areas in the upper part of the lithosphere at the onset of hot spot activity; (2) extensive lava flooding, interpreted as a consequence of the high thermal anomaly at the onset of plume (hot spot) activity; (3) wrinkle ridge geometry in the Tharsis hemisphere, the formation of which is interpreted as a consequence of buoyant subsidence of the brittle crust in response to the lava load; (4) Valles Marineris limited stretching by preliminary passive rifting, and uplift, viewed as a necessary consequence of adiabatic mantle decompression induced by stretching. The geometrical analysis of dyke swarms suggests the existence of a large, Tharsis-independent extensional state of stress during all the period of tectonic activity, in which the minimum compressive stress is roughly N---S oriented. Although magmatism must have loaded the lithosphere significantly after the plume activity ceased and be responsible for additional surface deformations, there is no requirement for further loading stress to explain surficial features. Comparison with succession of magmatic and tectonic events related to hot spots on the Earth suggests that the total time required to produce all the surface deformation observed in the Tharsis province over the last 3.8 Ga does probably not exceed 10 or 15 Ma.     

Budgets of Mn,Cd and Cu in the macrotidal Gironde estuary (SW France)

The behavior and budget of Mn, Cd and Cu in the Gironde estuary were investigated through data from both the water column (WC) and sediment depth profiles. In the estuarine freshwater reaches, Mn and Cd removal from and Cu addition to the dissolved phase occurs with a magnitude equivalent to 10%, 30% and 25% of their respective annual fluvial gross dissolved input, respectively. In the saline estuary, diffusive benthic outflow is the main source of dissolved Mn (74% of the total gross dissolved input within the estuary) to the WC. In contrast, Cd (96%) and Cu (89%) are mainly released into the dissolved phase of the WC from fluvial, estuarine and dredging-related particles through complexation (Cd) and organic carbon mineralization (Cu). Anthropogenic activities (sediment dredging) induce pore water inputs, particulate sulfide oxidation and sediment resuspension, significantly contributing to the metal budget of the WC. The related amounts of metals released could be equivalent to 20–50% (Cd) and up to 70% (Cu) of their respective net dissolved addition. Mass balances suggest that a large part of the metals previously released into the dissolved phase from processes within the estuary are removed by suspended particles due to (co-)precipitation of Fe/Mn (oxy)hydroxides and scavenging on autochthonous organic matter. On an annual basis, the Gironde estuary acts as a net sink of dissolved Mn, removing 60% of the dissolved fluvial inputs, and as a net source of dissolved Cd and Cu, contributing ∼ 85% and 20–45% to the dissolved Cd and Cu fluxes to the ocean.     

A linear Boltzmann equation to model wave scattering in the marginal ice zone

We present a linear Boltzmann equation to model wave scattering in the Marginal Ice Zone (the region of ocean which consists of broken ice floes). The equation is derived by two methods, the first based on Meylan et al. [Meylan, M.H., Squire, V.A., Fox, C., 1997. Towards realism in modeling ocean wave behavior in marginal ice zones. J. Geophys. Res. 102 (C10), 22981–22991] and second based on Masson and LeBlond [Masson, D., LeBlond, P., 1989. Spectral evolution of wind-generated surface gravity waves in a dispersed ice field. J. Fluid Mech. 202, 111–136]. This linear Boltzmann equation, we believe, is more suitable than the equation presented in Masson and LeBlond [Masson, D., LeBlond, P., 1989. Spectral evolution of wind-generated surface gravity waves in a dispersed ice field. J. Fluid Mech. 202, 111–136] because of its simpler form, because it is a differential rather than difference equation and because it does not depend on any assumptions about the ice floe geometry. However, the linear Boltzmann equation presented here is equivalent to the equation in Masson and LeBlond [Masson, D., LeBlond, P., 1989. Spectral evolution of wind-generated surface gravity waves in a dispersed ice field. J. Fluid Mech. 202, 111–136] since it is derived from their equation. Furthermore, the linear Boltzmann equation is also derived independently using the argument in Meylan et al. [Meylan, M.H., Squire, V.A., Fox, C., 1997. Towards realism in modeling ocean wave behavior in marginal ice zones. J. Geophys. Res. 102 (C10), 22981–22991]. We also present details of how the scattering kernel in the linear Boltzmann equation is found from the scattering by an individual ice floe and show how the linear Boltzmann equation can be solved straightforwardly in certain cases.     

The GeoForschungsZentrum Potsdam/Groupe de Recherche de Gèodésie Spatiale satellite-only and combined gravity field models: EIGEN-GL04S1 and EIGEN-GL04C

The recent improvements in the Gravity Recovery And Climate Experiment (GRACE) tracking data processing at GeoForschungsZentrum Potsdam (GFZ) and Groupe de Recherche de Géodésie Spatiale (GRGS) Toulouse, the availability of newer surface gravity data sets in the Arctic, Antarctica and North-America, and the availability of a new mean sea surface height model from altimetry processing at GFZ gave rise to the generation of two new global gravity field models. The first, EIGEN-GL04S1, a satellite-only model complete to degree and order 150 in terms of spherical harmonics, was derived by combination of the latest GFZ Potsdam GRACE-only (EIGEN-GRACE04S) and GRGS Toulouse GRACE/LAGEOS (EIGEN-GL04S) mean field solutions. The second, EIGEN-GL04S1 was combined with surface gravity data from altimetry over the oceans and gravimetry over the continents to derive a new high-resolution global gravity field model called EIGEN-GL04C. This model is complete to degree and order 360 and thus resolves geoid and gravity anomalies at half- wavelengths of 55 km at the equator. A degree-dependent combination method has been applied in order to preserve the high accuracy from the GRACE satellite data in the lower frequency band of the geopotential and to form a smooth transition to the high-frequency information coming from the surface data. Compared to pre-CHAMP global high-resolution models, the accuracy was improved at a spatial resolution of 200 km (half-wavelength) by one order of magnitude to 3 cm in terms of geoid heights. The accuracy of this model (i.e. the commission error) at its full spatial resolution is estimated to be 15 cm. The model shows a reduced artificial meridional striping and an increased correlation of EIGEN-GL04C-derived geostrophic meridional currents with World Ocean Atlas 2001 (WOA01) data. These improvements have led to select EIGEN-GL04C for JASON-1 satellite altimeter data reprocessing. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Ocean tide loading (OTL) displacements from global and local grids: comparisons to GPS estimates over the shelf of Brittany, France

In this paper we examine OTL displacements detected by GPS stations of a dedicated campaign and validate ocean tide models. Our area of study is the continental shelf of Brittany and Cotentin in France. Brittany is one of the few places in the world where tides provoke loading displacements of ∼10–12 cm vertically and a few cm horizontally. Ocean tide models suffer from important discrepancies in this region. Seven global and regional ocean tide models were tested: FES2004 corrected for K2, TPXO.7.0, TPXO.6.2, GOT00.2, CSR4.0, NAO.99b and the most recent regional grids of the North East Atlantic (NEA2004). These gridded amplitudes and phases of ocean tides were convolved in order to get the predicted OTL displacements using two different algorithms. Data over a period of 3.5 months of 8 GPS campaign stations located on the north coast of Brittany are used, in order to evaluate the geographical distribution of the OTL effect. We have modified and implemented new algorithms in our GPS software, GINS 7.1. GPS OTL constituents are estimated based on 1-day batch solutions. We compare the observed GPS OTL constituents of M₂, S₂, N₂ and K₁ waves with the selected ocean tide models on global and regional grids. Large phase-lag and amplitude discrepancies over 20° and 1.5 cm in the vertical direction in the semi-diurnal band of M₂ between predictions and GPS/models are detected in the Bay of Mont St-Michel. From a least squares spectral analysis of the GPS time-series, significant harmonic peaks in the integer multiples of the orbital periods of the GPS satellites are observed, indicating the existence of multipath effects in the GPS OTL constituents. The GPS OTL observations agree best with FES2004, NEA2004, GOT00.2 and CSR4.0 tide models.     

Numerical simulation of a stratigraphic model

A multi-lithology diffusive stratigraphic model is considered, which simulates at large scales in space and time the infill of sedimentary basins governed by the interaction between tectonics displacements, eustatic variations, sediment supply, and sediment transport laws. The model accounts for the mass conservation of each sediment lithology resulting in a mixed parabolic, hyperbolic system of partial differential equations (PDEs) for the lithology concentrations and the sediment thickness. It also takes into account a limit on the rock alteration velocity modeled as a unilaterality constraint. To obtain a robust, fast, and accurate simulation, fully and semi-implicit finite volume discre tization schemes are derived for which the existence of stable solutions is proved. Then, the set of nonlinear equations is solved using a Newton algorithm adapted to the unilaterality constraint, and preconditioning strategies are defined for the solution of the linear system at each Newton iteration. They are based on an algebraic approximate decoupling of the sediment thickness and the concentration variables as well as on a proper preconditioning of each variable. These algorithms are studied and compared in terms of robustness, scalability, and efficiency on two real basin test cases.     

Gravity field determination using boundary element methods

The Boundary Element Method (BEM), a numerical technique for solving boundary integral equations, is introduced to determine the earth's gravity field. After a short survey on its main principles, we apply this method to the fixed gravimetric boundary value problem (BVP), i.e. the determination of the earth's gravitational potential from measurements of the intensity of the gravity field in points on the earth's surface. We show how to linearize this nonlinear BVP using an implicit function theorem and how to transform the linearized BVP into a boundary integral equation using the single layer representation. A Galerkin method is used to transform the boundary integral equation using the single layer representation. A Galerkin method is used to transform the boundary integral equation into a linear system of equations. We discuss the major problems of this approach for setting up and solving the linear system. The BVP is numerically solved for a bounded part of the earth's surface using a high resolution reference gravity model, measured gravity values of high density, and a 50 50 m² digital terrain model to describe the earth's surface. We obtain a gravity field resolution of 1 1 km² with an accuracy of the order 10^–3 to 10^–4 in about 1 CPU-hour on a Siemens/Fujitsu SIMD vector pipeline machine using highly sophisticated numerical integration techniques and fast equation solvers. We conclude that BEM is a powerful numerical tool for solving boundary value problems and may be an alternative to classical geodetic techniques.     

Intercomparison of measurements of stratospheric hydrogen fluoride

Observations of the vertical profile of hydrogen fluoride (HF) vapor in the stratosphere and of the vertical column amounts of HF above certain altitudes were made using a variety of spectroscopic instruments in the 1982 and 1983 Balloon Intercomparison Campaigns. Both emission instruments working in the far infrared spectral region and absorption instruments using solar occultation in the 2.5m region were employed. No systematic differences were seen in results from the two spectral regions. A mean profile from 20–45 km is presented, with uncertainties ranging from 20% to 50%. Total columns measured from ground and from 12 km are consistent with the profile if the mixing ratio for HF is small in the tropophere and low stratosphere.     

Late Holocene sea-level changes and isostasy in western Denmark

Cores and exposed cliff sections in salt marshes around Ho Bugt, a tidal embayment in the northernmost part of the Danish Wadden Sea, were subjected to ¹⁴C dating and litho- and biostratigraphical analyses to reconstruct paleoenvironmental changes and to establish a late Holocene relative sea-level history. Four stages in the late Holocene development of Ho Bugt can be identified: (1) groundwater-table rise and growth of basal peat (from at least 2300 BC to AD 0); (2) salt-marsh formation (0 to AD 250); (3) a freshening phase (AD 250 to AD 1600?), culminating in the drying out of the marshes and producing a distinct black horizon followed by an aeolian phase with sand deposition; and (4) renewed salt-marsh deposition (AD 1600? to present). From 16 calibrated AMS radiocarbon ages on fossil plant fragments and 4 calibrated conventional radiocarbon ages on peat, we reconstructed a local relative sea-level history that shows a steady sea-level rise of 4 m since 4000 cal yr BP. Contrary to suggestions made in the literature, the relative sea-level record of Ho Bugt does not contain a late Holocene highstand. Relative sea-level changes at Ho Bugt are controlled by glacio-isostatic subsidence and can be duplicated by a glacial isostatic adjustment model in which no water is added to the world's oceans after ca. 5000 cal yr BP.