Swinging and tumbling of elastic capsules in shear flow

The deformation of an elastic micro-capsule in an infinite shear flow is studied numerically using a spectral method. The shape of the capsule and the hydrodynamic flow field are expanded into smooth basis functions. Analytic expressions for the derivative of the basis functions permit the evaluation of elastic and hydrodynamic stresses and bending forces at specified grid points in the membrane. Compared to methods employing a triangulation scheme, this method has the advantage that the resulting capsule shapes are automatically smooth, and few modes are needed to describe the deformation accurately. Computations are performed for capsules both with spherical and ellipsoidal unstressed reference shape. Results for small deformations of initially spherical capsules coincide with analytic predictions. For initially ellipsoidal capsules, recent approximative theories predict stable oscillations of the tank-treading inclination angle, and a transition to tumbling at low shear rate. Both phenomena have also been observed experimentally. Using our numerical approach we could reproduce both the oscillations and the transition to tumbling. The full phase diagram for varying shear rate and viscosity ratio is explored. While the numerically obtained phase diagram qualitatively agrees with the theory, intermittent behaviour could not be observed within our simulation time. Our results suggest that initial tumbling motion is only transient in this region of the phase diagram.Comment: 20 pages, 7 figure

A two-species continuum model for aeolian sand transport

Starting from the physics on the grain scale, we develop a simple continuum description of aeolian sand transport. Beyond popular mean-field models, but without sacrificing their computational efficiency, it accounts for both dominant grain populations, hopping (or "saltating") and creeping (or "reptating") grains. The predicted stationary sand transport rate is in excellent agreement with wind tunnel experiments simulating wind conditions ranging from the onset of saltation to storms. Our closed set of equations thus provides an analytically tractable, numerically precise, and computationally efficient starting point for applications addressing a wealth of phenomena from dune formation to dust emission.Comment: 23 pages, 9 figure

Diatom metabarcoding and microscopic analyses from sediment samples at Lake Nam Co, Tibet: The effect of sample-size and bioinformatics on the identified communities

Diatoms (Bacillariophyceae) are characterized by silicified cell walls that favor their long-term preservation in sediments, therefore widely used as bioindicators of present and past water conditions. Alongside with traditional morphological analyses, metabarcoding has become a valuable tool to study the community structures of various organisms, including diatoms. Here, we test whether the quantity of sediment sample used for DNA extraction affects the results obtained from high-throughput sequencing (metabarcoding) of the diatom rbcL region by isolating DNA from 10 g and 0.5 g (wet weight) of lake surface sediment samples. Because bioinformatics processing of metabarcoding data may affect the outcome, we also tested the consistency of the results from three different pipelines: 1) ESVs (exact sequence variants) pipeline; 2) clustering sequences at 95% sequence identity to form OTUs (operational taxonomic units; 95% OTUs); and 3) 97% OTUs pipeline. Additionally, the agreement between metabarcoding data and morphological inventories of corresponding samples were compared. Our results demonstrate highly uniform patterns between the diatom rbcL amplicons from 10 g and 0.5 g of sedimentary DNA (sedDNA) extracts (HTS 10 and HTS 0.5, respectively). Furthermore, after the careful curation of the sequencing data, metabarcoding results were highly consistent among the data sets produced by different bioinformatics pipelines. Comparing results from metabarcoding and microscopy, we identified some taxonomic mismatches: morphological analyses identified 59 diatom genera, whereas metabarcoding 49 to 54 genera. These mismatches are related to incompleteness of the sequence databases, but also to inconsistencies in diatom taxonomy in general and potential dissolution effects of diatom valves caused by high alkalinity of the investigated lake waters. Nevertheless, multivariate community analysis revealed consistent results between data sets identified by microscopy and metabarcoding – water depth and conductivity as the most significant variables in driving diatom communities in Lake Nam Co – further confirming that metabarcoding is a viable method for identifying diatom-environment relationships

Spin-resolved electron-impact ionization of lithium

Electron-impact ionization of lithium is studied using the convergent close-coupling (CCC) method at 25.4 and 54.4 eV. Particular attention is paid to the spin-dependence of the ionization cross sections. Convergence is found to be more rapid for the spin asymmetries, which are in good agreement with experiment, than for the underlying cross sections. Comparison with the recent measured and DS3C-calculated data of Streun et al (1999) is most intriguing. Excellent agreement is found with the measured and calculated spin asymmetries, yet the discrepancy between the CCC and DS3C cross sections is very large

Corridors of barchan dunes: stability and size selection

