Poisson -- Boltzmann Brownian Dynamics of Charged Colloids in Suspension

We describe a method to simulate the dynamics of charged colloidal particles suspended in a liquid containing dissociated ions and salt ions. Regimes of prime current interest are those of large volume fraction of colloids, highly charged particles and low salt concentrations. A description which is tractable under these conditions is obtained by treating the small dissociated and salt ions as continuous fields, while keeping the colloidal macroions as discrete particles. For each spatial configuration of the macroions, the electrostatic potential arising from all charges in the system is determined by solving the nonlinear Poisson--Boltzmann equation. From the electrostatic potential, the forces acting on the macroions are calculated and used in a Brownian dynamics simulation to obtain the motion of the latter. The method is validated by comparison to known results in a parameter regime where the effective interaction between the macroions is of a pairwise Yukawa form

Five decades of northern land carbon uptake revealed by the interhemispheric CO2 gradient

The global land and ocean carbon sinks have increased proportionally with increasing carbon dioxide emissions during the past decades 1 . It is thought that Northern Hemisphere lands make a dominant contribution to the global land carbon sink 2–7 ; however, the long-term trend of the northern land sink remains uncertain. Here, using measurements of the interhemispheric gradient of atmospheric carbon dioxide from 1958 to 2016, we show that the northern land sink remained stable between the 1960s and the late 1980s, then increased by 0.5 ± 0.4 petagrams of carbon per year during the 1990s and by 0.6 ± 0.5 petagrams of carbon per year during the 2000s. The increase of the northern land sink in the 1990s accounts for 65% of the increase in the global land carbon flux during that period. The subsequent increase in the 2000s is larger than the increase in the global land carbon flux, suggesting a coincident decrease of carbon uptake in the Southern Hemisphere. Comparison of our findings with the simulations of an ensemble of terrestrial carbon models 5,8 over the same period suggests that the decadal change in the northern land sink between the 1960s and the 1990s can be explained by a combination of increasing concentrations of atmospheric carbon dioxide, climate variability and changes in land cover. However, the increase during the 2000s is underestimated by all models, which suggests the need for improved consideration of changes in drivers such as nitrogen deposition, diffuse light and land-use change. Overall, our findings underscore the importance of Northern Hemispheric land as a carbon sink

Five decades of northern land carbon uptake revealed by the interhemispheric CO2 gradient

The global land and ocean carbon sinks have increased proportionally with increasing carbon dioxide emissions during the past decades 1 . It is thought that Northern Hemisphere lands make a dominant contribution to the global land carbon sink 2–7 ; however, the long-term trend of the northern land sink remains uncertain. Here, using measurements of the interhemispheric gradient of atmospheric carbon dioxide from 1958 to 2016, we show that the northern land sink remained stable between the 1960s and the late 1980s, then increased by 0.5 ± 0.4 petagrams of carbon per year during the 1990s and by 0.6 ± 0.5 petagrams of carbon per year during the 2000s. The increase of the northern land sink in the 1990s accounts for 65% of the increase in the global land carbon flux during that period. The subsequent increase in the 2000s is larger than the increase in the global land carbon flux, suggesting a coincident decrease of carbon uptake in the Southern Hemisphere. Comparison of our findings with the simulations of an ensemble of terrestrial carbon models 5,8 over the same period suggests that the decadal change in the northern land sink between the 1960s and the 1990s can be explained by a combination of increasing concentrations of atmospheric carbon dioxide, climate variability and changes in land cover. However, the increase during the 2000s is underestimated by all models, which suggests the need for improved consideration of changes in drivers such as nitrogen deposition, diffuse light and land-use change. Overall, our findings underscore the importance of Northern Hemispheric land as a carbon sink

The consolidated European synthesis of CO2 emissions and removals for the European Union and United Kingdom: 1990-2018

