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).
Appendix 2Table A2.1
Current and past insular area and distance to the mainland.