(Table 1) Lake elevation change derived from 4 or more years of measurements

In this study, ICESat altimetry data are used to provide precise lake elevations of the Tibetan Plateau (TP) during the period of 2003-2009. Among the 261 lakes examined ICESat data are available on 111 lakes: 74 lakes with ICESat footprints for 4-7 years and 37 lakes with footprints for 1 -3 years. This is the first time that precise lake elevation data are provided for the 111 lakes. Those ICESat elevation data can be used as baselines for future changes in lake levels as well as for changes during the 2003-2009 period. It is found that in the 74 lakes (56 salt lakes) examined, 62 (i.e. 84%) of all lakes and 50 (i.e. 89%) of the salt lakes show tendency of lake level increase. The mean lake water level increase rate is 0.23 m/year for the 56 salt lakes and 0.27 m/year for the 50 salt lakes of water level increase. The largest lake level increase rate (0.80 m/year) found in this study is the lake Cedo Caka. The 74 lakes are grouped into four subareas based on geographical locations and change tendencies in lake levels. Three of the four subareas show increased lake levels. The mean lake level change rates for subareas I, II, III, IV, and the entire TP are 0.12, 0.26, 0.19, -0.11, and 0.2 m/year, respectively. These recent increases in lake level, particularly for a high percentage of salt lakes, supports accelerated glacier melting due to global warming as the most likely cause.

Atmospheric, oceanic and hydrological polar motion excitation mechanisms based on a combination of geometric and gravimetric space observations

The goal of our study is to determine accurate time series of geophysical Earth rotation excitations to learn more about global dynamic processes in the Earth system. For this purpose, we developed an adjustment model which allows to combine precise observations from space geodetic observation systems, such as Satellite Laser Ranging (SLR), Global Navigation Satellite Systems (GNSS), Very Long Baseline Interferometry (VLBI), Doppler Orbit determination and Radiopositioning Integrated on Satellite (DORIS), satellite altimetry and satellite gravimetry in order to separate geophysical excitation mechanisms of Earth rotation. Three polar motion time series are applied to derive the polar motion excitation functions (integral effect). Furthermore we use five time variable gravity field solutions from Gravity Recovery and Climate Experiment (GRACE) to determine not only the integral mass effect but also the oceanic and hydrological mass effects by applying suitable filter techniques and a land-ocean mask. For comparison the integral mass effect is also derived from degree 2 potential coefficients that are estimated from SLR observations. The oceanic mass effect is also determined from sea level anomalies observed by satellite altimetry by reducing the steric sea level anomalies derived from temperature and salinity fields of the oceans. Due to the combination of all geodetic estimated excitations the weaknesses of the individual processing strategies can be reduced and the technique-specific strengths can be accounted for. The formal errors of the adjusted geodetic solutions are smaller than the RMS differences of the geophysical model solutions. The improved excitation time series can be used to improve the geophysical modeling.

Improved geoid solution for the Weddell Sea region

We present a geoid solution for the Weddell Sea and adjacent continental Antarctic regions. There, a refined geoid is of interest, especially for oceanographic and glaciological applications. For example, to investigate the Weddell Gyre as a part of the Antarctic Circumpolar Current and, thus, of the global ocean circulation, the mean dynamic topography (MDT) is needed. These days, the marine gravity field can be inferred with high and homogeneous resolution from altimetric height profiles of the mean sea surface. However, in areas permanently covered by sea ice as well as in coastal regions, satellite altimetry features deficiencies. Focussing on the Weddell Sea, these aspects are investigated in detail. In these areas, ground-based data that have not been used for geoid computation so far provide additional information in comparison with the existing high-resolution global gravity field models such as EGM2008. The geoid computation is based on the remove-compute-restore approach making use of least-squares collocation. The residual geoid with respect to a release 4 GOCE model adds up to two meters and more in the near-coastal and continental areas of the Weddell Sea region, also in comparison with EGM2008. Consequently, the thus refined geoid serves to compute new estimates of the regional MDT and geostrophic currents.

Geochemical and isotopic composition of basalts from the axial Mid-Atlantic Ridge, South Atlantic

