Transient Earth system responses to cumulative carbon dioxide emissions: linearities, uncertainties, and probabilities in an observation-constrained model ensemble

Information on the relationship between cumulative fossil CO2 emissions and multiple climate targets is essential to design emission mitigation and climate adaptation strategies. In this study, the transient response of a climate or environmental variable per trillion tonnes of CO2 emissions, termed TRE, is quantified for a set of impact-relevant climate variables and from a large set of multi-forcing scenarios extended to year 2300 towards stabilization. An  ∼ 1000-member ensemble of the Bern3D-LPJ carbon–climate model is applied and model outcomes are constrained by 26 physical and biogeochemical observational data sets in a Bayesian, Monte Carlo-type framework. Uncertainties in TRE estimates include both scenario uncertainty and model response uncertainty. Cumulative fossil emissions of 1000 Gt C result in a global mean surface air temperature change of 1.9 °C (68 % confidence interval (c.i.): 1.3 to 2.7 °C), a decrease in surface ocean pH of 0.19 (0.18 to 0.22), and a steric sea level rise of 20 cm (13 to 27 cm until 2300). Linearity between cumulative emissions and transient response is high for pH and reasonably high for surface air and sea surface temperatures, but less pronounced for changes in Atlantic meridional overturning, Southern Ocean and tropical surface water saturation with respect to biogenic structures of calcium carbonate, and carbon stocks in soils. The constrained model ensemble is also applied to determine the response to a pulse-like emission and in idealized CO2-only simulations. The transient climate response is constrained, primarily by long-term ocean heat observations, to 1.7 °C (68 % c.i.: 1.3 to 2.2 °C) and the equilibrium climate sensitivity to 2.9 °C (2.0 to 4.2 °C). This is consistent with results by CMIP5 models but inconsistent with recent studies that relied on short-term air temperature data affected by natural climate variability.

Estimating HIV Incidence, Time to Diagnosis, and the Undiagnosed HIV Epidemic Using Routine Surveillance Data.

BACKGROUND Estimates of the size of the undiagnosed HIV-infected population are important to understand the HIV epidemic and to plan interventions, including "test-and-treat" strategies. METHODS We developed a multi-state back-calculation model to estimate HIV incidence, time between infection and diagnosis, and the undiagnosed population by CD4 count strata, using surveillance data on new HIV and AIDS diagnoses. The HIV incidence curve was modelled using cubic splines. The model was tested on simulated data and applied to surveillance data on men who have sex with men in The Netherlands. RESULTS The number of HIV infections could be estimated accurately using simulated data, with most values within the 95% confidence intervals of model predictions. When applying the model to Dutch surveillance data, 15,400 (95% confidence interval [CI] = 15,000, 16,000) men who have sex with men were estimated to have been infected between 1980 and 2011. HIV incidence showed a bimodal distribution, with peaks around 1985 and 2005 and a decline in recent years. Mean time to diagnosis was 6.1 (95% CI = 5.8, 6.4) years between 1984 and 1995 and decreased to 2.6 (2.3, 3.0) years in 2011. By the end of 2011, 11,500 (11,000, 12,000) men who have sex with men in The Netherlands were estimated to be living with HIV, of whom 1,750 (1,450, 2,200) were still undiagnosed. Of the undiagnosed men who have sex with men, 29% (22, 37) were infected for less than 1 year, and 16% (13, 20) for more than 5 years. CONCLUSIONS This multi-state back-calculation model will be useful to estimate HIV incidence, time to diagnosis, and the undiagnosed HIV epidemic based on routine surveillance data.

Estimating Particle Shape and Orientation Using Volume Tensors

We develop statistical procedures for estimating shape and orientation of arbitrary three-dimensional particles. We focus on the case where particles cannot be observed directly, but only via sections. Volume tensors are used for describing particle shape and orientation, and we derive stereological estimators of the tensors. These estimators are combined to provide consistent estimators of the moments of the so-called particle cover density. The covariance structure associated with the particle cover density depends on the orientation and shape of the particles. For instance, if the distribution of the typical particle is invariant under rotations, then the covariance matrix is proportional to the identity matrix. We develop a non-parametric test for such isotropy. A flexible Lévy-based particle model is proposed, which may be analysed using a generalized method of moments in which the volume tensors enter. The developed methods are used to study the cell organization in the human brain cortex.

Estimating population and livestock density of mobile pastoralists and sedentary settlements in the southeastern Lake Chad area

