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.     

Numerical simulation of coronal heating by resonant absorption of Alfvén waves

The heating of coronal loops by resonant absorption of Alfvén waves is studied in compressible, resistive magnetohydrodynamics. The loops are approximated by straight cylindrical, axisymmetric plasma columns and the incident waves which excite the coronal loops are modelled by a periodic external driver. The stationary state of this system is determined with a numerical code based on the finite element method. Since the power spectrum of the incident waves is not well known, the intrinsic dissipation is computed. The intrinsic dissipation spectrum is independent of the external driver and reflects the intrinsic ability of the coronal loops to extract energy from incident waves by the mechanism of resonant absorption.The numerical results show that resonant absorption is very efficient for typical parameter values occurring in the loops of the solar corona. A considerable part of the energy supplied by the external driver, is actually dissipated Ohmically and converted into heat. The heating of the plasma is localized in a narrow resonant layer with a width proportional to ^1/3. The energy dissipation rate is almost independent of the resistivity for the relevant values of this parameter. The efficiency of the heating mechanism and the localization of the heating strongly depend on the frequency of the external driver. Resonant absorption is extremely efficient when the plasma is excited with a frequency near the frequency of a so-called collective mode.     

Application of FLATModel, a 2D finite volume code, to debris flows in the northeastern part of the Iberian Peninsula

FLATModel is a two-dimensional shallow-water approximation code with corrections and modifications that create a simulation tool adapted to debris-flows behaviour. FLATModel uses the finite volume method with the numerical implementation of the Godunov scheme and includes correction terms regarding the effect of flow over high slopes and curvature. Additionally, the stop-and-go phenomenon, the basal entrainment and a correction regarding the front inclination of the final deposit are incorporated into FLATModel. In addition, different flow resistance laws were integrated in the numerical code including Bingham, Herschel–Bulkley and Voellmy fluid model. Firstly, our numerical model was validated using analytical solutions of a dam-break scenario and published data on a laboratory experiment. Secondly, three real events, which occurred in the northeastern part of the Iberian Peninsula, were back-calculated. Although field observations of the three events are not very detailed, the back-analyses revealed interesting patterns on the flow dynamics, and the numerical results generally showed good agreement with field data. Comparing the different flow resistance laws, the Voellmy fluid model presents the best behaviour regarding both the flow behaviour and the deposit characteristics. Preliminary simulation runs incorporating the effect of basal entrainment offered satisfactory results, although the final volume is rather sensitive on the selected friction angle of channel-bed material. The outcomes regarding the correction of the calculated front inclination of the final deposit showed that this implementation strongly improves the simulation results and better represents steep fronts of final deposits.