The RICE discrete methane record was measured at Oregon State University (OSU) following methods described by with updates described in Appendix . A total of 702 samples were measured at 583 distinct depths between 60 and 753 m (Fig. a, f). Samples from 406 depths were measured between 60 and 670 m, dating from ∼1970 CE to 11.87 ka, with a mean sample spacing of 28.75 years. Between 670 and 718.13 m the record spans 11.87 to 29.9 ka; 96 samples in this interval provide age resolution of 189 years. Age resolution decreases significantly for deeper depths. The interval from 718.53 to 746.00 m corresponds to 30.1–64.6 ka with a mean resolution of 548 years. The deepest dated ice is at 752.95 m with an age of 83 ka (±2 ka).

Methane was also measured continuously with a laser spectroscopy technique during two separate continuous flow analysis (CFA) campaigns at GNS Science (Gracefield, New Zealand) in 2013 and 2014 . The CFA methane record was affected by variability in airflow to the measurement instrument and fractures that allowed drill fluid and modern air into the melt head. Exclusion of these artifacts caused significant gaps in the record, particularly at depths below 676 m (12.6 ka). We consider the CFA record to be a supplement to the more robust but lower-resolution discrete data. Between 29.9 and 59.1 ka, the CFA record is critical for establishing age control. Below ∼746 m (64.6 ka), the CFA data are difficult to interpret because of gaps in the record, uncertainty in measurement depth, and uncertainty in methane calibration (Appendix ).

Several anomalously high discrete methane measurements appear between 44.6 and 50.9 ka (729.05–736.05 m) and below 64.6 ka (746 m) (Fig. a). In the former interval, methane is enriched by ∼30 ppb compared with the WAIS Divide CFA data . This is likely due to healed fractures that include modern air . All RICE samples deeper than 500 m were visually inspected for fractures as they were prepared for measurement and for drill fluid during air extraction. In fractured ice it was common to see drill fluid in the flask as the samples were melted. However, neither observation was a strong indicator that a sample would have an elevated methane concentration. None of these high-concentration results were rejected.

Total air content (TAC) is defined as the amount trapped air in a sample, reported here in units of cubic centimeters of air at STP per gram of ice (cm3 g-1). TAC (Fig. e, j) was measured at OSU as part of the methane concentration measurement following methodology of . TAC is influenced by accumulation and temperature, total summer insolation , thermal gradients in the firn from multiannual climate trends, and surface air pressure .

Preparation of ice samples causes loss of air from bubbles which intersect the surface of the sample, an effect known as the cut-bubble effect . The amount of air lost due to the cut-bubble effect depends on the geometry of the sample and the bubble size. For bubble size, we assume that the bubble diameter at the bubble-close-off depth is the same as at the bubble-close-off depth at WAIS Divide and that the diameter shrinks exponentially to zero at the bubble–clathrate-transition depth (763 m). Ice samples were cut to uniform shapes to limit variability in TAC related to the cut-bubble effect so that TAC can be directly compared between samples. However, inclusion of some fractures was often unavoidable, and their contribution to gas loss is potentially large. TAC data were rejected when gas loss was believed to greatly impact the results, such as in samples with visible fractures or samples which consisted of multiple pieces. Of the 706 samples measured at OSU for TAC, 165 results were rejected based upon visual inspection of the sample. Many of these came from the 670 to 752.95 m (11.7–83 ka) interval where only 58 of 177 TAC measurements are considered reliable. Reproducibility of replicate TAC measurements is 0.7 % after application of a correction for the cut-bubble effect (reproducibility of 0.6 % without correction). Estimate of relative uncertainty in TAC measurements is 28 % with the largest contributions to the uncertainty from the shape of the sample and the assumed bubble diameter used for the cut-bubble correction.

The TAC record from the RICE ice core without the cut-bubble correction appears remarkably consistent (age weighted mean =0.1182 cm3 g-1; age weighted standard deviation =0.0023 cm3 g-1; n=410). After application of the cut-bubble correction, the TAC record shows a trend to higher values at more recent ages (Fig. e). This trend reflects the decrease in bubble size with depth and therefore a larger cut-bubble correction for shallower samples. Although this trend has a physical explanation, it is not statistically significant due to the uncertainty in the cut-bubble correction. We advise caution when interpreting the absolute values of TAC from the RICE ice core due to the potential artifacts caused by gas loss through fractures and the large uncertainty in the cut-bubble correction.

