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.     

Probing the physics of interacting galaxies

The morphological and velocity structures in the gaseous (HI and CO) and stellar components of two interacting systems are examined. Both Arp 140 and Arp 104 reveal extended tidal tails in the HI. The Hα and FIR fluxes of Arp 140 yield similar SFR of ∼ 0.8 M_⊙ yr^-1. In contrast the Hα flux of Arp 104 yields a SFR of ∼ 0.05 M_⊙ yr^-1, ∼ 20 times smaller than that obtained from the FIR flux. Spectra were used to examine the changing velocity of atomic and molecular gas in NGC 5218 (Arp 104). The atomic and molecular gas were found to be dynamically similar with comparable velocities and velocity widths across the galaxy; consistent with the two phases responding similarly to the interaction, or enhanced HI to CO conversion in the centre of the galaxy. This revised version was published online in August 2006 with corrections to the Cover Date.     

Resonant behaviour of MHD waves on magnetic flux tubes

Resonant absorption of MHD waves on a nonuniform flux tube is investigated as a driven problem for a 1D cylindrical equilibrium. The variation of the fractional absorption is studied as a function of the frequency and its relation to the eigenvalue problem of the MHD radiating eigenmodes of the nonuniform flux tube is established. The optimal frequencies producing maximal fractional absorption are determined and the condition for total absorption is obtained. This condition defines an impedance matching and is fulfilled for an equilibrium that is fine tuned with respect to the incoming wave. The variation of the spatial wave solutions with respect to the frequency is explained as due to the variation of the real and imaginary parts of the dispersion relation of the MHD radiating eigenmodes with respect to the real driving frequency.     

The carbon isotopic response of algae, (cyano)bacteria, archaea and higher plants to the late Cenomanian perturbation of the global carbon cycle: Insights from biomarkers in black shales from the Cape Verde Basin (DSDP Site 367)

The stable carbon isotopic compositions of free and sulfur (S)-bound biomarkers derived from algae, (cyano)bacteria, archaea and higher plants and total organic carbon (TOC) during the first phase of the late Cenomanian/Turonian oceanic anoxic event (OAE) were measured in black shales deposited in the southern proto-Atlantic Ocean in the Cape Verde basin (DSDP Site 367) to determine the response of these organisms to this major perturbation of the global carbon cycle resulting from widespread burial of marine organic matter. The average positive isotope excursions of TOC and biomarkers varied from 5.1‰ to 8.3‰. The δ¹³C values were cross correlated to infer potential common sources of biomarkers. This revealed common sources for C₃₁ and C₃₂ hopanes but no 1:1 relationship for pristane and phytane. The correlation of δ¹³C_TOC with the δ¹³C value of sulfur (S)-bound phytane is the strongest. This is because S-bound phytane is derived from phytol that originates from all marine primary producers (algae and cyanobacteria) and thus represents a weighted average of their carbon isotopic compositions. The δ¹³C values of S-bound phytane and C₃₅ hopane were also used to estimate pCO₂ levels. Before the OAE burial event, pCO₂ levels are estimated to be ca. 1300 ppmv using both biomarkers and the independent maximum Rubisco fractionation factors. At times of maximum organic carbon burial rates during the OAE, reconstructed pCO₂ levels are estimated to be ca. 700 ppmv. However, compared to other C/T OAE sections the positive isotope excursion of S-bound phytane is also affected by an increased production during the OAE. When we compensate for this, we arrive at pCO₂ levels around 1000 ppmv, a reduction of ca. 25%. This indicates that burial of organic matter can have a large effect on atmospheric CO₂ levels.     

Evaluation of approaches to calculate debris-flow parameters for hazard assessment

