Partial reconstruction of the trajectories of a discretely observed branching diffusion with immigration and an application to inference

In this treatise we consider finite systems of branching particles where the particles move independently of each other according to d-dimensional diffusions. Particles are killed at a position dependent rate, leaving at their death position a random number of descendants according to a position dependent reproduction law. In addition particles immigrate at constant rate (one immigrant per immigration time). A process with above properties is called a branching diffusion withimmigration (BDI). In the first part we present the model in detail and discuss the properties of the BDI under our basic assumptions. In the second part we consider the problem of reconstruction of the trajectory of a BDI from discrete observations. We observe positions of the particles at discrete times; in particular we assume that we have no information about the pedigree of the particles. A natural question arises if we want to apply statistical procedures on the discrete observations: How can we find couples of particle positions which belong to the same particle? We give an easy to implement 'reconstruction scheme' which allows us to redraw or 'reconstruct' parts of the trajectory of the BDI with high accuracy. Moreover asymptotically the whole path can be reconstructed. Further we present simulations which show that our partial reconstruction rule is tractable in practice. In the third part we study how the partial reconstruction rule fits into statistical applications. As an extensive example we present a nonparametric estimator for the diffusion coefficient of a BDI where the particles move according to one-dimensional diffusions. This estimator is based on the Nadaraya-Watson estimator for the diffusion coefficient of one-dimensional diffusions and it uses the partial reconstruction rule developed in the second part above. We are able to prove a rate of convergence of this estimator and finally we present simulations which show that the estimator works well even if we leave our set of assumptions.

High-resolution reconstruction of Holocene climate variability and environmental changes in the North Pacific using bivalve shells

Bivalve mollusk shells are useful tools for multi-species and multi-proxy paleoenvironmental reconstructions with a high temporal and spatial resolution. Past environmental conditions can be reconstructed from shell growth and stable oxygen and carbon isotope ratios, which present an archive for temperature, freshwater fluxes and primary productivity. The purpose of this thesis is the reconstruction of Holocene climate and environmental variations in the North Pacific with a high spatial and temporal resolution using marine bivalve shells. This thesis focuses on several different Holocene time periods and multiple regions in the North Pacific, including: Japan, Alaska (AK), British Columbia (BC) and Washington State, which are affected by the monsoon, Pacific Decadal Oscillation (PDO) and El Niño/Southern Oscillation (ENSO). Such high-resolution proxy data from the marine realm of mid- and high-latitudes are still rare. Therefore, this study contributes to the optimization and verification of climate models. However, before using bivalves for environmental reconstructions and seasonality studies, life history traits must be well studied to temporally align and interpret the geochemical record. These calibration studies are essential to ascertain the usefulness of selected bivalve species as paleoclimate proxy archives. This work focuses on two bivalve species, the short-lived Saxidomus gigantea and the long-lived Panopea abrupta. Sclerochronology and oxygen isotope ratios of different shell layers of P. abrupta were studied in order to test the reliability of this species as a climate archive. The annual increments are clearly discernable in umbonal shell portions and the increments widths should be measured in these shell portions. A reliable reconstruction of paleotemperatures may only be achieved by exclusively sampling the outer shell layer of multiple contemporaneous specimens. Life history traits (e.g., timing of growth line formation, duration of the growing season and growth rates) and stable isotope ratios of recent S. gigantea from AK and BC were analyzed in detail. Furthermore, a growth-temperature model based on S. gigantea shells from Alaska was established, which provides a better understanding of the hydrological changes related to the Alaska Coastal Current (ACC). This approach allows the independent measurement of water temperature and salinity from variations in the width of lunar daily growth increments of S. gigantea. Temperature explains 70% of the variability in shell growth. The model was calibrated and tested with modern shells and then applied to archaeological specimens. The time period between 988 and 1447 cal yrs BP was characterized by colder (~1-2°C) and much drier (2-5 PSU) summers, and a likely much slower flowing ACC than at present. In contrast, the summers during the time interval of 599-1014 cal yrs BP were colder (up to 3°C) and fresher (1-2 PSU) than today. The Aleutian Low may have been stronger and the ACC was probably flowing faster during this time.

Palaeo wind system reconstruction of the last glacial period over Europe, using high resolution proxy data and model-data-comparison

