Resolving cosmic structure formation with the Millennium-II Simulation

We present the Millennium-II Simulation (MS-II), a very large N -body simulation of dark matter evolution in the concordance Λ cold dark matter (ΛCDM) cosmology. The MS-II assumes the same cosmological parameters and uses the same particle number and output data structure as the original Millennium Simulation (MS), but was carried out in a periodic cube one-fifth the size  (100  h ⁻¹ Mpc)  with five times better spatial resolution (a Plummer equivalent softening of  1.0  h ⁻¹ kpc  ) and with 125 times better mass resolution (a particle mass of  6.9 × 10⁶  h ⁻¹ M_⊙  ). By comparing results at MS and MS-II resolution, we demonstrate excellent convergence in dark matter statistics such as the halo mass function, the subhalo abundance distribution, the mass dependence of halo formation times, the linear and non-linear autocorrelations and power spectra, and halo assembly bias. Together, the two simulations provide precise results for such statistics over an unprecedented range of scales, from haloes similar to those hosting Local Group dwarf spheroidal galaxies to haloes corresponding to the richest galaxy clusters. The 'Milky Way' haloes of the Aquarius Project were selected from a lower resolution version of the MS-II and were then resimulated at much higher resolution. As a result, they are present in the MS-II along with thousands of other similar mass haloes. A comparison of their assembly histories in the MS-II and in resimulations of 1000 times better resolution shows detailed agreement over a factor of 100 in mass growth. We publicly release halo catalogues and assembly trees for the MS-II in the same format within the same archive as those already released for the MS.     

Effects of seed particle size and composition on heterogeneous nucleation of n-nonane

Heterogeneous nucleation of supersaturated n-nonane vapour on seed particles of different size and composition has been investigated using a fast expansion chamber. Monodisperse seed particle sizes were ranging from about 4 nm up to about 24 nm in diameter. By using different types of particle generators WO_x, Ag and (NH₄)₂SO₄ particles were generated. For direct comparison between different particle compositions overlapping sizes have been generated for WO_x and Ag at about 7 nm particle diameter as well as for Ag and (NH₄)₂SO₄ at about 15 nm. Nucleation temperature was kept constant at about 278 K. Experimental data were compared to Kelvin equation and Fletcher theory including the effect of line tension. It was found that heterogeneous nucleation of n-nonane seems to be independent of seed particle composition and starts well below the Kelvin curve. Good agreement was achieved with Fletcher theory including the effect of line tension.     

The importance of three centuries of land-use change for the global and regional terrestrial carbon cycle

Large amounts of carbon (C) have been released into the atmosphere over the past centuries. Less than half of this C stays in the atmosphere. The remainder is taken up by the oceans and terrestrial ecosystems. Where does the C come from and where and when does this uptake occur? We address these questions by providing new estimates of regional land-use emissions and natural carbon fluxes for the 1700–2000 period, simultaneously considering multiple anthropogenic (e.g. land and energy demand) and biochemical factors in a geographically explicit manner. The observed historical atmospheric CO₂ concentration profile for the 1700 to 2000 period has been reproduced well. The terrestrial natural biosphere has been a major carbon sink, due to changes in climate, atmospheric CO₂, nitrogen and management. Due to land-use change large amounts of carbon have been emitted into the atmosphere. The net effect was an emission of 35 Pg C into the atmosphere for the 1700 to 2000 period. If land use had remained constant at its distribution in 1700, then the terrestrial C uptake would have increased by 142 Pg C. This overall difference of including or excluding land-use changes (i.e. 177 Pg C) comes to more than half of the historical fossil-fuel related emissions of 308 Pg C. Historically, global land-use emissions were predominantly caused by the expansion of cropland and pasture, while wood harvesting (for timber and fuel wood) only played a minor role. These findings are robust even when changing some of the important drivers like the extent of historical land-use changes. Under varying assumptions, land-use emissions over the past three centuries could have increased up to 20%, but remained significantly lower than from other sources. Combining the regional land-use and natural C fluxes, North America and Europe were net C sources before 1900, but turned into sinks during the twentieth century. Nowadays, these fluxes are a magnitude smaller than energy- and industry-related emissions. Tropical regions were C neutral prior to 1950, but then accelerated deforestation turned these regions into major C sources. The energy- and industry-related emissions are currently increasing in many tropical regions, but are still less than the land-use emissions. Based on the presented relevance of the land-use and natural fluxes for the historical C cycle and the significance of fossil-fuel emissions nowadays, there is a need for an integrated approach for energy, nature and land use in evaluating possible climate change mitigation policies.     

Changing hillslope and fluvial Holocene sediment dynamics in a Belgian loess catchment