Reliable quantification of the sources and sinks of atmospheric carbon dioxide (CO2), including that of their trends and uncertainties, is essential to monitoring the progress in mitigating anthropogenic emissions under the Kyoto Protocol and the Paris Agreement. This study provides a consolidated synthesis of estimates for all anthropogenic and natural sources and sinks of CO2 for the European Union and UK (EU27 + UK), derived from a combination of state-of-the-art bottom-up (BU) and top-down (TD) data sources and models. Given the wide scope of the work and the variety of datasets involved, this study focuses on identifying essential questions which need to be answered to properly understand the differences between various datasets, in particular with regards to the less-well-characterized fluxes from managed ecosystems. The work integrates recent emission inventory data, process-based ecosystem model results, data-driven sector model results and inverse modeling estimates over the period 1990-2018. BU and TD products are compared with European national greenhouse gas inventories (NGHGIs) reported under the UNFCCC in 2019, aiming to assess and understand the differences between approaches. For the uncertainties in NGHGIs, we used the standard deviation obtained by varying parameters of inventory calculations, reported by the member states following the IPCC Guidelines. Variation in estimates produced with other methods, like atmospheric inversion models (TD) or spatially disaggregated inventory datasets (BU), arises from diverse sources including within-model uncertainty related to parameterization as well as structural differences between models. In comparing NGHGIs with other approaches, a key source of uncertainty is that related to different system boundaries and emission categories (CO2 fossil) and the use of different land use definitions for reporting emissions from land use, land use change and forestry (LULUCF) activities (CO2 land). At the EU27 + UK level, the NGHGI (2019) fossil CO2 emissions (including cement production) account for 2624 Tg CO2 in 2014 while all the other seven bottom-up sources are consistent with the NGHGIs and report a mean of 2588 (± 463 Tg CO2). The inversion reports 2700 Tg CO2 (± 480 Tg CO2), which is well in line with the national inventories. Over 2011-2015, the CO2 land sources and sinks from NGHGI estimates report-90 Tg C yr-1 ± 30 Tg C yr-1 while all other BU approaches report a mean sink of-98 Tg C yr-1 (± 362 Tg of C from dynamic global vegetation models only). For the TD model ensemble results, we observe a much larger spread for regional inversions (i.e., mean of 253 Tg C yr-1 ± 400 Tg C yr-1). This concludes that (a) current independent approaches are consistent with NGHGIs and (b) their uncertainty is too large to allow a verification because of model differences and probably also because of the definition of "CO2 flux"obtained from different approaches. The referenced datasets related to figures are visualized. © 2021 Ana Maria Roxana Petrescu et al

Current systematic carbon-cycle observations and the need for implementing a policy-relevant carbon observing system

A globally integrated carbon observation and analysis system is needed to improve the fundamental understanding of the global carbon cycle, to improve our ability to project future changes, and to verify the effectiveness of policies aiming to reduce greenhouse gas emissions and increase carbon sequestration. Building an integrated carbon observation system requires transformational advances from the existing sparse, exploratory framework towards a dense, robust, and sustained system in all components: anthropogenic emissions, the atmosphere, the ocean, and the terrestrial biosphere. The paper is addressed to scientists, policymakers, and funding agencies who need to have a global picture of the current state of the (diverse) carbon observations. We identify the current state of carbon observations, and the needs and notional requirements for a global integrated carbon observation system that can be built in the next decade. A key conclusion is the substantial expansion of the ground-based observation networks required to reach the high spatial resolution for CO₂ and CH₄ fluxes, and for carbon stocks for addressing policy-relevant objectives, and attributing flux changes to underlying processes in each region. In order to establish flux and stock diagnostics over areas such as the southern oceans, tropical forests, and the Arctic, in situ observations will have to be complemented with remote-sensing measurements. Remote sensing offers the advantage of dense spatial coverage and frequent revisit. A key challenge is to bring remote-sensing measurements to a level of long-term consistency and accuracy so that they can be efficiently combined in models to reduce uncertainties, in synergy with ground-based data. Bringing tight observational constraints on fossil fuel and land use change emissions will be the biggest challenge for deployment of a policy-relevant integrated carbon observation system. This will require in situ and remotely sensed data at much higher resolution and density than currently achieved for natural fluxes, although over a small land area (cities, industrial sites, power plants), as well as the inclusion of fossil fuel CO₂ proxy measurements such as radiocarbon in CO₂ and carbon-fuel combustion tracers. Additionally, a policy-relevant carbon monitoring system should also provide mechanisms for reconciling regional top-down (atmosphere-based) and bottom-up (surface-based) flux estimates across the range of spatial and temporal scales relevant to mitigation policies. In addition, uncertainties for each observation data-stream should be assessed. The success of the system will rely on long-term commitments to monitoring, on improved international collaboration to fill gaps in the current observations, on sustained efforts to improve access to the different data streams and make databases interoperable, and on the calibration of each component of the system to agreed-upon international scales

Seismic Constraints on the Thickness and Structure of the Martian Crust from InSight

NASA¿s InSight mission [1] has for the first time placed a very broad-band seismometer on the surface of Mars. The Seismic Experiment for Interior Structure (SEIS) [2] has been collecting continuous data since early February 2019. The main focus of InSight is to enhance our understanding of the internal structure and dynamics of Mars, which includes the goal to better constrain the crustal thickness of the planet [3]. Knowing the present-day crustal thickness of Mars has important implications for its thermal evolution [4] as well as for the partitioning of silicates and heat-producing elements between the different layers of Mars. Current estimates for the crustal thickness of Mars are based on modeling the relationship between topography and gravity [5,6], but these studies rely on different assumptions, e.g. on the density of the crust and upper mantle, or the bulk silicate composition of the planet and the crust. The resulting values for the average crustal thickness differ by more than 100%, from 30 km to more than 100 km [7]. New independent constraints from InSight will be based on seismically determining the crustal thickness at the landing site. This single firm measurement of crustal thickness at one point on the planet will allow to constrain both the average crustal thickness of Mars as well as thickness variations across the planet when combined with constraints from gravity and topography [8]. Here we describe the determination of the crustal structure and thickness at the InSight landing site based on seismic receiver functions for three marsquakes compared with autocorrelations of InSight data [9].We acknowledge NASA, CNES, partner agencies and institutions (UKSA, SSO,DLR, JPL, IPGP-CNRS, ETHZ, IC, MPS-MPG) and the operators of JPL, SISMOC, MSDS, IRIS-DMC and PDS for providing SEED SEIS data. InSight data is archived in the PDS, and a full list of archives in the Geosciences, Atmospheres, and Imaging nodes is at https://pds-geosciences.wustl.edu/missions/insight/. This work was partially carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration. ©2021, California Institute of Technology. Government sponsorship acknowledge

