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.     

Growth of the Namche Barwa Syntaxis and associated evolution of the Tsangpo Gorge: Constraints from structural and thermochronological data

The eastern syntaxis of the Himalaya, Namche Barwa, is dominated by a north-plunging antiform which began to decompress/grow at approximately 4 Ma. New fission-track analyses on both apatite and zircon, combined with previous geochronological ages, indicate that the Namche Barwa Dome also extended laterally while growing vertically. Zircon fission-track ages range from 17.6 to 0.2 Ma and have a strong relationship to the main faults of the region, including the Tertiary Tsangpo Suture, with the younger ages inside the fault bounds towards the syntaxis core on the Indian Plate and the older ages away from the fault. Apatite ages reveal that the dome has grown laterally and now impinges over the older faulted margin onto the Asian Plate. The dome is traversed by the Tsangpo which has followed the trace of the Suture for over 1300 km from its source to the entrance of the dome near Dania. As the Tsangpo crosses the dome it departs from the Suture but rejoins it some 60 km northeastwards. We construe that the Suture has been displaced by the growing antiform and as a consequence, the antecedent river has been “dragged” in a left-lateral sense along the exhuming north-plunging dome. Restoring the Suture to its position prior to 4 Ma reveals a path of the Tsangpo eastwards across the present southwestern position of the Namche Barwa indentation. This geometric reconstrunction implies that the Tsangpo and the Brahmaputra were always one and the same river. In addition, the Tsangpo was tectonically forced into juxtaposition with a tributary of the Jiali-Parlung which it probably then captured. The capture was due to tectonic forcing, in the last 4 Ma, rather than headward retreat of the paleo-Brahmaputra as has been previously suggested.     

Partitioning of polycyclic aromatic hydrocarbons between water and particles compared to bioaccumulation in mussels: a harbour case

Water and particles from three sites located in Halifax Harbour, near and further away from combined municipal sewage effluents (CMSE) were analysed for parental and alkylated polycyclic aromatic hydrocarbons (par and alkPAH). Bioavailability of PAH was compared for inter-tidal mussels collected at the same sites in April and November. The PAH fingerprint determined over 9 sampling times covering a period of 19 months differed more between phases (water and particles) and seasons, than between sites. In the spring, more alkPAH associated with diesel and gasoline were detected in the soluble phase, along with more bioaccumulation of alkPAH in inter-tidal mussels. A broader number of parPAH were detected in mussels collected in the fall. The mean sum of dissolved alkPAH concentrations was higher in water at the site closer to raw CMSE than at the other two sites and particles of that site. However, lowest bioconcentration factors (BCF) were determined in mussels of this more contaminated site. Similar biota-particle accumulation factors (BPAF) were determined for parPAH in mussels from the three sites, all lower than the BCF of alkPAH. The study indicates that sewage treatment plants will reduce the amount of parPAH especially larger than fluoranthene and pyrene from being deposited in the harbour; that alk naphthalenes and fluorenes present in water will continue to be discharged; that the disturbance of sediments can make particle-bound PAH available to mussels.     

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.     

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.     

Direct observations of aluminosilicate weathering in the hyporheic zone of an Antarctic Dry Valley stream

This study focused on chemical weathering and bacterial ecology in the hyporheic zone of Green Creek, a McMurdo Dry Valley (Antarctica) stream. An in situ microcosm approach was used to observe dissolution features on the basal-plane surface of muscovite mica. Four mica chips were buried in December 1999 and dug up 39 d later. Atomic force microscopy (AFM) of the basal-plane surfaces revealed small, anhedral ∼10-Å-deep etch pits covering ∼4% of the surfaces, from which an approximate basal-plane dissolution rate of 8.3 × 10⁻¹⁸ mol muscovite cm⁻² s⁻¹ was calculated (on the basis of the geometric surface area) for the study period. This is an integrated initial dissolution rate on a fresh surface exposed for a relatively brief period over the austral summer and should not be compared directly to other long-term field rates. The observation of weathering features on mica agrees with previous stream- and watershed-scale studies in the Dry Valleys, which have demonstrated that weathering occurs where liquid water is present, despite the cold temperatures.AFM imaging of mica surfaces revealed biofilms including numerous small (<1 μm long), rounded, oblong bacteria. The AFM observations agreed well with X-ray photoelectron microscopy results showing increased organic C and N. Bacteriologic analysis of the hyporheic zone sediments also revealed <1-μm-long bacteria. α-Proteobacteria were observed, consistent with the oligotrophic conditions of the hyporheic zone. Nitrate-reducing bacteria were found, in agreement with a previous tracer test at Green Creek that suggested nitrate reduction occurs in the hyporheic zone. The results of this study thus provide direct evidence of dynamic geochemical and microbial processes in the hyporheic zone of a Dry Valley stream despite the extreme conditions; such processes were inferred previously from stream-scale hydrogeochemical studies.