Floodplain deposition is an essential part of the Holocene sediment dynamics of many catchments and a thorough dating control of these floodplain deposits is therefore essential to understand the driving forces of these sediment dynamics. In this paper we date floodplain and colluvial deposition in the Belgian Dijle catchment using accelerator mass spectrometric radiocarbon and optical stimulated luminescence dating. Relative mass accumulation curves for the Holocene were constructed for three colluvial sites and 12 alluvial sites. A database was constructed of all available radiocarbon ages of the catchment and this database was analysed using relative sediment mass accumulation rates and cumulative probability functions of ages and site‐specific sedimentation curves. Cumulative probability functions of ages were split into different depositional environments representing stable phases and phases of accelerated clastic deposition. The results indicate that there is an important variation between the different dated sites. After an initial stable early and middle Holocene phase with mainly peat growth in the floodplains, clastic sedimentation rates increased from 4000 BC on. This first phase was more pronounced and started somewhat earlier for colluvial deposits then for alluvial deposits. The main part of the Holocene deposits, both in colluvial and alluvial valleys, was deposited during the last 1 ka. The sedimentation pattern of the individual dated sites and the catchment‐wide pattern indicate that land use changes are responsible for the main variations in the Holocene sediment dynamics of this catchment, while the field data do not provide indications for a climatological influence on the sediment dynamics. Copyright © 2011 John Wiley & Sons, Ltd.     

Modelling mean annual sediment yield using a distributed approach

In this paper a spatially distributed model for the calculation of sediment delivery to river channels is presented (SEDEM: SEdiment DElivery Model). The model consists of two components: (1) the calculation of a spatial pattern of mean annual soil erosion rates in the catchment using a RUSLE (Revised Soil Erosion Equation) approach; and (2) the routing of the eroded sediment to the river channel network taking into account the transport capacity of each spatial unit. If the amount of routed sediment exceeds the local transport capacity, sediment deposition occurs. An existing dataset on sediment yield for 24 catchments in central Belgium was used to calibrate the transport capacity parameters of the model. A validation of the model results shows that the sediment yield for small and medium sized catchments (10–5000 ha) can be predicted with an average accuracy of 41 per cent. The predicted sediment yield values with SEDEM are significantly more accurate than the predictions using a lumped regression model. Moreover a spatially distributed approach allows simulation of the effect of different land use scenarios and soil conservation techniques. Copyright © 2001 John Wiley & Sons, Ltd.     

The formulation of an effective demulsifier for oil spill emulsions

The formation of stable sea water-in-crude oil emulsions after an oil spill was first reported after the Torrey Canyon wreck in 1967. The problem of handling this stable water (65%)-in-oil (35%) emulsion rather than oil is evidenced in the additional expensive on-site storage requirement as well as the handling difficulty with the viscous semi-solid emulsion. This paper details the results of an R&D program that has established the mechanism for the formation of this emulsion. This insight enabled a specific effective chemical demulsifier to be formulated. A treatment technique was devised that entails injecting the demulsifier into the oil spill emulsion as it is skimmed from the water surface. The emulsion is then rapidly separated into water-free oil and a clean-water phase that can be decanted back into the sea. Full-scale tests were successfully conducted by JBF Scientific Corp. using two types of commercial JBF DIP skimmers in a small concrete lined pond.     

A Note on Second Species Solutions Generated from p–q Resonant Orbits

The first image under the flow of the restricted three-body problem of the p—q resonant strips — that appear in the study of the p—q resonant orbits — do not have, in general, intersection with the strip. In this paper we show some particular situations in which the above intersections exist for some very simple p—q resonant orbits which, at the same time, are periodic second species solutions.This revised version was published online in October 2005 with corrections to the Cover Date.     

Three-Dimensional p–q Resonant Orbits Close to Second Species Solutions

The purpose of this paper is to study, for small values of , the three-dimensional p–q resonant orbits that are close to periodic second species solutions (SSS) of the restricted three-body problem. The work is based on an analytic study of the in- and out-maps. These maps are associated to follow, under the flow of the problem, initial conditions on a sphere of radius around the small primary, and consider the images of those initial points on the same sphere. The out-map is associated to follow the flow forward in time and the in-map backwards. For both mappings we give analytical expressions in powers of the mass parameter. Once these expressions are obtained, we proceed to the study of the matching equations between both, obtaining initial conditions of orbits that will be 'periodic' with an error of the order ^1–, for some (1/3,1/2). Since, as 0, the inner solution and the outer solution will collide with the small primary, these orbits will be close to SSS.     

3D wavepath migration

3D Kirchhoff migration (KM) smears a trace's time sample along a quasi-ellipsoid in the model space. This is a costly and sometimes noisy process as reflection energy is smeared far away from the actual reflector position, introducing far-field migration artefacts. As a reduced form of 3D KM, 3D wavepath migration (WM) smears a picked reflection arrival to a small Fresnel zone portion centred about the specular reflection point, leading to fewer migration artefacts and reduced computation time. Both the traveltime and the angle of incidence are required by WM for locating the specular reflection point. Our results with 3D prestack synthetic data show that WM generates fewer migration artefacts and can sometimes define complex structure better than KM. Our results with 3D prestack field data show that WM can mostly suppress migration artefacts and can sometimes resolve reflection interfaces better than KM. The CPU comparison shows that, for both the synthetic and field data examples, WM can be more than an order of magnitude faster than KM. The limitation with 3D WM is that the angle of incidence calculation is sensitive to the recording geometry and the signal-to-noise (S/N) ratio, which can lead to blurred images.