2000 years of annual ice core data from Law Dome, East Antarctica
Ice core records from Law Dome in East Antarctica collected over the last four decades provide high-resolution data for studies of the climate of Antarctica, Australia, and the Southern and Indo-Pacific oceans. Here, we present a set of annually dated records of trace chemistry, stable water isotopes and snow accumulation from Law Dome covering the period from −11 to 2017 CE (1961 to −66 BP 1950) and the level-1 chemistry data from which the annual chemistry records are derived. Law Dome ice core records have been used extensively in studies of the past climate of the Southern Hemisphere and in large-scale data syntheses and reconstructions in a region where few records exist, especially at high temporal resolution. This dataset provides an update and extensions both forward and back in time of previously published subsets of the data, bringing them together into a coherent set with improved dating to enable continued use of this record. The data are available for download from the Australian Antarctic Data Centre at https://doi.org/10.26179/5zm0-v192 (Curran et al., 2021).
Law Dome is a small ice cap on the coast of East Antarctica (Fig. 1), uniquely positioned to provide exceptionally high temporal resolution ice core records which preserve climate signals from a large sector of the Southern Hemisphere where there are limited historical and instrumental records available. It is situated on a promontory along the Budd Coast of East Antarctica which is bounded by the Totten and Vanderford glaciers, separating it from the main ice flow in the interior of the Antarctic ice sheet. Its climate is driven largely by its position in the path of the large cyclonic weather systems originating from the Southern Ocean (Bromwich, 1988; Udy et al., 2021), which in turn provide a link to lower latitude climate. These cyclonic weather systems result in relatively high annual snowfall with orographic influences resulting in a snowfall gradient across the dome from east to west with higher annual snowfall to the eastern side (Pfitzner and Australian National Antarctic Research Expeditions, 1980). The high snowfall accumulation rate at the summit is predominantly uniform throughout the year (van Ommen and Morgan, 1997) with well preserved annual layers, allowing for ice core records of annual to seasonal resolution and clear dating which have provided high-resolution proxy data of the Southern Hemisphere climate over the past 2000 years. In particular, it provides high-resolution records of the Pacific and Indian sectors of the Southern Ocean where there are limited historical or instrumental observations.
Law Dome's close proximity to Casey Station has made it a relatively easy site to access for ice coring activities and repeat site visits to update records with shallow cores. Ice cores have been collected from various sites on Law Dome and have been shown to record regional, hemispheric and global records of past climate over the common era. The ice core from the Dome Summit South (DSS) site provides a long record extending back to the last glacial maximum but with very high temporal resolution over the past 2000 years. Measurements in cores from elsewhere on Law Dome have included gases trapped within ice to investigate past atmospheric composition (Etheridge et al., 1998; Rubino et al., 2019). Stable isotopes of water from the DSS site show a regular seasonal cycle, helping to identify annual layers and enabling the estimation of the annual snow accumulation (Roberts et al., 2015) and have wider uses in southern Indian Ocean variability (Masson-Delmotte et al., 2003) and in large compilations of past global temperature variability reconstruction (e.g. Ahmed et al., 2013; Emile-Geay et al., 2017). Proxies derived from annually dated ice core records from Law Dome have been used to reconstruct climate modes such as El Niño–Southern Oscillation (ENSO) and the Interdecadal Pacific Oscillation (IPO) (Vance et al., 2015), have been used as records of Australian rainfall (van Ommen and Morgan, 2010; Vance et al., 2013) and streamflow (Tozer et al., 2018), and for evaluation of solar and meteorological influences on the 10Be solar activity proxy (Pedro et al., 2011a, b, 2012), volcanic activity (Plummer et al., 2012) and sea ice extent (Curran et al., 2003).
Palaeoclimate records from East Antarctica provide important longer term context of trends and variability for the shorter instrumental records that are available. The Law Dome record is significant as it is one of the few continually updated sources of palaeoclimate data in this region. Publishing the latest, version controlled record to ensure the most accurate data are available and citable is important for input into large community climate compilation and synthesis efforts such as the Antarctic Climate Change over the 21st Century (AntClim21) and Antarctic Climate Change over the 21st Century (AntClimNow) Scientific Committee on Antarctic Research (SCAR) initiatives and projects under the Past Global Changes (PAGES) 2k network. Examples of projects with particular reliance on ice core records include reconstructions of climate variability over the last 2000 years by Stenni et al. (2017), synthesis and reconstruction of snow accumulation patterns across Antarctica (Thomas et al., 2017) and global temperature trends (Emile-Geay et al., 2017).
This collection of records supersedes all previously released versions of data from the DSS site on Law Dome, with the composite record described here forming the best continuous record available with improved subannual dating. Repeated visits to the DSS site has allowed the record to be continuously updated, extending the overlap with the instrumental data period. Further analysis of previously collected ice core material at finer resolution than previously performed has enabled improvements to the subannual dating.
We provide trace chemistry, stable water isotopes and snow accumulation data with some higher-level products also available. Level-1 data consist of quality-controlled trace chemistry data for chloride (Cl−), sodium (Na+), magnesium (Mg2+), calcium (Ca2+), potassium (K+), sulfate () and nitrate () ions. At this time, level-1 stable water isotope data are not included. Level-2 datasets are derived from the level-1 products with annual averages provided for each chemical species, δ18O and snow accumulation and seasonal sea salts. These derived datasets serve as examples of what we believe to be the most suitable method for annualising each data stream. We also detail the limitations which should be considered in their application in further research.
Further visits to this site for collecting ice cores is expected as part of an ongoing monitoring project with subsequent updates to this record when ice cores become available for analysis. These regular visits to update the record by collecting shallow cores extend the period which overlaps with modern satellite observations, ensuring the record stays up to date and more data are available for calibration of palaeoclimate reconstructions.
This publication provides a reference point for the full suite of data with links to versioned and updated data publicly available for download from the Australian Antarctic Data Centre at https://doi.org/10.26179/5zm0-v192 (Curran et al., 2021). Note that in the text of this article, we use years CE for dates as the more natural dating scale and to be consistent with the naming of some cores using the year or season in which they were drilled. The datasets available for download provide years in both CE and BP 1950 scales and we use BP 1950 (i.e. years before 1950 CE) for plotting of data in figures here.
2.1 Drilling campaigns
The DSS site is located at 66∘46′11′′ S, 112∘48′25′′ E, approximately 4.7 km SSW from the dome summit (Morgan et al., 1997). The drill sites of the four cores included in this composite record are all located within approximately 1 km of each other with their relative locations shown in Fig. 1b. Future visits to the site are planned to ensure the record continues to be updated. A summary of the details of the drilling campaigns is found in Table 1 with a more detailed description of the original drilling site published by Morgan et al. (1997).
2.2 Dating and age horizon uncertainties
The DSS record has been dated using seasonal variations in the deposition of species to define calendar year boundaries, commonly known as annual layer counting. The stable water isotope record is the primary seasonal indicator with a well-defined maximum of January 10±2.2 d (van Ommen and Morgan, 1997). Confirmation of dating is provided by summer-peaking hydrogen peroxide where available and other chemical markers such as summer peaking methanesulfonic acid (MSA), non-sea salt sulfate (), sulfate/chloride () ratio, and the winter-peaking sea salt species (chloride (Cl−), sodium (Na+) and magnesium (Mg2+)). Layer counting of the DSS record between 1300 and 1995 CE was previously completed by Palmer et al. (2001b) and verified and extended to cover over 2000 years by Plummer et al. (2012). The DSS record currently spans −11 to 2017 CE with the addition of the DSS1617 core which is dated unambiguously from 1989 to 2017 CE. The timescale from Plummer et al. (2012) has been updated by refining the annual layer placement, assisted by new water isotope analysis and by applying a stricter error counting protocol. The number of years counted has not changed, however, the error estimate has increased in the oldest part of the record to +20/−7 at −11 CE where this date may be up to 20 years older or 7 years younger (−31 to −4 CE) as shown in Fig. 2. Where evidence for a year horizon was not clearly identifiable in the primary seasonal indicators or only weakly supported by confirmatory species, it was not counted. When the majority of (but not all) seasonal indicators show evidence of a year, these were counted. Both cases contribute to the total uncertainty estimate. The asymmetrical uncertainty is biased toward under-counting, reflecting increasing likelihood of encountering years with unclear signals with the decreasing layer thickness at depth. DSS is a long record from a high accumulation site (0.69 m yr−1 ice equivalent using an ice density of 917 kg m−3) and maintains six to eight chemistry samples per year at −11 CE, enabling development of a long, accurate annual layer-counted chronology. Intra-annual dating of individual sample depths is determined by interpolation between the layer counted annual horizon depths.
We take a conservative approach to calculating the age horizons and uncertainties. The depth of the top of each 1 m length of drilled core is locked to the top of the bag as recorded in the drilling logs so uncertainty from the sample resolution is limited only within each 1 m length and errors in depth arising from missing or badly fragmented core segments are not cumulative. The sample resolution of chemical measurements is generally much greater than the water isotope sample resolution, hence, the larger uncertainty value is used.
The DSS chronology presented here has been dated using annual layer counting only, without reference to externally dated reference events (e.g. volcanic events) and is independent of other ice core age scales. There are differences when compared to other ice core chronologies, however, those differences are within our error bounds. To assist users to include DSS data in large-scale reconstructions or regional synthesis products, we have provided a translation of the age of the volcanic events used in the PAGES Ant2k compilation (Ahmed et al., 2013) to the WD2014 (Sigl et al., 2016) chronology and our DSS age scale in Appendix A.
The depths of the year boundaries for the 2000 years of data, produced using annual layer counting methods as described above, are contained in the file titled “DSS_2k_age_horizons.csv”. The column headings for this file are described in Table 2.
The level-1 datasets have undergone calibration, quality control and post-processing of the raw instrument measurements. The datasets here are obtained from four cores (see Table 1), drilled at the Dome Summit South (DSS 66∘46′ S, 112∘48′ E) at different times. The main 1196 m (DSSMain) core was drilled between 1987 and 1993. The uppermost 117 m of DSSMain was thermally drilled and the presence of micro-fractures made it unsuitable for trace ion analysis. Two further shallow cores – DSS97 and DSS99 – were drilled to cover this period for improved data quality. The stable water isotope records from all three cores were used to establish and lock unambiguously the dating and overlap between them (Palmer et al., 2001b). The DSS site has been revisited subsequently with new short cores retrieved to update the DSS record. The most recent update to the DSS record was DSS1617, a 30 m core drilled during the 2016/17 CE austral summer season, covering 1989–2017 CE, providing 7 years of overlap with DSS97. Previous work from DSS has used other composite records using short cores overlapping with DSS97 (Plummer et al., 2012; Vance et al., 2013; Roberts et al., 2015), which are now no longer used as the periods they covered have been superseded by the longer DSS1617. The periods covered by each of these four cores are illustrated in Fig. 3. Where we have a choice of cores, we have typically selected records to avoid the uppermost annual cycle of any core (where the friability of surface snow and firn can result in core sections that are more easily compromised during core processing) and to use long, continuous sections from single cores where possible.
The stable water isotope (δ18O) record is sampled at 10–50 mm resolution with finer physical sampling as annual layers thin with depth. Level-2 annual mean data are provided (see Sect. 4.1), however, the level-1 datasets of stable water isotopes are not provided at this time but will be published in the future with a manuscript currently in preparation. To assist with dating, sections of cores were measured for hydrogen peroxide at an average resolution of 50 mm as described by van Ommen and Morgan (1996). Discrete chemistry samples were prepared using the clean techniques described by Curran and Palmer (2001). Cores DSS1617, DSS99, DSS97 and DSSMain to 402 m (∼1300 CE) were sampled at an average resolution of 50 mm, providing up to 60 samples per year in the near-surface firn. Due to a sampling error, a section of DSSMain from 251–273 m (1611–1568 CE) was sampled for chemistry at 100 mm resolution, however, isotope and peroxide measurements are available at 50 mm resolution to maintain dating integrity. Beyond 402 m, sample resolution changed to 30 mm and below 578 m (979 CE), resolution changed to 25 mm to offset the effects of layer thinning on seasonality.
3.1 Trace ion chemistry
Ice core samples have been analysed using ion chromatography for chloride (Cl−), sodium (Na+), magnesium (Mg2+), calcium (Ca2+), potassium (K+), sulfate (), nitrate () and methanesulfonic acid (MSA−). Analytical methods have been modified and updated over the course of the analysis and are summarised as follows (see also Table 3). DSSMain from 117–402 m was analysed using the methods of Buck et al. (1992) at the University of New Hampshire. DSS99, DSS97 and DSSMain from 402–748 m were analysed according to the methods of Curran and Palmer (2001), with a deeper section from 748–794 m using the methods of Plummer et al. (2012). The most recent analysis, DSS1617, was analysed according to Plummer et al. (2012) with the exception of the cation analytical column being changed to an Ionpac CS19 for improved detection and peak resolution of magnesium and calcium. Comparison of the analysis methods of Curran and Palmer (2001) and Buck et al. (1992) were discussed by Palmer et al. (2001b) and showed no significant difference. A comparison of the technique used by Plummer et al. (2012) and Curran and Palmer (2001) similarly found no significant differences.Plummer et al. (2012)Curran and Palmer (2001)Curran and Palmer (2001)Buck et al. (1992)Buck et al. (1992)Curran and Palmer (2001)Curran and Palmer (2001)Plummer et al. (2012)
∗ A sampling error resulted in the larger sample size in this section.
A number of species that were measured on parts of the core are not included in this dataset. MSA was not measured on a significant portion of the record. Additionally, deeper sections of the MSA record suffer from unquantified losses from storage prior to analysis (Roberts et al., 2009). The potassium and calcium records are incomplete and suffer issues with poor detection in some older analyses. Additionally, some sections of the DSSMain calcium record appear to have been affected by dust from storage in a concrete-floored freezer. Due to these quality concerns, the potassium and calcium records are not discussed further.
For the remaining datasets, we provide for each chemical species the concentration for each sample used in the composite record, an ID corresponding to the drilled core it was obtained from, the top and bottom depths of each sample and the corresponding age. Summary statistics for the level-1 trace chemistry species are included in Table 4 with the corresponding histograms shown for each species shown in Fig. 4. These show that the distributions of these species are generally non-normal. In further analysis (such as shown later in this paper in Sect. 4.3), the concentrations are log-transformed to partially compensate for the long-tail of the distributions.
The level-1 trace chemistry dataset provides a depth and age for the top and bottom of each sample and concentrations of the measured ions. These are contained in the file titled “DSS_2k_chemistry_level1.csv” with column headings provided for reference here in Table 5.
Further statistical analysis was performed to examine the consistency in the chemistry data obtained from the different cores within the composite even where there are little or no overlapping data. As the concentrations of trace ions in the ice may depend on the accumulation rate, we have analysed the cores for epochs of approximately constant accumulation rate and only compared sections which have similar (but not identical) accumulation rates. The detailed results of this analysis show that the distributions of the chemical species concentrations in the different cores are similar in epochs of comparable snow accumulation rate. The results and details of the analysis performed are found in Appendix B.
The firn densification in this region of Law Dome was established from density measurements taken from several cores in the vicinity of DSS. Samples were prepared from dry drilled cores or from melt-free interior of thermally drilled cores and were either machined to cylinders on a lathe or pressed (cookie-cutter) from softer firn to a precise volume.
An empirical fit to the firn density profile with depth was presented in van Ommen et al. (1999). Specifically, the density (ρ in kg m−3) is fitted as a compound function of depth (z, in m) as
The density observations and empirical fit as a function of depth (see Eq. 1) are plotted in Fig. 5, showing good agreement between the two.
Ice equivalent depths
For some derived data products (e.g. snow accumulation history), it is useful to be able to remove the effects of firn densification. One common method is to utilise “ice equivalent depths” that represent the depth of a column numerically compressed to the same density as glacial ice (density ρice=917 kg m−3) and the same cross-sectional area with the same mass as the firn column. Algebraically, the ice equivalent depth (𝒵) is given by
The level-2 datasets are derived from the level-1 datasets above. Those presented here have largely been previously published or included in large data compilations but are now updated with improved dating of the core. Where there are differences between the versions, these are remarked upon here. These level-2-derived datasets can be considered as examples of best practise for utilising the level-1 data. A plot of the time series for all level-2 datasets is shown in Fig. 6.
Two CSV format text files are provided for the level-2 datasets, “DSS_2k_winter_centred.csv” and “DSS_2k_summer_centred.csv”, which contain the annually averaged datasets as described in Sect. 4. Column headings are provided for reference in Tables 6 and 7. These datasets are provided together to ensure a consistent age scale across all the derived datasets.
4.1 Annual stable water isotopes
The DSS annual stable water isotope presented here updates and extends the stable isotope record included in Emile-Geay et al. (2017) and Stenni et al. (2017), using the DSS1617, DSS99, DSS97 and DSSMain cores (see Fig. 3) and currently spans the period −11 to 2016 CE. The previously published annual composite record for DSS, used in Antarctic temperature reconstructions (Ahmed et al., 2013; Emile-Geay et al., 2017; Stenni et al., 2017), included other short cores which are now no longer used as the time periods they covered and are now superseded by the longer DSS1617. Also, further analysis has extended the annual stable water isotope for the DSSMain record to cover −11 to 174 CE. The record presented here includes some changes from the previous records due to the new cores, minor changes in the year boundaries and the new analysis data. Annual averages are calculated using the year boundaries defined by the summer peak in δ18O. Changes from previously released records have occurred due to new data being obtained from the DSS1617 core and new analysis at depth from the original DSSMain core. Minor adjustments have been made to the year/depth horizons and corrections to core flips, sample mishandling and depth errors which have been identified since the last release and have been verified by the chemistry records. The time series for the annual-averaged δ18O data is shown in Fig. 7.
4.2 Annual trace ions
Annual average trace ion concentration values are presented here and the calculation method is species-dependent. A winter-centred average is used for winter-peaking species (e.g. sea salts) with the boundaries set at the beginning and end of each calendar year. The summer-centred value is calculated from mid-year to mid-year and used for the nominally summer-peaking species nitrate and non-sea-salt sulfate. This was done to reduce edge effects where small differences in year boundary placement could have strong effects on summer-peaking species. The non-sea-salt sulfate concentration is calculated according to the method of Palmer et al. (2001a) and has a fractionation correction applied to minimise non-physical negative values. The distributions of the annually averaged records are shown in Fig. 8.
Plots of the annual average time series for each of the measured chemistry species are given below along with the median value, 2σ threshold and the number of samples used for each annual average shown in the lower panel (see Figs. 9–14).
4.3 Seasonal sea salts
Characteristics of seasonal cycles of the different species of ions from the Law Dome ice core have been studied over short time periods (Curran et al., 1998), while annual sea salt data have been used previously for high-precision dating of volcanic events (Palmer et al., 2001b) and used for palaeoclimate studies of atmospheric circulation (Souney et al., 2002).
Three seasonal sea salt records were developed from aggregating the level-1 sodium record for climate proxy and reconstruction purposes in Vance et al. (2013, 2015) and further used by Crockart et al. (2021); Vance et al. (2021). Specifically, these three records include a prescribed summer season of December to March inclusive mean chloride concentrations, hereafter the “summer sea salt record” and a semi-annual aggregation into a defined warm (December to June inclusive) and cool (July to November inclusive) mean salt concentration. The summer sea salt record was developed as a proxy for rainfall variability in eastern Australia and a western tropical Pacific ENSO proxy over the Common Era, and is provided as log-transformed mean concentrations. Variability in summer sea salt concentrations is attributed to synoptic scale changes in the southern Indian Ocean (Udy et al., 2021). The warm and cool season mean sea salt concentrations were similarly log-transformed and provided input time series for reconstructing Pacific Ocean decadal variability, specifically, the IPO index (Parker et al., 2007; Vance et al., 2015, 2016). These datasets have been extended and undergone improved quality control procedures, as detailed in Vance et al. (2021) and adopted here, which extended the IPO reconstruction to span the Common Era.
These seasonal records were derived by binning the quality controlled, dated sea salt data from the level-1 datasets into 12 months per year. The salt concentrations were log-transformed to normalise (as sea salt data from Law Dome display long tails due to high frequency aerosol generation events from synoptic activity in the southern Indian Ocean). After binning, time series of seasonal averages (December–March, December–May and June–November) were produced with the following caveats. The level-1 sea salt records were examined in conjunction with the stable water isotope records and field and laboratory logbooks detailing core morphology and ice sample cutting. In some instances, the binned sea salt data were compromised due to either missing ice core material from core breaks or shattered sections, or analytical problems resulting in no data or suspected contamination. If this had occurred, the data in question were removed from the analysis. Where more than 1 month of data were compromised for the summer sea salt record or 2 months for the warm and cool season records, we elected to not include that season in our published time series. This results in missing values in the seasonal sea salt time series of summer, warm and cool season salt concentrations over the last 2000 years and are recorded as NaN in the dataset files. Time series for the summer, warm and cool season log-transformed sea salts are shown in Figs. 15, 16 and 17.
4.4 Annual snow accumulation
Snow accumulation is derived using the year horizons obtained through the layer-counted dating. This record has been previously published in Roberts et al. (2015) and included in PAGES 2k snow accumulation data compilation (Thomas et al., 2017). This update includes the newer DSS1617 core data and the improved dating. We assume steady state for depth profiles of both density (Sorge's law) and the vertical strain rate. To compensate for firn densification effects, the depths for year boundaries have been converted to ice equivalent depths using Eq. (2). Annual snow accumulation rates were then estimated using the power-law vertical strain rate method of Roberts et al. (2015), based on the assumption of no long-term trend in snow accumulation rate. The resulting snow accumulation rate is shown in Fig. 18 with the long-term mean snow accumulation rate of 0.691±0.004 IE m yr−1.
The dataset described in this paper can be accessed at the Australian Antarctic Data Centre under https://doi.org/10.26179/5zm0-v192 (Curran et al., 2021). Individual files may be accessed by clicking on “View Dataset Contents”.
In this paper, we have presented the suite of data covering the last 2000 years from the Law Dome ice cores. This set has been updated from previously published results with improved and extended annual dating. We provide quality-controlled trace ion chemistry datasets and derived, annualised products.
This dataset can be used for high-resolution studies of Southern Hemisphere climate drivers at annual to seasonal scales.
Our aim for this paper has been to provide the community with the most complete and consistent suite of Law Dome for the past 2000 years and examples of best-practise methods for generating derived datasets such as the annual means. The datasets here supersede all previous releases. This dataset is useful for future high-resolution studies of palaeoclimate proxies as has been useful in the past. Future users of the data should take the following summary of the points into consideration:
The dating presented has uncertainties due to occasionally unclear seasonal markers. The dating is biased towards undercounting of years with unclear boundaries.
As this dataset is dated through annual layer counting, independent of reference to volcanic or solar activity ties, there is an offset in the chronology compared to other ice core chronologies. A translation in ages of volcanic ties is provided to facilitate the inclusion of these data where a common age scale is required.
While we provide all the level-1 chemistry data, it is more useful to consider annualised values for the data, either summer- or winter-centred as we have described, depending on the particular analyte. Future users of the data remain free to define seasonal boundaries as best suits their particular needs.
Seasonal averages are possible as provided here in the recently published Vance et al. (2021) as an example. It should be noted that these data have undergone further quality controls resulting in no average value being recorded where there are an insufficient number of clear records within the seasonal bounds to calculate an average with confidence.
We encourage all potential users of the data to contact us if they have any further queries about how best to use the data.
The Law Dome record remains independently dated using annual layer counting techniques without reference to tie points in other ice core chronologies. As users of these data may wish to map the Law Dome data onto other commonly used ice core chronologies, we provide here a mapping of the age of volcanic tie points between the Law Dome age scale to the Ant2k (Ahmed et al., 2013), the AICC2012 (Veres et al., 2013) and WD2014 (Sigl et al., 2016) age scales. As the Law Dome record from DSS is a composite record consisting of cores drilled at different times over a period of 4 decades, it is difficult to provide a depth for each year using these chronologies so we provide a “year-to-year” mapping of the volcanic tie points across these age scales.
Due to only very short overlaps between the various ice cores, there are insufficient data for a reliable direct assessment of the consistency of chemical species between the cores. Therefore, as an alternative, we identify different epochs when we might expect similar distributions, specifically, periods of time covered by one of the ice cores where the average accumulation rate approximately the same in a different core making up the compilation but during a different period of time. Change-point analysis was performed on the cumulative summation (CUSUM) of the snow accumulation record to determine different epochs of approximately constant accumulation, a technique that has been used to detect changes in the Law Dome snowfall data by Zheng et al. (2021). The CUSUM simply sums the data anomalies and identifies a step change in the underlying data through changes in the gradient of the CUSUM. The change points in the CUSUM gradient were detected using the Ruptures Python package (Truong et al., 2020). We select epochs of the closest mean accumulation value and then compare between the different cores using a Kolmogorov–Smirnov (K–S) and a Welch's t test (W-T) on the log-transformed data for each of the chemical species. These two tests provide two different ways to test the populations, with the K–S test not assuming a normally distributed population, while log-transforming the data for the W-T test mitigates some of the long tail of the distribution. The epochs that were compared are shown in Fig. B1 with a summary of the statistical results in Table B1. Visual inspection was also performed using the empirical cumulative probability distributions of the log-transformed concentration data in Fig. B2. Overall, the results shown in Table B1 are broadly consistent across all the ice cores used in the composite record. We would not expect all of the tests to show statistical significance for several reasons: (a) the accumulation rates used for the comparisons are as similar as possible but not identical, (b) even if the accumulation rates where identical, we would expect variability in the species concentrations (i.e. the correlation with accumulation is not perfect) and (c) the presence of surface features such as sastrugi will result in different (but highly correlated) ice core records from ice cores drilled at different locations, even only a few metres apart.
All the authors were involved in the generation of the dataset through collection, analysis, quality control and dating. LMJ coordinated the compilation of the data and led the writing of the article with CP and contributions from all the other authors.
The contact author has declared that none of the authors has any competing interests.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The authors would like to acknowledge the contribution of the expeditioners involved in drilling the ice cores and laboratory support staff. This research has been supported by the Australian Government Department of Industry, Science, Energy and Resources (grant no. ASCI000002), the US National Science Foundation (OPP-9811857 and ATM-9808963), the Australian Research Council Discovery Project (DP180102522), Australian Research Council Special Research Initiative for Antarctic Gateway Partnership (SR140300001) and past funding from the Antarctic Climate & Ecosystems Cooperative Research Centre (ACE CRC, 2010-2019). The Australian Antarctic Division provided funding and logistical support (ASS project 757 and AAS projects, 4061, 4062, 4537).
This research has been supported by the Australian Antarctic Division (grant nos. AAS 757, 4061, 4062, and 4537), the Department of Industry, Science, Energy and Resources, Australian Government (grant no. ASCI000002), the National Science Foundation (grant nos. OPP-9811857 and ATM-9808963) and the Australian Research Council (grant nos. DP180102522 and SR140300001).
This paper was edited by David Carlson and reviewed by Elizabeth Thomas and two anonymous referees.
Ahmed, M., Anchukaitis, K. J., Asrat, A., Borgaonkar, H. P., Braida, M., Buckley, B. M., Büntgen, U., Chase, B. M., Christie, D. A., Cook, E. R., Curran, M. A., Diaz, H. F., Esper, J., Fan, Z. X., Gaire, N. P., Ge, Q., Gergis, J., González-Rouco, J. F., Goosse, H., Grab, S. W., Graham, N., Graham, R., Grosjean, M., Hanhijärvi, S. T., Kaufman, D. S., Kiefer, T., Kimura, K., Korhola, A. A., Krusic, P. J., Lara, A., Lézine, A. M., Ljungqvist, F. C., Lorrey, A. M., Luterbacher, J., Masson-Delmotte, V., McCarroll, D., McConnell, J. R., McKay, N. P., Morales, M. S., Moy, A. D., Mulvaney, R., Mundo, I. A., Nakatsuka, T., Nash, D. J., Neukom, R., Nicholson, S. E., Oerter, H., Palmer, J. G., Phipps, S. J., Prieto, M. R., Rivera, A., Sano, M., Severi, M., Shanahan, T. M., Shao, X., Shi, F., Sigl, M., Smerdon, J. E., Solomina, O. N., Steig, E. J., Stenni, B., Thamban, M., Trouet, V., Turney, C. S., Umer, M., van Ommen, T., Verschuren, D., Viau, A. E., Villalba, R., Vinther, B. M., Von Gunten, L., Wagner, S., Wahl, E. R., Wanner, H., Werner, J. P., White, J. W., Yasue, K., and Zorita, E.: Continental-scale temperature variability during the past two millennia, Nat. Geosci., 6, 339–346, 2013. a, b, c, d, e
Bromwich, D.: Snowfall in high southern latitudes, Reviews of Geophysics, 26, 149–168, 1988. a
Buck, C. F., Mayewski, P. A., Spencer, M. J., Whitlow, S., Twickler, M. S., and Barrett, D.: Determination of major ions in snow and ice cores by ion chromatography, J. Chromatogr. A, 594, 225–228, https://doi.org/10.1016/0021-9673(92)80334-Q, 1992. a, b, c, d
Crockart, C. K., Vance, T. R., Fraser, A. D., Abram, N. J., Criscitiello, A. S., Curran, M. A. J., Favier, V., Gallant, A. J. E., Kittel, C., Kjær, H. A., Klekociuk, A. R., Jong, L. M., Moy, A. D., Plummer, C. T., Vallelonga, P. T., Wille, J., and Zhang, L.: El Niño–Southern Oscillation signal in a new East Antarctic ice core, Mount Brown South, Clim. Past, 17, 1795–1818, https://doi.org/10.5194/cp-17-1795-2021, 2021. a
Curran, M., M., Jong, L., Moy, A., Plummer, C., Roberts, J., Van Ommen, T., and Vance, T.: The Law Dome ice core 2000 year dataset collection, Australian Antarctic Data Centre [data set], https://doi.org/10.26179/5zm0-v192, 2021. a, b, c
Curran, M. A. and Palmer, A. S.: Suppressed ion chromatography methods for the routine determination of ultra low level anions and cations in ice cores, J. Chromatogr. A, 919, 107–113, https://doi.org/10.1016/S0021-9673(01)00790-7, 2001. a, b, c, d, e, f, g, h
Curran, M. A., Van Ommen, T. D., and Morgan, V.: Seasonal characteristics of the major ions in the high-accumulation Dome Summit South ice core, Law Dome, Antarctica, Ann. Glaciol., 27, 385–390, https://doi.org/10.3189/1998AoG27-1-385-390, 1998. a
Curran, M. A. J., Van Ommen, T. D., Morgan, V. I., Phillips, K. L., and Palmer, A. S.: Ice Core Evidence for Antarctic Sea Ice Decline since the 1950s, Science, 302, 1203–1206, https://doi.org/10.1126/science.1087888, 2003. a
Emile-Geay, J., McKay, N. P., Kaufman, D. S., von Gunten, L., Wang, J., Anchukaitis, K. J., Abram, N. J., Addison, J. A., Curran, M. A., Evans, M. N., Henley, B. J., Hao, Z., Martrat, B., McGregor, H. V., Neukom, R., Pederson, G. T., Stenni, B., Thirumalai, K., Werner, J. P., Xu, C., Divine, D. V., Dixon, B. C., Gergis, J., Mundo, I. A., Nakatsuka, T., Phipps, S. J., Routson, C. C., Steig, E. J., Tierney, J. E., Tyler, J. J., Allen, K. J., Bertler, N. A., Björklund, J., Chase, B. M., Chen, M.-T., Cook, E., de Jong, R., DeLong, K. L., Dixon, D. A., Ekaykin, A. A., Ersek, V., Filipsson, H. L., Francus, P., Freund, M. B., Frezzotti, M., Gaire, N. P., Gajewski, K., Ge, Q., Goosse, H., Gornostaeva, A., Grosjean, M., Horiuchi, K., Hormes, A., Husum, K., Isaksson, E., Kandasamy, S., Kawamura, K., Kilbourne, K. H., Koç, N., Leduc, G., Linderholm, H. W., Lorrey, A. M., Mikhalenko, V., Mortyn, P. G., Motoyama, H., Moy, A. D., Mulvaney, R., Munz, P. M., Nash, D. J., Oerter, H., Opel, T., Orsi, A. J., Ovchinnikov, D. V., Porter, T. J., Roop, H. A., Saenger, C., Sano, M., Sauchyn, D., Saunders, K. M., Seidenkrantz, M.-S., Severi, M., Shao, X., Sicre, M.-A., Sigl, M., Sinclair, K., St. George, S., St. Jacques, J.-M., Thamban, M., Kuwar Thapa, U., Thomas, E. R., Turney, C., Uemura, R., Viau, A. E., Vladimirova, D. O., Wahl, E. R., White, J. W., Yu, Z., Zinke, J., and Zinke, J.: A global multiproxy database for temperature reconstructions of the Common Era, Sci. Data, 4, 170088, https://doi.org/10.1038/sdata.2017.88, 2017. a, b, c, d
Etheridge, D. M., Steele, L. P., Francey, R. J., and Langenfelds, R. L.: Atmospheric methane between 1000 A.D. and present: Evidence of anthropogenic emissions and climatic variability, J. Geophys. Res.-Atmos., 103, 15979–15993, https://doi.org/10.1029/98JD00923, 1998. a
Masson-Delmotte, V., Delmotte, M., Morgan, V., Etheridge, D., van Ommen, T., Tartarin, S., and Hoffmann, G.: Recent southern Indian Ocean climate variability inferred from a Law Dome ice core: New insights for the interpretation of coastal Antarctic isotopic records, Clim. Dynam., 21, 153–166, https://doi.org/10.1007/s00382-003-0321-9, 2003. a
Morgan, V., Wookey, C., Li, J., van Ommen, T., Skinner, W., and Fitzpatrick, M.: Site information and initial results from deep ice drilling on Law Dome, Antarctica, J. Glaciol., 43, 3–10, https://doi.org/10.3189/S0022143000002768, 1997. a, b
Palmer, A. S., Van Ommen, T. D., Curran, M. A. J., and Morgan, V.: Ice-core evidence for a small solar-source of atmospheric nitrate, Geophys. Res. Lett., 28, 1953–1956, https://doi.org/10.1029/2000GL012207, 2001a. a
Palmer, A. S., van Ommen, T. D., Curran, M. A. J., Morgan, V., Souney, J. M., and Mayewski, P. A.: High-precision dating of volcanic events (A.D. 1301–1995) using ice cores from Law Dome, Antarctica, J. Geophys. Res.-Atmos., 106, 28089–28095, https://doi.org/10.1029/2001JD000330, 2001b. a, b, c, d
Parker, D., Folland, C., Scaife, A., Knight, J., Colman, A., Baines, P., and Dong, B.: Decadal to multidecadal variability and the climate change background, J. Geophys. Res.-Atmos., 112, D18115, https://doi.org/10.1029/2007JD008411, 2007. a
Pedro, J. B., Heikkilä, U. E., Klekociuk, A., Smith, A. M., Van Ommen, T. D., and Curran, M. A.: Beryllium-10 transport to Antarctica: Results from seasonally resolved observations and modeling, J. Geophys. Res.-Atmos., 116, 1–14, https://doi.org/10.1029/2011JD016530, 2011a. a
Pedro, J. B., Smith, A. M., Simon, K. J., van Ommen, T. D., and Curran, M. A. J.: High-resolution records of the beryllium-10 solar activity proxy in ice from Law Dome, East Antarctica: measurement, reproducibility and principal trends, Clim. Past, 7, 707–721, https://doi.org/10.5194/cp-7-707-2011, 2011b. a
Pedro, J. B., McConnell, J. R., van Ommen, T. D., Fink, D., Curran, M. A., Smith, A. M., Simon, K. J., Moy, A. D., and Das, S. B.: Solar and climate influences on ice core 10Be records from Antarctica and Greenland during the neutron monitor era, Earth Planet. Sc. Lett., 355–356, 174–186, https://doi.org/10.1016/j.epsl.2012.08.038, 2012. a
Pfitzner, M. L. and Australian National Antarctic Research Expeditions: The Wilkes ice cap project, 1966, Australian Government Publishing Service, Canberra, https://nla.gov.au/nla.cat-vn8147198 (last access: 18 July 2022), 1980. a
Plummer, C. T., Curran, M. A. J., van Ommen, T. D., Rasmussen, S. O., Moy, A. D., Vance, T. R., Clausen, H. B., Vinther, B. M., and Mayewski, P. A.: An independently dated 2000-yr volcanic record from Law Dome, East Antarctica, including a new perspective on the dating of the 1450s CE eruption of Kuwae, Vanuatu, Clim. Past, 8, 1929–1940, https://doi.org/10.5194/cp-8-1929-2012, 2012. a, b, c, d, e, f, g, h, i
Roberts, J., Plummer, C., Vance, T., van Ommen, T., Moy, A., Poynter, S., Treverrow, A., Curran, M., and George, S.: A 2000-year annual record of snow accumulation rates for Law Dome, East Antarctica, Clim. Past, 11, 697–707, https://doi.org/10.5194/cp-11-697-2015, 2015. a, b, c, d
Roberts, J. L., Van Ommen, T. D., Curran, M. A. J., and Vance, T. R.: Methanesulphonic acid loss during ice-core storage: Recommendations based on a new diffusion coefficient, J. Glaciol., 55, 784–788, https://doi.org/10.3189/002214309790152474, 2009. a
Rubino, M., Etheridge, D. M., Thornton, D. P., Howden, R., Allison, C. E., Francey, R. J., Langenfelds, R. L., Steele, L. P., Trudinger, C. M., Spencer, D. A., Curran, M. A. J., van Ommen, T. D., and Smith, A. M.: Revised records of atmospheric trace gases CO2, CH4, N2O, and δ13C-CO2 over the last 2000 years from Law Dome, Antarctica, Earth Syst. Sci. Data, 11, 473–492, https://doi.org/10.5194/essd-11-473-2019, 2019. a
Scambos, T. A., Haran, T. M., Fahnestock, M. A., Painter, T. H., and Bohlander, J.: MODIS-based Mosaic of Antarctica (MOA) data sets: Continent-wide surface morphology and snow grain size, Remote Sens. Environ., 111, 242–257, https://doi.org/10.1016/j.rse.2006.12.020, 2007. a
Sigl, M., Fudge, T. J., Winstrup, M., Cole-Dai, J., Ferris, D., McConnell, J. R., Taylor, K. C., Welten, K. C., Woodruff, T. E., Adolphi, F., Bisiaux, M., Brook, E. J., Buizert, C., Caffee, M. W., Dunbar, N. W., Edwards, R., Geng, L., Iverson, N., Koffman, B., Layman, L., Maselli, O. J., McGwire, K., Muscheler, R., Nishiizumi, K., Pasteris, D. R., Rhodes, R. H., and Sowers, T. A.: The WAIS Divide deep ice core WD2014 chronology – Part 2: Annual-layer counting (0–31 ka BP), Clim. Past, 12, 769–786, https://doi.org/10.5194/cp-12-769-2016, 2016. a, b
Souney, J. M., Mayewski, P. A., Goodwin, I. D., Meeker, L. D., Morgan, V., Curran, M. A., Van Ommen, T. D., and Palmer, A. S.: A 700-year record of atmospheric circulation developed from the Law Dome ice core, East Antarctica, J. Geophys. Res.-Atmos., 107, ACL 1-1–ACL 1-9, https://doi.org/10.1029/2002JD002104, 2002. a
Stenni, B., Curran, M. A. J., Abram, N. J., Orsi, A., Goursaud, S., Masson-Delmotte, V., Neukom, R., Goosse, H., Divine, D., van Ommen, T., Steig, E. J., Dixon, D. A., Thomas, E. R., Bertler, N. A. N., Isaksson, E., Ekaykin, A., Werner, M., and Frezzotti, M.: Antarctic climate variability on regional and continental scales over the last 2000 years, Clim. Past, 13, 1609–1634, https://doi.org/10.5194/cp-13-1609-2017, 2017. a, b, c
Thomas, E. R., van Wessem, J. M., Roberts, J., Isaksson, E., Schlosser, E., Fudge, T. J., Vallelonga, P., Medley, B., Lenaerts, J., Bertler, N., van den Broeke, M. R., Dixon, D. A., Frezzotti, M., Stenni, B., Curran, M., and Ekaykin, A. A.: Regional Antarctic snow accumulation over the past 1000 years, Clim. Past, 13, 1491–1513, https://doi.org/10.5194/cp-13-1491-2017, 2017. a, b
Tozer, C. R., Kiem, A. S., Vance, T. R., Roberts, J. L., Curran, M. A., and Moy, A. D.: Reconstructing pre-instrumental streamflow in Eastern Australia using a water balance approach, J. Hydrol., 558, 632–646, https://doi.org/10.1016/j.jhydrol.2018.01.064, 2018. a
Truong, C., Oudre, L., and Vayatis, N.: Selective review of offline change point detection methods, Signal Processing, 167, 107299, https://doi.org/10.1016/j.sigpro.2019.107299, 2020. a
Udy, D. G., Vance, T. R., Kiem, A. S., Holbrook, N. J., and Curran, M. A. J.: Links between Large-Scale Modes of Climate Variability and Synoptic Weather Patterns in the Southern Indian Ocean, J. Climate, 34, 883–899, https://doi.org/10.1175/JCLI-D-20-0297.1, 2021. a, b
van Ommen, T. D. and Morgan, V.: Peroxide concentrations in the Dome Summit South ice core, Law Dome, Antarctica, J. Geophys. Res.-Atmos., 101, 15147–15152, https://doi.org/10.1029/96JD00838, 1996. a
van Ommen, T. D. and Morgan, V.: Calibrating the ice core paleothermometer using seasonality, J. Geophys. Res.-Atmos., 102, 9351–9357, https://doi.org/10.1029/96JD04014, 1997. a, b
van Ommen, T. D. and Morgan, V.: Snowfall increase in coastal East Antarctica linked with southwest Western Australian drought, Nat. Geosci., 3, 267–272, https://doi.org/10.1038/ngeo761, 2010. a
van Ommen, T. D., Morgan, V. I., Jacka, T. H., Woon, S., and Elcheikh, A.: Near-surface temperatures in the Dome Summit South (Law Dome, East Antarctica) borehole, Ann. Glaciol., 29, 141–144, https://doi.org/10.3189/172756499781821382, 1999. a
Vance, T. R., van Ommen, T. D., Curran, M. A. J., Plummer, C. T., and Moy, A. D.: A Millennial Proxy Record of ENSO and Eastern Australian Rainfall from the Law Dome Ice Core, East Antarctica, J. Climate, 26, 710–725, https://doi.org/10.1175/JCLI-D-12-00003.1, 2013. a, b, c
Vance, T. R., Roberts, J. L., Plummer, C. T., Kiem, A. S., and van Ommen, T. D.: Interdecadal Pacific variability and eastern Australian megadroughts over the last millennium, Geophys. Res. Lett., 42, 129–137, https://doi.org/10.1002/2014GL062447, 2015. a, b, c
Vance, T. R., Roberts, J. L., Moy, A. D., Curran, M. A. J., Tozer, C. R., Gallant, A. J. E., Abram, N. J., van Ommen, T. D., Young, D. A., Grima, C., Blankenship, D. D., and Siegert, M. J.: Optimal site selection for a high-resolution ice core record in East Antarctica, Clim. Past, 12, 595–610, https://doi.org/10.5194/cp-12-595-2016, 2016. a
Vance, T. R., Kiem, Anthony S., J., M., L., Roberts, J. L. R., Plummer, C. T., Moy, A. D., Curran, M. A., and van Ommen, T. D.: Pacific decadal variability over the last 2000 years and implications for climatic risk, Commun. Earth Environ., 3, 33, https://doi.org/10.1038/s43247-022-00359-z, 2022. a, b, c
Veres, D., Bazin, L., Landais, A., Toyé Mahamadou Kele, H., Lemieux-Dudon, B., Parrenin, F., Martinerie, P., Blayo, E., Blunier, T., Capron, E., Chappellaz, J., Rasmussen, S. O., Severi, M., Svensson, A., Vinther, B., and Wolff, E. W.: The Antarctic ice core chronology (AICC2012): an optimized multi-parameter and multi-site dating approach for the last 120 thousand years, Clim. Past, 9, 1733–1748, https://doi.org/10.5194/cp-9-1733-2013, 2013. a
Zheng, Y., Jong, L. M., Phipps, S. J., Roberts, J. L., Moy, A. D., Curran, M. A. J., and van Ommen, T. D.: Extending and understanding the South West Western Australian rainfall record using a snowfall reconstruction from Law Dome, East Antarctica, Clim. Past, 17, 1973–1987, https://doi.org/10.5194/cp-17-1973-2021, 2021. a