Data-model comparison of Holocene sea-level change in the circum-Caribbean region

Relative sea-level (RSL) reconstructions from the circum-Caribbean region were interpreted using a glacial isostatic adjustment (CIA) model with the aims to quantify the contribution of this process to both the temporal and spatial forms of the RSL observations and remove the GIA signal to estimate land ice volume change (eustasy) during the mid-to-late Holocene. To infer an optimal GIA parameter set, the RSL data were used to determine best-fitting Earth viscosity model parameters for two different global ice histories. The RSL data indicate a clear preference (95% confidence) for relatively high viscosity values: >0.8 x 10(21) Pas in the upper mantle and >3 x 10(22) Pas in the lower mantle. The data were not able to discriminate (at 95% confidence) between lithospheric thickness values ranging between 71 and 120 km, although the thickest value considered produced the best fits. RSL predictions based on the best-fitting model parameters indicate a spatial variability across the region of up to similar to 7 m during the early to mid-Holocene which is large enough to introduce significant error when using the entire dataset to produce a single regional RSL curve without correcting for CIA. When corrected for GIA, the most precise data from the region indicate about 3-4 m of land ice melt (eustasy) from similar to 7 cal kyr B.P. to similar to 3-2 cal kyr B.P., although we note that there is considerable scatter in these corrected data. This likely reflects CIA model limitations, such as the assumption of lateral homogeneity in a region that contains several plate boundaries, the influence of tectonic processes (which were not modelled), as well as errors and an underestimate of the uncertainty in the RSL reconstructions. (C) 2013 Elsevier B.V. All rights reserved.