Appendix 1
In this earlier work, the Aegean archipelago was reconstructed based on
a geophysical model of relative sea-level change that uses generalized
sea-level equations accounting for hydro isostatic adjustments and
applied on a topographic and bathymetric grid with a resolution of 30
arc-seconds (Simaiakis et al. 2017). For the median value we took a
global eustatic median sea-level stand of 65 m below present for last
nine glacial-interglacial cycles (Norder et al. 2019). We assumed that
the regional geophysical effects affecting for this eustatic median
sea-level stand were minimal and comparable by the post glacial
sea-level rise conditions at 11 ka BP when this sea level was achieved
lastly. Islands smaller than the surface of a grid cell
(~1 km2) were excluded from the
analysis, and islands separated by a distance smaller than the grid size
(~1 km2) were aggregated in the
paleogeographic reconstruction to obtain a conservative estimate of the
fragmentation dynamics occurring in the system (see Rijsdijk et al.
2014). All geographical data were processed in ArcGIS 10.2.2. Past and
present areas and distance to the mainland were used to calculate the
magnitude of geographical change caused by sea-level rise (difference
between current and past values). For land-bridge islands (distance = 0
at Median sea level or Late Glacial Maximum), the area considered is
that of the island at the time step (1 Kya) preceding the connection to
the mainland (see Appendix S1 in Kougioumoutzis et al. 2017).
Current and past insular area and distance to the mainland.