One of the essential problems of oceanic tectonics is estimation of the influence of plumes of the deep hot mantle on processes in the axial spreading zone. Areas of two giant (St. Helena and Tristan da Cunha) plumes in the Mid-Atlantic Ridge (MAR) rift zone (South Atlantic) are characterized by the effusion of basalts that differ from typical depleted riftogenic tholeiites by anomalously high contents of lithophile components and specific isotopic compositions. Moreover, the rift valley floor with basalt effusion is significantly uplifted above the adjacent sectors of the rift. The formation of the St. Helena Seamount located 400 km east of the MAR axis is related to magmatism that is active to this day. St. Helena Island is a member of the structural ensemble of large volcanic seamounts (Bonaparte, Bagration, and Kutuzov). Like St. Helena Island, each seamount incorporates a series of smaller rises of different morphologies and dimensions. Thus, a system of subparallel series of NE-trending (~45°) rises extend from the seamount ensemble to the African continent. According to the plate tectonics concept, the seamount series represent hotspots related to a deep mantle plume that can be projected onto the present-day St. Helena Island area (St. Helena plume). At the same time, the inferred topographic map based on satellite altimetry data shows that the seamount series also extend along the opposite southwestern direction (~225°) toward the axial MAR and even intersect the latter structure. This fact cannot be explained by the hotspot hypothesis, which suggests stationary positions of plumes relative to the mobile oceanic plate. In the course of Cruise 10 of the R/V Akademik Ioffe (2002), detailed geological and geophysical investigations were carried out at the junction of one structural series with the MAR rift zone located near the Martin Vaz Fracture Zone (Martin Vaz test area, 19°-20° S). The present communication is devoted to the study of lithology, geochemistry, and isotopy of basalts dredged at the test area.

(Table 1) Summary statistics of in situ sea-ice characteristics of Ice Station Belgica sites Fabra, Patria/Liège and Brussels in October 2007

Three ice type regimes at Ice Station Belgica (ISB), during the 2007 International Polar Year SIMBA (Sea Ice Mass Balance in Antarctica) expedition, were characterized and assessed for elevation, snow depth, ice freeboard and thickness. Analyses of the probability distribution functions showed great potential for satellite-based altimetry for estimating ice thickness. In question is the required altimeter sampling density for reasonably accurate estimation of snow surface elevation given inherent spatial averaging. This study assesses an effort to determine the number of laser altimeter 'hits' of the ISB floe, as a representative Antarctic floe of mixed first- and multi-year ice types, for the purpose of statistically recreating the in situ-determined ice-thickness and snow depth distribution based on the fractional coverage of each ice type. Estimates of the fractional coverage and spatial distribution of the ice types, referred to as ice 'towns', for the 5 km**2 floe were assessed by in situ mapping and photo-visual documentation. Simulated ICESat altimeter tracks, with spot size ~70 m and spacing ~170 m, sampled the floe's towns, generating a buoyancy-derived ice thickness distribution. 115 altimeter hits were required to statistically recreate the regional thickness mean and distribution for a three-town assemblage of mixed first- and multi-year ice, and 85 hits for a two-town assemblage of first-year ice only: equivalent to 19.5 and 14.5 km respectively of continuous altimeter track over a floe region of similar structure. Results have significant implications toward model development of sea-ice sampling performance of the ICESat laser altimeter record as well as maximizing sampling characteristics of satellite/airborne laser and radar altimetry missions for sea-ice thickness.

Topographic model of Greenland (elevation lines in PanMap and ASCII format)

In this paper, a new high-resolution elevation model of Greenland, including the ice sheet as well as the ice free regions, is presented. It is the first published full coverage model, computed with an average resolution of 2 km and providing an unprecedented degree of detail. The topography is modeled from a wide selection of data sources, including satellite radar altimetry from Geosat and ERS 1, airborne radar altimetry and airborne laser altimetry over the ice sheet, and photogrammetric and manual map scannings in the ice free region. The ice sheet model accuracy is evaluated by omitting airborne laser data from the analysis and treating them as ground truth observations. The mean accuracy of the ice sheet elevations is estimated to be 12-13 m, and it is found that on surfaces of a slope between 0.2° and 0.8°, corresponding to approximately 50% of the ice sheet, the model presents a 40% improvement over models based on satellite altimetry alone. On coastal bedrock, the model is compared with stereo triangulated reference points, and it is found that the model accuracy is of the order of 25-35 m in areas covered by stereo photogrammetry scannings and between 200 and 250 m elsewhere.

Hydrography and biogeochemistry during S.A. Agulhas II cruise Agu011, in the SW Indian Ocean, April 2014