Mobile pastoralists provide major contributions to the gross domestic product in Chad, but little information is available regarding their demography. The Lake Chad area population is increasing, resulting in competition for scarce land and water resources. For the first time, the density of people and animals from mobile and sedentary populations was assessed using randomly defined sampling areas. Four sampling rounds were conducted over two years in the same areas to show population density dynamics. We identified 42 villages of sedentary communities in the sampling zones; 11 (in 2010) and 16 (in 2011) mobile pastoralist camps at the beginning of the dry season and 34 (in 2011) and 30 (in 2012) camps at the end of the dry season. A mean of 64.0 people per km2 (95% confidence interval, 20.3-107.8) were estimated to live in sedentary villages. In the mobile communities, we found 5.9 people per km2 at the beginning and 17.5 people per km2 at the end of the dry season. We recorded per km2 on average 21.0 cattle and 31.6 small ruminants in the sedentary villages and 66.1 cattle and 102.5 small ruminants in the mobile communities, which amounts to a mean of 86.6 tropical livestock units during the dry season. These numbers exceed, by up to five times, the published carrying capacities for similar Sahelian zones. Our results underline the need for a new institutional framework. Improved land use management must equally consider the needs of mobile communities and sedentary populations.

Fast Update of Conditional Simulation Ensembles

Gaussian random field (GRF) conditional simulation is a key ingredient in many spatial statistics problems for computing Monte-Carlo estimators and quantifying uncertainties on non-linear functionals of GRFs conditional on data. Conditional simulations are known to often be computer intensive, especially when appealing to matrix decomposition approaches with a large number of simulation points. This work studies settings where conditioning observations are assimilated batch sequentially, with one point or a batch of points at each stage. Assuming that conditional simulations have been performed at a previous stage, the goal is to take advantage of already available sample paths and by-products to produce updated conditional simulations at mini- mal cost. Explicit formulae are provided, which allow updating an ensemble of sample paths conditioned on n ≥ 0 observations to an ensemble conditioned on n + q observations, for arbitrary q ≥ 1. Compared to direct approaches, the proposed formulae proveto substantially reduce computational complexity. Moreover, these formulae explicitly exhibit how the q new observations are updating the old sample paths. Detailed complexity calculations highlighting the benefits of this approach with respect to state-of-the-art algorithms are provided and are complemented by numerical experiments.

Estimating whole-lake fish catch per unit effort

The European standard for gillnetsampling to characterize lake fish communities stratifies sampling effort (i.e., number of nets) within depth strata. Nets to sample benthic habitats are randomly distributed throughout the lake within each depth strata. Pelagic nets are also stratified by depth, but are set only at the deepest point of the lake. Multiple authors have suggested that this design under-represents pelagic habitats, resulting in estimates of whole-lake CPUE and community composition which are disproportionately influenced by ecological conditions of littoral and benthic habitats. To address this issue, researchers have proposed estimating whole-lake CPUE by weighting the catch rate in each depth-compartment by the proportion of the volume of the lake contributed by the compartment. Our study aimed to assess the effectiveness of volume-weighting by applying it to fish communities sampled according to the European standard (CEN), and by a second whole-lake gillnetting protocol (VERT), which prescribes additional fishing effort in pelagic habitats. We assume that convergence between the protocols indicates that volume-weighting provides a more accurate estimate of whole-lake catch rate and community composition. Our results indicate that volume-weighting improves agreement between the protocols for whole-lake total CPUE, estimated proportion of perch and roach and the overall fish community composition. Discrepancies between the protocols remaining after volume-weighting maybe because sampling under the CEN protocol overlooks horizontal variation in pelagic fish communities. Analyses based on multiple pelagic-set VERT nets identified gradients in the density and biomass of pelagic fish communities in almost half the lakes that corresponded with the depth of water at net-setting location and distance along the length of a lake. Additional CEN pelagic sampling effort allocated across water depths and distributed throughout the lake would therefore help to reconcile differences between the sampling protocols and, in combination with volume-weighting, converge on a more accurate estimate of whole-lake fish communities.

Associating optical measurements and estimating orbits of geocentric objects through population-based meta-heuristic methods