Many different runout prediction methods can be applied to estimate the mobility of future debris flows during hazard assessment. The present article reviews the empirical, analytical, simple flow routing and numerical techniques. All these techniques were applied to back-calculate a debris flow, which occurred in 1982 at La Guingueta catchment, in the Eastern Pyrenees. A sensitivity analysis of input parameters was carried out, while special attention was paid to the influence of rheological parameters. We used the Voellmy fluid rheology for our analytical and numerical modelling, since this flow resistance law coincided best with field observations. The simulation results indicated that the “basal” friction coefficients rather affect the runout distance, while the “turbulence” terms mainly influence flow velocity. A comparison of the velocity computed on the fan showed that the analytical model calculated values similar to the numerical ones. The values of our rheological parameters calibrated at La Guingueta agree with data back-calculated for other debris flows. Empirical relationships represent another method to estimate total runout distance. The results confirmed that they contain an important uncertainty and they are strictly valid only for the conditions, which were the basis for their development. With regards to the simple flow routing algorithm, this methods could satisfactorily simulate the total area affected by the 1982 debris flow, but it was not able to directly calculate total runout distance and velocity. Finally, a suggestion on how different runout prediction methods can be applied to generate debris-flow hazard maps is presented. Taking into account the definition of hazard and intensity, the best choice would be to divide the resulting hazard maps into two types: “final hazard maps” and “preliminary hazard maps”. Only the use of numerical models provided final hazard maps, because they could incorporate different event magnitudes and they supplied output-values for intensity calculation. In contrast, empirical relationships and flow routing algorithms, or a combination of both, could be applied to create preliminary hazard maps. The present study only focussed on runout prediction methods. Other necessary tasks to complete the hazard assessment can be looked up in the “Guidelines for landslide susceptibility, hazard and risk zoning” included in this Special Issue.     

Determination of fluid/melt partition coefficients by LA-ICPMS analysis of co-existing fluid and silicate melt inclusions: Controls on element partitioning

Analyses of co-existing silicate melt and fluid inclusions, entrapped in quartz crystals in volatile saturated magmatic systems, allowed direct quantitative determination of fluid/melt partition coefficients. Investigations of various granitic systems (peralkaline to peraluminous in composition, log fO₂ = NNO−1.7 to NNO+4.5) exsolving fluids with various chlorinities (1-14 mol/kg) allowed us to assess the effect of these variables on the fluid/melt partition coefficients (D). Partition coefficients for Pb, Zn, Ag and Fe show a nearly linear increase with the chlorinity of these fluid (D_Pb ∼ 6 ∗ m_Cl, D_Zn ∼ 8 ∗ m_Cl, D_Ag ∼ 4 ∗ m_Cl, D_Fe ∼ 1.4 ∗ m_Cl, where m_Cl is the molinity of Cl). This suggests that these metals are dissolved primarily as Cl-complexes and neither oxygen fugacity nor the composition of the melt affects significantly their fluid/melt partitioning. By contrast, partition coefficients for Mo, B, As, Sb and Bi are highest in low salinity (1-2 mol/kg Cl) fluids with maximum values of D_Mo ∼ 20, D_B ∼ 15, D_As ∼ 13, D_Sb ∼ 8, D_Bi ∼ 15 indicating dissolution as non-chloride (e.g., hydroxy) complexes. Fluid/melt partition coefficients of copper are highly variable, but highest between vapor like fluids and silicate melt (D_Cu ? 2700), indicating an important role for ligands other than Cl. Partition coefficients for W generally increase with increasing chlorinity, but are exceptionally low in some of the studied brines which may indicate an effect of other parameters. Fluid/melt partition coefficients of Sn show a high variability but likely increase with the chlorinity of the fluid (D_Sn = 0.3-42, D_W = 0.8-60), and decrease with decreasing oxygen fugacity or melt peraluminosity.     

Ionospheric conditions

A general analysis of ionospheric conditions has been made in the light of possible ionic reactions occurring in the upper atmosphere. Data obtained on various parameters, such as ionic production and recombination, show that precise knowledge of the spectral distribution of solar radiation is needed and that other experimental determinations on dissociative recombinations are required.

The ionic complexity of the ionosphere is underlined by describing how the atomic ions O⁺ and N⁺ react with N₂, O₂ and NO molecules. The behavior of the molecular ions N⁺₂, O⁺₂and NO⁺depends on a group of simultaneous processes involving charge transfers and ionatom interchanges which are more important than dissociative recombinations. The altitude distribution of ions is exemplified by discussing the relative importance of various loss coefficients in the D-, E- and F-regions. It is seen that molecular nitrogen ions are subject to important charge transfer processes, that nitric oxide ions are always final products destroyed only by dissociative recombination. Additionally, the entire production of atomic oxygen ions is related to the photoionization of molecular nitrogen. Some information is also given on possible anomalies in the ratio of O⁺₂ and NO⁺ densities in the lower ionosphere. From the lack of sufficient experimental information on ionic processes it is shown that a precise analysis of ionospheric behavior remains highly speculative.     

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.