δ18Oatm and δ15N-N2 were measured on samples adjacent to the discrete methane samples (Fig. c, h and d, i, respectively). Analysis was conducted at Scripps Institution of Oceanography following and . Pooled standard deviation of replicate measurements is 0.006 ‰ for δ18Oatm and 0.0027 ‰ for δ15N-N2 (both scales are relative to modern atmospheric composition).

Variations in δ18Oatm are primarily caused by changes in location and intensity of low-latitude rainfall that affect the δ18O of leaf water used in photosynthesis and changes in seawater δ18O caused by ice sheets on glacial cycles . Importantly, these variations are well known in independently dated ice cores, and the atmosphere is well mixed on the relevant timescales, so δ18Oatm variability may be used as a chronostratigraphic marker . These variations are well sampled in the RICE record until ∼64.6 ka (746.00 m), beyond which the chronology is no longer continuous.

Atmospheric N2 is isotopically stable over very long timescales . Variability in δ15N-N2 in the ice core records primarily reflects changes in gravitational fractionation and thermal fractionation within the firn due to changes in surface temperature and accumulation rate . In Antarctica, temperature changes gradually, and the effect of thermal separation is minimal.

The RICE δD record (Fig. b, g) was measured using off-axis integrated cavity output spectroscopy (OA-ICOS) using a Los Gatos Research (LGR) isotope water analyzer (IWA-35EP) in combination with a modified LGR water vapor isotope standard source (WVISS) evaporation unit . The data were calibrated using a time-weighted, two-point linear calibration with two standards. The cumulative uncertainty is calculated as the quadrature sum of Allan variance error, scatter over the averaging (15 s) interval (error of the variance) and calibration error (error of the mean). Final mean total uncertainty for the RICE δD record is 0.74 ‰ . The record shown in Fig. b has been smoothed with a 500-year low-pass filter (LOESS filter) and sampled at 50-year intervals. The smoothed δD record is provided in the Supplement to this paper.

An automated matching algorithm was adapted from to synchronize between gas records. Prior ACPs, shown as open triangles in Figs.  and , all correspond to well-defined increases or decreases in either methane or δ18Oatm. We assume that the uncertainty (2 standard deviations) for an ACP corresponds to the time elapsed between 25 % and 75 % of the change in either methane or δ18Oatm (Table ). The goal of the routine is to iteratively adjust the interpolation between ACPs to improve a goodness of fit parameter by following steps 1–9 described below.

“Goodness of fit” is calculated as the χ2 value comparing the normalized methane record from the RICE ice core to expectant values interpolated from the WAIS Divide ice core plus the analogous χ2 value comparing the normalized δ18Oatm records. The algorithm can accept additional parameters for synchronization, such as CO2 or N2O, if available. In this analysis we normalized both the methane and δ18Oatm records to have the same mean (5) and variance (1) in order to equally weight their χ2 values and their importance to the synchronization.

A realization starts by randomly perturbing the age of ACPs within their prescribed uncertainty (Table ) to define an initial depth–age scale. The perturbed ACPs remain fixed throughout the realization. ACPs are given the subscript k, where k=1 is the shallowest (youngest) ACP. The records are then broken into N subsections which are distributed between ACPs. Subsections are designated with a subscript i, where i=1,2,3,…N. Floating control points (FCPs) are defined as the bounding depths (ages) of these subsections, which will include the initial ACPs. Prior to any perturbations, the durations of the subsections (Δt) are approximately the same. We have chosen Δt to roughly match the recurrence interval of variations in methane in the reference record and so that each subsection contains about five methane and five δ18Oatm samples. The following steps are repeated to optimize goodness of fit:

Random scaling factors (pi), which perturb the durations of the subsections, are drawn for each subsection from a normal distribution of μ=1 and σ1=0.25.1Δti′=Δti⋅pi.

A Monte Carlo analysis repeats these steps (Steps 1–9) 1000 times, initiated from a different prior depth–age relationship by randomly perturbing ACPs within their age uncertainty (Table ). Parameters used for the synchronization are given in Table  for the 60.05–670.13 m and 670.13–719.30 m intervals. Code for the synchronization routine is provided in the Supplement.

