Holocene record of Tuggerah Lake estuary development on the Australian east coast: sedimentary responses to sea-level fluctuations and climate variability

We investigated the Holocene palaeo-environmental record of the Tuggerah Lake barrier estuary on the south-east coast of Australia to determine the influence of local, regional and global environmental changes on estuary development. Using multi-proxy approaches, we identified significant down-core variation in sediment cores relating to sea-level rise and regional climate change. Following erosion of the antecedent land surface during the post-glacial marine transgression, sediment began to accumulate at the more seaward location at ~8500. years before present, some 1500. years prior to barrier emplacement and ~4000. years earlier than at the landward site. The delay in sediment accumulation at the landward site was a consequence of exposure to wave action prior to barrier emplacement, and due to high river flows of the mid-Holocene post-barrier emplacement. As a consequence of the mid-Holocene reduction in river flows, coupled with a moderate decline in sea-level, the lake experienced major changes in conditions at ~4000. years before present. The entrance channel connecting the lake with the ocean became periodically constricted, producing cyclic alternation between intervals of fluvial- and marine-dominated conditions. Overall, this study provides a detailed, multi-proxy investigation of the physical evolution of Tuggerah Lake with causative environmental processes that have influenced development of the estuary.

The stratigraphic record of the quaternary sea level fluctuations and the impact of the post-glacial sea level rise (Termination I) in the Adriatic basin (Mediterranean sea)

The modern stratigraphy of clastic continental margins is the result of the interaction between several geological processes acting on different time scales, among which sea level oscillations, sediment supply fluctuations and local tectonics are the main mechanisms. During the past three years my PhD was focused on understanding the impact of each of these process in the deposition of the central and northern Adriatic sedimentary successions, with the aim of reconstructing and quantifying the Late Quaternary eustatic fluctuations. In the last few decades, several Authors tried to quantify past eustatic fluctuations through the analysis of direct sea level indicators, among which drowned barrier-island deposits or coral reefs, or indirect methods, such as Oxygen isotope ratios (δ18O) or modeling simulations. Sea level curves, obtained from direct sea level indicators, record a composite signal, formed by the contribution of the global eustatic change and regional factors, as tectonic processes or glacial-isostatic rebound effects: the eustatic signal has to be obtained by removing the contribution of these other mechanisms. To obtain the most realistic sea level reconstructions it is important to quantify the tectonic regime of the central Adriatic margin. This result has been achieved integrating a numerical approach with the analysis of high-resolution seismic profiles. In detail, the subsidence trend obtained from the geohistory analysis and the backstripping of the borehole PRAD1.2 (the borehole PRAD1.2 is a 71 m continuous borehole drilled in -185 m of water depth, south of the Mid Adriatic Deep - MAD - during the European Project PROMESS 1, Profile Across Mediterranean Sedimentary Systems, Part 1), has been confirmed by the analysis of lowstand paleoshorelines and by benthic foraminifera associations investigated through the borehole. This work showed an evolution from inner-shelf environment, during Marine Isotopic Stage (MIS) 10, to upper-slope conditions, during MIS 2. Once the tectonic regime of the central Adriatic margin has been constrained, it is possible to investigate the impact of sea level and sediment supply fluctuations on the deposition of the Late Pleistocene-Holocene transgressive deposits. The Adriatic transgressive record (TST - Transgressive Systems Tract) is formed by three correlative sedimentary bodies, deposited in less then 14 kyr since the Last Glacial Maximum (LGM); in particular: along the central Adriatic shelf and in the adjacent slope basin the TST is formed by marine units, while along the northern Adriatic shelf the TST is represented by costal deposits in a backstepping configuration. The central Adriatic margin, characterized by a thick transgressive sedimentary succession, is the ideal site to investigate the impact of late Pleistocene climatic and eustatic fluctuations, among which Meltwater Pulses 1A and 1B and the Younger Dryas cold event. The central Adriatic TST is formed by a tripartite deposit bounded by two regional unconformities. In particular, the middle TST unit includes two prograding wedges, deposited in the interval between the two Meltwater Pulse events, as highlighted by several 14C age estimates, and likely recorded the Younger Dryas cold interval. Modeling simulations, obtained with the two coupled models HydroTrend 3.0 and 2D-Sedflux 1.0C (developed by the Community Surface Dynamics Modeling System - CSDMS), integrated by the analysis of high resolution seismic profiles and core samples, indicate that: 1 - the prograding middle TST unit, deposited during the Younger Dryas, was formed as a consequence of an increase in sediment flux, likely connected to a decline in vegetation cover in the catchment area due to the establishment of sub glacial arid conditions; 2 - the two-stage prograding geometry was the consequence of a sea level still-stand (or possibly a fall) during the Younger Dryas event. The northern Adriatic margin, characterized by a broad and gentle shelf (350 km wide with a low angle plunge of 0.02° to the SE), is the ideal site to quantify the timing of each steps of the post LGM sea level rise. The modern shelf is characterized by sandy deposits of barrier-island systems in a backstepping configuration, showing younger ages at progressively shallower depths, which recorded the step-wise nature of the last sea level rise. The age-depth model, obtained by dated samples of basal peat layers, is in good agreement with previous published sea level curves, and highlights the post-glacial eustatic trend. The interval corresponding to the Younger Dyas cold reversal, instead, is more complex: two coeval coastal deposits characterize the northern Adriatic shelf at very different water depths. Several explanations and different models can be attempted to explain this conundrum, but the problem remains still unsolved.