The atmosphere is a global influence on the movement of heat and humidity between the continents, and thus significantly affects climate variability. Information about atmospheric circulation are of major importance for the understanding of different climatic conditions. Dust deposits from maar lakes and dry maars from the Eifel Volcanic Field (Germany) are therefore used as proxy data for the reconstruction of past aeolian dynamics.rnrnIn this thesis past two sediment cores from the Eifel region are examined: the core SM3 from Lake Schalkenmehren and the core DE3 from the Dehner dry maar. Both cores contain the tephra of the Laacher See eruption, which is dated to 12,900 before present. Taken together the cores cover the last 60,000 years: SM3 the Holocene and DE3 the marine isotope stages MIS-3 and MIS-2, respectively. The frequencies of glacial dust storm events and their paleo wind direction are detected by high resolution grain size and provenance analysis of the lake sediments. Therefore two different methods are applied: geochemical measurements of the sediment using µXRF-scanning and the particle analysis method RADIUS (rapid particle analysis of digital images by ultra-high-resolution scanning of thin sections).rnIt is shown that single dust layers in the lake sediment are characterized by an increased content of aeolian transported carbonate particles. The limestone-bearing Eifel-North-South zone is the most likely source for the carbonate rich aeolian dust in the lake sediments of the Dehner dry maar. The dry maar is located on the western side of the Eifel-North-South zone. Thus, carbonate rich aeolian sediment is most likely to be transported towards the Dehner dry maar within easterly winds. A methodology is developed which limits the detection to the aeolian transported carbonate particles in the sediment, the RADIUS-carbonate module.rnrnIn summary, during the marine isotope stage MIS-3 the storm frequency and the east wind frequency are both increased in comparison to MIS-2. These results leads to the suggestion that atmospheric circulation was affected by more turbulent conditions during MIS-3 in comparison to the more stable atmospheric circulation during the full glacial conditions of MIS-2.rnThe results of the investigations of the dust records are finally evaluated in relation a study of atmospheric general circulation models for a comprehensive interpretation. Here, AGCM experiments (ECHAM3 and ECHAM4) with different prescribed SST patterns are used to develop a synoptic interpretation of long-persisting east wind conditions and of east wind storm events, which are suggested to lead to an enhanced accumulation of sediment being transported by easterly winds to the proxy site of the Dehner dry maar.rnrnThe basic observations made on the proxy record are also illustrated in the 10 m-wind vectors in the different model experiments under glacial conditions with different prescribed sea surface temperature patterns. Furthermore, the analysis of long-persisting east wind conditions in the AGCM data shows a stronger seasonality under glacial conditions: all the different experiments are characterized by an increase of the relative importance of the LEWIC during spring and summer. The different glacial experiments consistently show a shift from a long-lasting high over the Baltic Sea towards the NW, directly above the Scandinavian Ice Sheet, together with contemporary enhanced westerly circulation over the North Atlantic.rnrnThis thesis is a comprehensive analysis of atmospheric circulation patterns during the last glacial period. It has been possible to reconstruct important elements of the glacial paleo climate in Central Europe. While the proxy data from sediment cores lead to a binary signal of the wind direction changes (east versus west wind), a synoptic interpretation using atmospheric circulation models is successful. This shows a possible distribution of high and low pressure areas and thus the direction and strength of wind fields which have the capacity to transport dust. In conclusion, the combination of numerical models, to enhance understanding of processes in the climate system, with proxy data from the environmental record is the key to a comprehensive approach to paleo climatic reconstruction.rn

The reconstruction of paleoclimate variability in Northern Germany, during marine isotope stage 5

In order to obtain a better understanding about the influence of post-depositional diagenesis on speleothem 230Th/U-ages and paleoclimate variability during Marine Isotope Stage (MIS) 5 in northern Germany, four stalagmites from the Riesenberghöhle (RBH) were investigated by thin section analysis, 230Th/U-dating as well as stable oxygen and carbon isotope and laser ablation inductively coupled mass spectrometry (LA-ICPMS) trace element analysis. The RBH is located in the Weser Hills and is one of the northernmost limestone caves in Germany.rnMulti collector (MC) ICPMS 230Th/U-ages and thin section analysis of the RBH stalagmites shows that some growth phases of the stalagmites were diagenetically altered after their deposition. The impact of post-depositional diagenesis (PDD) on the 230Th/U-ages is modeled, and potential processes leading to PDD are discussed. In this context, it is suggested that PDD may be induced by rapid climate change at the inception of the GIS.rnDespite of the dating uncertainties resulting from PDD, 230Th/U-dating shows that the RBH stalagmites grew during the Eemian and most of the Greenland Interstadials (GIS) during MIS 5. Thus, the growth phases of the RBH stalagmites might be related to a reorganization of the Atlantic Meridional Overturning Circulation (AMOC). The stable isotope (δ13C and δ18O) and the trace element variability of the stalagmites reflects rapid changes of past temperature and precipitation on millennial and sub-millennial timescales. These past climate changes can be amplified by orbitally forced variations of the July solar insolation at 65°N.

