Search

We describe the compilation of two early instrumental daily temperature series from Bern and Zurich, Switzerland, starting from 1760 and 1756, respectively. The series are a combination of numerous small segments from different observers at different locations within and outside the two cities that are converted to modern units and homogenized. In addition, we introduce a methodology to estimate the errors affecting daily and monthly mean values derived from early instrumental observations. Given the frequent small data gaps, we merge the two daily series into a more complete series representing the central Swiss Plateau. We finally compare the homogenized monthly series with other temperature reconstructions for Switzerland. We find significant differences before 1860, pointing to biases that might affect some of the most widely used instrumental data sets. In general, the homogenization of temperature measurements at the transition between the early instrumental and national weather service eras remains a problematic issue in historical climatology and has significant implications for other fields of climate research.

European flood frequency and intensity change on a multidecadal scale. Floods were more frequent in the 19th (central Europe) and early 20th century (western Europe) than during the mid-20th century and again more frequent since the 1970s. The causes of this variability are not well understood and the relation to climate change is unclear. Palaeoclimate studies from the northern Alps suggest that past flood-rich periods coincided with cold periods. In contrast, some studies suggest that more floods might occur in a future, warming world. Here we address the contribution of atmospheric circulation and of warming to multidecadal flood variability. For this, we use long series of annual peak streamflow, daily weather data, reanalyses, and reconstructions. We show that both changes in atmospheric circulation and moisture content affected multidecadal changes of annual peak streamflow in central and western Europe over the past two centuries. We find that during the 19th and early 20th century, atmospheric circulation changes led to high peak values of moisture flux convergence. The circulation was more conducive to strong and long-lasting precipitation events than in the mid-20th century. These changes are also partly reflected in the seasonal mean circulation and reproduced in atmospheric model simulations, pointing to a possible role of oceanic variability. For the period after 1980, increasing moisture content in a warming atmosphere led to extremely high moisture flux convergence. Thus, the main atmospheric driver of flood variability changed from atmospheric circulation variability to water vapour increase.

Data assimilation approaches such as the ensemble Kalman filter method have become an important technique for paleoclimatological reconstructions and reanalysis. Different sources of information, from proxy records and documentary data to instrumental measurements, were assimilated in previous studies to reconstruct past climate fields. However, precipitation reconstructions are often based on indirect sources (e.g., proxy records). Assimilating precipitation measurements is a challenging task because they have high uncertainties, often represent only a small region, and generally do not follow a Gaussian distribution. In this paper, experiments are conducted to test the possibility of using information about precipitation in climate reconstruction with monthly resolution by assimilating monthly instrumental precipitation amounts or the number of wet days per month, solely or in addition to other climate variables such as temperature and sea-level pressure, into an ensemble of climate model simulations. The skill of all variables (temperature, precipitation, sea-level pressure) improved over the pure model simulations when only monthly precipitation amounts were assimilated. Assimilating the number of wet days resulted in similar or better skill compared to assimilating the precipitation amount. The experiments with different types of instrumental observations being assimilated indicate that precipitation data can be useful, particularly if no other variable is available from a given region. Overall the experiments show promising results because with the assimilation of precipitation information a new data source can be exploited for climate reconstructions. The wet day records can become an especially important data source in future climate reconstructions because many existing records date several centuries back in time and are not limited by the availability of meteorological instruments.

Differences between paleoclimatic reconstructions are caused by two factors: the method and the input data. While many studies compare methods, we will focus in this study on the consequences of the input data choice in a state-of-the-art Kalman-filter paleoclimate data assimilation approach. We evaluate reconstruction quality in the 20th century based on three collections of tree-ring records: (1) 54 of the best temperature-sensitive tree-ring chronologies chosen by experts; (2) 415 temperature-sensitive tree-ring records chosen less strictly by regional working groups and statistical screening; (3) 2287 tree-ring series that are not screened for climate sensitivity. The three data sets cover the range from small sample size, small spatial coverage and strict screening for temperature sensitivity to large sample size and spatial coverage but no screening. Additionally, we explore a combination of these data sets plus screening methods to improve the reconstruction quality. A large, unscreened collection generally leads to a poor reconstruction skill. A small expert selection of extratropical Northern Hemisphere records allows for a skillful high-latitude temperature reconstruction but cannot be expected to provide information for other regions and other variables. We achieve the best reconstruction skill across all variables and regions by combining all available input data but rejecting records with insignificant climatic information (p value of regression model >0.05) and removing duplicate records. It is important to use a tree-ring proxy system model that includes both major growth limitations, temperature and moisture.