Several open-ocean mesoscale features, a "young" warm-core (or anti-cyclonic) eddy at 52°S, an "older" warm-core eddy at 57.5°S, as well as an adjacent cold-core (or cyclonic) eddy at 56°S, were surveyed during a R/V S.A. Agulhas II cruise in April 2014. The main aim of the survey was to obtain hydrographical and biogeochemical profile data for contrasting open-ocean eddies in the Southern Ocean, that will be suitable for comparison and modelling of their heat, salt and nutrient characteristics, and the changes that occur in these properties as warm-core eddies migrate from the polar front southwards into the Southern Ocean. A total of 18 CTD stations were occupied in a sector south of the South-West Indian Ridge, along three transects crossing several mesoscale features identified from satellite altimetry data prior to the cruise.

Comparación de la precisión de los Modelos Geopotenciales EGM96 y EGM08 en la zona del Caribe.

Recientemente fue publicado un nuevo modelo geopotencial global, el EGM08. Este modelo ha mostrado una notable mejoría en la calidad de sus tres fuentes de datos; las observaciones del movimiento perturbado de los satélites artificiales,la altimetría por satélite y la gravimetría terrestre, por lo que se ha conseguido mejorar su precisión. En Puerto Rico, nuestra área de estudio, encontramos que al comprar las diferencias de los incrementos de la ondulación del geoide geométrico (calculado con medidas de campo) con los valores de los incrementos de la ondulación del geoide obtenidos utilizando estos modelos geopotenciales globales, la precisión del EGM08 fue ± 0,029 metros mientras que la precisión del EGM96 fue ± 0,055 metros. Estos resultados demuestran que en nuestra región, el modelo EGM08 ha presentado una mejoría considerable sobre su predecesor el EGM96 al momento de determinar los valores de los incrementos de la ondulación del geoide. Abstract: Recently, the new global geopotential model, the EGM08 was published. This model has shown a marked improvement in the quality of its three sources of data; the observations of the disturbed motion of artificial satellites, satellite altimetry and terrestrial gravity, so it has improved its precision. In our study area, Puerto Rico, we found that when we compare the differences of the increments of the geometric geoid undulation (computed with field data) with the values of the increments of the geoid undulation obtained using these models, the EGM08 accuracy was ± 0,029 meters, while the EGM96 accuracy was ± 0,055 meters. These results confirm that in our region, the EGM08 model has presented a significant improvement over its predecessor the EGM96 when determining the values of the increments of the geoid undulation.

Determinación de un modelo de geoide gravimétrico para Puerto Rico como sistema de referencia para las altitudes ortométricas