Currently several thousands of objects are being tracked in the MEO and GEO regions through optical means. The problem faced in this framework is that of Multiple Target Tracking (MTT). In this context both, the correct associations among the observations and the orbits of the objects have to be determined. The complexity of the MTT problem is defined by its dimension S. The number S corresponds to the number of fences involved in the problem. Each fence consists of a set of observations where each observation belongs to a different object. The S ≥ 3 MTT problem is an NP-hard combinatorial optimization problem. There are two general ways to solve this. One way is to seek the optimum solution, this can be achieved by applying a branch-and- bound algorithm. When using these algorithms the problem has to be greatly simplified to keep the computational cost at a reasonable level. Another option is to approximate the solution by using meta-heuristic methods. These methods aim to efficiently explore the different possible combinations so that a reasonable result can be obtained with a reasonable computational effort. To this end several population-based meta-heuristic methods are implemented and tested on simulated optical measurements. With the advent of improved sensors and a heightened interest in the problem of space debris, it is expected that the number of tracked objects will grow by an order of magnitude in the near future. This research aims to provide a method that can treat the correlation and orbit determination problems simultaneously, and is able to efficiently process large data sets with minimal manual intervention.

Estimating Intra-Party Preferences: Comparing Speeches to Votes

Well-established methods exist for measuring party positions, but reliable means for estimating intra-party preferences remain underdeveloped. While most efforts focus on estimating the ideal points of individual legislators based on inductive scaling of roll call votes, this data suffers from two problems: selection bias due to unrecorded votes and strong party discipline, which tends to make voting a strategic rather than a sincere indication of preferences. By contrast, legislative speeches are relatively unconstrained, as party leaders are less likely to punish MPs for speaking freely as long as they vote with the party line. Yet, the differences between roll call estimations and text scalings remain essentially unexplored, despite the growing application of statistical analysis of textual data to measure policy preferences. Our paper addresses this lacuna by exploiting a rich feature of the Swiss legislature: on most bills, legislators both vote and speak many times. Using this data, we compare text-based scaling of ideal points to vote-based scaling from a crucial piece of energy legislation. Our findings confirm that text scalings reveal larger intra-party differences than roll calls. Using regression models, we further explain the differences between roll call and text scalings by attributing differences to constituency-level preferences for energy policy.

Loss on ignition as a method for estimating organic and carbonate content in sediments: reproducibility and comparability of results

Five test runs were performed to assess possible bias when performing the loss on ignition (LOI) method to estimate organic matter and carbonate content of lake sediments. An accurate and stable weight loss was achieved after 2 h of burning pure CaCO3 at 950 °C, whereas LOI of pure graphite at 530 °C showed a direct relation to sample size and exposure time, with only 40-70% of the possible weight loss reached after 2 h of exposure and smaller samples losing weight faster than larger ones. Experiments with a standardised lake sediment revealed a strong initial weight loss at 550 °C, but samples continued to lose weight at a slow rate at exposure of up to 64 h, which was likely the effect of loss of volatile salts, structural water of clay minerals or metal oxides, or of inorganic carbon after the initial burning of organic matter. A further test-run revealed that at 550 °C samples in the centre of the furnace lost more weight than marginal samples. At 950 °C this pattern was still apparent but the differences became negligible. Again, LOI was dependent on sample size. An analytical LOI quality control experiment including ten different laboratories was carried out using each laboratory's own LOI procedure as well as a standardised LOI procedure to analyse three different sediments. The range of LOI values between laboratories measured at 550 °C was generally larger when each laboratory used its own method than when using the standard method. This was similar for 950 °C, although the range of values tended to be smaller. The within-laboratory range of LOI measurements for a given sediment was generally small. Comparisons of the results of the individual and the standardised method suggest that there is a laboratory-specific pattern in the results, probably due to differences in laboratory equipment and/or handling that could not be eliminated by standardising the LOI procedure. Factors such as sample size, exposure time, position of samples in the furnace and the laboratory measuring affected LOI results, with LOI at 550 °C being more susceptible to these factors than LOI at 950 °C. We, therefore, recommend analysts to be consistent in the LOI method used in relation to the ignition temperatures, exposure times, and the sample size and to include information on these three parameters when referring to the method.

Estimating the benefit of a second bone anchored hearing implant in unilaterally implanted users with a testband.

Conclusion Using a second bone anchored hearing implant (BAHI) mounted on a testband in unilaterally implanted BAHI users to test its potential advantage pre-operatively under-estimates the advantage of two BAHIs placed on two implants. Objectives To investigate how well speech understanding with a second BAHI mounted on a testband approaches the benefit of bilaterally implanted BAHIs. Method Prospective study with 16 BAHI users. Eight were implanted unilaterally (group A) and eight were implanted bilaterally (group B). Aided speech understanding was measured. Speech was presented from the front and noise came either from the left, right, or from the front in two conditions for group A (with one BAHI, and with two BAHIs, where the second device was mounted on a testband) and in three conditions for group B (same two conditions as group A, and in addition with both BAHIs mounted on implants). Results Speech understanding in noise improved with the additional device for noise from the side of the first BAHI (+0.7 to +2.1 dB) and decreased for noise from the other side (-1.8 dB to -3.9 dB). Improvements were highest (+2.1 dB, p = 0.016) and disadvantages were smallest (-1.8 dB, p = 0.047) with both BAHIs mounted on implants. Testbands yielded smaller advantages and higher disadvantages of the additional BAHI (average difference = -0.9 dB).

