Historical and recent sea level rise and land subsidence in Marina di Ravenna, northern Italy ANNALS OF GEOPHYSICS, 59, 5, 2016, P0546; doi:10.4401/ag-7022 P0546 Historical and recent sea level rise and land subsidence in Marina di Ravenna, northern Italy Ines Cerenzia1,2, Davide Putero3,2, Flavio Bonsignore1, Gaia Galassi4, Marco Olivieri5,*, Giorgio Spada4,5 1 ARPAE Emilia-Romagna, Bologna, Italy 2 Università di Bologna, Dipartimento di Fisica e Astronomia (DIFA), Bologna, Italy 3 Istituto di Scienze dell’Atmosfera e del Clima (ISAC-CNR), Bologna, Italy 4 Università di Urbino, Dipartimento di Scienze Pure e Applicate (DiSPeA), Urbino, Italy 5 Istituto Nazionale di Geofisica e Vulcanologia (INGV), Sezione di Bologna, Italy ABSTRACT The regions facing the northern Adriatic Sea are particularly vulnerable to sea-level rise. Several trade ports are located there, and the area is im- portant from social and economical viewpoints. Since tourism and cul- tural heritage are a significant source of income, an increase in sea-level could hinder the development of these regions. One of the longest sea- level time series in the northern Adriatic, which goes back to the late 1880s, has been recorded at Marina di Ravenna, in Emilia-Romagna re- gion. The record is anomalous, showing a rate of increase that largely ex- ceeds that observed in nearby stations. During the last few decades, geodetic campaigns based on geometric high precision leveling, SAR in- terferometry, and GPS have monitored the Ravenna area. In this work, tide gauge observations are merged with yet unpublished geodetic data, aiming at a coherent interpretation of vertical land movements. We con- firm that land subsidence is the major cause of relative sea-level change at Marina di Ravenna, at least during the period allowing for a quan- titative analysis (1990-2011). The rate of absolute sea-level change (2.2±1.3 mm yr−1 during the same time period), given by the difference between the rate of relative sea-level change and the rate of subsidence, is consistent with the rate of absolute sea-level change observed by al- timetry in the northern Adriatic Sea. 1. Introduction Sea-level change is one of the key indicators of global warming. According to the latest report of the Intergovernmental Panel on Climate Change (IPCC), global sea level rose at a rate in the range between 1.3 and 1.7 mm yr−1 between 1901 and 1990 [Church et al. 2013], with an acceleration of ~0.01 mm yr−2 [Spada et al. 2015]. However, satellite data [Wunsch and Stammer 1995] and hydrographic observations [Roemmich and Owens 2000], in agreement with climate models [Spada et al. 2013], show that sea level has not been rising uni- formly across the oceans. In some regions, rates are even several times larger than the global mean sea-level rise, while in others sea level is falling [Cazenave and Llovel 2010]. Globally, sea-level rise is mainly driven by the thermal expansion of water masses and by the melting of glaciers and ice sheets [Lombard et al. 2005, Plag 2006, Cazenave and Remy 2011], but both factors also imply significant local and regional variations [Spada and Galassi 2016]. These are also controlled by changes of atmospheric pressure and, in the case of semi-enclosed basins like the Mediterranean Sea, by the forcing exerted from nearby ocean masses driving vari- ations of relative sea level [Umgiesser et al. 2011]. Other causes are significantly contributing, such as the depo- sition and compaction of sediments, tectonic move- ments, the soil reaction to fluid or gas withdrawal, and the response of the Earth to changes in surface loading [Milne et al. 2009]. Understanding the mechanisms driving sea-level changes is important to promote mit- igation measures, when these are possible [McGrana- han et al. 2007, Nicholls and Cazenave 2010]. It has long been recognised that the northern Adri- atic Sea is particularly vulnerable to sea-level rise. This stems from the combined effect of low elevation coast- lines, subsidence and intense urbanization [Armaroli et al. 2012, Perini et al. 2015]. The area is of great impor- tance from the social and economical standpoints, mainly for several ports located here (e.g., Venice, Trieste, Ravenna, Dubrovnik, Split; see Figure 1), as well as for Article history Received March 30, 2016; accepted September 14, 2016. Subject classification: Sea-level change, Tide gauges, Land subsidence. tourism and cultural heritage [e.g. Cori 1999, Preti and De Nigris 2009]. An important cause of relative sea- level change in some parts of this area facing the north- ern Adriatic Sea is land subsidence [Bitelli et al. 2010], of both natural and anthropogenic origins [Brunetti et al. 1998, Carminati et al. 2003, Stright Jr et al. 2008, Baldi et al. 2009]. The effects of subsidence are particularly important in the area of Ravenna in Emilia-Romagna region, which belongs to the delta portion of the Po River plain (Figure 1). Due to its sedimentary origin, the Po valley is subject to natural subsidence, which is partly compensated by the transportation of sediments and by the advance of the coastline [Teatini et al. 2005]. Over the last century, the observed subsidence showed a significant increase in the area of Ravenna, with rates as large as tens of millimetres per year between 1972 and 1977 [Bertoni et al. 1973, Carbognin et al. 1978, Carbognin et al. 1984, Teatini et al. 2005]. A first cause of subsidence is the ground water withdrawal from subsurface aquifers, promoted by the economic growth after the end of World War II and by the development of tourism along the coast [Preti and De Nigris 2009, Bitelli et al. 2010]. Another source of subsidence is the gas extraction from deep on-shore and off-shore reser- voirs, with the first producing a diffuse sinking along the coasts and the latter having much more localised effects [Gambolati et al. 1991, 1999]. In addition to these causes, sea-side processes as the steep drop in the sand transportation by rivers due to the overbuilding of the river banks, or the reclamation of wetlands and lagoons areas on the north, are playing a significant role [Teatini et al. 2005, Preti and De Nigris 2009]. For its vulnerability and for the high subsidence rates, the area of Ravenna has been intensively studied. On one hand, efforts from integrated modelling ap- proaches [Gambolati et al. 1991, 1999] have allowed a quantitative evaluation of the main causes impacting the stability of the whole area. On the other hand, var- ious geodetic surveys have provided new insights on the recent rates of subsidence. Over time, the surveys have been planned by the Italian Geographic Military Insti- tute (IGM), the regional environmental service (ARPAE, formerly ARPA), the Consorzio di Bonifica di Ravenna, the Geological Service of the Ravenna Municipality, and the Italian oil company ENI-E&P. Furthermore, in the northern Adriatic Sea, the Porto Corsini tide gauge sited in Marina di Ravenna (formerly known as Porto Corsini) has collected a long record (1897-2014) of rel- ative sea level, in a location that during last century has been particularly affected by land subsidence [Teatini et al. 2005]. No other comparatively long sea-level time series is available from instruments deployed along the coasts of Emilia-Romagna region. It has long been known (see e.g., Caputo [1971]) that long-term relative sea-level variations at the Porto Corsini tide gauge are strongly affected by the vertical land movements and possibly by pier instability. Consequently, the correla- tion with sea-level signals from other tide gauges de- ployed along the Adriatic and northern Mediterranean coast (Venice, Trieste and Genoa) is poor, as also pointed out by Raicich [2011] and Tsimplis et al. [2012]. Indeed, the mean rate of relative sea-level change shown CERENZIA ET AL. 2 Figure 1. Map of the northern Adriatic Sea and surrounding regions, showing the toponyms of interest for this work. Italy and Emilia- Romagna region are coloured in light grey and yellow, respectively. A blue line marks the path of the Po River. 3 by the Porto Corsini tide gauge record exceeds by 3-6 times the typical Mediterranean rates during the period 1897-2011 [Tsimplis et al. 2012]. Since they can help to assess the long-term subsidence, the relative sea-level data from the tide gauge of Porto Corsini are of re- markable importance. In the framework outlined above, the present study aims at investigating how the sea-level rise ob- served at the Porto Corsini tide gauge is reflecting the observed local land subsidence. During the last few decades, this is accomplished by means of yet unpub- lished geodetic data, which allow for estimating the rate of change of absolute sea level. We believe that this study may improve the reliability of the Porto Corsini tide gauge observations, which are now included in the quality-controlled section of the Permanent Service for Mean Sea Level database (PSMSL, see Holgate et al. [2012]) only for a very limited time span (1969-1972). The paper is organised as follows. In Section 2, we describe the data and the methods employed in our study, in Section 3 we first estimate the long-term rela- tive sea level and subsidence at the Porto Corsini tide gauge, and then we focus on the last few decades. The results are discussed in Section 4 and the conclusions are drawn in Section 5. 2. Data and methods 2.1. Sea level measurements The Porto Corsini tide gauge station (44.49° N, 12.28° E) is placed on the Italian Adriatic coastline in the municipality of Marina di Ravenna, east of the town of Ravenna in Emilia-Romagna region (see Figure 1). Here, instrumental sea-level observations are available since 1897, when the station became part of the newly born Italian tide gauge network. The annual mean sea-level data of Porto Corsini used in this work come from several sources: (i) the Permanent Service for Mean Sea Level (PSMSL) archive (data retrieved on January 1, 2016, from http://www. psmsl.org/), for the time period 1897-1922 [Wood- worth and Player 2003, Holgate et al. 2012], (ii) Annali Idrologici di Bologna (Hydrology Annals of Bologna) for the period 1934-1979 [Annali 1934-1979], and (iii) DSTN/APAT/ISPRA (Institute for Environmental Pro- tection and Research) for the period 1980-2014 (data are available from http://www.mareografico.it). A plot of the annual mean sea-level record for Porto Corsini is shown in Figure 2. In order to obtain a coherent time series, the different datasets have been homogenised re- lying upon the information on benchmarks and data overlaps (see also Mosetti [1969]). A relocation of the tide gauge was, for example, necessary due to works for the re-construction of the dock next to the Marina di Ravenna lighthouse. The sea-level measurements at the Porto Corsini tide gauge are historically referred to the IGM bench- mark located in the Marina di Ravenna lighthouse, sited 150 m away from the tide gauge. The lighthouse was connected in 1885 to the first IGM Italian leveling net- work [Salvioni 1957], gathering at the time the main Italian harbours, with the establishment of the Porto Corsini tide gauge. However, this initial network re- sulted inadequate and after the end of World War II it was completely re-established. The first leveling of the lighthouse benchmark with respect to the national al- timetric datum, defined by the tide gauge in Genoa in 1942, was accomplished in 1970. Subsequently, other SUBSIDENCE IN MARINA DI RAVENNA Figure 2. Annual values of mean sea-level at the Porto Corsini tide gauge in Marina di Ravenna during 1897-2013. The blue line represents the long-term trend obtained by applying the least squares regression method over the entire time series, while the red one shows the trend during 1970-2013. level surveys were performed by IGM in 1977 and 1990, which provided a reference for the sea-level data col- lected by the Porto Corsini tide gauge since 1990. In 2012, the lighthouse benchmark was connected to the IGM network and since September 2013 the height ob- served in that occasion became the reference value used by the tide gauge. 2.2. Land subsidence measurements In this work, different measurements of land alti- tude and of its rate of change are considered and com- pared, during different time periods. As mentioned in the Introduction, land subsidence in the Ravenna area has been monitored by several agencies. Here, we con- sider only locations geographically or historically re- lated to the Porto Corsini tide gauge site. A summary of the monitoring campaigns is presented in Table 1. Concerning the IGM measurements, points prior to World War II resulted inconsistent with respect to the rest of the period. Salvioni [1957] homogenized the 1897-1957 dataset by using a compensation approach and the same IGM datum. Other surveys were performed by the Consorzio di Bonifica di Ravenna in 1949 and 1972, with the establishment of a local leveling network in Ravenna Municipality. The elaboration in contour maps of these IGM and Consorzio di Bonifica di Ra- venna surveys for the periods 1897-1957 and 1949-1972 [Salvioni 1957, Carbognin et al. 1984, Arca and Beretta 1985] was later georeferenced and gridded in GIS maps by Teatini et al. [2005], to which we refer to qualita- tively characterise the land subsidence at Marina di Ravenna prior to 1970. First IDROSER SpA and after- words the ARPAE agency, on behalf of Emilia-Romagna region, monitored the subsidence by using the high precision geometric leveling method, assuming as a sta- ble reference the IGM benchmark at Sasso Marconi in the Apennine Mountains [Regione Emilia-Romagna - Arpa 2006]. Many sites along the Ravenna coast were monitored, including the Marina di Ravenna lighthouse and the tide gauge station on the dock (Table 1). The re- sults from the 1993 survey campaign are not considered here, since doubts on their validity have been raised [Regione Emilia-Romagna - Arpa 2001]. Since 2005, ARPAE continued the land subsidence monitoring by the support of SAR Interferometry (PSInSARTM tech- nique), to provide rates of subsidence [Regione Emilia- Romagna - Arpa 2007, 2012]. Two different calibration methods were used, pro- viding rates of vertical movement for periods 1992- 2000, 2003-2006, and 2006-2011 (see Table 1). For the first two periods SAR Interferometry was calibrated by 2005 leveling surveys [Regione Emilia-Romagna - Arpa 2007], while, for the last period, a network of 16 GPS- GNSS permanent stations located in the Emilia-Ro- magna plain [Bitelli et al. 2014, 2015] was used [Regione Emilia-Romagna - Arpa 2012]. Since none of these GPS stations is located in or near Marina di Ravenna, we did not include any GPS data in the present analysis. We consider the scatter points available in the vicinity of the tide gauge site, i.e., on the dock structure and on top of the lighthouse or in its neighbourhood (an example of these points is shown in Figure 3). When several scatter points were available on the same structure (this occurred during period 2003-2006), averaged values have been considered. Table 2 reports raw data for the height measurement and for the rate of vertical movement, CERENZIA ET AL. 4 Agency Date Type of survey Notes IGM and Consorzio di Bonifica di Ravenna 1897, 1949, 1972 Geometric leveling GIS maps in Teatini et al. [2005] IGM 1970, 1977, 1990 Geometric leveling Lighthouse benchmark IDROSER 1984, 1987 Geometric leveling Lighthouse benchmark and tide gauge benchmark on the dock ARPAE 1999, 2005 Geometric leveling Lighthouse benchmark and tide gauge benchmark on the dock ARPAE 1992-2000 InSAR calibrated on 2005 geometric leveling Two scatter points: one on the dock and one over the lighthouse ARPAE 2003-2006 InSAR calibrated on 2005 geometric leveling Four scatter points: two on the dock and two over the lighthouse ARPAE 2006-2011 InSAR calibrated on GPS permament network Two scatter points: one on the dock and one over the closest building to the lighthouse Table 1. Main surveys carried out in the Marina di Ravenna area by various agencies, with data used in the present work. 5 collected on the dock structure and on to the light- house or in its neighbourhood from 1970 to 2011. 3. Results 3.1. Long-term relative sea level and land subsidence in Marina di Ravenna The historical time series of relative sea level at the Porto Corsini tide gauge is shown in Figure 2 (the data sources are summarized in Section 2.1). Here, annual mean relative sea-level data are presented, over a period exceeding one century (1897-2014). By considering the whole time series, the total variation of relative sea- level is ~1.0 m. To assess the rate of relative sea-level change, we applied the least squares regression, which provided a value of 8.5±0.2 mm yr−1. This trend, dis- played as a blue line in Figure 2, is comparable to the one that can be obtained by applying the Theil-Sen es- timator [Theil 1950, Sen 1968], i.e. 8.7 mm yr−1 (the 95% confidence interval for the rate ranges between 8.1 and 9.5 mm yr−1). To obtain this statistically significant (p < 0.001) trend, we adopted the Theil-Sen implemen- SUBSIDENCE IN MARINA DI RAVENNA Date/Interval Lighthouse Dock Level (m a.s.l.) Rate (mm yr−1) Level (m a.s.l.) Rate (mm yr−1) 1970 0.9325 1977 0.7896 1984 0.6866 1.0688 1987 0.6315 1.0187 1990 0.6273 1999 0.5205 0.9298 2005 0.4639 0.8848 1992-2000 −6.7±1.0 −5.9±1.0 2003-2006 −7.2±2.6 −5.7±2.6 2006-2011 −5.3±2.0 −3.1±2.0 Figure 3. Aerial view of Marina di Ravenna showing the Porto Corsini tide gauge station and the surroundings. The ARPAE benchmark is sited in the same place of the tide gauge. The altimetric reference is the benchmark of the lighthouse. Red dots mark the scatter points from the interferometric satellite considered in this work, for the period 2006-2011 (units are mm yr−1). Table 2. Direct and undirect height measurements at the lighthouse and at the dock of Marina di Ravenna from the subsidence surveys pre- sented in Table 1. A distinction is made between data sampled as height a.s.l. (e.g. levelling data, column “Level”) and as average rates dur- ing specific time intervals (e.g. InSAR, column “Rate”). tation in the openair software [Carslaw and Ropkins 2012], after having filled the gaps in the time series by applying a linear interpolation. Such rate of relative sea- level rise largely exceeds that obtained for other stations along the northern Adriatic and Mediterranean coasts over comparable time periods. According to Raicich [2011] and Tsimplis et al. [2012], the rates of relative sea level observed in Venice, Trieste and Genoa are 2.6, 1.3 and 1.2 mm yr−1, respectively. By applying a quad- ratic (three-parameter) regression, we find that the value of the sea-level acceleration at Marina di Ravenna is 0.08±0.02 mm yr−2, which provides evidence for a sig- nificant increase in the long-term rate of relative sea- level change. Although existing estimates of global long-term sea-level acceleration are scattered (see re- view in Spada et al. [2015], and references therein), the value for Marina di Ravenna largely exceeds all of these. The historical subsidence in the area of Marina di Ravenna is qualitatively inferred from the GIS slicing of the leveling surveys along the shoreline of Ravenna Municipality from the work of Teatini et al. [2005], by associating the subsidence between the intervals 1897- 1949 and 1949-1972 with the benchmark elevation measured by IGM in 1970 at the Marina di Ravenna lighthouse (see Figure 14 in Teatini et al. [2005]). The time history of subsidence is plotted in Figure 4, where the observations are referred to the lighthouse of Ma- rina di Ravenna. Coherently with the general trend in the Ravenna area (see the case study of Porta Adriana in the town of Ravenna [Carbognin et al. 1978]), the subsidence in the area of Marina di Ravenna experi- enced a marked increase during 1950-1970, with rates close to 24 mm yr−1. These exceed the “natural trend” during 1897-1950 by about 6 mm yr−1. In the ensuing years, the rate decreased, but still exceeding the natural trend. The last values obtained from the interferomet- ric satellite analysis during 2006-2011 are close to 5.5 mm yr−1. Overall, these data are consistent with the in- crease in the groundwater withdrawal between the 50s and late 70s and its subsequent decline. While the main- land appears to be substantially stable, many sites along the coast are still subsiding at rates exceeding the natu- ral trend [Teatini et al. 2005]. Marina di Ravenna is co- herent with this general behaviour, although the most recent interferometric analysis (2006-2011) has evi- denced a recovering of the slower natural trend of sub- sidence. New data from the interferometric analysis scheduled for year 2016 will be useful to confirm or deny this trend. 3.2. Sea level at Marina di Ravenna since 1970 During the time period 1970-2013, monthly aver- aged values of relative sea level are available for Marina di Ravenna, characterised by a high degree of com- pleteness. We have applied the least squares regression to this relatively short time series, obtaining a trend of 7.7±0.3 mm yr−1, displayed as the red line in Figure 2. Again, this trend resulted comparable to the significant (p < 0.001) trend obtained by the Theil-Sen method, i.e. 7.5 mm yr−1 (the 95% confidence interval for the rate ranges between 7.0 and 8.1 mm yr−1). To analyse the sea-level variability coherently, seasonal cycles should CERENZIA ET AL. 6 Figure 4. Altitude measurements in Marina di Ravenna at various epochs. Given the various elaborations applied to data before 1970 (see Section 2.2), the association with the 1970 IGM leveling is approximate (dashed line). Data after 1970 refer to different sources and survey methods, marked by various symbols (see legend). Further details are given in Tables 1 and 2. 7 be removed, since these derive from regional coastal ef- fects. In many coastal areas, a large fraction of the sea- sonal cycle is associated with the steric component of sea-level change, induced by density variations in re- sponse to temperature and salinity fluctuations [Tsim- plis et al. 2012]. However, due to the influence of regional hydrology, the non-steric component of the seasonal cycle of sea level was found to be important in some cases. The seasonal cycle and the amplitude of coastal sea-level height thus represent the local response to a seasonal pattern of large-scale sea-level variability and hydrology (see Bergant et al. [2005], and references therein). In Marina di Ravenna, we estimated the “cli- matological” annual variation considering monthly val- ues, although the details of this analysis are not shown here. The minimum sea-level height was found in Feb- ruary (1.51±0.26 m), while two maxima were found in November (1.64±0.25 m) and June (1.55±0.22 m). This is consistent with with the results of Bergant et al. [2005], showing a minimum during the first months of the year, and two maxima in November and April. However, by a Mann-Whitney “U test” [Mann and Whitney 1947], we have found that the difference be- tween the amplitude of the two peaks is not statistically significant (p = 0.764). To detect possible cyclic components in the tide gauge time series of Porto Corsini since 1970, we have employed the empirical mode decomposition (EMD) method of Huang et al. [1998], in the improved version of Torres et al. [2011] (the Complete ensemble EMD with adaptive noise). The EEMD has been adopted by Breaker and Ruzmaikin [2011] to analyse the San Fran- cisco tide gauge record and, more recently, by Spada et al. [2014] to study the Nuuk/Godthab (SW Greenland) time series, by Galassi and Spada [2014] to characterise the cyclic sea-level variations across the Adriatic Sea and by Vecchio et al. [2014] to study the transient effect of a tsunami in the Mediterranean Sea. The EEMD allows to split non-linear and non-stationary time series into a limited sequence of empirically orthogonal “intrinsic mode functions” (IMFs) that describe cyclic modes not necessarily characterised by a constant amplitude nor phase. Differently from traditional spectral methods, the EEMD is not requiring assumptions on the func- tional expression of the regression model. Figure 5 shows the seven IMFs that compose the non monoto- nic portion of the Porto Corsini tide gauge record. IMF2 and IMF3 capture the semi-annual and annual pe- SUBSIDENCE IN MARINA DI RAVENNA Figure 5. Results of the EMD analysis of the Marina di Ravenna relative sea-level time series since 1970, with the IMFs in the top frames and the residual (RES) in the bottom, respectively. Units are cm throughout. riodicities, respectively. IMF4 corresponds to an oscilla- tion at about 1.7 years whereas the IMF5 reveals a pe- riod of about 3.6 years; this latter has been detected, adopting a similar approach, on the whole Adriatic Sea scale (see Galassi and Spada [2014]). IMF6 and IMF7 show longer periods, with 8.4 and 14.2 years, respec- tively. RES represents the non-cyclic residual, revealing the long-term natural trend of the time series [Huang et al. 1998] that remains once all the cyclic components have been sifted. The linear trend of RES since 1970 is 6.8 mm yr−1, somewhat smaller than the one obtained above by the Theil-Sen method. 3.3. Land subsidence at Marina di Ravenna since 1970 To estimate the rate of subsidence in Marina di Ravenna since 1970, we have considered the full dataset from the IGM and ARPAE leveling surveys and the re- sults of the interferometric analysis (see Tables 1 and 2). Values for years 1992, 2000, 2003, 2006 and 2011 have been derived from the three interferometric analyses, by computing the height variation corresponding to the interferometric velocity multiplied by the time dif- ference from the previous or the following leveling datum. This has allowed to compute regression curves and to estimate the rates of subsidence of the ground settlement at the lighthouse and at the tide gauge scat- ter points located over the dock, separately. The data shown in Figure 6, clearly indicate that the rate of land subsidence at the lighthouse benchmark (vl ) has been significantly larger than the rate at the tide gauge point over the dock (vd ). Indeed, by two linear regressions we obtain vl = −13.0±0.7 mm yr −1 (1970 - 2011) (1) and vd = −8.3±0.4 mm yr −1 (1984 - 2011), (2) for the two points, respectively. The different rates of subsidence may be explained as a local effect due to the weight of the structures on which the lighthouse and the dock are located, being the embankment lighter compared to the base of the building. The same con- clusion can be reached only considering the level sur- veys data, available for the period 1984-2005, which provide the values vl = −10.3±0.7 mm yr −1, and vd = −8.3±0.7 mm yr−1, for the lighthouse and for the dock, respectively. This further motivates the choice of con- sidering separately the rates observed on the dock structure and those observed close to the lighthouse. An acceleration in the rate of land settlement is found at both points, with comparable amplitudes. From a quadratic regression we obtain al = 0.19±0.03 mm yr −2 for period 1970-2011 and ad = 0.18±0.05 mm yr −2 dur- ing 1984-2011. 4. Discussion A simultaneous evaluation of the sea-level rise and land subsidence in Marina di Ravenna is possible dur- ing the time interval 1990-2011. In fact, the elevation of the dock structure was measured for the last time in 2011 and the reference altitude value used by the tide gauge was not recoverable before 1990, when the IGM level of 1990 was adopted. During 1990-2011, the dock CERENZIA ET AL. 8 Figure 6. Elevation measurements at the Porto Corsini lighthouse (red) and dock (green) benchmarks at various epochs. The various data sources and sampling methods are denoted by different symbols, as indicated in the legend (see Tables 1 and 2 for details). 9 structure subsided by 167±2 mm, where this value is derived from the linear regression of the dock time se- ries that indicates a rate of subsidence of 7.6±0.1 mm yr−1. During the same time span, in Marina di Ravenna, relative sea-level rose by 216±29 mm, a value obtained by a linear regression of the annual averaged time series between 1990 and 2013 that provides a rate of change of 9.8±1.3 mm yr−1. On longer time periods, the analysis becomes more problematic, since the land movement is measured only at the Marina di Ravenna lighthouse since 1970 and in the previous years only one estimate across the whole Marina di Ravenna area was available. It is likely that the vertical movement at this site had a dominant role on the observed relative sea level, although the available data do not allow to obtain any quantitative estimate. Our results firmly indicate that, during the period 1990-2011, local land subsidence was the major re- sponsible for relative sea level rise observed at the site of Porto Corsini. Anthropogenic factors, as water extrac- tion and soil compaction, are known to be the major causes of land subsidence, contributing to the relative sea-level variations observed, as explained in Teatini et al. [2005]. Other possible factors contributing to sea- level change in the Adriatic Sea are worth to be dis- cussed, although their role appears to be less important. The modelled contribution of glacial isostatic adjust- ment (GIA) to the rate of relative sea-level change at the Porto Corsini tide gauge location only amounts to 0.1 mm yr−1, according to the ICE-6G model (http://www.atmosp.physics.utoronto.ca/~peltier/ data.php) [Peltier et al. 2015]. This implies that for the time period 1990-2011, only ~2 mm of the observed relative sea-level change can be attributed to the ef- fects of GIA. The contribution of co-seismic and post- seismic vertical deformation to relative sea-level change at tide gauges has been assessed at global level by Melini et al. [2004] and Melini and Piersanti [2006]. Using data from the Seismic Moment Tensor cata- logue, they have clearly shown that earthquakes have cumulatively the tendency to produce a positive sea- level trend of ~ 0.2 mm yr−1 at the tide gauges loca- tions. The low level of seismicity along the coasts of the northern Adriatic Sea may suggest that this global value may represent an upper bound for the site of Marina di Ravenna. In addition, in their work on the effects of seismicity on subsidence in foreland basins with an application to the case study of Venice, Carminati et al. [2007] have suggested that accelera- tions in the rate of subsidence due to co-seismic dis- placements associated to earthquakes that occurred along active faults around and within the Po Plain are negligible. According to the “sea-level equation” (see, e.g., Spada [2016]), the rate of absolute sea-level change at the Porto Corsini tide gauge can be estimated by sub- tracting the rate of subsidence from the rate of rela- tive sea-level change. During the time period 1990-2011, we obtain a trend of 2.2±1.3 mm yr−1. Therefore, de- spite the large influence of land vertical movement at Marina di Ravenna, the absolute sea-level increase can be recognised combining the rate of subsidence with the tidal record. Since the Porto Corsini tide gauge is not facing the open sea (see Figure 3), a direct com- parison with the rate of absolute sea-level change ob- served by altimetry is not feasible. However, the satellite altimetry trend during 1993-2012 in the north- ern Adriatic Sea has been close to 3 mm yr−1 [Bona- duce et al. 2016], in fair agreement with the value obtained from in situ observations at the Porto Corsini tide gauge. 5. Conclusions The analysis performed in our work has permit- ted, for the first time, the simultaneous evaluation of the sea-level rise and of land subsidence in Marina di Ravenna, along the North-Adriatic coast of Italy. Al- though long records exist for both the relative sea- level variations at the Porto Corsini tide gauge (1897-2014) and for the historical land movements in the Marina di Ravenna area (1897-2011), a quantita- tive analysis has been possible only for the period 1990-2011, due to missing information about tide gauge setup before 1990 and since the last useful data of height measurement dates 2011. We have enlight- ened the importance of studying land movements at the Porto Corsini tide gauge by relying upon meas- urements collected on the dock structure above which the tide gauge is located. Indeed, the analysis of the available geodetic data indicates that the reference point at which the tide gauge is usually connected (i.e., the lighthouse) has been subject to a rate of sub- sidence significantly exceeding that of the dock struc- ture. In the time interval 1990-2011, the land subsided at a rate of 7.6±0.1 mm yr−1 while relative sea-level rose at a rate of 9.8±1.3 mm yr−1. Therefore, our study confirms that land subsidence has been the major cause of relative sea-level change at Marina di Ravenna, at least during the period permitting a quantitative comparison (1990-2011) between relative sea level and vertical land motion signals. The resid- ual (rate of relative sea-level change minus rate of land subsidence) of 2.2±1.3 mm yr−1 is found to be coherent with the altimetric observations in the northern Adriatic for the rate of change of absolute sea-level during the same time period. SUBSIDENCE IN MARINA DI RAVENNA Acknowledgements. Some of the Porto Corsini 1993-2014 data were obtained by the DSTN/APAT/ISPRA website. The sea- coast unit of Regione Emilia-Romagna is also acknowledged for useful discussion. This work has been partly funded by a research grant of the Department of Pure and Applied Sciences (DiSPeA) of the Urbino University “Carlo Bo” (CUP H32I160000000005 and H32I15000160001). Fabio Raicich and two anonymous reviewers have greatly contributed to improve the manuscript. References Annali (1934-1979). Annali Idrologici, Tech. rep., Uffi- cio Idrografico del Po, Parma; http://www.arpa.emr. it/cms3/documenti/subsidenza/Relfin 2001_PDF. Arca, S., and G. Beretta (1985). Prima sintesi geodetico- geologica sui movimenti verticali del suolo nell’Ita- lia Settentrionale (1897-1957), Bollettino di Geode- sia e Scienze Affini, 44 (2), 125-156. Armaroli, C., P. Ciavola, L. Perini, L. Calabrese, S. Lo- rito, A. Valentini and M. Masina (2012). Critical storm thresholds for significant morphological changes and damage along the Emilia-Romagna coastline, Italy, Geomorphology, 143, 34-51. Baldi, P., G. Casula, N. Cenni, F. Loddo and A. Pesci (2009). GPS-based monitoring of land subsidence in the Po Plain (northern Italy), Earth and Planetary Science Letters, 288 (1), 204-212. Bergant, K., M. Sušnik, I. Strojan and A.G.P. Shaw (2005). Sea level variability at Adriatic coast and its relationship to atmospheric forcing, Annales Geo- physicae, 23 (6), 1997-2010. Bertoni, W., L. Carbognin, P. Gatto and G. Mozzi (1973). Note interpretative preliminari sulle cause della sub- sidenza in atto a Ravenna, Tech. Rep. 65, CNR - La- boratorio Dinamica Grandi Masse, Venezia (Italy). Bitelli, G., F. Bonsignore, L. Carbognin, A. Ferretti, T. Strozzi, P. Teatini, L. Tosi and L. Vittuari (2010). Radar interferometry-based mapping of the present land subsidence along the low-lying northern Adri- atic coast of Italy. Land Subsidence, Associated Haz- ards and the Role of Natural Resources Development, In: Proceedings of EISOLS 2010, vol. 339 (IAHS, Queretaro, Mexico, October 17-22, 2010). Bitelli, G., F. Bonsignore, S. Del Conte, F. Novali, I. Pellegrino and L. Vittuari (2014). Integrated use of Advanced InSAR and GPS data for subsidence mon- itoring, In: G. Lollino et al. (eds.), Engineering Ge- ology for Society and Territory, vol. 5, III, Springer, 147-150. Bitelli, G., F. Bonsignore, I. Pellegrino and L. Vittuari (2015). Evolution of the techniques for subsidence monitoring at regional scale: the case of Emilia-Ro- magna region, (Italy), In: Proceedings of the Ninth International Symposium on Land Subsidence, vol. 92 (IAHS, Nagoya, Japan, November 15-19), 1-7. Bonaduce, A., N. Pinardi, P. Oddo, G. Spada and G. Larnicol (2016). Sea-level variability in the Mediter- ranean Sea from altimetry and tide gauges, Climate Dynamics, 1-16. Breaker, L., and A. Ruzmaikin (2011). The 154-year record of sea level at San Francisco: extracting the long-term trend, recent changes, and other tidbits, Climate Dynamics, 36 (3/4), 545-559. Brunetti, A., M. Denèfle, M. Fontugne, C. Hatté and P. Pirazzoli (1998). Sea-level and subsidence data from a Late Holocene back-barrier lagoon (Valle Standia- na, Ravenna, Italy), Marine Geology, 150 (1), 29-37. Caputo, M. (1971). Intervento, In: I Movimenti del Suo- lo nel Ravennate. Atti della tavola rotonda tenutasi il 3 aprile 1971, Lions Club / Rotary Club Ravenna, 45-57. Carbognin, L., P. Gatto, G. Mozzi and G. Gambolati (1978). Land subsidence of Ravenna and its similar- ities with the Venice case, In: Evaluation and pre- diction of subsidence, ASCE, 254-266. Carbognin, L., P. Gatto and F. Marabini (1984). Guide- book of the eastern Po plain (Italy): A short illus- tration about the environment and land subsidence, Ufficio Stampa, Comunicazione e Informazione. Carminati, E., C. Doglioni and D. Scrocca (2003). Apen- nines subduction-related subsidence of Venice (Italy), Geophysical Research Letters, 30 (13). Carminati, E., S. Enzi and D. Camuffo (2007). A study on the effects of seismicity on subsidence in fore- land basins: An application to the Venice area, Global and Planetary Change, 55 (4), 237-250. Carslaw, D.C., and K. Ropkins (2012). openair – An R package for air quality data analysis, Environmen- tal Modelling & Software, 27, 52-61. Cazenave, A., and W. Llovel (2010). Contemporary sea level rise, Annual Review of Marine Science, 2, 145- 173. Cazenave, A., and F. Remy (2011). Sea level and climate: measurements and causes of changes, Wiley Inter- disciplinary Reviews: Climate Change, 2 (5), 647-662. Church, J., P. Clark, A. Cazenave, J. Gregory, S. Jevre- jeva, A. Levermann, M. Merrifield, G. Milne, R. Nerem, P. Nunn, A. Payne, W. Pfeffer, D. Stammer and A. Unnikrishnan (2013). Sea Level Change, In: T. Stocker, D. Qin, G.-K. Plattner, M. Tignor, S. Allen, J. Boschung, A. Nauels, Y. Xia, V. Bex and P. Midgley (eds.), Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovern- mental Panel on Climate Change, Cambridge Uni- versity Press, Cambridge, 1138-1191. Cori, B. (1999). Spatial dynamics of Mediterranean coastal regions, Journal of Coastal Conservation, 5, CERENZIA ET AL. 10 11 105-112. Galassi, G., and G. Spada (2014). Linear and non-linear sea-level variations in the Adriatic Sea from tide gauge records (1872-2012), Annals of Geophysics, 57 (6), P0658. Gambolati, G., G. Ricceri, W. Bertoni, G. Brighenti and E. Vuillermin (1991). Mathematical simulation of the subsidence of Ravenna, Water Resources Research, 27 (11), 2899-2918. Gambolati, G., P. Teatini and L. Tomasi (1999). Coast- line regression of the Romagna region, Italy, due to natural and anthropogenic land subsidence and sea level rise, Water Resources Research, 35 (1), 163-184. Holgate, S., A. Matthews, P. Woodworth, L. Rickards, M. Tamisiea, E. Bradshaw, P. Foden, K. Gordon, S. Jevrejeva and J. Pugh (2012). New data systems and products at the permanent service for mean sea level, Journal of Coastal Research, 29 (3), 493-504. Huang, N., Z. Shen, S. Long, M. Wu, H. Shih, Q. Zheng, N.-C. Yen, C. Tung and H. Liu (1998). The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time se- ries analysis, In: Proceedings of the Royal Society of London, Series A: Mathematical, Physical and En- gineering Sciences, 454 (1971), 903-995. Lombard, A., A. Cazenave, P. Le Traon and M. Ishii (2005). Contribution of thermal expansion to pres- ent-day sea-level change revisited, Global and Plan- etary Change, 47, 1-16. Mann, H.B., and D.R. Whitney (1947). On a test of whether one of two random variables is stochasti- cally larger than the other, The Annals of Mathe- matical Statistics, 50-60. McGranahan, G., D. Balk and B. Anderson (2007). The rising tide: assessing the risk of climate change and human settlements in low elevation coastal zones, Environment and Urbanization, 19, 17-37. Melini, D., A. Piersanti, G. Spada, G. Soldati, E. Casa- rotti and E. Boschi (2004). Earthquakes and relative sealevel changes, Geophysical Research Letters, 31, L09601. Melini, D., and A. Piersanti (2006). Impact of global seismicity on sea level change assessment, Journal of Geophysical Research, 111, B03406. Milne, G.A., W.R. Gehrels, C.W. Hughes and M.E. Tamisiea (2009). Identifying the causes of sea-level change, Nature Geoscience, 2 (7), 471-478. Mosetti, F. (1969). Le variazioni relative del livello ma- rino nell’Adriatico dal 1896 al 1967 e il problema dello sprofondamento di Venezia, Bollettino di Geo- fisica Teorica ed Applicata, 11 (43/44), 243. Nicholls, J., and A. Cazenave (2010). Sea-level rise and its impact on coastal zones, Science, 328 (5985), 1517-1520. Peltier, W., D. Argus and R. Drummond (2015). Space geodesy constrains ice age terminal deglaciation: The global ICE-6G C (VM5a) model, Journal of Geophysical Research: Solid Earth, 120 (1), 450-487. Perini, L., L. Calabrese, G. Salerno, P. Ciavola and C. Armaroli (2015). Evaluation of coastal vulnerability to flooding: comparison of two different method- ologies adopted by the Emilia-Romagna Region (Italy), Natural Hazards and Earth System Sciences, 3, 4315- 4352. Plag, H.-P. (2006). Recent relative sea-level trends: an at- tempt to quantify the forcing factors, Philosophical Transactions of the Royal Society A, 364, 821-44. Preti, M., and N. De Nigris (2009). One meter of sub- sidence along Emilia Romagna coast in the last 55 years: effect and defence strategies, In: Proc. of the tenth international conference on the mediter- ranean coastal environmenit, MEDCOAST09 (Erdal O zhan, Sochi, Russia, November 10-14). Raicich, F. (2011). On Sea-Level changes in the Mediter- ranean Sea with focus on the Adriatic Sea, In: E. Brugnoli, G. Cavarretta, S. Mazzola, F. Trincardi, M. Ravaioli and S. Santoleri (eds.), Marine Research at CNR, Consiglio Nazionale delle Ricerche, Rome, Italy, 1435-1443. Regione Emilia-Romagna - Arpa (2001). Misura della rete regionale di controllo della subsidenza, misura di linee della rete costiera non comprese nella rete regionale, rilievi batimetrici, Tech. rep., Bologna, Arpa, Direzione Tecnica, F. Bonsignore (ed.); http:// www.arpae.it/cms3/documenti/subsidenza/Relfin _2001.PDF. Regione Emilia-Romagna - Arpa (2006). Rilievo della subsidenza nella pianura emiliano-romagnola. Mi- sura della rete costiera di controllo della subsidenza, Tech. rep., Bologna, Arpa, Direzione Tecnica, F. Bonsignore (ed.); http://www.arpae.it/cms3/docu menti/subsidenza/Relfin rete costiera.pdf. Regione Emilia-Romagna - Arpa (2007). Rilievo della subsidenza nella pianura emiliano-romagnola. Ana- lisi interferometrica, Tech. rep., Bologna, Arpa, Di- rezione Tecnica, F. Bonsignore (ed.); URL:http:// www.arpae.it/cms3/documenti/subsidenza/Relfin _2007 rid.pdf. Regione Emilia-Romagna - Arpa (2012). Rilievo della subsidenza nella pianura emiliano-romagnola. Se- conda fase, Tech. rep., Bologna, Arpa, Direzione Tecnica, F. Bonsignore (ed.); http://www.arpae.it/ cms3/documenti/subsidenza/Relfin_2012.pdf. Roemmich, D., and W. Owens (2000). The ARGO proj- ect: global ocean observations for understanding for understanding and prediction of climate variability, SUBSIDENCE IN MARINA DI RAVENNA Oceanography, 13 (2), 45-50. Salvioni, G. (1957). I movimenti del suolo nell’Italia cen- tro-settentrionale, Bollettino di Geodesia e Scienze Affini, 16, 325-366. Sen, P.K. (1968). Estimates of the regression coefficient based on Kendall’s tau, Journal of the American Sta- tistical Association, 63 (324), 1379-1389. Spada, G., J. Bamber and R. Hurkmans (2013). The gravitationally consistent sea-level fingerprint of fu- ture terrestrial ice loss, Geophysical Research Let- ters, 40 (3), 482-486. Spada, G., G. Galassi and M. Olivieri (2014). A study of the longest tide gauge sea-level record in Greenland (Nuuk/Godthab, 1958-2002), Global and Planetary Change, 118, 42-51. Spada, G., M. Olivieri and G. Galassi (2015). A heuristic evaluation of long-term global sea level accelera- tion, Geophysical Research Letters, 42, 4166-4172. Spada, G. (2016). Glacial Isostatic Adjustment and Con- temporary Sea Level Rise: An Overview, Surveys in Geophysics, in press. Spada, G., and G. Galassi (2016). Spectral analysis of sea level during the altimetry era, and evidence for GIA and glacial melting fingerprints, Global and Plane- tary Change, 143, 34-49. Stright Jr, D., A. Settari, D.A. Walters and K. Aziz (2008). Characterization of the Pliocene gas reservoir aquifers for predicting subsidence on the Ravenna coast, Petroleum Science and Technology, 26 (10/11), 1267-1281. Teatini, P., M. Ferronato, G. Gambolati, W. Bertoni and M. Gonella (2005). A century of land subsidence in Ravenna, Italy, Environmental Geology, 47 (6), 831- 846. Theil, H. (1950). A rank invariant method of linear and polynomial regression analysis, i, ii, iii, In: Proceed- ings of the Koninklijke Nederlandse Akademie Wetenschappen, Series A - Mathematical Sciences, vol. 53, Stichting Mathematisch Centrum, 386-392, 521-525, 1397-1412. Torres, M., M. Colominas, G. Schlotthauer and P. Flan- drin (2011). A complete ensemble empirical mode decomposition with adaptive noise, In: Acoustics, Speech and Signal Processing (ICASSP), Interna- tional Conference, IEEE, 4144-4147. Tsimplis, M.N., F. Raicich, L. Fenoglio-Marc, A.G. Shaw, M. Marcos, S. Somot and A. Bergamasco (2012). Re- cent developments in understanding sea level rise at the Adriatic coasts, Physics and Chemistry of the Earth, Parts A/B/C, 40, 59-71. Umgiesser, G., J. Anderson, V. Artale, M. Breil, S. Gualdi, P. Lionello, N. Marinova, M. Orlic, P. Piraz- zoli, S. Rahmstorf et al. (2011). From Global to re- gional: Local Sea Level Rise Scenarios. Focus on the Mediterranean Sea and the Adriatic Sea, In: Work- shop Report No. 1. Vecchio, A., M. Anzidei and V. Carbone (2014). New in- sights on the tsunami recording of the May, 21, 2003, Mw 6.9 Boumerdès earthquake from tidal data analysis, Journal of Geodynamics, 79, 39-49. Woodworth, P.L., and R. Player (2003). The Permanent Service for Mean Sea Level: an update to the 21st century, Journal of Coastal Research, 19, 287-295. Wunsch, C., and D. Stammer (1995). The global fre- quency-wavenumber spectrum of oceanic variabil- ity estimated from TOPEX/POSEIDON altimetric measurements, Journal of Geophysical Research: Oceans (1978-2012), 100 (C12), 24895-24910. *Corresponding author: Marco Olivieri, Istituto Nazionale di Geofisica e Vulcanologia (INGV), Sezione di Bologna, Italy; email: marco.olivieri@ingv.it. © 2016 by the Istituto Nazionale di Geofisica e Vulcanologia. All rights reserved. CERENZIA ET AL. 12 << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles false /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Warning /CompatibilityLevel 1.3 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.1000 /ColorConversionStrategy /LeaveColorUnchanged /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams true /MaxSubsetPct 100 /Optimize false /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments false /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile (None) /AlwaysEmbed [ true /AndaleMono /Apple-Chancery /Arial-Black /Arial-BoldItalicMT /Arial-BoldMT /Arial-ItalicMT /ArialMT /CapitalsRegular /Charcoal /Chicago /ComicSansMS /ComicSansMS-Bold /Courier /Courier-Bold /CourierNewPS-BoldItalicMT /CourierNewPS-BoldMT /CourierNewPS-ItalicMT /CourierNewPSMT /GadgetRegular /Geneva /Georgia /Georgia-Bold /Georgia-BoldItalic /Georgia-Italic /Helvetica /Helvetica-Bold /HelveticaInserat-Roman /HoeflerText-Black /HoeflerText-BlackItalic /HoeflerText-Italic /HoeflerText-Ornaments /HoeflerText-Regular /Impact /Monaco /NewYork /Palatino-Bold /Palatino-BoldItalic /Palatino-Italic /Palatino-Roman /SandRegular /Skia-Regular /Symbol /TechnoRegular /TextileRegular /Times-Bold /Times-BoldItalic /Times-Italic /Times-Roman /TimesNewRomanPS-BoldItalicMT /TimesNewRomanPS-BoldMT /TimesNewRomanPS-ItalicMT /TimesNewRomanPSMT /Trebuchet-BoldItalic /TrebuchetMS /TrebuchetMS-Bold /TrebuchetMS-Italic /Verdana /Verdana-Bold /Verdana-BoldItalic /Verdana-Italic /Webdings ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 150 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.10000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 150 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.10000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.08250 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile (None) /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName (http://www.color.org) /PDFXTrapped /Unknown /CreateJDFFile false /SyntheticBoldness 1.000000 /Description << /ENU (Use these settings to create PDF documents with higher image resolution for high quality pre-press printing. The PDF documents can be opened with Acrobat and Reader 5.0 and later. These settings require font embedding.) /JPN /FRA /DEU /PTB /DAN /NLD /ESP /SUO /NOR /SVE /KOR /CHS /CHT /ITA >> >> setdistillerparams << /HWResolution [2400 2400] /PageSize [595.000 842.000] >> setpagedevice