Middle to late Miocene oxygen isotope stratigraphy of ODP site 1085 (SE Atlantic): new constrains on Miocene climate variability and sea-level fluctuations

The middle Miocene delta18O increase represents a fundamental change in earth's climate system due to a major expansion and permanent establishment of the East Antarctic Ice Sheet accompanied by some effect of deepwater cooling. The long-term cooling trend in the middle to late Miocene was superimposed by several punctuated periods of glaciations (Mi-Events) characterized by oxygen isotopic shifts that have been related to the waxing and waning of the Antarctic ice-sheet and bottom water cooling. Here, we present a high-resolution benthic stable oxygen isotope record from ODP Site 1085 located at the southwestern African continental margin that provides a detailed chronology for the middle to late Miocene (13.9-7.3 Ma) climate transition in the eastern South Atlantic. A composite Fe intensity record obtained by XRF core scanning ODP Sites 1085 and 1087 was used to construct an astronomically calibrated chronology based on orbital tuning. The oxygen isotope data exhibit four distinct delta18O excursions, which have astronomical ages of 13.8, 13.2, 11.7, and 10.4 Ma and correspond to the Mi3, Mi4, Mi5, and Mi6 events. A global climate record was extracted from the oxygen isotopic composition. Both long- and short-term variabilities in the climate record are discussed in terms of sea-level and deep-water temperature changes. The oxygen isotope data support a causal link between sequence boundaries traced from the shelf and glacioeustatic changes due to ice-sheet growth. Spectral analysis of the benthic delta18O record shows strong power in the 400-kyr and 100-kyr bands documenting a paleoceanographic response to eccentricity-modulated variations in precession. A spectral peak around 180-kyr might be related to the asymmetry of the obliquity cycle indicating that the response of the dominantly unipolar Antarctic ice-sheet to obliquityinduced variations probably controlled the middle to late Miocene climate system. Maxima in the delta18O record, interpreted as glacial periods, correspond to minima in 100-kyr eccentricity cycle and minima in the 174-kyr obliquity modulation. Strong middle to late Miocene glacial events are associated with 400-kyr eccentricity minima and obliquity modulation minima. Thus, fluctuations in the amplitude of obliquity and eccentricity seem to be the driving force for the middle to late Miocene climate variability.

Granulometric characterization and evaluation of annually banded mid-Holocene estuarine silts, Welsh Severn Estuary (UK): coastal change, sea level and climate

Holocene silts (salt marshes) and highest intertidal-supratidal peats are superbly exposed on a 15 kin coastal transect which reveals two laterally extensive units of annually banded silts (Beds 3, 7) associated with three transgressive-regressive silt-peat cycles (early sixth-early fourth millennium BC). Bed 3 in places is concordantly and gradationally related to peats above and below, but in others transgresses older strata. Bed 7 also grades up into peat, but everywhere overlies a discordance. The banding in Bed 3 at three main and two minor sites was resolved and characterized texturally at high-resolution (2.5/5 mm contiguous slices) using laser granulometry (LS230 with PIDS) and a comprehensive scheme of data-assessment. Most of Bed 3 formed very rapidly, at peak values of several tens of millimetres annually, in accordance with modelled effects of sea-level fluctuations on mature marshes (bed concordant and gradational) and on marshes growing up after coastal erosion and retreat (bed with discordant base). Using data from the modern Severn Estuary, the textural contrast within bands, and its variation between bands, points to a variable but overall milder mid-Holocene climate than today. The inter-annual variability affected marsh dynamics, as shown by the behaviour of the finely divided plant tissues present. Given local calibration, the methodology is applicable to other tidal systems with banded silts in Britain and mainland northwest Europe. (c) 2006 Elsevier Ltd. All rights reserved.