Whole-genome phylogeography of the intertidal snail Littorina saxatilis

Understanding the factors that have shaped the current distributions and diversity of species is a central and longstanding aim of evolutionary biology. The recent inclusion of genomic data into phylogeographic studies has dramatically improved our understanding in organisms where evolutionary relationships have been challenging to infer. We used whole-genome sequences to study the phylogeography of the intertidal snail Littorina saxatilis, which has successfully colonized and diversified across a broad range of coastal environments in the Northern Hemisphere amid repeated cycles of glaciation. Building on past studies based on short DNA sequences, we used genome-wide data to provide a clearer picture of the relationships among samples spanning most of the species natural range. Our results confirm the trans-Atlantic colonization of North America from Europe, and have allowed us to identify rough locations of glacial refugia and to infer likely routes of colonization within Europe. We also investigated the signals in different datasets to account for the effects of genomic architecture and non-neutral evolution, which provides new insights about diversification of four ecotypes of L. saxatilis (the crab, wave, barnacle, and brackish ecotypes) at different spatial scales. Overall, we provide a much clearer picture of the biogeography of L. saxatilis, providing a foundation for more detailed phylogenomic and demographic studies

Retracing the history and planning the future of the red squirrel (Sciurus vulgaris) in Ireland using non-invasive genetics

The Eurasian red squirrel’s (Sciurus vulgaris) history in Ireland is largely unknown, but the original population is thought to have been driven to extinction by humans in the 17th Century, and multiple records exist for its subsequent reintroduction in the 19th 4 Century. However, it is currently unknown how these reintroductions affect the red squirrel population today, or may do so in the future. In this study, we report on the development of a DNA toolkit for the non-invasive genetic study of the red squirrel. Non-invasively collected red squirrel samples were combined with other samples collected throughout Ireland and previously published mitochondrial DNA (mtDNA) data from Ireland, Great Britain and continental Europe to give an insight into population genetics and historical introductions of the red squirrel in Ireland. Our findings demonstrate that the Irish red squirrel population is on a national scale quite genetically diverse, but at a local level contains relatively low levels of genetic diversity and evidence of genetic structure. This is likely an artefact of the introduction of a small number of genetically similar animals to specific sites. A lack of continuous woodland cover in Ireland has prevented further mixing with animals of different origins that may have been introduced even to neighbouring sites. Consequently, some of these genetically isolated populations are or may in the future be at risk of extinction. The Irish red squirrel population contains mtDNA haplotypes of both a British and Continental European origin, the former of which are now extinct or simply not recorded in contemporary Great Britain. The Irish population is therefore important in terms of red squirrel conservation not only in Ireland, but also for Great Britain, and should be appropriately managed

Inferring dispersal across a fragmented landscape using reconstructed families in the Glanville fritillary butterfly

Dispersal is important for determining both species ecological processes, such as population viability, and its evolutionary processes, like gene flow and local adaptation. Yet obtaining accurate estimates in the wild through direct observation can be challenging or even impossible, particularly over large spatial and temporal scales. Genotyping many individuals from wild populations can provide detailed inferences about dispersal. We therefore utilized genomewide marker data to estimate dispersal in the classic metapopulation of the Glanville fritillary butterfly (Melitaea cinxia L.), in the Aland Islands in SW Finland. This is an ideal system to test the effectiveness of this approach due to the wealth of information already available covering dispersal across small spatial and temporal scales, but lack of information at larger spatial and temporal scales. We sampled three larvae per larval family group from 3732 groups over a six-year period and genotyped for 272 SNPs across the genome. We used this empirical data set to reconstruct cases where full-sibs were detected in different local populations to infer female effective dispersal distance, that is, dispersal events directly contributing to gene flow. On average this was one kilometre, closely matching previous dispersal estimates made using direct observation. To evaluate our power to detect full-sib families, we performed forward simulations using an individual-based model constructed and parameterized for the Glanville fritillary metapopulation. Using these simulations, 100% of predicted full-sibs were correct and over 98% of all true full-sib pairs were detected. We therefore demonstrate that even in a highly dynamic system with a relatively small number of markers, we can accurately reconstruct full-sib families and for the first time make inferences on female effective dispersal. This highlights the utility of this approach in systems where it has previously been impossible to obtain accurate estimates of dispersal over both ecological and evolutionary scales.Peer reviewe