Source author record

Christian Ritter

Christian Ritter appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.SR astro-ph.GA Applications

Catalog footprint

What is connected

5works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2021arXiv

Nonparametric monitoring of sunspot number observations: a case study

Solar activity is an important driver of long-term climate trends and must be accounted for in climate models. Unfortunately, direct measurements of this quantity over long periods do not exist. The only observation related to solar activity whose records reach back to the seventeenth century are sunspots. Surprisingly, determining the number of sunspots consistently over time has remained until today a challenging statistical problem. It arises from the need of consolidating data from multiple observing stations around the world in a context of low signal-to-noise ratios, non-stationarity, missing data, non-standard distributions and many kinds of errors. The data from some stations experience therefore severe and various deviations over time. In this paper, we propose the first systematic and thorough statistical approach for monitoring these complex and important series. It consists of three steps essential for successful treatment of the data: smoothing on multiple timescales, monitoring using block bootstrap calibrated CUSUM charts and classifying of out-of-control situations by support vector techniques. This approach allows us to detect a wide range of anomalies (such as sudden jumps or more progressive drifts), unseen in previous analyses. It helps us to identify the causes of major deviations, which are often observer or equipment related. Their detection and identification will contribute to improve future observations. Their elimination or correction in past data will lead to a more precise reconstruction of the world reference index for solar activity: the International Sunspot Number.

preprint2020arXiv

Uncertainty quantification in sunspot counts

Observing and counting sunspots constitutes one of the longest-running scientific experiment, with first observations dating back to Galileo and the invention of the telescope around 1610. Today the sunspot number (SN) time series acts as a benchmark of solar activity in a large range of physical models. An appropriate statistical modelling, adapted to the time series' complex nature, is however still lacking. In this work, we provide the first comprehensive uncertainty quantification analysis of sunspot counts. Our interest lies in the following three components: the number of spots ($N_s$), the number of sunspot groups ($N_g$), and the composite $N_c$, defined as $N_c:=N_s+10N_g$. Those are reported by a network of observatories around the world, and are corrupted by errors of various types. We use a multiplicative framework to provide, for each of the three components, an estimation of their error distribution in various regimes (short-term, long-term, minima of solar activity). We also propose a robust estimator for the underlying solar signal and fit a density distribution that takes into account intrinsic characteristics such as over-dispersion, excess of zeros, and multiple modes. The estimation of the solar signal underlying the composite $N_c$ may be seen as a robust version of the International Sunspot Number (ISN), a quantity widely used as a proxy of solar activity. Therefore our results on $N_c$ may serve to characterize the uncertainty on ISN as well. Our results paves the way for a future monitoring of the observatories in quasi-real time, with the aim to alert the observers when they start deviating from the network and prevent large drifts from occurring in the network.

preprint2016arXiv

Mass and Metallicity Requirement in Stellar Models for Galactic Chemical Evolution Applications

We used a one-zone chemical evolution model to address the question of how many masses and metallicities are required in grids of massive stellar models in order to ensure reliable galactic chemical evolution predictions. We used a set of yields that includes seven masses between 13 and 30 Msun, 15 metallicities between 0 and 0.03 in mass fraction, and two different remnant mass prescriptions. We ran several simulations where we sampled subsets of stellar models to explore the impact of different grid resolutions. Stellar yields from low- and intermediate-mass stars and from Type Ia supernovae have been included in our simulations, but with a fixed grid resolution. We compared our results with the stellar abundances observed in the Milky Way for O, Na, Mg, Si, Ca, Ti, and Mn. Our results suggest that the range of metallicity considered is more important than the number of metallicities within that range, which only affects our numerical predictions by about 0.1 dex. We found that our predictions at [Fe/H] < -2 are very sensitive to the metallicity range and the mass sampling used for the lowest metallicity included in the set of yields. Variations between results can be as high as 0.8 dex, for any remnant mass prescription. At higher [Fe/H], we found that the required number of masses depends on the element of interest and on the remnant mass prescription. With a monotonic remnant mass prescription where every model explodes as a core-collapse supernova, the mass resolution induces variations of 0.2 dex on average. But with a remnant mass prescription that includes islands of non-explodability, the mass resolution can cause variations of about 0.2 to 0.7 dex depending on the choice of metallicity range. With such a prescription, explosive or non-explosive models can be missed if not enough masses are selected, resulting in over- or under-estimations of the mass ejected by massive stars.

preprint2016arXiv

Uncertainties in Galactic Chemical Evolution Models

We use a simple one-zone galactic chemical evolution model to quantify the uncertainties generated by the input parameters in numerical predictions, for a galaxy with properties similar to those of the Milky Way. We compiled several studies from the literature to gather the current constraints for our simulations regarding the typical value and uncertainty of seven basic parameters, which are: the lower and upper mass limits of the stellar initial mass function (IMF), the slope of the high-mass end of the stellar IMF, the slope of the delay-time distribution function of Type Ia supernovae (SNe Ia), the number of SNeIa per Msun formed, the total stellar mass formed, and the final mass of gas. We derived a probability distribution function to express the range of likely values for every parameter, which were then included in a Monte Carlo code to run several hundred simulations with randomly selected input parameters. This approach enables us to analyze the predicted chemical evolution of 16 elements in a statistical way by identifying the most probable solutions along with their 68% and 95% confidence levels. Our results show that the overall uncertainties are shaped by several input parameters that individually contribute at different metallicities, and thus at different galactic ages. The level of uncertainty then depends on the metallicity and is different from one element to another. Among the seven input parameters considered in this work, the slope of the IMF and the number of SNe Ia are currently the two main sources of uncertainty.

preprint2014arXiv

MESA and NuGrid simulations of classical novae: CO and ONe nova nucleosynthesis

Classical novae are the result of thermonuclear flashes of hydrogen accreted by CO or ONe white dwarfs, leading eventually to the dynamic ejection of the surface layers. These are observationally known to be enriched in heavy elements, such as C, O and Ne that must originate in layers below the H-flash convection zone. Building on our previous work, we now present stellar evolution simulations of ONe novae and provide a comprehensive comparison of our models with published ones. Some of our models include exponential convective boundary mixing to account for the observed enrichment of the nova ejecta even when accreted material has a solar abundance distribution. Our models produce maximum temperature evolution profiles and nucleosynthesis yields in good agreement with models that generate enriched ejecta by assuming that the accreted material was pre-mixed. We confirm for ONe novae the result we reported previously, i.e.\ we found that $^3$He could be produced {\it in situ} in solar-composition envelopes accreted with slow rates ($\dot{M} < 10^{-10}\,M_\odot/\mbox{yr}$) by cold ($T_{\rm WD} < 10^7$ K) CO WDs, and that convection was triggered by $^3$He burning before the nova outburst in that case. In addition, we now find that the interplay between the $^3$He production and destruction in the solar-composition envelope accreted with an intermediate rate, e.g.\ $\dot{M} = 10^{-10}\,M_\odot/\mbox{yr}$, by the $1.15\,M_\odot$ ONe WD with a relatively high initial central temperature, e.g.\ $T_{\rm WD} = 15\times 10^6$ K, leads to the formation of a thick radiative buffer zone that separates the bottom of the convective envelope from the WD surface. (Abridged)