Identification of Potential Paleoislands in the Mediterranean Sea During the Last Glacial Cycle

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.

Here we demonstrate the applicability of using altimetry data and Landsat imagery to provide the most accurate digital elevation model (DEM) of Australia's largest playa lake - Lake Eyre. We demonstrate through the use of geospatial techniques a robust assessment of lake area and volume of recent lake-filling episodes whilst also providing the most accurate estimates of area and volume for larger lake filling episodes that occurred throughout the last glacial cycle. We highlight that at a depth of 25 m Lake Mega-Eyre would merge with the adjacent Lake Mega-Frome to form an immense waterbody with a combined area of almost 35,000 km**2 and a combined volume of ~520 km**3. This would represent a vast water body in what is now the arid interior of the Australian continent. The improved DEM is more reliable from a geomorphological and hydrological perspective and allows a more accurate assessment of water balance under the modern hydrological regime. The results presented using GLAS/ICESat data suggest that earlier historical soundings were correct and the actual lowest topographic point in Australia is -15.6 m below sea level. The results also contrast nicely the different basin characteristics of two adjacent lake systems; Lake Eyre and Lake Frome.

We have developed a web-based processing system that can simulate positive and negative sea level changes globally by selecting the best Digital Elevation Model (DEM) for a target region from multiple DEMs. A PNG elevation tile format is used as the DEM format, which reduces the DEM data size. The PNG tile format implements client-based processing, and the DEM data are provided from different websites. In addition, the smart tile architecture is adopted, which enables on-demand simulation by adding a tile conversion process (i.e., a DEM selection process) during image drawing by using JavaScript. To demonstrate the system, we have employed three DEMs, i.e., the Geospatial Information Authority of Japan (GSI) map (~10-m resolution), the ASTER Global Digital Elevation Models (ASTER GDEM version 3) as global land area (~30-m resolution), and the General Bathymetric Chart of the Oceans as bathymetric data (~1000-m resolution). The ASTER Global Water Bodies Database is also used in the data selection process. The GSI provides their DEM in a PNG elevation tile format, and the other data are provided by the Geological Survey of Japan in PNG elevation tile format. We assume the current DEM sea level as 0 m, and the sea level can be changed to an arbitrary integer value (−10,000 to 10,000 m). Combining ASTER GDEM for land and GEBCO for sea makes it possible to target DEM of the whole earth. Moreover, it was shown that if a higher resolution DEM is available, it is possible to combine the higher resolution DEM in that area. The combining the PNG elevation tile format with the smart tile architecture demonstrates the possibilities of a client-based web processing service like that of the server-based OGC Web Processing Service.

Based on the glacial to postglacial delta13C differences between endobenthic Uvigerina peregrina species from the Alboran basin and from other mediterranean basins, changes in the fertility of the western part of this basin during the last deglaciation are reconstructed. As a result of particulate organic carbon (POC) rain from the highly productive upwelling cell along the northwestern margin of the Alboran basin, U. peregrina is presently depleted by about 1.6per mil with respect to the measured delta13C values of bottom water SumCO2 and by about 0.9per mil with respect to specimens from other areas of the western Mediterranean or from the Gulf of Cadiz within the Mediterranean Outflow Water. The Uvigerina delta13C difference between the Alboran Sea and the Gulf of Cadiz (Delta delta13C), was close to 0per mil at the beginning of the last deglaciation and during the late glacial time. This suggests that highly fertile systems set in the Alboran Sea near 16 kyr B.P. Two rapid increases in the Delta delta13C offset are recorded near 15 kyr and 11 kyr B.P. Fluctuations around 1.1 to 1.2per mil occurred during the early Holocene, and a maximum was reached near 9 kyr B.P. After 4 kyr the Delta delta13C offset decreased to its present-day average value of 0.9per mil. Changes in the intensity of surficial production cannot account for all the observed fluctuations, especially in the early Holocene time. A strong decrease in the intermediate and deep water ventilation of the Alboran basin may have occurred near 8-9 kyr, in phase with the last stagnant phase in the eastern Mediterranean and the deposition of Sapropel S1. As a result, the redistribution and remineralization at depth of the produced organic matter was incomplete. The POC rain reaching the sediment was locally intensified and caused the lowering of the delta13C values of endobenthic foraminifers such as U. peregrina. The benthic 13C signal suggests that the difference between the Alboran Sea and the Gulf of Cadiz was at its maximum. At the same time, an important modification in the water masses structure may have occurred near 9-8 kyr B.P. The deepening of the permanent pycnocline probably related to a thicker Atlantic jet at a stage of high sea level stand is recorded by the replacement of the right coiling N. pachyderma dominance (coincident with a shallow pycnocline) by the G. inflata dominance (coincident with a deep pycnocline). Diatom abundances were strongly reduced indicating an important modification of the productive system. The glacial-postglacial evolution of productivity within the Alboran basin was therefore more complex than in the adjacent Atlantic Ocean and opposite to the global one which displays a general increase in productivity during glacial time. Although it is the global budget of paleoproductivity that would drive the partitioning of carbon within the ocean, local or regional discrepancies with the global glacial-interglacial model must be addressed. Local winds and regional atmospheric pressure systems, which are the forcing factors for circulation and exchange with the Atlantic, control the fertile systems of the Alboran basin.

