Pleistocene terrace formation, Quaternary rock uplift rates and geodynamics of the Hellenic Subduction Zone revealed from dating of paleoshorelines on Crete, Greece

Quaternary paleoshorelines are common landforms on the island of Crete, a forearc high above the Hellenic Subduction Zone. These geomorphic markers are useful on Crete and elsewhere in determining coastal uplift rates, the identification of active geologic structures, and to constrain geodynamic models and seismic hazards. Controversy exists in the literature regarding the formation mechanisms and age of late Pleistocene paleoshorelines on Crete that has led to competing models of the uplift history, tectonic evolution, and seismic hazards of the Hellenic forearc. We present new mapping and results from luminescence and radiocarbon geochronology of paleoshoreline deposits that constrain the spatial and temporal pattern of rock uplift around the Cretan coastline. Existing and new radiocarbon data are variable and show no obvious age-elevation trends within individual terrace sequences. By contrast, nearly all luminescence ages, some from shorelines dated with radiocarbon, show positive age-elevation trends and range from 60–220 ka suggesting that all dated paleoshorelines are beyond the limits of radiocarbon. We propose that the inconsistencies between the different geochronological methods are the result of secondary contamination of young carbonate, possibly from meteoric waters, that bias radiocarbon in Cretan Pleistocene marine fossils. Most luminescence ages closely correlate with the timing of mid-to-late Pleistocene relative sea level highstands, consistent with stratigraphic observations. Calculated coastal uplift rates using a Monte-Carlo error analysis range from ∼0–1.2 mm/yr; the lowest uplift rates are found along the northern and eastern coasts of the island, while the most rapid are focused along the southern and western coasts where active normal faults are observed offsetting paleoshoreline sequences. Based on this new data, we favor a tectonic model where slip along upper crustal normal faults acts to locally augment a steady regional signal of uplift along the south and west coast, interpreted to result from the deep underplating of rock at the base of the subduction wedge beneath Crete. Arcward of the contact between the upper plate Moho and the top of the subducting slab, crustal thinning will occur in the orogenic wedge resulting in subsidence along the north coast of Crete