Binary recombinase systems for high-resolution conditional mutagenesis.

Conditional mutagenesis using Cre recombinase expressed from tissue specific promoters facilitates analyses of gene function and cell lineage tracing. Here, we describe two novel dual-promoter-driven conditional mutagenesis systems designed for greater accuracy and optimal efficiency of recombination. Co-Driver employs a recombinase cascade of Dre and Dre-respondent Cre, which processes loxP-flanked alleles only when both recombinases are expressed in a predetermined temporal sequence. This unique property makes Co-Driver ideal for sequential lineage tracing studies aimed at unraveling the relationships between cellular precursors and mature cell types. Co-InCre was designed for highly efficient intersectional conditional transgenesis. It relies on highly active trans-splicing inteins and promoters with simultaneous transcriptional activity to reconstitute Cre recombinase from two inactive precursor fragments. By generating native Cre, Co-InCre attains recombination rates that exceed all other binary SSR systems evaluated in this study. Both Co-Driver and Co-InCre significantly extend the utility of existing Cre-responsive alleles.

Estimating the uptake of traffic-derived NO2 from 15N abundance in Norway spruce needles

The 15N ratio of nitrogen oxides (NOx) emitted from vehicles, measured in the air adjacent to a highway in the Swiss Middle Land, was very high [δ15N(NO2) = +5.7‰]. This high 15N abundance was used to estimate long-term NO2 dry deposition into a forest ecosystem by measuring δ15N in the needles and the soil of potted and autochthonous spruce trees [Picea abies (L.) Karst] exposed to NO2 in a transect orthogonal to the highway. δ15N in the current-year needles of potted trees was 2.0‰ higher than that of the control after 4 months of exposure close to the highway, suggesting a 25% contribution to the N-nutrition of these needles. Needle fall into the pots was prevented by grids placed above the soil, while the continuous decomposition of needle litter below the autochthonous trees over previous years has increased δ15N values in the soil, resulting in parallel gradients of δ15N in soil and needles with distance from the highway. Estimates of NO2 uptake into needles obtained from the δ15N data were significantly correlated with the inputs calculated with a shoot gas exchange model based on a parameterisation widely used in deposition modelling. Therefore, we provide an indication of estimated N inputs to forest ecosystems via dry deposition of NO2 at the receptor level under field conditions.

Analysis of BMP signalling through the receptor type IA in embryogenesis using conditional gene inactivation in mice

Although bone morphogenetic proteins (BMPs) were initially identified for their potent bone-inducing activity, their precise roles in processes of endochondral and intramembranous bone formation are far from being clear. Tissue-specific loss-of-function experiments using the BMP receptor type IA (BMPR-IA) are particularly attractive since this receptor is thought to be essential for signaling by the closely related BMPs -2, 4, and 7. To ablate signaling through this receptor during chondrogenesis, we have generated transgenic mice expressing Cre recombinase under the control of the collagen type II (Col2a1) gene regulatory sequences. Mice lacking BMPR-IA function in chondrocytes display a number of skeletal abnormalities, including defects in bones of the chondrocranium, abnormal dorsal vertebral processes, scapulae with severe hypoplasia of dorsal elements, and shortening of the long bones. Alterations in the growth plate of long bones in mutants suggest that BMPR-IA is not required for early steps of the chondrocyte specification, but is rather important in regulation of terminal differentiation. Molecular analysis revealed noticeable downregulation of the Ihh/Ptch signalling pathway, decreased chondrocyte proliferation rate and deregulation of hypertrophy. ^ In order to elucidate the role of BMP signalling in development of the limb and intramembranous ossification, we have used mice expressing Cre recombinase under control of the Prx1 (MHox) regulatory elements (M. Logan, pers comm.). Cre activity was found in those mice in the developing limb bud mesenchyme, as well as in a subset of cranial neural crest cells. Prx1-Cre-induced conditional mutants display prominent defects in distal limb outgrowth, as well as ossification defects in a number of neural crest-derived calvarial bones. Intriguingly, mutant limbs displayed alterations in patterning along all three axes. Molecular analysis revealed ectopic anterior Shh/Ptch signalling pathway activation and expression of some Hox genes. Observed loss of Msx1 and Msx2 expression in the progress zone correlates with downregulation of Cyclin D1 and decreased distal outgrowth. Abnormal ventral localization of Lmx1b-expressing cells along with observed later morphological abnormalities suggest a novel role for BMP signalling in establishment or maintaining of the dorso-ventral polarity in the limb mesoderm. ^