Many beaches in Penang island were severely inundated by the 26 December 2004 Indian Ocean mega tsunami with 57 deaths recorded. It is anticipated that the next big tsunami will cause even more damages to beaches in Penang. Hence, developing community resilience against the risks of the next tsunami is essential. Resilience entails many interlinked components, beginning with a good understanding of the inundation scenarios critical to community evacuation and resilience preparation. Inundation scenarios are developed from tsunami simulations involving all three phases of tsunami generation, propagation and run-up. Accurate and high-resolution bathymetric–topographic maps are essential for simulations of tsunami wave inundation along beaches. Bathymetric maps contain information on the depths of landforms below sea level while topographic maps reveal the elevation of landforms above sea level. Bathymetric and topographic datasets for Malaysia are, however, currently not integrated and are available separately and in different formats, not suitable for inundation simulations. Bathymetric data are controlled by the National Hydrographic Centre (NHC) of the Royal Malaysian Navy while topographic data are serviced by the Department of Survey and Mapping Malaysia (JUPEM). It is highly desirable to have seamless integration of high-resolution bathymetric and topographic data for tsunami simulations and for other scientific studies. In this paper, we develop a robust method for integrating the NHC bathymetric and JUPEM topographic data into a regularly-spaced grid system essential for tsunami simulation. A primary objective of this paper is to develop the best Digital Elevation and Bathymetry Model (DEBM) for Penang based upon the most suitable and accurate interpolation method for integrating bathymetric and topographic data with minimal interpolation errors. We analyze four commonly used interpolation methods for generating gridded topographic and bathymetric surfaces, namely (i) Kriging, (ii) Multiquadric (MQ), (iii) Thin Plate Spline (TPS) and (iv) Inverse Distance to Power (IDP). The study illustrated that the Kriging interpolation method produces an integrated bathymetric and topographic surface that best approximates the admiralty nautical chart of Penang essential for tsunami run-up and inundation simulations. Tsunami inundation scenarios critical to risk analysis and mitigation could then be developed using this DEBM for various earthquake scenarios, as presented in this paper for the 2004 Indian Ocean Tsunami.

The scarcity of high-resolution urban digital elevation model (DEM) datasets, particularly in certain developing countries, has posed a challenge for many water-related applications such as flood risk management. A solution to address this is to develop effective approaches to reconstruct high-resolution DEMs from their low-resolution equivalents that are more widely available. However, the current high-resolution DEM reconstruction approaches mainly focus on natural topography. Few attempts have been made for urban topography, which is typically an integration of complex artificial and natural features. This study proposed a novel multi-scale mapping approach based on convolutional neural network (CNN) to deal with the complex features of urban topography and to reconstruct high-resolution urban DEMs. The proposed multi-scale CNN model was firstly trained using urban DEMs that contained topographic features at different resolutions, and then used to reconstruct the urban DEM at a specified (high) resolution from a low-resolution equivalent. A two-level accuracy assessment approach was also designed to evaluate the performance of the proposed urban DEM reconstruction method, in terms of numerical accuracy and morphological accuracy. The proposed DEM reconstruction approach was applied to a 121 km2 urbanized area in London, United Kingdom. Compared with other commonly used methods, the current CNN-based approach produced superior results, providing a cost-effective innovative method to acquire high-resolution DEMs in other data-scarce regions.

