Parametric uncertainty in complex environmental models: a cheap emulation approach for models with high-dimensional output

In order to understand underlying processes governing environmental and physical processes, and predict future outcomes, a complex computer model is frequently required to simulate these dynamics. However there is inevitably uncertainty related to the exact parametric form or the values of such parameters to be used when developing these simulators, with \emph{ranges} of plausible values prevalent in the literature. Systematic errors introduced by failing to account for these uncertainties have the potential to have a large effect on resulting estimates in unknown quantities of interest. Due to the complexity of these types of models, it is often unfeasible to run large numbers of training runs that are usually required for full statistical emulators of the environmental processes. We therefore present a method for accounting for uncertainties in complex environmental simulators without the need for very large numbers of training runs and illustrate the method through an application to the Met Office's atmospheric transport model NAME. We conclude that there are two principle parameters that are linked with variability in NAME outputs, namely the free tropospheric turbulence parameter and particle release height. Our results suggest the former should be significantly larger than is currently implemented as a default in NAME, whilst changes in the latter most likely stem from inconsistencies between the model specified ground height at the observation locations and the true height at this location. Estimated discrepancies from independent data are consistent with the discrepancy between modelled and true ground height

The position of graptolites within Lower Palaeozoic planktic ecosystems.

An integrated approach has been used to assess the palaeoecology of graptolites both as a discrete group and also as a part of the biota present within Ordovician and Silurian planktic realms. Study of the functional morphology of graptolites and comparisons with recent ecological analogues demonstrates that graptolites most probably filled a variety of niches as primary consumers, with modes of life related to the colony morphotype. Graptolite coloniality was extremely ordered, lacking any close morphological analogues in Recent faunas. To obtain maximum functional efficiency, graptolites would have needed varying degrees of coordinated automobility. A change in lifestyle related to ontogenetic changes was prevalent within many graptolite groups. Differing lifestyle was reflected by differing reproductive strategies, with synrhabdosomes most likely being a method for rapid asexual reproduction. Direct evidence in the form of graptolithophage 'coprolitic' bodies, as well as indirect evidence in the form of probable defensive adaptations, indicate that graptolites comprised a food item for a variety of predators. Graptolites were also hosts to a variety of parasitic organisms and provided an important nutrient source for scavenging organisms

The dynamic centres of infrared-dark clouds and the formation of cores

High-mass stars have an enormous influence on the evolution of the interstellar medium in galaxies, so it is important that we understand how they form. We examine the central clumps within a sample of seven infrared-dark clouds (IRDCs) with a range of masses and morphologies. We use 1-pc-scale observations from the Northern Extended Millimeter Array (NOEMA) and the IRAM 30m telescope to trace dense cores with 2.8-mm continuum, and gas kinematics in C18O, HCO+, HNC, and N2H+ (J = 1–0). We supplement our continuum sample with six IRDCs observed at 2.9 mm with the Atacama Large Millimeter/submillimeter Array (ALMA), and examine the relationships between core- and clump-scale properties. We have developed a fully automated multiple-velocity component hyperfine line-fitting code called MWYDYN which we employ to trace the dense gas kinematics in N2H+ (1–0), revealing highly complex and dynamic clump interiors. We find that parsec-scale clump mass is the most important factor driving the evolution; more massive clumps are able to concentrate more mass into their most massive cores – with a log-normally distributed efficiency of around 9 per cent – in addition to containing the most dynamic gas. Distributions of linewidths within the most massive cores are similar to the ambient gas, suggesting that they are not dynamically decoupled, but are similarly chaotic. A number of studies have previously suggested that clumps are globally collapsing; in such a scenario, the observed kinematics of clump centres would be the direct result of gravity-driven mass inflows that become ever more complex as the clumps evolve, which in turn leads to the chaotic mass growth of their core populations

CHIMPS: the 13CO/C18O (J = 3 to 2) Heterodyne Inner Milky Way Plane Survey

