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.     

Structured ionospheric outflow during the Cassini T55-T59 Titan flybys

During the final three of the five consecutive and similar Cassini Titan flybys T55-T59 we observe a region characterized by high plasma densities (electron densities of 1-8 cm⁻³) in the tail/nightside of Titan. This region is observed progressively farther downtail from pass to pass and is interpreted as a plume of ionospheric plasma escaping Titan, which appears steady in both location and time. The ions in this plasma plume are moving in the direction away from Titan and are a mixture of both light and heavy ions with composition revealing that their origin are in Titan's ionosphere, while the electrons are more isotropically distributed. Magnetic field measurements indicate the presence of a current sheet at the inner edge of this region. We discuss the mechanisms behind this outflow, and suggest that it could be caused by ambipolar diffusion, magnetic moment pumping or dispersive Alfvén waves.     

Characteristics and regional classification of the copepod community in Ariake Bay with note on comparison with 3 decades ago

Ariake Bay, Kyushu Island, has recently exhibited environmental degradation in the form of red tides and anoxic bottom water. To determine the characteristics and regional classification of the copepod community, zooplankton surveys were made throughout the bay in three cruises each in October 2004, January 2005 and March 2009 by vertical tows of a 0.1-mm-mesh plankton net. Oithona davisae was the most abundant in January and March, and Microsetella norvegica in October. Cluster analysis revealed that the copepod community from each cruise was generally separated into the inner to eastern-middle region and the western-middle to mouth region. A SIMPER analysis revealed that dissimilarity between the groups and similarity within each group were mainly due to the dominant species, but the similarities within the outer-region group in January and March were mostly contributed by Paracalanus parvus s.l. A non-metric multidimensional scaling with bubble plots of environmental variables and a BIOENV analysis showed that transparency was more different between the two groups than temperature and salinity. As compared with two similarly sized, eutrophic bays (Tokyo and Osaka Bays), Ariake Bay differs from Tokyo Bay in the dominance of M. norvegica and from Osaka Bay in the high abundance of O. davisae even in the colder season. The dominant species and the regional classification of the copepod community in the bay were similar to those in the studies in the 1970s, suggesting that the copepod community has not notably changed during the last 3 decades.     

Effect of temperature on the methyl chloride production rate in a marine phytoplankton, <Emphasis Type="Italic">Phaeodactylum tricornutum</Emphasis>

Methyl halides such as methyl chloride (CH₃Cl) are known to be important carriers of halogen from the ocean to the atmosphere, and the halogens they release into the stratosphere by photolysis catalyze ozone depletion. Marine phytoplankton have been reported as a source of CH₃Cl, but the effects of environmental temperature on the CH₃Cl production by phytoplankton have not been investigated. In this study, we investigated the effects of temperature on the production of CH₃Cl in the culture of a marine diatom, Phaeodactylum tricornutum CCMP 630, incubated at 10, 15, 20, 25, and 30 °C. CH₃Cl concentrations in cultured samples were determined using purge and trap gas chromatograph–mass spectrometry. Phytoplankton growth was monitored by measuring the chlorophyll a concentrations. CH₃Cl production was observed for several weeks at four different temperatures ranging from 10 to 25 °C. The CH₃Cl production from P. tricornutum was increased with increasing temperature from 10 to 25 °C, and the maximum production rate for CH₃Cl was 0.21~0.26 μmol (g chlorophyll a)^?1 d^?1 at 25 °C, which was several times higher than that at 10 °C (~0.03 μmol (g chlorophyll a)^?1 d^?1). The Arrhenius equation was successfully used to characterize the effects of temperature on the production rates of CH₃Cl in the culture of P. tricornutum. Our results suggest that water temperature directly affects CH₃Cl production derived from P. tricornutum and that water temperature would be a significant factor for estimating the emissions of CH₃Cl from marine environments.     

Magmatic fluids of Tatun volcanic group,Taiwan

Two distinctive magmatic fluids were recognized in the Tatun volcanic group (TVG), Taiwan. One is a relatively reduced fluid represented by the fumarolic gases at Hsiao-you-ken (HYK) geothermal field. Another is an oxidized fluid containing high concentrations of HCl represented by the fumarolic gases at Da-you-ken (DYK). An intermediate gas was recognized at Gung-tze-ping (GTP) and She-hung-ping (SHP). The fumarolic gases at HYK and GTP possess the features of so-called primary steam generated on mixing of magmatic gas and meteoric groundwater. The fumarolic gases at DYK are a simple mixture between magmatic gas and water vapor of meteoric origin. The CO₂/H₂O molar ratio of the magmatic component in the fumarolic gases at DYK was estimated to be 0.018, meanwhile it was estimated to be 0.027 for the fumarolic gases at HYK and GTP, suggesting the magma beneath DYK is depleted in volatiles relative to the magma beneath HYK and GTP. The estimated CO₂/H₂O ratio for the magmatic component is comparable to that of some active volcanoes in Japan, suggesting the enrichment of volatiles in the magmas beneath TVG.     