Geomorphological and geological evidence for former Quaternary glaciation has been mapped in the Pindus Mountains of northwest Greece. The dynamics and chronology of glaciation in this area has been established through sedimentological analysis, soil analysis and Uranium-series dating. Four glacial events are recorded in the sedimentological and geomorphological records. The most extensive recorded glaciation pre-dates 350,000 years BP and was characterised by extensive valley glaciers and ice-fields. A second glaciation occurred prior to the last interglacial, before ca. 127,000 years BP, and was characterised by glaciers that reached mid-valley positions. The height of the last glacial stage in Greece (30-20,000 14C years BP) is recorded by small cirque glacier moraines and relict periglacial rock glaciers. Evidence for a fourth glacial phase is recorded only in the highest cirques of Mount Smolikas (2637 m a.s.l.), the highest peak in the Pindus Mountains. This phase of glaciation is likely to have occurred during the Late-glacial Substage (14-10,000 14C years BP). All of the glaciers during the different glacial stages were reconstructed and used alongside periglacial rock glaciers to determine palaeoclimate. During the glacial maximum of the last glacial stage, mean annual temperatures were ca. 8-9 degrees C lower than at present and mean annual precipitation greater than 2000 mm - similar to modern values. Earlier glacial maxima are likely to have been colder but with mean annual precipitation still greater than 2000 mm. Maximum glacier extent in the Pindus Mountains is likely to have preceded the most severe arid phase of glacial cycles indicated in the pollen record and also global glacial maxima. This was because of the small size of the former Pindus glaciers and their rapid response to climate change, as well as the increased prevalence of aridity around the global glacial maxima. The glacial sequence in the Pindus Mountains represents the best-dated and longest recognised record of glaciation in the Mediterranean region and provides a stratigraphical framework for Quaternary cold-stage climates in Greece.

Anatolia is situated in the Eastern Mediterranean region between 36 – 42N and 26 – 45E. The geological records of paleoglaciations in the high terrains of Anatolia are key archives to quantify paleoclimate change in the Eastern Mediterranean area. The climate of the Eastern Mediterranean region is influenced by three main atmospheric systems: the main middle to high latitude westerlies, the mid-latitude subtropical high-pressure systems, and the monsoon climate. Glacial geological studies in Turkey have started in the late 19th century. Glacial deposits are found mainly in the eastern, northeastern and southern part of the Anatolian Peninsula. Anatolia is the fundamental element to understand the interactions between paleoenvironment, climatic variations, and development of the human societies. As the Taurus and Black Sea Mountains are sensitively situated for the paleoclimatic reconstructions, a chronostratigraphic framework on the paleoglaciation should be elaborated. The timing of the Last Glacial Maximum (LGM) in Anatolia is still unknown. Our first results from Kavron Valley (Kackar Mountains, NE Turkey) are encouraging for the reconstruction of paleoglaciations in Turkey and related paleoclimatological interpretations although it is presently difficult to pinpoint the classical Last Glacial Maximum – Younger Dryas – Little Ice Age moraine sequences in the field.

Climate change and the associated sea level rise (SLR) are presenting newfound challenges to the port systems and coastal transportation infrastructure of southeast Texas. This paper introduces a Geographic Information Systems (GIS) based model designed to simulate inundation scenarios under various sea-level projections, aiming to assess the vulnerabilities of both port facilities and road networks. The study area encompasses a specific region within Jefferson County, southeast Texas, encompassing three major ports: Port Arthur, Beaumont, and Orange. Utilizing a high-resolution (1-m) Digital Elevation Model (DEM) derived from the 2017 LiDAR dataset, this model is integrated with NASA’s sea-level rise projections to compute the extent and volume of inundation across low, medium, and high SLR scenarios. Drawing from monthly mean sea level data spanning from 1958 to 2020, the lowest SLR projections, derived from the relative sea-level trend measured at the Sabine Pass, TX gauge station, indicate a yearly increase of 6.16 mm, with a 95% confidence interval of +/- 0.74 mm. Projections for 2050 and 2,100 show the lowest SLR at 0.17 m and 0.48 m, respectively. In contrast, the medium to high RSLR projections under the IPCC SSP3-7.0 scenario for 2050 and 2,100 stand at 0.54 m and 1.34 m, respectively. The findings reveal that, under medium to high SLR scenarios, the extent of inundated areas in the study region is expected to expand by 12.4% in 2050 and 19.9% in 2,100, compared to the lowest SLR projection. Additionally, the length of submerged roadways is predicted to increase by 6.9% in 2050 and 13.3% in 2,100, in comparison to the lowest SLR projection. It is worth noting that some margin of error may be introduced due to factors such as the width of the port area and access roads, the high-resolution DEM, and the alignment of computed inundated areas with the existing topography. Overall, the manuscript highlights the urgency of proactive planning and underscores the importance of safeguarding critical infrastructure in the context of climate change and SLR.