Data assimilation has been adapted in paleoclimatology to reconstruct past climate states. A key component of some assimilation systems is the background-error covariance matrix, which controls how the information from observations spreads into the model space. In ensemble-based approaches, the background-error covariance matrix can be estimated from the ensemble. Due to the usually limited ensemble size, the background-error covariance matrix is subject to the so-called sampling error. We test different methods to reduce the effect of sampling error in a published paleoclimate data assimilation setup. For this purpose, we conduct a set of experiments, where we assimilate early instrumental data and proxy records stored in trees, to investigate the effect of (1) the applied localization function and localization length scale; (2) multiplicative and additive inflation techniques; (3) temporal localization of monthly data, which applies if several time steps are estimated together in the same assimilation window. We find that the estimation of the background-error covariance matrix can be improved by additive inflation where the background-error covariance matrix is not only calculated from the sample covariance but blended with a climatological covariance matrix. Implementing a temporal localization for monthly resolved data also led to a better reconstruction.

For an assessment of the roles of soil and vegetation in the climate system, a further understanding of the flux components of H2O and CO2 (e.g., transpiration, soil respiration) and their interaction with physical conditions and physiological functioning of plants and ecosystems is necessary. To obtain magnitudes of these flux components, we applied source partitioning approaches after Scanlon and Kustas (2010; SK10) and after Thomas et al. (2008; TH08) to high-frequency eddy covariance measurements of 12 study sites covering different ecosystems (croplands, grasslands, and forests) in different climatic regions. Both partitioning methods are based on higher-order statistics of the H2O and CO2 fluctuations, but proceed differently to estimate transpiration, evaporation, net primary production, and soil respiration. We compared and evaluated the partitioning results obtained with SK10 and TH08, including slight modifications of both approaches. Further, we analyzed the interrelations among the performance of the partitioning methods, turbulence characteristics, and site characteristics (such as plant cover type, canopy height, canopy density, and measurement height). We were able to identify characteristics of a data set that are prerequisites for adequate performance of the partitioning methods. SK10 had the tendency to overestimate and TH08 to underestimate soil flux components. For both methods, the partitioning of CO2 fluxes was less robust than for H2O fluxes. Results derived with SK10 showed relatively large dependencies on estimated water use efficiency (WUE) at the leaf level, which is a required input. Measurements of outgoing longwave radiation used for the estimation of foliage temperature (used in WUE) could slightly increase the quality of the partitioning results. A modification of the TH08 approach, by applying a cluster analysis for the conditional sampling of respiration–evaporation events, performed satisfactorily, but did not result in significant advantages compared to the original method versions developed by Thomas et al. (2008). The performance of each partitioning approach was dependent on meteorological conditions, plant development, canopy height, canopy density, and measurement height. Foremost, the performance of SK10 correlated negatively with the ratio between measurement height and canopy height. The performance of TH08 was more dependent on canopy height and leaf area index. In general, all site characteristics that increase dissimilarities between scalars appeared to enhance partitioning performance for SK10 and TH08.

The Modern Era Reanalysis (ModE-RA) is a global monthly paleo-reanalysis covering the period between 1421 and 2008. To reconstruct past climate fields an offline data assimilation approach is used, blending together information from an ensemble of transient atmospheric model simulations and observations. In the early period, ModE-RA utilizes natural proxies and documentary data, while from the 17 th century onward instrumental measurements are also assimilated. The impact of each observation on the reconstruction is stored in the observation feedback archive, which provides additional information on the input data such as preprocessing steps and the regression-based forward models. The monthly resolved reconstructions include estimates of the most important climate fields. Furthermore, we provide a reconstruction, ModE-RAclim, which together with ModE-RA and the model simulations allows to disentangle the role of observations and model forcings. ModE-RA is best suited to study intra-annual to multi-decadal climate variability and to analyze the causes and mechanisms of past extreme climate events., (© 2024. The Author(s).)

