A multiphysics and multiscale software environment for modeling astrophysical systems

We present MUSE, a software framework for combining existing computational tools for different astrophysical domains into a single multiphysics, multiscale application. MUSE facilitates the coupling of existing codes written in different languages by providing inter-language tools and by specifying an interface between each module and the framework that represents a balance between generality and computational efficiency. This approach allows scientists to use combinations of codes to solve highly coupled problems without the need to write new codes for other domains or significantly alter their existing codes. MUSE currently incorporates the domains of stellar dynamics, stellar evolution and stellar hydrodynamics for studying generalized stellar systems. We have now reached a “Noah’s Ark” milestone, with (at least) two available numerical solvers for each domain. MUSE can treat multiscale and multiphysics systems in which the time- and size-scales are well separated, like simulating the evolution of planetary systems, small stellar associations, dense stellar clusters, galaxies and galactic nuclei. In this paper we describe three examples calculated using MUSE: the merger of two galaxies, the merger of two evolving stars, and a hybrid N-body simulation. In addition, we demonstrate an implementation of MUSE on a distributed computer which may also include special-purpose hardware, such as GRAPEs or GPUs, to accelerate computations. The current MUSE code base is publicly available as open source at http://muse.li.     

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.     

The High-Resolution Solar Reference Spectrum between 250 and 550 nm and its Application to Measurements with the Ozone Monitoring Instrument

We have constructed a new high resolution solar reference spectrum in the spectral range between 250 and 550 nm. The primary use of this spectrum is for the calibration of the Dutch – Finnish Ozone Monitoring Instrument (OMI), but other applications are mentioned. The incentive for deriving a new high resolution solar reference spectrum is that available spectra do not meet our requirements on radiometric accuracy or spectral resolution. In this paper we explain the steps involved in constructing the new spectrum, based on available low and high resolution spectra and discuss the main sources of uncertainty. We compare the result with solar measurements obtained with the OMI as well as with other UV-VIS space-borne spectrometers with a similar spectral resolution. We obtain excellent agreement with the OMI measurements, which indicates that both the newly derived solar reference spectrum and our characterization of the OMI instrument are well understood. We also find good agreement with previously published low resolution spectra. The absolute intensity scale, wavelength calibration and representation of the strength of the Fraunhofer lines have been investigated and optimized to obtain the resulting high resolution solar reference spectrum.     

The solar spectral irradiance and its action in the atmospheric photodissociation processes

A critical analysis has been made of solar irradiance in the spectral region covering wavelengths from 100 nm upwards; the absorption characteristics of molecules of oxygen and ozone have been taken into account with a view to the direct application of the results to atmospheric photochemistry. The absorption of radiation by these molecules results in the photodissociation of both of them in the homosphere, and it also makes possible the penetration of solar radiation from the thermosphere, through the mesosphere and the stratosphere, down to the troposphere.Special attention has been given to each of the following spectral regions: Lyman-alpha radiation at 121.6 nm, the O₂ Schumann-Runge continuum at wavelengths less than 175 nm, the O₂ Schumann-Runge band system from 200 to 175 nm, and the O₂ Herzberg continuum at 242.4 nm. For absorption by ozone, the solar spectrum has been analysed in the following regions: the Hartley band at wavelengths less than 310 nm, the Huggins bands at wavelengths above 310 nm and the visible Chappuis bands. Finally, for the photodissociation of O₃, particular attention has been given to the transition region (300–320 nm) in which there is a change-over from the production of the excited atom O(¹D) to that of the atom in its ground state O(³p).     

Earthquake awareness and perception of risk among the residents of Istanbul

Awareness and perception of risk are among the most crucial steps in the process of taking precautions at individual level for various hazards. In this study, we investigated the factors affecting better knowledge and greater risk perception about earthquakes among residents of Istanbul. A field survey was carried out, and a total of 1,123 people were interviewed in two districts of Istanbul with different seismic risk levels and from three (low, moderate and high) socio-economic levels (SEL). The findings showed that although the level of knowledge regarding earthquakes and preparedness for them was promising, it could be improved. The results indicated that future preparedness programmes should target people with lower educational and socio-economic levels. The media were the leading source of information among the respondents. Location of the home was a strong influence on individuals having above average earthquake knowledge and even more on high risk perception. Socio-economic parameters (educational level, economic status, SEL of the sub-district and tenure of the home), gender and attitude score were other factors influencing greater risk perception with regard to earthquakes.     

Tracing of variabilities within a geological barrier by molecular organic geochemistry.: Case of the Callovo-Oxfordian sedimentary series in the East of the Paris Basin (France)

The Callovo-Oxfordian claystones located at 500 m depth at Bure (Meuse, France) are currently being investigated by Andra (the French National Radioactive Waste Management Agency) for testing the feasibility of long-term and deep geological nuclear waste disposal. In order to evaluate its potential as a geological barrier, it is very important to study, assess and describe its physico-chemical variability. The molecular biomarker composition of 150 samples of these claystones and their surrounding limestones carry diverse information on the sources of the sedimentary organic matter, the chemistry of the depositional environment, the preservation and diagenesis conditions. It also allows assessing the degree of lateral and vertical variability of the organic matter within these sedimentary series. The abundance of unsaturated biomarkers, the distribution of steroids and hopanoids and CPI values >2 prove the thermal immaturity of the organic matter. The co-occurrence of plankton, bacteria and land plant biomarkers indicate that the organic matter is a mixture of marine and continental contributions. The data also reveal that the organic matter was deposited under oxic and open-sea conditions except for a brief event of photic zone anoxia at the beginning of the Middle Callovian. In the claystones, the geosynthesis of diasterenes is favored to the detriment of the formation of steranes, especially in smectite-rich levels, and the organic matter is rapidly isolated from oxidizing then reducing conditions after the deposition due to the protective effect of clays. On the scale investigated, the claystones are characterized by a unique molecular facies and are thus homogenous from their organic content point of view. Yet, detailed investigation of specific molecular families indicates changes related to major claystone–limestone transitions. The homogeneity of these claystones can be explained by the paleogeographic position of their depositional setting and the plane and sub-horizontal paleotopography on which they were deposited. This study demonstrates the efficiency of organic geochemistry in the assessment of the variability of geological barriers.     

Buchbesprechungen

Ohne Zusammenfassung     

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.     

An ecosystem-based composite spatial model for floodplain vulnerability assessment: a case study of Artigas,Uruguay

Floods are natural processes that constitute a hazard to society when associated to improper land use. Anthropic activities in floodplains are a factor of vulnerability that converts a natural hazard into a threat factor, eventually leading to disaster. Nowadays, natural and social complex processes demand integrated assessments in order to improve their understanding, helping decision making over sustainable use of territory, as well as integrating society’s activity in ecosystems and potentials, restrictions and benefits that society obtain from them. In this context, the objective of this work was to build a composite vulnerability model for a floodplain under urban influence, using an integrated assessment approach. This model was based on three dimensions; threat, fragility and an ecosystem services provision. These dimensions were calculated using both primary and secondary information, and weights by specialists. Main results show that the area presents high vulnerability with an increasing gradient towards high and urbanized areas, associated with an important number and relevant ecosystem services. Also, a spatial heterogeneity of the three dimensions emerged, making evident this area’s complexity and the need of integrated assessments to approach it. The composite vulnerability model proposed presents an elevated potential for natural and social processes analysis in floodplains, which is crucial for these territory management. Moreover, these integrated dimensions could contribute to decision making in different levels, as well as generating important supplies for environmental management and land planning.