This paper presents the results of acoustic surveys and video imaging of the seabed off Waterford Harbour, beyond a previously identified submerged palaeochannel. This system extends south-west into waters at c. -56 m OD and terminates in an area of possible glacigenic sediments. The aim of the new surveys was to obtain multibeam sonar imagery and to correlate it with seismic profiles in order to understand the genetic relationships between the morpho-sedimentological elements of the seabed in this area. The multibeam sonar imagery revealed four arcuate, morphologically complex features that on-lap the regional bedrock and confine the previously identified palaeochannel. Both the morphology of these forms and their sedimentology, as revealed by video imaging, imply formation in an ice-marginal environment. These forms are hypothesised to reflect several stages of sedimentation at the margins of ice progressively retreating from the nearshore shelf of the Celtic Sea during the Last Glacial Maximum.

The spatial variability of input parameters plays a crucial role in the interpretation of geomorphic indices, with digital elevation models (DEMs) being the primary data source. However, the influence of DEM resolution on these indices has rarely been investigated. This study investigated the influence of DEM resolution on the assessment of tectonic activity using the normalized stream length–gradient (SLk) index, which reflects variations along river profiles. The SLk index is sensitive to changes in river gradients that may result from active faulting or differential uplift, making it a valuable tool for identifying zones of active tectonic deformation. Therefore, understanding the impact of DEM resolution on SLk analysis is critical for accurately detecting and interpreting subtle tectonic signals, particularly in intraplate regions where deformation is slow and geomorphic expressions are faint and discontinuous. By comparing high-resolution LiDAR-derived DEMs (L-DEMs) and low-resolution topographic map-derived DEMs (T-DEMs), we analyzed the SLk index distributions along the Yangsan Fault, Korean Peninsula, an intraplate setting with Quaternary activity. According to the results, SLk anomalies derived from L-DEMs had a continuous distribution along the fault, closely aligning with known surface ruptures and indicating active tectonic deformation. In contrast, SLk anomalies derived from T-DEMs were sporadic and less continuous, especially in low-relief landscapes such as alluvial fans and floodplains, highlighting the limitations of T-DEMs in detecting fault-related features. High-resolution DEMs were better able to capture finer-scale geomorphic features, such as fault scarps, deflected streams, and lineaments associated with active tectonics, providing a more comprehensive view of fault-related deformation. This discrepancy highlights the importance of resolution choice in tectonic assessments, as low-resolution DEMs may underestimate the tectonic activities of intraplate faults by missing subtle topographic variations. While the choice of DEM resolution may depend on study area, scope, and data availability, high-resolution DEMs are critical for identifying tectonic activity in intraplate regions where geomorphic features of faulting due to slow deformation are subtle and dispersed.