Barchans are crescentic dunes propagating on a solid ground. They form dune fields in the shape of elongated corridors in which the size and spacing between dunes are rather well selected. We show that even very realistic models for solitary dunes do not reproduce these corridors. Instead, two instabilities take place. First, barchans receive a sand flux at their back proportional to their width while the sand escapes only from their horns. Large dunes proportionally capture more than they loose sand, while the situation is reversed for small ones: therefore, solitary dunes cannot remain in a steady state. Second, the propagation speed of dunes decreases with the size of the dune: this leads -- through the collision process -- to a coarsening of barchan fields. We show that these phenomena are not specific to the model, but result from general and robust mechanisms. The length scales needed for these instabilities to develop are derived and discussed. They turn out to be much smaller than the dune field length. As a conclusion, there should exist further - yet unknown - mechanisms regulating and selecting the size of dunes.Comment: 13 pages, 13 figures. New version resubmitted to Phys. Rev. E. Pictures of better quality available on reques

A Continuum Saltation Model for Sand Dunes

We derive a phenomenological continuum saltation model for aeolian sand transport that can serve as an efficient tool for geomorphological applications. The coupled differential equations for the average density and velocity of sand in the saltation layer reproduce both known equilibrium relations for the sand flux and the time evolution of the sand flux as predicted by microscopic saltation models. The three phenomenological parameters of the model are a reference height for the grain-air interaction, an effective restitution coefficient for the grain-bed interaction, and a multiplication factor characterizing the chain reaction caused by the impacts leading to a typical time or length scale of the saturation transients. We determine the values of these parameters by comparing our model with wind tunnel measurements. Our main interest are out of equilibrium situations where saturation transients are important, for instance at phase boundaries (ground/sand) or under unsteady wind conditions. We point out that saturation transients are indispensable for a proper description of sand flux over structured terrain, by applying the model to the windward side of an isolated dune, thereby resolving recently reported discrepancies between field measurements and theoretical predictions.Comment: 11 pages, 7 figure

Implementation of U.K. Earth system models for CMIP6

We describe the scientific and technical implementation of two models for a core set of experiments contributing to the sixth phase of the Coupled Model Intercomparison Project (CMIP6). The models used are the physical atmosphere-land-ocean-sea ice model HadGEM3-GC3.1 and the Earth system model UKESM1 which adds a carbon-nitrogen cycle and atmospheric chemistry to HadGEM3-GC3.1. The model results are constrained by the external boundary conditions (forcing data) and initial conditions.We outline the scientific rationale and assumptions made in specifying these. Notable details of the implementation include an ozone redistribution scheme for prescribed ozone simulations (HadGEM3-GC3.1) to avoid inconsistencies with the model's thermal tropopause, and land use change in dynamic vegetation simulations (UKESM1) whose influence will be subject to potential biases in the simulation of background natural vegetation.We discuss the implications of these decisions for interpretation of the simulation results. These simulations are expensive in terms of human and CPU resources and will underpin many further experiments; we describe some of the technical steps taken to ensure their scientific robustness and reproducibility

Comparative study of argon 3p electron-impact ionization at low energies

(c) 2001 The American Physical SocietyAn experimental and theoretical study of electron-impact ionization of the 3p orbital in argon is presented. The (e,2e) technique was used to measure the relative triple-differential cross section for this process in the coplanar asymmetric geometry. The experimental results were obtained at an incident electron energy of 113.5 eV, a scattering angle of 15°, and ejected electron energies of 10, 7.5, 5, and 2 eV. The experimental data are compared with a distorted-wave Born approximation (DWBA) calculation, and also with previous results for argon 3s ionization obtained under identical kinematic conditions. Discrepancies between the experimental and theoretical data are attributed to the effects of charge-cloud polarization and higher-order scattering processes, which are not incorporated in the DWBA calculation.Haynes, Matthew A. and Lohmann, Birgi

Evidence for early warming and cooling in North Atlantic surface waters during the last interglacial

In-depth analysis of planktic foraminiferal census data paired with δ18O records of specific indicator species provides new insight into the surface ocean evolution of the northeast Atlantic during the previous interglacial warm period (oxygen isotope stage (OIS) 5e). Full interglacial conditions existed at the study site for a maximum of only 8 kyr, between 125 and 117 ka. Highest sea surface temperatures (SSTs) occurred during early OIS 5e concomitant with high summer insolation but after the main phase of ice sheet melting of the preceding glaciation (Saalian). This early peak SST interval is marked by the appearance of tropical-subtropical species and lasted for 4 kyr until 121 ka, as corroborated by a major change in planktic δ18O. Relative stability in global ice volume continued for another 3–4 kyr before SSTs dropped further toward the next stadial. During early OIS 5e the situation of the surface water vertical structure appears to have been different from the early Holocene. For OIS 5e it is therefore suggested that the particular melting history of late Saalian ice had a long-lasting and profound effect on both postdeglacial surface water mass configuration in the North Atlantic and heat-moisture transfer into Europe