Abrupt climate change as an important agent of ecological change in the Northeast U.S. throughout the past 15,000 years

We use a series of tests to evaluate two competing hypotheses about the association of climate and vegetation trends in the northeastern United States over the past 15 kyrs. First, that abrupt climate changes on the scale of centuries had little influence on long-term vegetation trends, and second, that abrupt climate changes interacted with slower climate trends to determine the regional sequence of vegetation phases. Our results support the second. Large dissimilarity between temporally close fossil pollen samples indicates large vegetation changes within 500 years across >4° of latitude at ca 13.25–12.75, 12.0–11.5, 10.5, 8.25, and 5.25 ka. The evidence of vegetation change coincides with independent isotopic and sedimentary indicators of rapid shifts in temperature and moisture balance. In several cases, abrupt changes reversed long-term vegetation trends, such as when spruce (Picea) and pine (Pinus) pollen percentages rapidly declined to the north and increased to the south at ca 13.25–12.75 and 8.25 ka respectively. Abrupt events accelerated other long-term trends, such as a regional increase in beech (Fagus) pollen percentages at 8.5–8.0 ka. The regional hemlock (Tsuga) decline at ca 5.25 ka is unique among the abrupt events, and may have been induced by high climatic variability (i.e., repeated severe droughts from 5.7 to 2.0 ka); autoregressive ecological and evolutionary processes could have maintained low hemlock abundance until ca 2.0 ka. Delayed increases in chestnut (Castanea) pollen abundance after 5.8 and 2.5 ka also illustrate the potential for multi-century climate variability to influence species' recruitment as well as mortality. Future climate changes will probably also rapidly initiate persistent vegetation change, particularly by acting as broad, regional-scale disturbances.     

Coastal marsh response to historical and future sea-level acceleration

We consider the response of marshland to accelerations in the rate of sea-level rise by utilizing two previously described numerical models of marsh elevation. In a model designed for the Scheldt Estuary (Belgium–SW Netherlands), a feedback between inundation depth and suspended sediment concentrations allows marshes to quickly adjust their elevation to a change in sea-level rise rate. In a model designed for the North Inlet Estuary (South Carolina), a feedback between inundation and vegetation growth allows similar adjustment. Although the models differ in their approach, we find that they predict surprisingly similar responses to sea-level change. Marsh elevations adjust to a step change in the rate of sea-level rise in about 100 years. In the case of a continuous acceleration in the rate of sea-level rise, modeled accretion rates lag behind sea-level rise rates by about 20 years, and never obtain equilibrium. Regardless of the style of acceleration, the models predict approximately 6–14 cm of marsh submergence in response to historical sea-level acceleration, and 3–4 cm of marsh submergence in response to a projected scenario of sea-level rise over the next century. While marshes already low in the tidal frame would be susceptible to these depth changes, our modeling results suggest that factors other than historical sea-level acceleration are more important for observations of degradation in most marshes today.     

Cassini imaging of Saturn's rings: II. A wavelet technique for analysis of density waves and other radial structure in the rings

We describe a powerful signal processing method, the continuous wavelet transform, and use it to analyze radial structure in Cassini ISS images of Saturn's rings. Wavelet analysis locally separates signal components in frequency space, causing many structures to become evident that are difficult to observe with the naked eye. Density waves, generated at resonances with saturnian satellites orbiting outside (or within) the rings, are particularly amenable to such analysis. We identify a number of previously unobserved weak waves, and demonstrate the wavelet transform's ability to isolate multiple waves superimposed on top of one another. We also present two wave-like structures that we are unable to conclusively identify. In a multi-step semi-automated process, we recover four parameters from clearly observed weak spiral density waves: the local ring surface density, the local ring viscosity, the precise resonance location (useful for pointing images, and potentially for refining saturnian astrometry), and the wave amplitude (potentially providing new constraints upon the masses of the perturbing moons). Our derived surface densities have less scatter than previous measurements that were derived from stronger non-linear waves, and suggest a gentle linear increase in surface density from the inner to the mid-A Ring. We show that ring viscosity consistently increases from the Cassini Division outward to the Encke Gap. Meaningful upper limits on ring thickness can be placed on the Cassini Division (3.0 m at r∼118,800 km, 4.5 m at r∼120,700 km) and the inner A Ring (10-15 m for r<127,000 km).     

Climatic and environmental controls on speleothem oxygen-isotope values