We present the 13CO/C18O (J = 3 → 2) Heterodyne Inner Milky Way Plane Survey (CHIMPS) which has been carried out using the Heterodyne Array Receiver Program on the 15 m James Clerk Maxwell Telescope (JCMT) in Hawaii. The high-resolution spectral survey currently covers |b| ≤ 0.5° and 28° ≲ l ≲ 46°, with an angular resolution of 15 arcsec in 0.5 km s-1 velocity channels. The spectra have a median rms of ˜0.6 K at this resolution, and for optically thin gas at an excitation temperature of 10 K, this sensitivity corresponds to column densities of NH2 ˜ 3 × 1020 cm-2 and NH2 ˜ 4 × 1021 cm-2 for 13CO and C18O, respectively. The molecular gas that CHIMPS traces is at higher column densities and is also more optically thin than in other publicly available CO surveys due to its rarer isotopologues, and thus more representative of the three-dimensional structure of the clouds. The critical density of the J = 3 → 2 transition of CO is ≳104 cm-3 at temperatures of ≤20 K, and so the higher density gas associated with star formation is well traced. These data complement other existing Galactic plane surveys, especially the JCMT Galactic Plane Survey which has similar spatial resolution and column density sensitivity, and the Herschel infrared Galactic Plane Survey. In this paper, we discuss the observations, data reduction and characteristics of the survey, presenting integrated-emission maps for the region covered. Position-velocity diagrams allow comparison with Galactic structure models of the Milky Way, and while we find good agreement with a particular four-arm model, there are some significant deviations.Peer reviewedFinal Accepted Versio

Micronutrient and amino acid losses during renal replacement therapy for acute kidney injury

© 2019 International Society of Nephrology Introduction: Malnutrition is common in patients with acute kidney injury (AKI), particularly in those requiring renal replacement therapy (RRT). Use of RRT removes metabolic waste products and toxins, but it will inevitably also remove useful molecules such as micronutrients, which might aggravate malnutrition. The RRT modalities vary in mechanism of solute removal; for example, intermittent hemodialysis (IHD) uses diffusion, continuous veno-venous hemofiltration (CVVH) uses convection, and sustained low-efficiency diafiltration (SLEDf) uses a combination of these. Methods: We assessed micronutrient and amino acid losses in 3 different RRT modalities in patients with AKI (IHD, n = 27; SLEDf, n = 12; CVVH, n = 21) after correction for dialysis dose and plasma concentrations. Results: Total losses were affected by modality; generally CVVH >> SLEDf > IHD (e.g., amino acid loss was 18.69 ± 3.04, 8.21 ± 4.07, and 5.13 ± 3.1 g, respectively; P < 0.001). Loss of specific trace elements (e.g., copper and zinc) during RRT was marked, with considerable heterogeneity between RRT types (e.g., +849 and +2325 μg/l lost during SLEDf vs. IHD, respectively), whereas effluent losses of copper and zinc decreased during CVVH (effect size relative to IHD, −3167 and −1442 μg/l, respectively). B vitamins were undetectable in effluent, but experimental modeling estimated 40% to 60% loss within the first 15 minutes of RRT. Conclusion: Micronutrient and amino acid losses are marked during RRT in patients with AKI, with variation between RRT modalities and micronutrients

SEDIGISM-ATLASGAL: Dense Gas Fraction and Star Formation Efficiency Across the Galactic Disk

By combining two surveys covering a large fraction of the molecular material in the Galactic disk we investigate the role the spiral arms play in the star formation process. We have matched clumps identified by ATLASGAL with their parental GMCs as identified by SEDIGISM, and use these giant molecular cloud (GMC) masses, the bolometric luminosities, and integrated clump masses obtained in a concurrent paper to estimate the dense gas fractions ($DGF_{gmc} = ∑M_{clump}/M_{gmc}$) and the instantaneous star forming efficiencies (i.e., $SFE_{gmc} = ∑L_{clump}/M_{gmc}$). We find that the molecular material associated with ATLASGAL clumps is concentrated in the spiral arms (∼60 per cent found within ±10 km s$^{−1}$ of an arm). We have searched for variations in the values of these physical parameters with respect to their proximity to the spiral arms, but find no evidence for any enhancement that might be attributable to the spiral arms. The combined results from a number of similar studies based on different surveys indicate that, while spiral-arm location plays a role in cloud formation and HI to H$_2$ conversion, the subsequent star formation processes appear to depend more on local environment effects. This leads us to conclude that the enhanced star formation activity seen towards the spiral arms is the result of source crowding rather than the consequence of a any physical process