We choose the best realization for gas age constraints for the RICE17 gas age-scale. The best age estimate (realization) is not necessarily the most frequent age estimate. Figure d shows an example from sample depth 621.28 m where there is a large difference between these two age estimates. Large differences can occur because the prior age estimate (i.e., the age estimate based only on prior ACPs) differs by a large amount from the “true” age of that sample and because the goodness-of-fit parameter considers the fit over the whole record. For depth 621.28 m most realizations resulted in an age estimate of 9200 years BP, similar to its prior age estimate of 9240 years BP, but the best realization estimated the age to be 9012 years BP. There are two possible reasons for this type of result: (1) this realization is the best because it managed to push the age of this sample by >200 years towards younger ages while not significantly changing the ages of nearby sections which already fit well or (2) no significant improvement in the goodness of fit was found by adjusting the age of this depth, and the goodness of fit was dictated by other sections of the record.

Uncertainty in gas age constraints is calculated as the root-mean-square error of FCP ages from the 1000 Monte Carlo realizations (Fig. d), although we acknowledge that this can overly simplify the empirical distribution of possible ages for a given depth (Fig. ). Assessing uncertainty in this way integrates two types of error. The first is the ability to exactly match timing of two features. This uncertainty is determined by how well features are resolved in the records, measurement error, and the limited degrees of freedom in synchronization (i.e., the synchronization routine assumes constant accumulation during each subsections in a piecewise manner, although the true accumulation history is more likely to smoothly vary). The second type of error is analogous to deciding which feature in the reference record is the correct match. For example, the methane peak centered at 459.05 m (4675 years BP) is, in some realizations, matched to a peak at ∼4550 years BP in the WAIS Divide record, providing two distinct age possibilities (Fig. c).

The more limited resolution of the RICE methane record below 719.3 m (30.66 ka) prevented use of the automated synchronization routine. Gas age control points were visually chosen between 719.30 m (30.66 ka) and 746.00 m (65.6 ka) and between 746.00 and 752.95 m (∼83 ka) (Table ). Between 719.30 (30.66 ka) and 746.00 m (65.6 ka), ACPs were visually chosen by comparing the RICE discrete and CFA methane records and the δ18Oatm record to the WAIS Divide methane and δ18Oatm records .

The WAIS Divide records and the WD2014 chronology end at 67.8 ka. For the deeper (older) section of the RICE ice core, between 746.00 and 752.95 m, we build target records from records of methane and δ18Oatm from the NGRIP ice core . Using the NGRIP ice core is internally consistent with the WD2014 chronology, which is tied to a modified GICC05modelext age-scale between 30 and 67.8 ka . Because NGRIP is in the Northern Hemisphere the methane concentration is higher than that in the RICE record. To account for this interpolar difference, we scale the NGRIP target methane record to Antarctic values using an empirical least-squares fit between WAIS and NGRIP records between 55 and 67.8 ka.

found that the annual-layer-counted portion of the GICC05modelext chronology (0–60 ka) is systematically younger than ages of corresponding features found in the U–Th absolute dated Hulu speleothem record . A fit to Hulu ages was optimized by scaling the GICC05modelext ages linearly by 1.0063. This suggests that on average 6.3 out of every 1000 annual layers were not counted. For our NGRIP-based target records we adjust the NGRIP age-scale by adopting the scaling of . Older ages in the GICC05modelext chronology are derived from the ss09sea Dansgaard–Johnsen model . We have added a constant 378 years (0.0063×60000) for depths older than 60 ka in the target GICC05modelext ages to make this section continuous with the adjusted annual-layer-counted section.

Synchronization requires well-sampled identifiable features in both records. We estimate age uncertainty for visually matched ACPs (ACPs older than 29.9 ka) as the larger of the sample spacing of the two records being synchronized, following methods of . Gas age uncertainty is plotted in Fig. d. The largest uncertainty is ∼1500–1700 years and occurs between 41.7 and 47.6 ka, where the age-scale is very compressed, and the RICE records are coarsely sampled.

Continuity of the RICE ice core appears to end at 746.00 m below surface, where a discontinuity of 20 ‰ is observed in the δD record (746.00–746.10 m) (Fig. b), and a 0.35 ‰ change is observed between sequential δ18Oatm samples (745.81 and 746.20 m) (Fig. c). The continuous chronology described in the previous section dates this ice to the end of DO-18 (∼64.6 ka). This age is supported by very negative δD values and a trend to more enriched values at shallower depths, which is consistent with warming trends observed in other Antarctic ice cores at this time period (Fig. f).