Variations in speleothem oxygen-isotope values (δ¹⁸O) result from a complicated interplay of environmental controls and processes in the ocean, atmosphere, soil zone, epikarst, and cave system. As such, the controls on speleothem δ¹⁸O values are extremely complex. An understanding of the processes that control equilibrium and kinetic fractionation of oxygen isotopes in water and carbonate species is essential for the proper interpretation of speleothem δ¹⁸O as paleoclimate and paleoenvironmental proxies, and is best complemented by study of site-specific cave processes such as infiltration, flow routing, drip seasonality and saturation state, and cave microclimate, among others. This review is a process-based summary of the multiple controls on δ¹⁸O in the atmosphere, soil, epikarst, and speleothem calcite, illustrated with case studies. Primary controls of δ¹⁸O in the atmosphere include temperature and relative humidity through their role in the multiple isotope “effects”. Variability and modifications of water δ¹⁸O values in the soil and epikarst zones are dominated by evaporation, mixing, and infiltration of source waters. The isotopically effective recharge into a cave system consists of those waters that participate in precipitation of CaCO₃, resulting in calcite deposition rates which may be biased to time periods with optimal dripwater saturation state. Recent modeling, experimental, and observational data yield insight into the significance of kinetic fractionation between dissolved carbonate phases and solid CaCO₃, and have implications for the ‘Hendy’ test. To assist interpretation of speleothem δ¹⁸O time series, quantitative and semi-quantitative δ¹⁸O-climate calibrations are discussed with an emphasis on some of the difficulties inherent in using modern spatial and temporal isotope gradients to interpret speleothems as paleoclimate proxy records. Finally, several case studies of globally significant speleothem paleoclimate records are discussed that show the utility of δ¹⁸O to reconstruct past climate changes in regions that have been typically poorly represented in paleoclimate records, such as tropical and subtropical terrestrial locations. The new approach to speleothem paleoclimatology emphasizes climate teleconnections between regions and attribution of forcing mechanisms. Such investigations allow paleoclimatologists to infer regional to global-scale climate dynamics.     

Late Quaternary moisture export across Central America and to Greenland: evidence for tropical rainfall variability from Costa Rican stalagmites

We present a high-resolution terrestrial archive of Central American rainfall over the period 100–24 and 8.1–6.5 ka, based on δ¹⁸O time series from U-series dated stalagmites collected from a cave on the Pacific Coast of Costa Rica. Our results indicate substantial δ¹⁸O variability on millennial to orbital time scales that is interpreted to reflect rainfall variations over the cave site. Correlations with other paleoclimate proxy records suggest that the rainfall variations are forced by sea surface temperatures (SST) in the Atlantic and Pacific Oceans in a fashion analogous to the modern climate cycle. Higher rainfall is associated with periods of a warm tropical North Atlantic Ocean and large SST gradients between the Atlantic and Pacific Oceans. Rainfall variability is likely linked to the intensity and/or latitudinal position of the intertropical convergence zone (ITCZ). Periods of higher rainfall in Costa Rica are also associated with an enhanced sea surface salinity gradient on either side of the isthmus, suggesting greater freshwater export from the Atlantic Basin when the ITCZ is stronger and/or in a more northerly position. Further, wet periods in Central America coincide with high deuterium excess values in Greenland ice, suggesting a direct link between low latitude SSTs, tropical rainfall, and moisture delivery to Greenland. Our results indicate that a stronger tropical hydrological cycle during warm periods and large inter-ocean SST gradients enhanced the delivery of low latitude moisture to Greenland.     

Calibrating a glaciological model of the Greenland ice sheet from the Last Glacial Maximum to present-day using field observations of relative sea level and ice extent

We constrain a three-dimensional thermomechanical model of Greenland ice sheet (GrIS) evolution from the Last Glacial Maximum (LGM, 21 ka BP) to the present-day using, primarily, observations of relative sea level (RSL) as well as field data on past ice extent. Our new model (Huy2) fits a majority of the observations and is characterised by a number of key features: (i) the ice sheet had an excess volume (relative to present) of 4.1 m ice-equivalent sea level at the LGM, which increased to reach a maximum value of 4.6 m at 16.5 ka BP; (ii) retreat from the continental shelf was not continuous around the entire margin, as there was a Younger Dryas readvance in some areas. The final episode of marine retreat was rapid and relatively late (c. 12 ka BP), leaving the ice sheet land based by 10 ka BP; (iii) in response to the Holocene Thermal Maximum (HTM) the ice margin retreated behind its present-day position by up to 80 km in the southwest, 20 km in the south and 80 km in a small area of the northeast. As a result of this retreat the modelled ice sheet reaches a minimum extent between 5 and 4 ka BP, which corresponds to a deficit volume (relative to present) of 0.17 m ice-equivalent sea level. Our results suggest that remaining discrepancies between the model and the observations are likely associated with non-Greenland ice load, differences between modelled and observed present-day ice elevation around the margin, lateral variations in Earth structure and/or the pattern of ice margin retreat.     

Predicting suitability of cockle Austrovenus stutchburyi restoration sites using hydrodynamic models of larval dispersal

An important aspect of shellfish restoration projects is to evaluate whether potential sites are likely to be recolonised naturally once disturbances are removed and habitats are restored. Similarly, if active translocation or reseeding of juveniles or adults is undertaken, it is important to understand whether these introduced populations will be self‐maintaining in future, and whether any seed they produce are likely to be retained within the restored area. We used a combined hydrodynamic and particle tracking model to predict larval dispersion patterns for the common cockle Austrovenus stutchburyi under different hydrodynamic scenarios for seven release locations in Whangarei Harbour, New Zealand. Our results implied that sites varied substantially in their potential for self‐seeding and for exporting seed to other locations. For sites with more restricted dispersal, the model predicted that most larvae originating at these sites would settle inside the release region (68–94% for passive particle simulations), whereas relatively few larvae originating from the other release sites settled at these sites. In contrast, model larvae released from sites exhibiting high connectivity dispersed to all sub‐regions in the harbour, and export outside of the model region was high. Forthcoming field validation of these model predictions will result in better integration of hydrodynamic connectivity in whole estuary restoration programs.