El geoide, definido como la superficie equipotencial que mejor se ajusta (en el sentido de los mínimos cuadrados) al nivel medio del mar en una determinada época, es la superficie que utilizamos como referencia para determinar las altitudes ortométricas. Si disponemos de una superficie equipotencial de referencia como dátum altimétrico preciso o geoide local, podemos entonces determinar las altitudes ortométricas de forma eficiente a partir de las altitudes elipsoidales proporcionadas por el Sistema Global de Navegación por Satélite (Global Navigation Satellite System, GNSS ). Como es sabido uno de los problemas no resueltos de la geodesia (quizás el más importante de los mismos en la actualidad) es la carencia de un dátum altimétrico global (Sjoberg, 2011) con las precisiones adecuadas. Al no existir un dátum altimétrico global que nos permita obtener los valores absolutos de la ondulación del geoide con la precisión requerida, es necesario emplear modelos geopotenciales como alternativa. Recientemente fue publicado el modelo EGM2008 en el que ha habido una notable mejoría de sus tres fuentes de datos, por lo que este modelo contiene coeficientes adicionales hasta el grado 2190 y orden 2159 y supone una sustancial mejora en la precisión (Pavlis et al., 2008). Cuando en una región determinada se dispone de valores de gravedad y Modelos Digitales del Terreno (MDT) de calidad, es posible obtener modelos de superficies geopotenciales más precisos y de mayor resolución que los modelos globales. Si bien es cierto que el Servicio Nacional Geodésico de los Estados Unidos de América (National Geodetic Survey, NGS) ha estado desarrollando modelos del geoide para la región de los Estados Unidos de América continentales y todos sus territorios desde la década de los noventa, también es cierto que las zonas de Puerto Rico y las Islas Vírgenes Estadounidenses han quedado un poco rezagadas al momento de poder aplicar y obtener resultados de mayor precisión con estos modelos regionales del geoide. En la actualidad, el modelo geopotencial regional vigente para la zona de Puerto Rico y las Islas Vírgenes Estadounidenses es el GEOID12A (Roman y Weston, 2012). Dada la necesidad y ante la incertidumbre de saber cuál sería el comportamiento de un modelo del geoide desarrollado única y exclusivamente con datos de gravedad locales, nos hemos dado a la tarea de desarrollar un modelo de geoide gravimétrico como sistema de referencia para las altitudes ortométricas. Para desarrollar un modelo del geoide gravimétrico en la isla de Puerto Rico, fue necesario implementar una metodología que nos permitiera analizar y validar los datos de gravedad terrestre existentes. Utilizando validación por altimetría con sistemas de información geográfica y validación matemática por colocación con el programa Gravsoft (Tscherning et al., 1994) en su modalidad en Python (Nielsen et al., 2012), fue posible validar 1673 datos de anomalías aire libre de un total de 1894 observaciones obtenidas de la base de datos del Bureau Gravimétrico Internacional (BGI). El aplicar estas metodologías nos permitió obtener una base de datos anomalías de la gravedad fiable la cual puede ser utilizada para una gran cantidad de aplicaciones en ciencia e ingeniería. Ante la poca densidad de datos de gravedad existentes, fue necesario emplear un método alternativo para densificar los valores de anomalías aire libre existentes. Empleando una metodología propuesta por Jekeli et al. (2009b) se procedió a determinar anomalías aire libre a partir de los datos de un MDT. Estas anomalías fueron ajustadas utilizando las anomalías aire libre validadas y tras aplicar un ajuste de mínimos cuadrados por zonas geográficas, fue posible obtener una malla de datos de anomalías aire libre uniforme a partir de un MDT. Tras realizar las correcciones topográficas, determinar el efecto indirecto de la topografía del terreno y la contribución del modelo geopotencial EGM2008, se obtuvo una malla de anomalías residuales. Estas anomalías residuales fueron utilizadas para determinar el geoide gravimétrico utilizando varias técnicas entre las que se encuentran la aproximación plana de la función de Stokes y las modificaciones al núcleo de Stokes, propuestas por Wong y Gore (1969), Vanicek y Kleusberg (1987) y Featherstone et al. (1998). Ya determinados los distintos modelos del geoide gravimétrico, fue necesario validar los mismos y para eso se utilizaron una serie de estaciones permanentes de la red de nivelación del Datum Vertical de Puerto Rico de 2002 (Puerto Rico Vertical Datum 2002, PRVD02 ), las cuales tenían publicados sus valores de altitud elipsoidal y elevación. Ante la ausencia de altitudes ortométricas en las estaciones permanentes de la red de nivelación, se utilizaron las elevaciones obtenidas a partir de nivelación de primer orden para determinar los valores de la ondulación del geoide geométrico (Roman et al., 2013). Tras establecer un total de 990 líneas base, se realizaron dos análisis para determinar la 'precisión' de los modelos del geoide. En el primer análisis, que consistió en analizar las diferencias entre los incrementos de la ondulación del geoide geométrico y los incrementos de la ondulación del geoide de los distintos modelos (modelos gravimétricos, EGM2008 y GEOID12A) en función de las distancias entre las estaciones de validación, se encontró que el modelo con la modificación del núcleo de Stokes propuesta por Wong y Gore presentó la mejor 'precisión' en un 91,1% de los tramos analizados. En un segundo análisis, en el que se consideraron las 990 líneas base, se determinaron las diferencias entre los incrementos de la ondulación del geoide geométrico y los incrementos de la ondulación del geoide de los distintos modelos (modelos gravimétricos, EGM2008 y GEOID12A), encontrando que el modelo que presenta la mayor 'precisión' también era el geoide con la modificación del núcleo de Stokes propuesta por Wong y Gore. En este análisis, el modelo del geoide gravimétrico de Wong y Gore presento una 'precisión' de 0,027 metros en comparación con la 'precisión' del modelo EGM2008 que fue de 0,031 metros mientras que la 'precisión' del modelo regional GEOID12A fue de 0,057 metros. Finalmente podemos decir que la metodología aquí presentada es una adecuada ya que fue posible obtener un modelo del geoide gravimétrico que presenta una mayor 'precisión' que los modelos geopotenciales disponibles, incluso superando la precisión del modelo geopotencial global EGM2008. ABSTRACT The geoid, defined as the equipotential surface that best fits (in the least squares sense) to the mean sea level at a particular time, is the surface used as a reference to determine the orthometric heights. If we have an equipotential reference surface or a precise local geoid, we can then determine the orthometric heights efficiently from the ellipsoidal heights, provided by the Global Navigation Satellite System (GNSS). One of the most common and important an unsolved problem in geodesy is the lack of a global altimetric datum (Sjoberg, 2011)) with the appropriate precision. In the absence of one which allows us to obtain the absolute values of the geoid undulation with the required precision, it is necessary to use alternative geopotential models. The EGM2008 was recently published, in which there has been a marked improvement of its three data sources, so this model contains additional coefficients of degree up to 2190 and order 2159, and there is a substantial improvement in accuracy (Pavlis et al., 2008). When a given region has gravity values and high quality digital terrain models (DTM), it is possible to obtain more accurate regional geopotential models, with a higher resolution and precision, than global geopotential models. It is true that the National Geodetic Survey of the United States of America (NGS) has been developing geoid models for the region of the continental United States of America and its territories from the nineties, but which is also true is that areas such as Puerto Rico and the U.S. Virgin Islands have lagged behind when to apply and get more accurate results with these regional geopotential models. Right now, the available geopotential model for Puerto Rico and the U.S. Virgin Islands is the GEOID12A (Roman y Weston, 2012). Given this need and given the uncertainty of knowing the behavior of a regional geoid model developed exclusively with data from local gravity, we have taken on the task of developing a gravimetric geoid model to use as a reference system for orthometric heights. To develop a gravimetric geoid model in the island of Puerto Rico, implementing a methodology that allows us to analyze and validate the existing terrestrial gravity data is a must. Using altimetry validation with GIS and mathematical validation by collocation with the Gravsoft suite programs (Tscherning et al., 1994) in its Python version (Nielsen et al., 2012), it was possible to validate 1673 observations with gravity anomalies values out of a total of 1894 observations obtained from the International Bureau Gravimetric (BGI ) database. Applying these methodologies allowed us to obtain a database of reliable gravity anomalies, which can be used for many applications in science and engineering. Given the low density of existing gravity data, it was necessary to employ an alternative method for densifying the existing gravity anomalies set. Employing the methodology proposed by Jekeli et al. (2009b) we proceeded to determine gravity anomaly data from a DTM. These anomalies were adjusted by using the validated free-air gravity anomalies and, after that, applying the best fit in the least-square sense by geographical area, it was possible to obtain a uniform grid of free-air anomalies obtained from a DTM. After applying the topographic corrections, determining the indirect effect of topography and the contribution of the global geopotential model EGM2008, a grid of residual anomalies was obtained. These residual anomalies were used to determine the gravimetric geoid by using various techniques, among which are the planar approximation of the Stokes function and the modifications of the Stokes kernel, proposed by Wong y Gore (1969), Vanicek y Kleusberg (1987) and Featherstone et al. (1998). After determining the different gravimetric geoid models, it was necessary to validate them by using a series of stations of the Puerto Rico Vertical Datum of 2002 (PRVD02) leveling network. These stations had published its values of ellipsoidal height and elevation, and in the absence of orthometric heights, we use the elevations obtained from first - order leveling to determine the geometric geoid undulation (Roman et al., 2013). After determine a total of 990 baselines, two analyzes were performed to determine the ' accuracy ' of the geoid models. The first analysis was to analyze the differences between the increments of the geometric geoid undulation with the increments of the geoid undulation of the different geoid models (gravimetric models, EGM2008 and GEOID12A) in function of the distance between the validation stations. Through this analysis, it was determined that the model with the modified Stokes kernel given by Wong and Gore had the best 'accuracy' in 91,1% for the analyzed baselines. In the second analysis, in which we considered the 990 baselines, we analyze the differences between the increments of the geometric geoid undulation with the increments of the geoid undulation of the different geoid models (gravimetric models, EGM2008 and GEOID12A) finding that the model with the highest 'accuracy' was also the model with modifying Stokes kernel given by Wong and Gore. In this analysis, the Wong and Gore gravimetric geoid model presented an 'accuracy' of 0,027 meters in comparison with the 'accuracy' of global geopotential model EGM2008, which gave us an 'accuracy' of 0,031 meters, while the 'accuracy ' of the GEOID12A regional model was 0,057 meters. Finally we can say that the methodology presented here is adequate as it was possible to obtain a gravimetric geoid model that has a greater 'accuracy' than the geopotential models available, even surpassing the accuracy of global geopotential model EGM2008.