Figure  presents our best effort to date the ice in the disturbed section, between 746.00 and 752.95 m. The left panel of Fig.  shows the data as a function of depth, color-coded by age. The right panel shows the reconstructed time history of these variables, color-coded by depth. Below, we discuss how the ages are assigned to different depth sections (divided by vertical red lines in the left side of Fig. ), starting with the shallowest portion. Dating folded ice remains challenging, and the solution we found may not be unique.

Immediately below the discontinuity, between 746.10 and 747.85 m, methane and δ18Oatm values best match our target records during DO-20 between 74.5 and 77.3 ka (Fig. ), although the RICE methane record appears ∼30 ppb too high. Dating of this section to DO-20 implies either that 9.4 ka of climate history is absent from the RICE ice core because climate was not recorded or that 9.4 ka is compressed into ∼10 cm of ice. Extremely thin layers can be explained by ice flow or by an absence of accumulation, but the latter scenario would cause a collapse of the firn and δ15N-N2 values approaching 0 ‰, which is not observed.

A cluster of samples between 747.85 and 750.46 m are ∼ 0.16 ‰ more enriched in δ18Oatm than was observed in the shallower section dated to DO-20 (Fig. c). In the NGRIP, EDML, and Siple Dome ice cores, a long-term trend in δ18Oatm towards more enriched values is observed from 80 to 65 ka (Fig. g). The enriched values between 748.34 and 750.46 m most likely indicate that this ice is younger than the adjacent shallower depths and that the stratigraphic order is reversed. This depth range is best matched to DO-19 between 71.8 and 74.3 ka (Fig. ). Below 750.46 m, methane and δ18Oatm return to values that best match DO-20.

We note that two significant measurement gaps exist in the δD record, at 746.89–747.07 and at 750.00–750.25 m (vertical bars in Fig. b). A 12 ‰ shift in δD accompanies both of these gaps. The second measurement gap (750.46–750.56 m) is at nearly the same depth which separates the samples dating to DO-19 from the grouping immediately below dated to DO-20. The discontinuity of the δD record at these gaps may signify another hiatus in the climate record or highly contorted layers that are typically found around folds .

Below 750.46 meters, trends in methane and δ18Oatm indicate that the age of ice increases with depth until at least 753 m (∼83 ka). Age of this depth is constrained by a measured δ18Oatm value of -0.175 ‰. Such a depleted value is rare and only occurs during periods of high sea level and small ice sheets, the most recent time period prior to the Holocene being MIS-5a (80–85 ka) . Negative values were observed at two other depths below 753 m (Fig. h), but the stratigraphic order of these depths is difficult to assess.

Air transport in the firn causes an age difference between ice and air trapped in the ice, commonly referred to as Δage. Firn densification models are typically used to simulate past Δage . Input parameters (for example, temperature, accumulation rate, and surface density) are normally the largest sources of uncertainty. Using δ15N-N2 as a proxy for past firn column thickness and assuming δD records past site temperature faithfully, firn densification models can be run in an inverse mode to estimate both past Δage and accumulation rates , which is the approach we employ here.

We use a dynamic version of the Herron–Langway model , which was also used for construction of the WD2014 chronology . The model simulates firn compaction rates and vertical heat diffusion and advection. The model domain is the full 763 m ice column at 0.5 m resolution; a time step of 1 year is used. The model simulates both gravitational enrichment of δ15N-N2 and fractionation in the presence of thermal gradients. The model is forced with a temperature history derived from the δD record assuming a constant isotope sensitivity of 6 ‰ K-1 . A smoothed (500 years LOESS filter) version of the δD data is shown in Fig. b and provided in the Supplement. The Δage model uses a heat diffusion model to determine firn temperatures. Differences between patterns of modeled temperature variation and that assumed from the smoothed δD record are minor (equivalent to Δage variations of decades, within model uncertainty). The assumed isotope-temperature sensitivity is similar to that observed for West Antarctica but somewhat higher than that obtained by comparing the RICE isotope record to ERA interim data . The dependency of Δage on the assumed isotope sensitivity was explored in a model sensitivity test (Appendix ) and is incorporated into the Δage uncertainty estimates. The model further assumes a constant ice thickness, a constant 2 m convective zone similar to other high-accumulation dome sites with low mean wind speeds , surface firn density of 400 kg m-3 to match the modern surface firn density, and geothermal heat flux of 78 mW m-2.