Conditional expression of Sry, and Wnt4 by gene-targeting: Studies in sexual and skeletal development

Sry and Wnt4 cDNAs were individually introduced into the ubiquitously-expressed Rosa26 ( R26) locus by gene targeting in embryonic stem (ES) cells to create a conditional gene expression system in mice. In the targeted alleles, expression of these cDNAs should be blocked by a neomycin resistance selection cassette that is flanked by loxP sites. Transgene expression should be activated after the blocking cassette is deleted by Cre recombinase. ^ To test this conditional expression system, I have bred R26-stop- Sry and R26-stop-Wnt4 heterozygotes with a MisRII-Cre mouse line that expresses Cre in the gonads of both sexes. Analysis of these two types of bigenic heterozygotes indicated that their gonads developed normally like those of wild types. However, one XX R26-Sry/R26-Sry; MisR2-Cre/+ showed epididymis-like structures resembling those of males. In contrast, only normal phenotypes were observed in XY R26-Wnt4/R26-Wnt4; MisR2-Cre /+ mice. To interpret these results, I have tested for Cre recombinase activity by Southern blot and transcription of the Sry and Wnt4 transgenes by RT-PCR. Results showed that bigenic mutants had insufficient activation of the transgenes in their gonads at E12.5 and E13.5. Therefore, the failure to observe mutant phenotypes may have resulted from low activity of MisR2-Cre recombination at the appropriate time. ^ Col2a1-Cre transgenic mice express Cre in differentiating chondrocytes. R26-Wnt4; Col2a1-Cre bigenic heterozygous mice were found to exhibit a dramatic alteration in growth presumably caused by Wnt4 overexpression during chondrogenesis. R26-Wnt4; Col2a1-Cre mice exhibited dwarfism beginning approximately 10 days after birth. In addition, they also had craniofacial abnormalities, and had delayed ossification of the lumbar vertebrate and pelvic bones. Histological analysis of the growth plates of R26-Wnt4; Col2a1-Cre mice revealed less structural organization and a delay in onset of the primary and secondary ossification centers. Molecular studies confirmed that overexpression of Wnt4 causes decreased proliferation and early maturation of chondrocytes. In addition, R26-Wnt4; Col2a1-Cre mice had decreased expression of vascular endothelial growth factor (VEGF), suggesting that defects in vascularization may contribute to the dwarf phenotype. Finally, 9-month-old R26-Wnt4; Col2a1-Cre mice had significantly more fat cells in the marrow cavities of their metaphysis long bones, implying that long-term overexpression of Wnt4may cause bone marrow pathologies. In conclusion, Wnt4 was activated by Col2a1-Cre recombinase and was overexpressed in the growth plate, resulting in aberrant proliferation and differentiation of chondrocytes, and ultimately leads to dwarfism in mice. ^

Empirical Analysis of Credit Risk Regime Switching and Temporal Conditional Default Correlation in Credit Default Swap Valuation: The Market liquidity effect

In this paper, we extend the debate concerning Credit Default Swap valuation to include time varying correlation and co-variances. Traditional multi-variate techniques treat the correlations between covariates as constant over time; however, this view is not supported by the data. Secondly, since financial data does not follow a normal distribution because of its heavy tails, modeling the data using a Generalized Linear model (GLM) incorporating copulas emerge as a more robust technique over traditional approaches. This paper also includes an empirical analysis of the regime switching dynamics of credit risk in the presence of liquidity by following the general practice of assuming that credit and market risk follow a Markov process. The study was based on Credit Default Swap data obtained from Bloomberg that spanned the period January 1st 2004 to August 08th 2006. The empirical examination of the regime switching tendencies provided quantitative support to the anecdotal view that liquidity decreases as credit quality deteriorates. The analysis also examined the joint probability distribution of the credit risk determinants across credit quality through the use of a copula function which disaggregates the behavior embedded in the marginal gamma distributions, so as to isolate the level of dependence which is captured in the copula function. The results suggest that the time varying joint correlation matrix performed far superior as compared to the constant correlation matrix; the centerpiece of linear regression models.