Relation of Lagragian structures and drifter dynamics in the Gulf of Mexico

We use a Lagrangian descriptor (the so called function M) which measures the length of particle trajectories on the ocean surface over a given interval of time. With this tool we identify the Lagrangian skeleton of the flow and compare it on three datasets over the Gulf of Mexico during the year 2010. The satellite altimetry data used come from AVISO and simulations from HYCOM GOMl0.04 experiments 30.1 and 31.0. We contrast the Lagrangian structure and transport using the evolution of several surface drifters. We show that the agreement in relevant cases between Lagrangian structures and dynamics of drifters depends on the quality of the data on the studied area.

Steric and mass-induced sea level variations in the Mediterranean Sea revisited

The total sea level variation (SLV) is the combination of steric and mass␣induced SLV, whose exact shares are key to understanding the oceanic response to climate system changes. Total SLV can be observed by radar altimetry satellites such as TOPEX/POSEIDON and Jason 1/2. The steric SLV can be computed through temperature and salinity profiles from in situ measurements or from ocean general circulation models (OGCM), which can assimilate the said observations. The mass-induced SLV can be estimated from its time-variable gravity (TVG) signals. We revisit this problem in the Mediterranean Sea estimating the observed, steric, and mass-induced SLV, for the latter we analyze the latest TVG data set from the GRACE (Gravity Recovery and Climate Experiment) satellite mission launched in 2002, which is 3.5 times longer than in previous studies, with the application of a two-stage anisotropic filter to reduce the noise in high-degree and -order spherical harmonic coefficients. We confirm that the intra-annual total SLV are only produced by water mass changes, a fact explained in the literature as a result of the wind field around the Gibraltar Strait. The steric SLV estimated from the residual of “altimetry minus GRACE” agrees in phase with that estimated from OGCMs and in situ measurements, although showing a higher amplitude. The net water fluxes through both the straits of Gibraltar and Sicily have also been estimated accordingly.