Relative quadrupole moments of exotic shapes at ultrahigh spin in 154Er: Calibrating the TSD/SD puzzle

Transition quadrupole moments, Qt, of two ultrahigh-spin, collective structures in 154Er have been measured for the first time using the Doppler Shift Attenuation Method (DSAM). Data were acquired at the ATLAS accelerator facility of Argonne National Laboratory, using the Gammasphere detector array. A thick, gold-backed 110Pd foil was bombarded by a beam of 48Ti ions at 215 MeV. The Qt for each band was determined from the Doppler shift of gamma rays emitted by the resulting recoil nuclei. The extracted transition quadrupole moments are significantly different in magnitude, suggesting the two structures in 154Er represent distinct exotic nuclear shapes, namely axial superdeformed (SD) with Q t 20 eb, and triaxial strongly deformed (TSD) with Qt ≈ 11 eb. Indeed, the results calibrate the quadrupole moments of TSD bands recently measured in light erbium nuclei, 157,158Er

Quadrupole moments of coexisting collective shapes at high spin in 154Er

Four high-spin collective bands have been populated in 68154Er86 via the 110Pd(48Ti,4nγ)154Er reaction. Average transition quadrupole moments Qt have been measured for three of the bands by using the Doppler-shift attenuation method. The strongest band has a value of Q t=11.0±1.0eb, similar to values found recently for four triaxial strongly deformed (TSD) bands in 157,158Er. The second band has a value of Qt=19.5±3.2eb, consistent with a predicted axially symmetric superdeformed (SD) shape, similar in deformation to the 152Dy isotone, and is used as a calibration point. The third, new band has a value of Qt=9.9±2.2eb. The results confirm the unexpectedly large Qt moments for the favored TSD bands in light erbium isotopes

Large-scale Velocity-coherent Filaments in the SEDIGISM Survey: Association with Spiral Arms and Fraction of Dense Gas

Context. Filamentary structures in the interstellar medium are closely related to star formation. Dense gas mass fraction (DGMF) or clump formation efficiency in large-scale filaments possibly determine their hosting star formation activities. Aims. We aim to automatically identify large-scale filaments, characterize them, investigate their association with Galactic structures, and study their DGMFs. Methods. We use a modified minimum spanning tree (MST) algorithm to chain parsec-scale 13CO clumps previously extracted from the SEDIGISM (Structure, Excitation, and Dynamics of the Inner Galactic InterStellar Medium) survey. The MST connects nodes in a graph such that the sum of edge lengths is minimum. Modified MST also ensures velocity coherence between nodes, so the identified filaments are coherent in position-position-velocity (PPV) space. Results. We generate a catalog of 88 large-scale (> 10 pc) filaments in the inner Galactic plane (with −60◦ < l < 18◦ and |b| < 0.5◦). These SEDIGISM filaments are larger and less dense than MST filaments previously identified from the BGPS and ATLASGAL surveys. We find that eight of the filaments run along spiral arms and can be regarded as “bones” of the Milky Way. We also find three bones associated with the Local Spur in PPV space. By compiling 168 large-scale filaments with available DGMF across the Galaxy, an order of magnitude more than previously investigated, we find that DGMFs do not correlate with Galactic location, but bones have higher DGMFs than other filaments

Exploring the interstellar medium of NGC 891 at millimeter wavelengths using the NIKA2 camera

In the framework of the IMEGIN Large Program, we used the NIKA2 camera on the IRAM 30-m telescope to observe the edge-on galaxy NGC 891 at 1.15 mm and 2 mm and at a FWHM of 11.1" and 17.6", respectively. Multiwavelength data enriched with the new NIKA2 observations fitted by the HerBIE SED code (coupled with the THEMIS dust model) were used to constrain the physical properties of the ISM. Emission originating from the diffuse dust disk is detected at all wavelengths from mid-IR to mm, while mid-IR observations reveal warm dust emission from compact HII regions. Indications of mm excess emission have also been found in the outer parts of the galactic disk. Furthermore, our SED fitting analysis constrained the mass fraction of the small (< 15 Angstrom) dust grains. We found that small grains constitute 9.5% of the total dust mass in the galactic plane, but this fraction increases up to ~ 20% at large distances (|z| > 3 kpc) from the galactic plane.Comment: To appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Science