The electromagnetic pickup of submicron-sized dust above Enceladus’s northern hemisphere

As the saturnian magnetoplasma sweeps past Enceladus, it experiences both a decrease in electron content and sharp slowdown in the northern hemisphere region within ～5 Enceladus Radii (R_e). This slowdown is observed by Cassini in regions not obviously associated with the southern directed plume-originating ions. We suggest herein that the decrease in northern hemisphere electron content and plasma slowdown could both be related to the presence of fine dust grains that are being accelerated by the Lorentz force created within the saturnian magnetic field system.     

Estimation of vertical profiles of leaf drying times after daytime rainfall within a Bornean tropical rainforest

South‐east Asian Bornean tropical rainforests have large and complex canopy structures. To clarify how forest structure affects wet‐canopy evaporation, it is necessary to consider wet‐canopy evaporation processes within the forest canopy, such as vertical profiles of canopy drying time. In a previous study a method was proposed that utilizes sap flow measurements to estimate canopy drying time after rainfall (CDTobs) during daytime. In this present study, the method was applied to estimations of the vertical variations in CDTobs in 11 individual trees of various heights, ranging from 2·7 m to 53 m. The study derived vertical profiles and showed that the lengths of CDTobs in lower canopy trees were 2–4 h longer than those in the upper canopy trees. The new method for CDTobs profiles presented in this study, which is available for validation of multi‐layer biosphere‐atmospheric models, is a useful illustration for clarifying wet‐canopy evaporation processes in tropical rainforests. Copyright © 2008 John Wiley & Sons, Ltd.     

The calibration of the Cassini-Huygens CAPS Electron Spectrometer

We present the two-stage method used to calibrate the electron spectrometer (ELS), part of the plasma spectrometer (CAPS) on board the Cassini spacecraft currently in orbit around Saturn. The CAPS-ELS is a top-hat electrostatic analyser designed to measure electron fluxes between 0.5 eV and 26 keV. The on-ground calibration method described here includes the production of photoelectrons, which are energised and passed into the CAPS-ELS in a purpose designed calibration facility. Knowledge of the intensity of these incident electrons and the subsequent instrument output provides an on-ground calibrated geometric factor. Comparative studies of physical quantities such as plasma density and electron differential flux calculated using on-ground calibration factor with the quantities deduced from the wave experiment and high energy electron detector provide in-flight calibration. The results of this are presented together with a comparison of the experimentally calibrated values with simulated calibration values.     

Late cretaceous to early paleogene paleomagnetic results from Sikhote Alin, far eastern Russia: implications for deformation of East Asia

Welded tuffs in the Bogopol and Sijanov groups were sampled at 27 sites from 12 caldera formations in the Sikhote Alin mountain range around Kavalerovo (44.3°N, 135.0°E) for chronological and paleomagnetic studies. KAr age dates show that the welded tuffs erupted between 66 Ma and 46 Ma. All sites yield reliable paleomagnetic directions, with unblocking temperatures higher than 560°C. The high-temperature component at 12 sites and the medium-temperature component at 3 sites in the Bogopol Group show reversed polarity (D = 193.7°, I = −57.6°,₉₅ = 8.1°). The high-temperature component at 11 sites in the Sijanov Group showed both reversed and normal polarities and its mean direction reveals no detectable deflection from north (D = −2.9°, I = 59.6°,₉₅ = 11.2°). The combined paleomagnetic direction of the two groups yields a paleomagnetic pole of 250.5°E, 84.1°N (A₉₅ = 8.8°), which falls near Cretaceous paleomagnetic poles from Outer Mongolia, Inner Mongolia, the North China Block and the South China Block. The Sikhote Alin area appears not to have been subjected to detectable motion with respect to East Asia since about 50 Ma. This implies that the Sikhote Alin area behaved as an integral part of East Asia during the opening of the Japan Sea at about 15 Ma. However, significant separation between the paleomagnetic poles of East Asia and Europe during the Jurassic-Paleogene implies a major relative movement between these two blocks since the Paleogene.