Benthic foraminifera and sea level analyses from IODP Hole 310-M0005D

The course of sea-level fluctuations during Termination II (TII; the penultimate deglaciation), which is critical for understanding ice-sheet dynamics and suborbital climate variability, has yet to be established. This is partly because most shallow-water sequences encompassing TII were eroded during sea-level lowstands of the last glacial period or were deposited below the present sea level. Here we report a new sequence recording sea-level changes during TII in the Pleistocene sequence at Hole M0005D (water depth: 59.63 m below sea level [mbsl]) off Tahiti, French Polynesia, which was drilled during Integrated Ocean Drilling Program Expedition 310. Lithofacies variations and stratigraphic changes in the taxonomic composition, preservation states, and intraspecific test morphology of large benthic foraminifers indicate a deepening-upward sequence in the interval from Core 310-M0005D-26R (core depth: 134 mbsl) through -16R (core depth: 106 mbsl). Reconstruction of relative sea levels, based on paleodepth estimations using large benthic foraminifers, indicated a rise in sea level of about 90 m during this interval, suggesting its correlation with one of the terminations. Assuming that this rise in sea level corresponds to that during TII, after correcting for subsidence since the time of deposition, a highstand sea-level position would be 2 ± 15 m above present sea level (masl), which is generally consistent with highstand sea-level positions in MIS 5e (4 ± 2 masl). If this rise in sea level corresponds to that during older terminations, the subsidence-corrected highstand sea-level positions (30 ± 15 masl for Termination III and 54 ± 15 masl for Termination IV) are not consistent with reported ranges of interglacial sea-level highstands (-18 to 15 masl). Therefore, the studied interval likely records the rise in sea level and associated environmental changes during TII. In particular, the intervening cored materials between the two episodes of sea-level rise found in the studied interval might record the sea-level reversal event during TII. This conclusion is consistent with U/Th ages of around 133 ka, which were obtained from slightly diagenetically altered (i.e., < 1% calcite) in situ corals in the studied interval (Core 310-M0005D-20R [core depth: 118 mbsl]). This study also suggests that our inverse approach to correlate a stratigraphic interval with an approximate time frame could be useful as an independent check on the accuracy of uranium-series dating, which has been applied extensively to fossil corals in late Quaternary sea-level studies.

Orbital-resolution Pliocene-Pleistocene simulations of CO2, sea level, global temperature and benthic d18O, and d11B-based proxy-CO2 data

During the past five million yrs, benthic d18O records indicate a large range of climates, from warmer than today during the Pliocene Warm Period to considerably colder during glacials. Antarctic ice cores have revealed Pleistocene glacial-interglacial CO2 variability of 60-100 ppm, while sea level fluctuations of typically 125 m are documented by proxy data. However, in the pre-ice core period, CO2 and sea level proxy data are scarce and there is disagreement between different proxies and different records of the same proxy. This hampers comprehensive understanding of the long-term relations between CO2, sea level and climate. Here, we drive a coupled climate-ice sheet model over the past five million years, inversely forced by a stacked benthic d18O record. We obtain continuous simulations of benthic d18O, sea level and CO2 that are mutually consistent. Our model shows CO2 concentrations of 300 to 470 ppm during the Early Pliocene. Furthermore, we simulate strong CO2 variability during the Pliocene and Early Pleistocene. These features are broadly supported by existing and new d11B-based proxy CO2 data, but less by alkenone-based records. The simulated concentrations and variations therein are larger than expected from global mean temperature changes. Our findings thus suggest a smaller Earth System Sensitivity than previously thought. This is explained by a more restricted role of land ice variability in the Pliocene. The largest uncertainty in our simulation arises from the mass balance formulation of East Antarctica, which governs the variability in sea level, but only modestly affects the modeled CO2 concentrations.

Alpha-cellulose δ13C variation in mangrove tree rings correlates well with annual sea level trend between 1982 and 1999