In a first iteration, we assume a prior ice age-scale for the temperature history by adding a constant 150 years to the gas chronology (Sect.  and ). A new Δage solution is then calculated using the dynamic firn densification model. The Δage solution (and thereby the ice age-scale) is refined iteratively until it no longer changes appreciably (consistent within ∼1 year).

The climate at Roosevelt Island, with high accumulation and relatively warm temperatures, results in small Δage values and consequently relatively low absolute uncertainty in Δage compared to most Antarctic sites. Modern Δage (estimated at 60 m depth, the shallowest measurement of δ15N-N2) is 146 years with a reconstructed lock-in depth (LID) of 48 m, consistent with the modern density profile . Holocene Δage is small, ranging between 140 and 182 years. During the last glacial period, simulated Δage values fluctuate between ∼150 and 350 years.

Uncertainties in past Δage include the uncertainty in model inputs as well as the model itself. The uncertainty in the Herron–Langway model is conservatively estimated to be 20 %, based on differences between firn models . We assessed the uncertainty due to model inputs using a steady-state Herron–Langway model that approximates the dynamic version but requires less computational time (Appendix ). In a sensitivity test, we randomly perturb the parameterizations of temperature and LID and assumptions of convective-zone thickness, surface firn density, and geothermal heat flux and recalculate the Δage solution (model parameters and their base values and ranges used in the sensitivity test are provided in Table ). A total of 6000 iterations were run. The sensitivity tests include a wide range of temperature histories to account for the possibility that some variations in δD were caused by nonthermal effects such as variability in precipitation seasonality, moisture sources and pathways, and post-depositional vapor exchange. Model sensitivity is reported as the root-mean-square error of Δage calculations for each depth. Combined Δage uncertainty, provided in the included age-scale, is the root sum squares of the model uncertainty and the model sensitivity.

The gas-derived ice age-scale provides a chronology from 60 to 753 m depth (Fig. ) that is independent of but overlaps the annual-layer-counted section (0–343.7 m; ). Figure  shows both chronologies and differences between them (positive values indicate that the annual layer counts are older than the gas-derived ages). The gas-derived ice age-scale agrees well with the annual-layer-counted age-scale, within 33 years, with similar trends in the implied annual layer thickness (Fig. c). Differences between chronologies can result from error in the synchronization of gas records, calculation of Δage for either the RICE or WAIS Divide ice cores, interpolation between ACPs, or annual layer counts in either core (because the WAIS Divide age-scale is used as a reference for the gas age-scale).

The average age difference at depths of FCPs is -1 years (n=18, implied age from gas-derived ice age-scale being older than the layer counted chronology). Root mean square of the age difference is 17.3 years. Discrete points in the gas-derived ice age-scale can also be compared to the age of 67 volcanic peaks identified in the RICE ice core and correlated to peaks identified in the WAIS Divide ice core (open red circles in Fig. ; ). Compared to these volcanic peaks, the root-mean-square ice age difference in the gas-derived ice age-scale is 13.6 years from their WD2014 ages. Good agreement between the two approaches gives confidence in the methodology used for the deeper section of the RICE17 chronology. The two largest differences occur at 89.72 m (243 years BP) and at 169.11 m (757 years BP) and are +30 and -33 years, respectively (Fig. b). These offsets are similar in magnitude to the individual uncertainties in calculating Δage or in synchronizing the gas records. The small age differences between the two ice chronologies also indicate that our approach to calculating uncertainty is likely conservative.

The RICE17 timescale transitions between the annual-layer-counted and gas-derived age-scales at 343.7 m, the deepest (oldest) FCP for which annual layers were identifiable (Sect. ). The age difference at this depth is 3 years, with the gas-derived ice age implying an older age than the annual-layer-counted chronology. This age difference is much less than the respective age uncertainties of 45 and 111 years (2σ) for the annual layer counts and gas-derived ice age-scale, respectively. Because of the good agreement between the two ice age-scales, no correction is made for the 3-year offset, and the annual-layer-counted age is used.

From 65 to 32 ka, the low variability in δ15N-N2, ranging from 0.20 ‰ to 0.24 ‰, and the cooling trend observed for Antarctica for this time period suggest a decreasing trend in accumulation (Fig. ).