Calibration study of bivalve shells : implications for paleoenvironmental reconstruction

The accretionary shells of bivalve mollusks can provide environmental information, such as water temperature, precipitation, freshwater fluxes, primary productivity and anthropogenic activities in the form of variable growth rates and variable geochemical properties, such as stable oxygen and carbon isotopes. However, paleoenvironmental reconstructions are constrained by uncertainties about isotopic equilibrium fractionation during shell formation, which is generally acknowledged as a reasonable assumption for bivalves, but it has been disputed in several species. Furthermore, the variation in shell growth rates is accepted to rely on multiple environmental variables, such as temperature, food availability and salinity, but can differ from species to species. Therefore, it is necessary to perform species-specific calibration studies for both isotope proxies and shell growth rates before they can be used with confidence for environmental interpretations of the past. Accordingly, the principal objective of this Ph.D research is to examine the reliability of selected bivalve species, the long-lived Eurhomalea exalbida (Dillwyn), the short-lived and fast growing species Paphia undulata (Born 1778), and the freshwater mussel Margaritifera falcata (Gould 1850), as paleoenvironmental proxy archives.rnThe first part is focused on δ18Oshell and shell growth history of live-collected E. exalbida from the Falkland Islands. The most remarkable finding, however, is that E. exalbida formed its shell with an offset of -0.48‰ to -1.91‰ from the expected oxygen isotopic equilibrium with the ambient water. If this remained unnoticed, paleotemperature estimates would overestimate actual water temperatures by 2.1-8.3°C. With increasing ontogenetic age, the discrepancy between measured and reconstructed temperatures increased exponentially, irrespective of the seasonally varying shell growth rates. This study clearly demonstrates that, when the disequilibrium fractionation effect is taken into account, E. exalbida can serve as a high-resolution paleoclimate archive for the southern South America. The species therefore provides quantifiable temperature estimates, which yields new insights into long-term paleoclimate dynamics for mid to high latitudes on the southern hemisphere.rnThe stable carbon isotope of biogenic carbonates is generally considered to be useful for reconstruction of seawater dissolved inorganic carbon. The δ13Cshell composition of E. exalbida was therefore, investigated in the second part of this study. This chapter focuses on inter-annual and intra-annual variations in δ13Cshell. Environmental records in δ13Cshell are found to be strongly obscured by changes in shell growth rates, even if removing the ontogenetic decreasing trend. This suggests that δ13Cshell in E. exalbida may not be useful as an environmental proxy, but a potential tool for ecological investigations. rnIn addition to long-lived bivalve species, short-lived species that secrete their shells extremely fast, can also be useful for environmental reconstructions, especially as a high-resolution recorder. Therefore, P. undulata from Daya Bay, South China Sea was utilized in Chapter 4 to evaluate and establish a potential proxy archive for past variations of the East Asian monsoon on shorter time-scales. The δ18Oshell can provide qualitative estimates of the amount of monsoonal rain and terrestrial runoff and the δ13Cshell likely reflect the relative amount of isotopically light terrestrial carbon that reaches the ocean during the summer monsoon season. Therefore, shells of P. undulata can provide serviceable proxy archives to reconstruct the frequency of exceptional summer monsoons in the past. The relative strength of monsoon-related precipitation and associated changes in ocean salinity and the δ13C ratios of the dissolved inorganic carbon signature (δ13CDIC) can be estimated from the δ18Oshell and δ13Cshell values as well as shell growth patterns. rnIn the final part, the freshwater pearl shell M. falcata from four rivers in British Columbia, Canada was preliminarily studied concerning the lifespans and the shell growth rates. Two groups separated by the Georgia Strait can be clearly distinguished. Specimens from the western group exhibit a shorter lifespan, while the eastern group live longer. Moreover, the average lifespan seems to decrease from south to north. The computed growth equations from the eastern and western groups differ as well. The western group exhibits a lower growth rate, while bivalves from the eastern group grow faster. The land use history seems to be responsible for the differences in lifespans of the specimens from the two groups. Differences in growth rate may be induced by differences in water temperature or nutrient input also related to the land use activities.

Paleoclimate reconstruction using LA-ICP-MS analyses of ostracods and stalagmite