A pilot study of tree rings in a modern mangrove tree (Rhizophora apiculata) from Leizhou Peninsula, northern South China Sea shows that ( 1) the tree-rings are annual; ( 2) the ring widths decrease; and ( 3) their alpha-cellulose delta(13)C values increase from 1982 to 1999 AD, consistent with the trends of annual sea level, salinity and sea surface temperatures in the same period. We propose that such changes were caused by increasingly longer duration of waterlogging in response to sea-level rise. If this is the case, alpha-cellulose delta(13)C in mangrove tree rings can be used as a potential indicator of past sea level fluctuations.

The seasonal variations of sea level due to density variations in the Arabian Gulf and Gulf of Oman

From the distribution of oceanographic data (temperature and salinity) in both Arabian Gulf and Gulf of Oman, the steric components (thermal, haline and steric heights) are calculated for the upper 50m layer during different seasons. The analysis reveals relevant evidence, that temperature variations (thermal component) play a role in the fluctuations of sea level within the investigated area. The salinity variations (haline component) is only significant near the entrance. The sea level variations due to density (steric component) is low during winter and spring and high during summer and autumn. The steric height is always lower in the northern and central regions of Arabian Gulf and higher in eastern region of Arabian Gulf and in the Gulf of Oman, i.e. the surface water must flow from the Gulf of Oman to the Arabian Gulf. The steric sea level gradient around the Strait of Hormuz are 0.04 cm/km in winter, 0.04 cm/km in spring, and 0.025 cm/km in summer and 0.014 cm/km in autumn.

Onset of recent rapid sea-level rise in the western Atlantic Ocean

A high-resolution record of sea-level change spanning the past 1000 years is derived from foraminiferal and chronological analyses of a 2m thick salt-marsh peat sequence at Chezzetcook, Nova Scotia, Canada. Former mean tide level positions are reconstructed with a precision of +/- 0.055 in using a transfer function derived from distributions of modern salt-marsh foraminifera. Our age model for the core section older than 300 years is based on 19 AMS C-14 ages and takes into account the individual probability distributions of calibrated radiocarbon ages. The past 300 years is dated by pollen and the isotopes Pb-206, Pb-207, Pb-210, Cs-137 and Am-241. Between AD 1000 and AD 1800, relative sea level rose at a mean rate of 17cm per century. Apparent pre-industrial rises of sea level dated at AD 1500-1550 and AD 1700-1800 cannot be clearly distinguished when radiocarbon age errors are taken into account. Furthermore, they may be an artefact of fluctuations in atmospheric C-14 production. In the 19th century sea level rose at a mean rate of 1.6mm/yr. Between AD 1900 and AD 1920, sea-level rise accelerated to the modern mean rate of 3.2mm/yr. This acceleration corresponds in time with global temperature rise and may therefore be associated with recent global warming. (c) 2005 Elsevier Ltd. All rights reserved.

Grain size and heavy minerals of the Late Quaternary eolian sediments from the Imbituba-Jaguaruna coast, Southern Brazil: Depositional controls linked to relative sea-level changes