Island nations may have potential long-term survival value for humanity in global catastrophes such as sun-blocking catastrophes from nuclear winter and large magnitude volcanic eruptions. One way to explore this issue further is to understand the impact on islands after the largest historically observed volcanic eruption: that of Mt Tambora in 1815. For each of the 31 large, populated islands selected, we conducted literature searches for relevant historical and palaeoclimate studies. We also analysed results from a reconstruction (EKF400v2), which uses atmospheric-only general circulation model simulations with assimilated observational and proxy data. From the literature review, there was widespread evidence for weather/climate anomalies in 1815-1817 for these islands (29/29 for those with data). But missing data was an issue for other dimensions such as impaired food production (seen in 8 islands out of only 12 with data). Based on the EKF400v2 reconstruction for temperature anomalies (compared to the relatively "non-volcanic" reference period of 1779 to 1808), the islands had lower temperature anomalies in the 1815-1818 period than latitudinally equivalent continental sites (at 100 km and 1000 km inland). This was statistically significant for the great majority of the comparisons for group analyses by hemisphere, oceans, and temperate/tropical zone. When considering just the islands, all but four showed statistically anomalous temperature reductions in the 1816-1817 period (for most p < 0.00001). In the peak impact year of 1816, the lowest anomalies were seen for islands in the Southern Hemisphere (p < 0.0001), the Indian Ocean (p < 0.0001), and in the tropics and subtropics of the Southern Hemisphere (p = 0.0057). In conclusion, the findings of both the literature review and reconstruction simulations suggest climatic impacts of the Tambora eruption for nearly all these 31 large islands, albeit less than for continental sites. Islands with the smallest temperature anomalies were in the Southern Hemisphere, in particular the Indian Ocean and the tropics and subtropics of the Southern Hemisphere., (© 2023. The Author(s).)

Data assimilation techniques are becoming increasingly popular for climate reconstruction. They benefit from estimating past climate states from both observation information and from model simulations. The first monthly global paleo-reanalysis (EKF400) was generated over the 1600 and 2005 time period, and it provides estimates of several atmospheric fields. Here we present a new, considerably improved version of EKF400 (EKF400v2). EKF400v2 uses atmospheric-only general circulation model simulations with a greatly extended observational network of early instrumental temperature and pressure data, documentary evidences and tree-ring width and density proxy records. Furthermore, new observation types such as monthly precipitation amounts, number of wet days and coral proxy records were also included in the assimilation. In the version 2 system, the assimilation process has undergone methodological improvements such as the background-error covariance matrix is estimated with a blending technique of a time-dependent and a climatological covariance matrices. In general, the applied modifications resulted in enhanced reconstruction skill compared to version 1, especially in precipitation, sea-level pressure and other variables beside the mostly assimilated temperature data, which already had high quality in the previous version. Additionally, two case studies are presented to demonstrate the applicability of EKF400v2 to analyse past climate variations and extreme events, as well as to investigate large-scale climate dynamics., Competing Interests: The authors declare that they have no conflict of interest., (© 2021 The Authors. Geoscience Data Journal published by Royal Meteorological Society and John Wiley & Sons Ltd.)

Annual-to-decadal variability in northern midlatitude temperature is dominated by the cold season. However, climate field reconstructions are often based on tree rings that represent the growing season. Here we present cold-season (October-to-May average) temperature field reconstructions for the northern midlatitudes, 1701-1905, based on extensive phenological data (freezing and thawing dates of rivers, plant observations). Northern midlatitude land temperatures exceeded the variability range of the 18th and 19th centuries by the 1940s, to which recent warming has added another 1.5 °C. A sequences of cold winters 1808/9-1815/6 can be explained by two volcanic eruptions and unusual atmospheric flow. Weak southwesterlies over Western Europe in early winter caused low Eurasian temperatures, which persisted into spring even though the flow pattern did not. Twentieth century data and model simulations confirm this persistence and point to increased snow cover as a cause, consistent with sparse information on Eurasian snow in the early 19th century., (© 2022. The Author(s).)

This paper describes a global monthly gridded Sea Surface Temperature (SST) and Sea Ice Concentration (SIC) dataset for the period 1000-1849, which can be used as boundary conditions for atmospheric model simulations. The reconstruction is based on existing coarse-resolution annual temperature ensemble reconstructions, which are then augmented with intra-annual and sub-grid scale variability. The intra-annual component of HadISST.2.0 and oceanic indices estimated from the reconstructed annual mean are used to develop grid-based linear regressions in a monthly stratified approach. Similarly, we reconstruct SIC using analog resampling of HadISST.2.0 SIC (1941-2000), for both hemispheres. Analogs are pooled in four seasons, comprising of 3-months each. The best analogs are selected based on the correlation between each member of the reconstructed SST and its target. For the period 1780 to 1849, We assimilate historical observations of SST and night-time marine air temperature from the ICOADS dataset into our reconstruction using an offline Ensemble Kalman Filter approach. The resulting dataset is physically consistent with information from models, proxies, and observations., (© 2021. The Author(s).)