Als paläoklimatische Archive können unter anderem Stalagmiten und Klappen von Ostrakoden herangezogen werden. Bisher gab es noch keine veröffentlichten Ergebnisse räumlich hochaufgelöster Spurenelementverteilungen in einzelnen Ostrakodenklappen. Das gleiche gilt für die Bestimmung radiogener Blei-Isotope in Stalagmiten. Um klimatische Prozesse vergangener Zeiten zu untersuchen, wurde eine neue LA-ICP-MS (Laserablations – Massenspektrometrie mit induktiv gekoppeltem Plasma) Technik für hochaufgelöste in-situ Messungen von Spurenelementen und Bleiisotopen entwickelt. rnrnZunächst wurden geeignete Materialien für die Kalibrierung der Technik untersucht; als Proben dienten die Silikatreferenzgläser BAM-S005-A und BAM-S005B der Bundesanstalt für Materialforschung und –prüfung (BAM). Die Homogenität dieser Referenzgläser wurde mit LA-ICP-MS und den anderen mikroanalytischen Methoden EPMA (Elektronenmikrosonde) und SIMS (Ionensonde) überprüft. Die Ergebnisse zeigten, dass alle Haupt- und die meisten Spurenelemente in beiden Gläsern selbst im Mikrometerbereich homogen verteilt sind. Ausnahmen sind einige Spurenelemente, wie Cs, Cl, Cr, Mo und Ni, die heterogen in den Gläsern verteilt sind. Die Hauptelementzusammensetzung von BAM-S005-A und BAM-S005-B wurde mit Hilfe der EPMA bestimmt, wobei die Ergebnisse die Referenzwerte des BAM-Zertifikats bestätigten. Mit Ausnahme von Sr, Ba, Ce und Pb, stimmten die LA-ICP-MS-Spurenelementwerte mit den zertifizierten Werten innerhalb der angegebenen Fehlergrenzen überein. Gründe für die Diskrepanz der vier oben erwähnten Elemente sind noch unklar, aber sind möglicherweise durch fehlerhafte Referenzwerte zu erklären. Zusätzlich wurden 22 Spurenelemente gemessen, deren Gehalte von BAM nicht zertifiziert wurden. Aufgrund dieser Untersuchungen konnte festgestellt werden, dass beide BAM-Gläser für mikroanalytische Anwendungen geeignet sind.rnrnUm neuartige paläoklimatische Proxies im Calcit von Ostrakoden aus tibetanischen Seesedimenten zu untersuchen, wurde die Spurenelementvariabilität in einzelnen Ostrakodenklappen durch eine neue LA-ICP-MS-Technik bestimmt. Klappen von Ostrakoden der drei Arten (Leucocytherella sinensis Huang, 1982, ?Leucocythere dorsotuberosa Huang, 1982 und ?L. dorsotuberosa f. postilirata sensu Pang, 1985) wurden aus zwei Sedimentkernen des Nam Co Sees auf dem Hochplateau von Tibet gewonnen. Zwei LA-ICP-MS-Varianten, Spot- bzw. Linienanalyse, wurden verwendet, um die Elementkonzentrationen der Spurenelemente Mg, Sr, Ba, U und die der Seltenen Erdelemente (SEE) in den Klappen einzelner Ostrakoden zu bestimmen. Die Ergebnisse zeigten, dass die Linienanalyse präzisere Daten als die Spotanalyse liefert und sie wurde daher vorgezogen. Signifikante Unterschiede in der Spurenelementzusammensetzung zwischen den verschiedenen Arten der Ostrakoden wurden nicht gefunden. Variationen der Elementverhältnisse Mg/Ca und Sr/Ca in den Klappen stimmen mit veröffentlichten Seenspiegelschwankungen während des Holozäns überein, was zeigt, dass Mg- und Sr-Messungen in den Ostrakoden zur Untersuchung paläohydrochemischer Prozesse in diesem Gebiet herangezogen werden kann. Die gute Korrelation, die in dieser Arbeit zwischen Ba/Ca und Sr/Ca gefunden wurde, ist ein Hinweis darauf, dass der Einbau von Ba und Sr in die Klappen durch den gleichen Mechanismus erfolgte. Eine mögliche Beziehung zwischen dem U/Ca-Verhältnis in den Ostrakoden und den Redoxbedingungen auf dem Boden des Sees in der Vergangenheit wird diskutiert. Relativ geringe und konstante La/Ca-Verhältnisse wurden festgestellt, deren Ursache möglicherweise entweder auf der SEE-Charakteristik des Seewassers, auf biologischen Prozessen in den Ostrakoden oder auf Kontamination von Fe-Mn und/oder organischen Substanzen beruhen. Weitere Untersuchungen an Proben aus diesem Gebiet, speziell Klappen von lebenden Ostrakoden, sind notwendig, um den Gehalt von Ba, U und den SEE in Ostrakoden als paläoklimatische Proxies von Umweltbedingungen zu verwenden.