Glacial-interglacial cycles, occurring at a period of approximately 100,000 years, have dominated Earth's climate over the past 800,000 years. These cycles involve major changes in land ice, global sea level, ocean circulation, and the carbon cycle. While it is generally agreed that the ultimate driver of global climate is changes in insolation, glacial cycles do not look like insolation forcing. Notably, there is a highly non-linear warming response at 100,000 years to a relatively small forcing, implicating a more complicated system of biogeochemical and physical drivers. The ocean plays a pivotal role in glacial-interglacial climate through direct equator-to-pole transport of heat and its role in the carbon cycle. The deep ocean contains 60 times more carbon than the atmosphere, and therefore even small changes in ocean circulation can have a large impact on atmospheric CO2, a crucial amplifier in the climate system. In order to better understand the role that ocean circulation plays in glacial-interglacial climate we focus on the last glacial-interglacial transition. In this thesis, we present reconstructions of changes in intermediate water circulation and explore a new time-dependent dynamical box model. We reconstruct circulation using radiocarbon and clumped isotope measurements on U/Th dated deep-sea corals from the New England and Corner Rise Seamounts in the western basin of the North Atlantic and from south of Tasmania in the Indo-Pacific sector of the Southern Ocean. Our new time-dependent model contains key aspects of ocean physics, including Southern Ocean Residual Mean theory, and allows us to explore dynamical mechanisms which drive abrupt climate transitions during the last glacial period. In Chapter 2 we present a compilation of reconnaissance dated deep-sea corals from the Caltech collection. Reconnaissance dating facilitates sample selection for our high-precision radiocarbon and temperature time series and patterns in the depth distribution of deep-sea corals over time contain additional relevant climate information. In Chapter 3, we present a high-resolution radiocarbon record from south of Tasmania which highlights variability in Southern Ocean Intermediate Water radiocarbon during the deglaciation, particularly during the Antarctic Cold Reversal. We use our radiocarbon data, in combination with other deglacial climate records, to infer changes in overturning circulation configuration across this time interval. In Chapter 4 we present our time-dependent dynamical box model. Our model displays hysteresis in basin stratification and Southern Ocean isopycnal outcrop position as a function of North Atlantic Deep Water formation rate. In a dynamical system, hysteresis implies that there are multiple stable states, and switches between these states can lead to abrupt transitions, such as those observed during the middle of the last glacial period. In Chapter 5 we present paired radiocarbon and temperature time series from the North Atlantic and Southern Ocean spanning the late part of the last glacial. We explore the mechanisms driving trends in radiocarbon and temperature by looking at cross-plots of the data, and we make inferences about changes in circulation configuration using insight gained from our dynamical box model.

Making post-glacial uplift visible: A model based high-resolution animation of shore displacement

<p>Tsunami simulations using high resolution datasets would always resemble the results that are closer to the reality. However, high resolution airborne or spaceborne local datasets have not yet been available for many regions and acquisition of this data is costly or might not even be possible for some locations. This hard-to-reach situation of high resolution datasets obliged researchers to work with open source datasets in their studies, which forces them to cope with the uncertainties of low spatial resolution.</p><p>Tsunami numerical models require both bathymetric and topographic data in order to calculate wave propagation in the water and inundation on the land. Leaving aside the availability of reliable bathymetric data, there are different open source global Digital Elevation Model (DEM) datasets, which are freely available. ASTER GDEM, SRTM and ALOS World 3D are present global open source DEMs that have highest spatial resolution of 30 meters. These three different sources of DEMs are generated using different technologies during data acquisition and different methodologies while processing. Even if they are the best available open source datasets, they all include variable sources of differences and errors.</p><p>This study aims evaluation of the sensitivity of open source DEM datasets against high resolution DEM datasets for tsunami hazard assessment and examination of accuracy of the simulations’ results. A small area in Silivri district of Istanbul is selected as study area, where 1m resolution of topographic data is available. Tsunami simulations are performed using NAMI DANCE GPU with topography data of 1m resolution based on LiDAR measurements and topography data of 30m resolution based on ASTER GDEM, SRTM and ALOS World 3D datasets. The resulted inundation on land and flow depth distributions are plotted and discussed with comparisons.</p><p> </p><p>Acknowledgement: MSc. Bora Yalciner and Assoc. Prof. Dr. Andrey Zaytsev are acknowledged for their contributions in developing tsunami numerical model NAMI DANCE GPU used in this study. The authors also thank Istanbul Metropolitan Municipality, Directorate of Earthquake and Ground Investigation for providing high quality data and close cooperation.</p><p>Keywords: tsunami, hazard assessment, numerical modeling, open source DEM, high-resolution DEM</p>

New York City Panel on Climate Change 2019 Report Chapter 5: Mapping Climate Risk