The stratigraphic subdivision and correlation of dune deposits is difficult, especially when age datings are not available. A better understanding of the controls on texture and composition of eolian sands is necessary to interpret ancient eolian sediments. The Imbituba-Jaguaruna coastal zone (Southern Brazil, 28 degrees-29 degrees S) stands out due to its four well-preserved Late Pleistocene (eolian generation 1) to Holocene eolian units (eolian generations 2, 3, and 4). In this study, we evaluate the grain-size and heavy-mineral characteristics of the Imbituba-Jaguartma eolian units through statistical analysis of hundreds of sediment samples. Grain-size parameters and heavy-mineral content allow us to distinguish the Pleistocene from the Holocene units. The grain size displays a pattern of fining and better sorting from generation 1 (older) to 4 (younger), whereas the content of mechanically stable (dense and hard) heavy minerals decreases from eolian generation 1 to 4. The variation in grain size and heavy-mineral content records shifts in the origin and balance (input versus output) of eolian sediment supply attributable mainly to relative sea-level changes. Dunefields submitted to relative sea-level lowstand conditions (eolian generation 1) are characterized by lower accumulation rates and intense post-depositional dissection by fluvial incision. Low accumulation rates favor deflation in the eolian system, which promotes concentration of denser and stable heavy minerals (increase of ZTR index) as well as coarsening of eolian sands. Dissection involves the selective removal of finer sediments and less dense heavy minerals to the coastal source area. Under a high rate of relative sea-level rise and transgression (eolian generation 2), coastal erosion prevents deflation through high input of sediments to the coastal eolian source. This condition favors dunefield growth. Coastal erosion feeds sand from local sources to the eolian system. including sands from previous dunefields (eolian generation 1) and from drowned incised valleys. Therefore, dunefields corresponding to transgressive phases inherit the grain-size and heavy-mineral characteristics of previous dunefields, leading to selective enrichment of finer sands and lighter minerals. Eolian generations 3 and 4 developed during a regressive-progradational phase (Holocene relative sea level highstand). The high rate of sediment supply during the highstand phase prevents deflation. The lack of coastal erosion favors sediment supply from distal sources (fluvial sediments rich in unstable heavy minerals). Thus, dunefields of transgressive and highstand systems tracts may be distinguished from dunefields of the lowstand systems tract through high rates of accumulation (low deflation) in the former. The sediment source of the transgressive dunefields (high input of previously deposited coastal sands) differs from that of the highstand dunefields (high input of fluvial distal sands). Based on this case study, we propose a general framework for the relation between relative sea level, sediment supply and the texture and mineralogy of eolian sediments deposited in siliciclastic wet coastal zones similar to the Imbituba-Jaguaruna coast (C) 2009 Elsevier B.V. All rights reserved.

The last mangroves of Marajo Island - Eastern Amazon: Impact of climate and/or relative sea-level changes

The dynamics, over the last 7500 years, of a mangrove at Marajo Island in northern Brazil were studied by pollen and sedimentary facies analyses using sediment cores. This island, located at the mouth of the Amazon River. is influenced by riverine inflow combined with tidal fluctuations of the equatorial Atlantic Ocean. Herbaceous vegetation intermingled with rainforest dominates the central area of the island, while varzea is the main vegetation type along the littoral. In particular, the modem northeastern coastal zone is covered by a mosaic of dense rainforest, herbaceous vegetation, mangroves, varzea, and restinga. The integration of pollen data and fades descriptions indicates a tidal mud flat colonized by mangroves in the interior of Marajo Island between similar to 7500 cal yr BP and similar to 3200 cal yr BP. During the late Holocene, mangroves retracted to a small area (100-700 m in width) along the northeastern coastal plain. Mangrove expansion during the early and mid Holocene was likely caused by the post-glacial sea-level rise which, combined with tectonic subsidence, led to a rise in tidal water salinity. Salinity must have further increased due to low river discharge resulting from increased aridity during the early and mid Holocene. The shrinking of the area covered by mangrove vegetation during the late Holocene was likely caused by the increase in river discharge during the late Holocene, which has maintained relatively low tidal water salinity in Marajo Island. Tidal water salinity is relatively higher in the northeastern part of the island than in others, due to the southeast-northwest trending current along the littoral. The mixing of marine and riverine freshwater inflows has provided a refuge for mangroves in this area. The increase in flow energy during the last century is related to landward sand migration, which explains the current retraction of mangroves. These changes may indicate an increased exposure to tidal influence driven by the relative sea-level rise, either associated with global fluctuations or tectonic subsidence, and/or by an increase in river water discharge. (C) 2012 Elsevier B.V. All rights reserved.

Does sea-level rise have an impact on saltwater intrusion?

Climate change effects are expected to substantially raise the average sea level. It is widely assumed that this raise will have a severe adverse impact on saltwater intrusion processes in coastal aquifers. In this study we hypothesize that a natural mechanism, identified as the “lifting process” has the potential to mitigate or in some cases completely reverse the adverse intrusion effects induced by sea-level rise. A detailed numerical study using the MODFLOW-family computer code SEAWAT, was completed to test this hypothesis and to understand the effects of this lifting process in both confined and unconfined systems. Our conceptual simulation results show that if the ambient recharge remains constant, the sea-level rise will have no long-term impact (i.e., it will not affect the steady-state salt wedge) on confined aquifers. Our transient confined flow simulations show a self-reversal mechanism where the wedge which will initially intrude into the formation due to the sea-level rise would be naturally driven back to the original position. In unconfined systems, the lifting process would have a lesser influence due to changes in the value of effective transmissivity. A detailed sensitivity analysis was also completed to understand the sensitivity of this self-reversal effect to various aquifer parameters.