While δ15N-N2 values appear steady from 65 to 32 ka and during the Holocene (Fig. d), they exhibit large variability between 32 and 10 ka (Fig. a). δ15N-N2 generally trends to heavier values beginning at 32 until 14.71 ka, indicating a growing firn column. An inflection point is observed at 25.3 ka, which is interpreted as an acceleration of firn thickening concluding in a peak of 0.293 ‰ at 21.8 ka. After a short depletion the steep trend resumed with a second δ15N-N2 maximum of 0.326 ‰ reached at 15.7 ka. Following this maximum, δ15N-N2 abruptly decreased from 0.308 ‰ to 0.220 ‰ at 14.71 ka, which corresponds with the beginning of the Antarctic Cold Reversal (ACR). At 12.38 ka, after the ACR, δ15N-N2 partially recovers to pre-ACR values with an abrupt increase to 0.260 ‰.

Between 25.3 and 32 ka, accumulation is estimated to be ∼ 10 cm ice equivalent per year (Fig. e), and the increasing firn thickness is largely attributable to decreasing temperature. After 25.3 ka accumulation starts to increase, and by the first δ15N-N2 maximum (21.8 ka), accumulation had increased to ∼17 cm ice equivalent per year. This feature is not apparent in other ice cores from the Ross Sea region, but those records tend to be difficult to interpret because of chronological uncertainties, such as in the case of Taylor Dome , or because of unexplained jumps in δ15N-N2, such as in the case of Siple Dome . By the second maximum (15.7 ka), accumulation increased to ∼25 cm ice equivalent per year, similar to accumulation rates observed during the Holocene. An accumulation peak at the end of glacial terminations is consistent with evidence from trimlines in interior WAIS and is also observed in the accumulation histories from WAIS Divide and Siple Dome but to a smaller degree. The early accumulation peak at 21.8 ka is unique to the RICE ice core.

A large, rapid depletion in δ15N-N2 is observed at 14.71 ka with low values lasting until 12.38 ka. Analysis of the lead–lag relationship between δ15N-N2 and methane has been used to infer near-synchronous climate changes throughout the tropics and Northern Hemisphere . These climate events are believed to originate in the Northern Hemisphere and propagate to the Southern Hemisphere . Curiously, the abrupt decrease in δ15N-N2 at 14.71 ka (0.088 ‰, 683.70–681.80 m) is observed in samples deeper than the increase in methane marking the onset of the Bølling-Allerød (14.66–14.52 ka, 683.13–682.52 m), meaning that this climate event unambiguously precedes the Northern Hemisphere event (Fig. ). Where as specifically considered the thermal component of δ15N-N2, the thermal component is considered to be negligible in Antarctic cores because air temperature only changes gradually. At Roosevelt Island, this period of low δ15N-N2 values is interpreted to represent shallow firn thickness which in the firn model is caused by low snow accumulation (< 10 cm yr-1; Fig. e). Low accumulation during this period is consistent with thin annual layers interpreted from the age–depth relationship (Fig. f), i.e., 0.3–0.6 cm yr-1 annual layer thickness compared with 1.6–3.2 cm yr-1 between 10.09 and 11.01 ka (642.75–661.07 m). Following the ACR, the modeled accumulation rate fully recovers to the ∼25 cm yr-1 observed before the ACR.

We propose three potential explanations for the thin annual layers and low accumulation rate observed during the ACR, with the last explanation being our preferred. (1) A large accumulation gradient is observed across the Roosevelt Island ice dome , which implies that a period of low accumulation at the RICE site may be the result of changes to the geometry of the Roosevelt Island ice dome. However, ice divide migration typically occurs over long timescales. (2) Accelerated ice flow may also cause thin annual layers and could potentially even affect layers within the firn. This flow could be the result of the ice streams which surrounded and buttressed Roosevelt Island ice dome suddenly being removed . The timing of the low-accumulation interval is similar to when dust records from the Taylor Dome ice core, a site located in the Transantarctic Mountains near the Ross Sea (Fig. ), show an increase in dust diameter and a change in radiogenic isotope composition, which implies that the provenance of dust changed to a local source . interpreted this source as newly exposed land created by the withdrawal of Ross Ice Shelf. The timing of the low-accumulation period is also similar to some estimates of the timing of retreat of the WAIS based on sediment facies succession and radiocarbon dating (minimum date of 8.6 ka; ) and ice-sheet modeling . Although annual layer thickness immediately above the ACR section is observed to be nearly 5 times as thick as the ACR section, this large change in thickness, if solely the result of changes in ice flow, would require an unrealistic change in the thickness of the dome or its ice flow which is not supported by TAC measurements or ice-flow modeling.