On the steric and mass-induced contributions to the annual sea level variations in the Mediterranean Sea

The sea level variation (SLVtotal) is the sum of two major contributions: steric and mass-induced. The steric SLVsteric is that resulting from the thermal and salinity changes in a given water column. It only involves volume change, hence has no gravitational effect. The mass-induced SLVmass, on the other hand, arises from adding or subtracting water mass to or from the water column and has direct gravitational signature. We examine the closure of the seasonal SLV budget and estimate the relative importance of the two contributions in the Mediterranean Sea as a function of time. We use ocean altimetry data (from TOPEX/Poseidon, Jason 1, ERS, and ENVISAT missions) to estimate SLVtotal, temperature, and salinity data (from the Estimating the Circulation and Climate of the Ocean ocean model) to estimate SLVsteric, and time variable gravity data (from Gravity Recovery and Climate Experiment (GRACE) Project, April 2002 to July 2004) to estimate SLVmass. We find that the annual cycle of SLVtotal in the Mediterranean is mainly driven by SLVsteric but moderately offset by SLVmass. The agreement between the seasonal SLVmass estimations from SLVtotal – SLVsteric and from GRACE is quite remarkable; the annual cycle reaches the maximum value in mid-February, almost half a cycle later than SLVtotal or SLVsteric, which peak by mid-October and mid-September, respectively. Thus, when sea level is rising (falling), the Mediterranean Sea is actually losing (gaining) mass. Furthermore, as SLVmass is balanced by vertical (precipitation minus evaporation, P–E) and horizontal (exchange of water with the Atlantic, Black Sea, and river runoff) mass fluxes, we compared it with the P–E determined from meteorological data to estimate the annual cycle of the horizontal flux.

Reply to comment by L. Fenoglio-Marc et al. on “On the steric and mass-induced contributions to the annual sea level variations in the Mediterranean Sea”

Reply to comment by L. Fenoglio-Marc et al. on “On the steric and mass-induced contributions to the annual sea level variations in the Mediterranean Sea”.