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.     

A continental analysis of correlations between tree patterns in African savannas and human and environmental variables

This study analyses possible relationships between natural processes taking place in savannas and the tree patterns found in savannas. This can lead to new hypotheses about which processes are driving savanna physiognomy. To do so tree patterns were quantified for African savannas from historical aerial photographs applying frequently used landscape metrics. Also, additional data for these areas were collected to quantify the processes taking place at these locations. Correlations between tree pattern indices and explaining factors were analysed. We found a negative trend between tree cover and density of sheep and goats, but no relationship between tree cover and density of cattle, suggesting that small livestock have an effect on tree cover, but that larger livestock (or obligate grazers) do not. Also, a positive correlation between human population density and tree cover was found. Possible explanations for the found relations are discussed. Subsequent ways to analyse the latter correlation are discussed, and the potential of the presented historical database of aerial photographs is highlighted.     

Origin and Extent of Fresh Paleowaters on the Atlantic Continental Shelf, USA

While the existence of relatively fresh groundwater sequestered within permeable, porous sediments beneath the Atlantic continental shelf of North and South America has been known for some time, these waters have never been assessed as a potential resource. This fresh water was likely emplaced during Pleistocene sea-level low stands when the shelf was exposed to meteoric recharge and by elevated recharge in areas overrun by the Laurentide ice sheet at high latitudes. To test this hypothesis, we present results from a high-resolution paleohydrologic model of groundwater flow, heat and solute transport, ice sheet loading, and sea level fluctuations for the continental shelf from New Jersey to Maine over the last 2 million years. Our analysis suggests that the presence of fresh to brackish water within shallow Miocene sands more than 100 km offshore of New Jersey was facilitated by discharge of submarine springs along Baltimore and Hudson Canyons where these shallow aquifers crop out. Recharge rates four times modern levels were computed for portions of New England's continental shelf that were overrun by the Laurentide ice sheet during the last glacial maximum. We estimate the volume of emplaced Pleistocene continental shelf fresh water (less than 1 ppt) to be 1300 km³ in New England. We also present estimates of continental shelf fresh water resources for the U.S. Atlantic eastern seaboard (10⁴ km³) and passive margins globally (3 × 10⁵ km³). The simulation results support the hypothesis that offshore fresh water is a potentially valuable, albeit nonrenewable resource for coastal megacities faced with growing water shortages.     

Stormflow generation in a small rainforest catchment in the Luquillo Experimental Forest,Puerto Rico

Various complementary techniques were used to investigate the stormflow generating processes in a small headwater catchment in northeastern Puerto Rico. Over 100 samples were taken of soil matrix water, macropore flow, streamflow and precipitation, mainly during two storms of contrasting magnitude, for the analysis of calcium, magnesium, silicon, potassium, sodium and chloride. These were combined with hydrometric information on streamflow, return flow, precipitation, throughfall and soil moisture to distinguish water following different flow paths. Geo‐electric sounding was used to survey the subsurface structure of the catchment, revealing a weathering front that coincided with the elevation of the stream channel instead of running parallel to surface topography. The hydrometric data were used in combination with soil physical data, a one‐dimensional soil water model (VAMPS ) and a three‐component chemical mass‐balance mixing model to describe the stormflow response of the catchment. It is inferred that most stormflow travelled through macropores in the top 20 cm of the soil profile. During a large event, saturation overland flow also accounted for a considerable portion of the stormflow, although it was not possible to quantify the associated volume fully. Although the mass‐balance mixing model approach gave valuable information about the various flow paths within the catchment, it was not possible to distill the full picture from the model alone; additional hydrometric and soil physical evidence was needed to aid in the interpretation of the model results. Copyright © 2004 John Wiley & Sons, Ltd.     

Climatic niche breadth can explain variation in geographical range size of alpine and subalpine plants

Understanding the environmental factors determining the distribution of species with different range sizes can provide valuable insights for evolutionary ecology and conservation biology in the face of expected climate change. However, little is known about what determines the variation in geographical and elevational ranges of alpine and subalpine plant species. Here, we examined the relationship between geographical and elevational range sizes for 80 endemic rhododendron species in China using Spearman’s rank-order correlation. We ran the species distribution model – maximum entropy modelling (MaxEnt) – with 27 environmental variables. The importance of each variable to the model prediction was compared for species groups with different geographical and elevational range sizes. Our results showed that the correlation between geographical and elevational range sizes of rhododendron species was not significant. Climate-related variables were found to be the most important factors in shaping the distributional ranges of alpine and subalpine plant species across China. Species with geographically and elevationally narrow ranges had distinct niche requirements. For geographical ranges, the narrow-ranged species showed less tolerance to niche conditions than the wide-ranged species. For elevational ranges, compared with the wide-ranged species, the narrow-ranged species showed an equivalent niche breadth, but occurred at different niche position along the environmental gradient. Our findings suggest that over large spatial extents the elevational range size can be a complementary trait of alpine and subalpine plant species to geographical range size. Climatic niche breadth, especially the range of seasonal variability, can explain species’ geographical range sizes. Changes in climate may influence the distribution of rhododendrons, with the effects likely being felt most by species with either a narrow geographical or narrow elevational range.     

Adaptive port-starboard beamforming of triplet sonar arrays

For a low-frequency active sonar (LFAS) with a triplet receiver array, it is not clear in advance which signal processing techniques optimize its performance. Here, several advanced beamformers are analyzed theoretically, and the results are compared to experimental data obtained in sea trials. Triplet arrays are single line arrays with three hydrophones on a circular section of the array. The triplet structure provides the ability to solve the notorious port-starboard (PS) ambiguity problem of ordinary single-array receivers. More importantly, the PS rejection can be so strong that it allows to unmask targets in the presence of strong coastal reverberation or traffic noise. The theoretical and experimental performance of triplet array beamformers is determined in terms of two performance indicators: array gain and PS rejection. Results are obtained under several typical acoustic environments: sea noise, flow noise, coastal reverberation, and mixtures of these. A new algorithm for (beam space) adaptive triplet beamforming is implemented and tuned. Its results are compared to those of other triplet beamforming techniques (optimum and cardioid beamforming). These beamformers optimize for only one performance indicator, whereas in theory, the adaptive beamformer gives the best overall performance (in any given environment). The different beamformers are applied to data obtained with an LFAS at sea. Analysis shows that adaptive triplet beamforming outperforms conventional beamforming algorithms. Adaptive triplet beamforming provides strong PS rejection, allowing the unmasking of targets in the presence of strong directional reverberation (e.g., from a coast) and at the same time provides positive array gain in most environments.