(3) Our preferred explanation is that an interval of low accumulation between 12.38 and 14.71 ka resulted in a shallow firn structure and depleted δ15N-N2 values. In recent times, moisture arriving at Roosevelt Island is frequently related to enhanced cyclonic airflow in the Ross Sea and a strong Amundsen Sea Low . This period of low accumulation may indicate a changed atmospheric structure in the South Pacific where southward airflow is blocked by a persistent low-pressure zone north of the Ross Sea such as observed in more recent periods in the accumulation record from RICE and potentially related to the past ice shelf extent. Such a pronounced minimum in accumulation is not observed in the Siple Dome ice core, where annual layers are observed to thicken during the ACR . If nonthermal effects influenced the RICE δD record, which is interpreted as temperature in our firn model, then the magnitude of accumulation change during this period may not be as large as reconstructed from the firn model. While not currently available, measurements of dexcess, dust particle size distributions, and dust geochemistry may be helpful in explaining the temperature and accumulation history of RICE.

Early reconstructions of the Ross Ice sheet during the LGM were based on glacial geological constraints from the western margin of the embayment. describe a maximum scenario in which thick ice in the Ross Embayment overrode both Siple Dome and Roosevelt Island. However, more recent observations and model experiments indicate a “fast and thin” scenario in which Siple Dome was not overrun by the interior ice sheet during the LGM , and although the Ross ice streams likely slowed down during the LGM, they remained active, maintaining a low elevation profile of the ice sheet in the Ross Sea .

Our results further support the fast and thin scenario and add a key new constraint on ice thickness and thinning in the eastern Ross Sea. Specifically, our results suggest that ice deposited on Roosevelt Island originated as accumulation local to the drilling site, which may not be true if WAIS was thick during the LGM and overrode Roosevelt Island. If remnants of WAIS were stranded on Roosevelt Island, this would likely result in a hiatus in the gas and ice chronology, in values of δD or TAC characteristic of more continental precipitation or much higher elevations, or in discontinuities in the δ15N-N2 record, indicating a much different firn structure.

While the RICE ice core chronology does exhibit an abrupt shift in δ15N-N2 and TAC during the ACR (Fig. ), no discontinuity in δD was observed, meaning that it is unlikely that any of this ice originated as part of WAIS. At least one age discontinuity (at 64.6 ka) as well as an age reversal was observed deeper in the core, but these depths do not coincide with the timing of the retreat of WAIS in the Ross Sea . The largest discontinuity, at 746.00–746.10 m, is accompanied a 20 ‰ change in δD. Methane and δ18Oatm indicate that this discontinuity represents a 9.4 ka age gap. Over the same age range, the EDML ice core records a similar change in δD (Fig. f), indicating that this change in δD is explained by Antarctic climate patterns alone and without invoking large changes in ice sheet configuration. A possible second discontinuity is observed at 747.00 m depth at a small gap in measurements of 2.7 cm. Ice immediately below 747.00 m is interpreted as being warmer, meaning that this is probably not derived from somewhere upstream in WAIS or from a higher elevation (Fig. b). Additionally, no dramatic or sudden change in δ15N-N2 or TAC was observed in association with either of these discontinuities (Fig. ). We conclude that it is highly unlikely that the accumulation site of the RICE ice core changed during the deglaciation and that Roosevelt Island ice dome probably remained independent of an advanced WAIS during the LGM. A similar conclusion was reached about Siple Dome during the LGM .

Geomorphological features on the Ross Sea bed do provide evidence of grounded ice north of Roosevelt Island during the LGM . If these features were formed by an extended WAIS, it would imply that ice flowed around Roosevelt Island and Siple Dome and therefore was limited in thickness. These conditions would indicate that ice streams were active throughout the last glacial period. Alternatively, these geomorphic features may be the result of ice from a different origin. proposed that during the LGM, an ice dome may have existed on Mary Byrd Land. In this scenario, thick, grounded ice could exist north of Roosevelt Island without flowing over or around the Roosevelt Island sea rise. The RICE records can not distinguish between these scenarios.