Monitoring of long-term land subsidence from 2003 to 2017 in coastal area of Semarang,Indonesia by SBAS DInSAR analyses using Envisat-ASAR,ALOS-PALSAR,and Sentinel-1A SAR data

Land subsidence is a critical issue that large cities located in coastal areas, such as Semarang, Indonesia, must address. The monitoring of land subsidence is vital for predicting and mitigating the disasters that such subsidence may cause. Therefore, an economical and effective monitoring method, which can continuously provide accurate measurements over extensive areas, is highly required. Differential Interferometry Synthetic Aperture Radar (DInSAR) has the potential to be a powerful technique that can meet the above demands. Actually, DInSAR has been applied to monitor the subsidence in Semarang, but it was for a limited period before 2012.In order to clarify the transition of the long-term subsidence behavior in Semarang, the Small Baseline Subset (SBAS) method, which is one type of time-series DInSAR, is employed in this research. The sets of data of Envisat-ASAR (2003–2007), ALOS-PALSAR (2007–2011), and Sentinel-1A (2015–2017) are employed for the analyses. Then, the validity of the SBAS results is discussed from the viewpoints of both spatial distribution and temporal transition using GPS displacement measurement results and the geological conditions of the ground.On the other hand, as the lifespan of SAR satellites is commonly designed to be around 5–7?years, an appropriate method, which can connect the subsidence provided independently by the unlinked time-series data sets of the three different SAR satellite data, is required. This study uses the Hyperbolic Method (HM) to connect the above unlinked SBAS results. The HM is often used to fit the monitored subsidence in practice as a geotechnical engineering tool. Using this method, 14?years of the temporal behavior of the subsidence in Semarang is evaluated.It is found that the transition of the subsidence is different depending on the location, and that the subsidence rate is still increasing in the north and northeast parts of the coastal area. This study shows that SBAS DInSAR can be a useful tool for long-term continuous subsidence monitoring.     

RS-GIS based morphometrical and geological multi-criteria approach to the landslide susceptibility mapping in Gish River Basin,West Bengal,India

Darjeeling Himalaya is one of the several mountainous areas of India which is often suffered from landslide hazards. In this paper, a multi criteria evaluation is applied using 16 morphometric indicators, geology and lineaments to identify the areas vulnerable in respect to drainage and relief conditions. As both drainage and relief parameters exert strong influences on landslide intensity, both the diversity maps are integrated for final landslide susceptibility mapping. The obtained results show that 20.17?sq.?km (7.61%) area within the basin is highly susceptible for landslides, where average drainage density is 3.78?km/sq.?km, relative relief is greater than 408?m and slope is greater than 12°. The validation result shows that very high landslide susceptible zone is associated with very high frequency of landslide occurrence. Beside this, ROC curve also suggests good predicted rate (86.60%) for the model. So, the proposed method can be applied for predicting landslide susceptible zone.     

Estimates of vertical land motion along the southwestern coasts of Turkey from coastal altimetry and tide gauge data

The differences between coastal altimetry and sea level time series of tide gauges in between March 1993 and December 2009 are used to estimate the rates of vertical land motion at three tide gauge locations along the southwestern coasts of Turkey. The CTOH/LEGOS along-track coastal altimetry retrieves altimetric sea level anomalies closer to the coast than the standard along-track altimetry products. However, the use of altimetry very close to the coast is not found to improve the results. On the contrary, the gridded and interpolated AVISO merged product exhibits the best agreement with tide gauge data as it provides the smoothest variability both in space and time compared with along track altimetry data. The Antalya gauge to the south (in the Mediterranean Sea) and the Mentes/Izmir gauge to the west (in the Aegean Sea) both show subsidence while the Bodrum tide gauge to the south (in the Aegean Sea) shows no significant vertical land motion. The results are compared and assessed with three independent geophysical vertical land motion estimates like from GPS. The GIA effect in the region is negligible. The VLM estimates from altimetry and tide gauge data are in good agreement both with GPS derived vertical velocity estimates and those inferred from geological and archaeological investigations.     

Sea level trend over Malaysian seas from multi-mission satellite altimetry and vertical land motion corrected tidal data

Rise in sea levels is one of the disastrous effects of climate change. A relatively small increase in sea level could affect natural coastal systems. In a study of long-term changes in sea level and measurements of postglacial rebound, monitoring vertical land motion (VLM) is of crucial interest. This study presents an approach to estimate precise sea level trends based on a combination of multi-sensor techniques in the Malaysian region over 19?years. In this study, satellite altimeters (SALT) were used to derive absolute sea levels (ASLs). Tide gauge (TG) stations along the coast of Malaysia were utilised to derive the rate of relative sea levels using sea level changes and VLMs. To obtain ASL at TGs, VLM at these stations were computed using Global Positioning System (GPS), Persistent Scatterer Interferometric Synthetic Aperture Radar (PS InSAR), and SALT minus TG. The computed VLMs mostly show similarities in signs rather than magnitude. The findings from the multi-sensor techniques showed that regional sea level trends ranged from 2.65?±?0.86?mm/yr to 6.03?±?0.79?mm/yr for chosen sub-areas, with an overall mean of 4.47?±?0.71?mm/yr and overall subsidence. This information is expected to be valuable for a wide variety of climatic applications and for studying environmental issues related to flooding and global warming in Malaysia.     

Influence of land surface parameters on the spatio-seasonal land surface temperature regime in rural West Bengal,India

The magnitude and causes of changes in the land surface temperature of rural areas have not been extensively studied. The thermal band of Landsat imagery is taken to extract winter, summer, and monsoon season land surface temperature (LST) and relate it to surface parameters over a 30-year period. From the extracted parameters constructed a prospective surface temperature (PST) model using Multivariate Adaptive Regression Splines. The Chandrabhaga river basin in West Bengal of the lateritic Rarh Tract at the Chota Nagpur Plateau fringe was chosen as the study area because it is far from urban influences, to avoid the well-known heat island effect. Over the study period, summer and winter average LST increased linearly by 0.085?°C/y and 0.016?°C/y respectively. These results were validated with air temperature (RMSE?=?x and y, respectively). Over time more of the area is in the higher temperature zones, e.g., in April 2011, 4% area exceeded >32°, whereas in 2015 this proportion reached 52%. PST models of all the seasons were moderate to highly correlate (0.57–0.87) with actual LST, showing the value of this model. It also revealed the relative importance of the regional factors. Based on this information factor management is a scientific step to restrict or minimize the temperature rise effect.     

Evaluation and improvement of coastal GNSS reflectometry sea level variations from existing GNSS stations in Taiwan

Global sea level rise due to an increasingly warmer climate has begun to induce hazards, adversely affecting the lives and properties of people residing in low-lying coastal regions and islands. Therefore, it is important to monitor and understand variations in coastal sea level covering offshore regions. Signal-to-noise ratio (SNR) data of Global Navigation Satellite System (GNSS) have been successfully used to robustly derive sea level heights (SLHs). In Taiwan, there are a number of continuously operating GNSS stations, not originally installed for sea level monitoring. They were established in harbors or near coastal regions for monitoring land motion. This study utilizes existing SNR data from three GNSS stations (Kaohsiung, Suao, and TaiCOAST) in Taiwan to compute SLHs with two methods, namely, Lomb–Scargle Periodogram (LSP)-only, and LSP aided with tidal harmonic analysis developed in this study. The results of both methods are compared with co-located or nearby tide gauge records. Due to the poor quality of SNR data, the worst accuracy of SLHs derived from traditional LSP-only method exceeds 1?m at the TaiCOAST station. With our procedure, the standard deviations (STDs) of difference between GNSS-derived SLHs and tide gauge records in Kaohsiung and Suao stations decreased to 10?cm and the results show excellent agreement with tide gauge derived relative sea level records, with STD of differences of 7?cm and correlation coefficient of 0.96. In addition, the absolute GNSS-R sea level trend in Kaohsiung during 2006–2011 agrees well with that derived from satellite altimetry. We conclude that the coastal GNSS stations in Taiwan have the potential of monitoring absolute coastal sea level change accurately when our proposed methodology is used.     

Seasonal,annual, and interannual variability in MLT quasi-two-day waves over the low-latitude region Kolhapur (16.8°N; 74.2°E)

This study presents the quasi-two-day wave (Q2DW) characteristics of the mesosphere and lower thermosphere (MLT) region obtained by taking hourly mean values of horizontal wind velocities for 4? years (August 2013–July 2017) through continuous measurements using a medium-frequency (MF) radar (operating frequency – 1.98?MHz) located at the low-latitude Indian station Kolhapur (16.8°N; 74.2°E). The MF radar located at Kolhapur was upgraded in 2013, and these results of Q2DW have been reported for the first time after upgrading. The present study investigated variability in seasonal, annual, interannual, and solar indices of Q2DWs traveling in zonal (EW) and meridional (NS) components in the MLT region. The Q2DW activity is observed to be stronger during austral summer (January–February) (EW?=?～5?m/s and NS?=?～8–10?m/s) than during boreal summer (June–July) (EW = ～5 m/s and NS = ～6–8?m/s). The Q2DW amplitudes are larger in the meridional component than in the zonal one. A strong semiannual oscillation (SAO) has been observed in Q2DWs, with peak during January–February and June–July. In addition, small enhancement is seen in meridional Q2DW in October (～5–6?m/s). It is observed that the entire spectrum (40–60?h) measured between 86 and 94?km contributes to the SAO amplitudes during January–February and June–July, whereas the waves measured between 42?h and 52?h contribute to enhancement in October similar to that reported elsewhere. In general, the Q2DW amplitude shows large interannual variability. The easterlies developed in the global circulation model in Northern hemisphere during May intensify up to around summer solstice. Q2DW activity peaks during westerly shear zone and intensifies with time at a lower thermospheric altitude (above 90?km). Small positive correlations (r?=?0.2 for sunspot number and r?=?0.1 for 10.7?cm solar flux) have been observed between Q2DW amplitudes and solar activity.     

Mesosphere/lower thermosphere wind regime parameters using a newly installed SKiYMET meteor radar at Kazan (56°N, 49°E)

New meteor radar (MR) horizontal wind data obtained during 2015–2018 at Kazan (56°N, 49°E) are presented. The measurements were carried out with a state-of-the-art SKiYMET meteor radar. Monthly mean vertical profiles of zonal and meridional components of the prevailing wind speeds, also amplitudes and phases of the components of diurnal (DT) and semidiurnal tide (SDT) winds are displayed as contour plots for a mean calendar year over the four recent years and compared with distributions of these parameters provided by the previous multiyear (1986–2002) meteor radar (MR) measurements at Kazan and by the recent HWM07 empirical model. The analysis shows that the SKiYMET zonal and meridional prevailing wind speeds are generally in good agreement, sharing the same seasonal features, with the earlier MR seasonal winds. Comparisons with the HWM07 model are not favourable: eastward solstitial cells as modelled are significantly larger, >30?m/s compared to 15–20?m/s. Also, reversal lines are too variable with height, and the positions of modelled cells (positive and negative) are unlike those of either MRs at Kazan or other MLT radars. Both MR systems provide the large SDT amplitudes, approximately 30?m/s and vertical wavelengths, approximately 55?km, for both components at middle latitudes in winter. They also show the well known strong SDT September feature (heights 85–100?km, the vertical wavelength ～55–60?km), and the weak summer SDT for 80–91?km. HWM07 shows unrealistic amplitudes and phases above 90?km by height and month: minimal amplitudes in equinoxes and no September feature.The weak DT of middle to high latitudes provide similar amplitude and phase structures from both MRs, 1986–2002 and 2015–2017: largest amplitudes (10–12 or 8–10?m/s) for the evanescent meridional tide in summer, peaking in late July; weakest (0–2, 2–4?m/s) at 80 to 92–96?km, when the tide is vertically propagating (January, February, November, December) with a vertical wavelength near 40?km. Again, HWM07 differs in amplitude and phase structures: showing peak amplitudes in equinoxes: April, 15?m/s at 88?km; October, 21?m/s at 89?km.Coupling of the MR wind parameters with the ERA5 wind parameters is studied for a case in 2016. It is shown that the prevailing winds and DT amplitudes and phases of both datasets can be simply linked together, but that the ERA5 SDT amplitudes are significantly underestimated at the top model levels of the ERA5 reanalysis project.     

Monitoring coastal sea level using reflected GNSS signals

A continuous monitoring of coastal sea level changes is important for human society since it is predicted that up to 332 million people in coastal and low-lying areas will be directly affected by flooding from sea level rise by the end of the 21st century. The traditional way to observe sea level is using tide gauges that give measurements relative to the Earth’s crust. However, in order to improve the understanding of the sea level change processes it is necessary to separate the measurements into land surface height changes and sea surface height changes. These measurements should then be relative to a global reference frame. This can be done with satellite techniques, and thus a GNSS-based tide gauge is proposed. The GNSS-based tide gauge makes use of both GNSS signals that are directly received and GNSS signals that are reflected from the sea surface. An experimental installation at the Onsala Space Observatory (OSO) shows that the reflected GNSS signals have only about 3 dB less signal-to-noise-ratio than the directly received GNSS signals. Furthermore, a comparison of local sea level observations from the GNSS-based tide gauge with two stilling well gauges, located approximately 18 and 33 km away from OSO, gives a pairwise root-mean-square agreement on the order of 4 cm. This indicates that the GNSS-based tide gauge gives valuable results for sea level monitoring.     

Determination of the mean dynamic ocean topography model through combining multi-source gravity data and DTU15 MSS around China's coast

Mean dynamic ocean topography (or MDT) is closely related to ocean circulation and global climate change. It has important scientific significance and application value for the development and utilization of marine resources in China's coastal areas. Based on the terrain gravity, marine gravity, and SRTM 3?s data, an algorithm to reduce the problem of gravity data gaps between land and sea is proposed. A consistent land-sea gravity model is established based on point-mass fusion method. Then geoid model, which accuracy was estimated to be 8.5?cm through the verification of 348 GNSS/level data from the coastal provinces, of China's coastal areas was calculated through remove-restore technique. Connecting the above geoid model with DTU15 MSS model to establish a MDT model in China's coastal areas using the direct method in space domain. The effect of gravity field model, dominant factors of sea surface topography, and low pass filter are analyzed. Taking Bohai Sea and Yellow Sea as an example, and comparing MDT with the two international models CNES_2013_MDT and DTU15_MDT. The results show that the MDT has the potential to construct a vertical datum of the ocean and carry out related scientific research and application.     

Response of the mid-latitude D-region ionosphere to the total solar eclipse of 22 July 2009 studied using VLF signals in South Korean peninsula

In this paper, we analyze VLF signals received at Busan to study the the D-region changes linked with the solar eclipse event of 22 July 2009 for very short (∼390 km) transmitter–receiver great circle path (TRGCP) during local noon time 00:36–03:13 UT (09:36–12:13 KST). The eclipse crossed south of Busan with a maximum obscuration of ∼84%. Observations clearly show a reduction of ∼6.2 dB in the VLF signal strength at the time of maximum solar obscuration (84% at 01:53 UT) as compared to those observed on the control days. Estimated values of change in Wait ionospheric parameters: reflection height (h′) in km and inverse scale height parameter (β) in km⁻¹ from Long Wave Propagation Capability (LWPC) model during the maximum eclipse phase as compared to unperturbed ionosphere are 7 km and 0.055 km⁻¹, respectively. Moreover, the D-region electron density estimated from model computation shows 95% depletion in electron density at the height of ∼71 km. The reflection height is found to increase by ∼7 km in the D-region during the eclipse as compared to those on the control days, implying a depletion in the Lyman-α flux by a factor of ∼7. The present observations are discussed in the light of current understanding on the solar eclipse induced D-region dynamics.     

Altimetry,bathymetry and geoid variations at the Gavdos permanent Cal/Val facility

The aim of this work has been to examine the relationship of steep bathymetry in the coastal areas around the permanent Cal/Val facility of Gavdos, and their influence on the produced calibration values for the Jason-2 satellite altimeter. The paper describes how changes in seafloor topography (from 200 to 3500 m depth over a distance of 10 km) are reflected on the determined altimeter parameters using different reference surfaces for satellite calibration. Finally, it describes the relation between these parameter trends and the region’s local characteristics.     

The X-TRACK/ALES multi-mission processing system: New advances in altimetry towards the coast

In the context of the ESA Climate Change Initiative project, a new coastal sea level altimetry product has been developed in order to support advances in coastal sea level variability studies. Measurements from Jason-1,2&3 missions have been retracked with the Adaptive Leading Edge Subwaveform (ALES) Retracker and then ingested in the X-TRACK software with the best possible set of altimetry corrections. These two coastal altimetry processing approaches, previously successfully validated and applied to coastal sea level research, are combined here for the first time in order to derive a 16-year-long (June 2002 to May 2018), high-resolution (20-Hz), along-track sea level dataset in six regions: Northeast Atlantic, Mediterranean Sea, West Africa, North Indian Ocean, Southeast Asia and Australia. The study demonstrates that this new coastal sea level product called X-TRACK/ALES is able to extend the spatial coverage of sea level altimetry data ~3.5 km in the land direction, when compared to the X-TRACK 1-Hz dataset. We also observe a large improvement in coastal sea level data availability from Jason-1 to Jason-3, with data at 3.6 km, 1.9 km and 0.9 km to the coast on average, for Jason-1, Jason-2 and Jason-3, respectively. When combining measurements from Jason-1 to Jason-3, we reach a distance of 1.2–4 km to the coast. When compared to tide gauge data, the accuracy of the new altimetry near-shore sea level estimations also improves. In terms of correlations with a large set of independent tide gauge observations selected in the six regions, we obtain an average value of 0.77. We also show that it is now possible to derive from the X-TRACK/ALES product an estimation of the ocean current variability up to 5 km to the coast. This new altimetry dataset, freely available, will provide a valuable contribution of altimetry in coastal marine research community.     

Seismogenic ionospheric anomalies associated with the strong Indonesian earthquake occurred on 11 April 2012 (M = 8.5)

Ionospheric perturbations in possible association with a major earthquake (EQ) (M?=?8.5) which occurred in India-Oceania region are investigated by monitoring subionospheric propagation of VLF signals transmitted from the NWC transmitter (F?=?19.8?kHz), Australia to a receiving station at Varanasi (geographic lat. 25.3°N, long 82.99°E), India. The EQ occurred on 11 April 2012 at 08:38:35?h UT (magnitude?≈?8.5, depth?=?10?km, and lat.?=?2.3°N, long.?=?93.0°E). A significant increase of few days before the EQ has been observed by using the VLF nighttime amplitude fluctuation method (fixed frequency transmitter signal). The analysis of total electron contents (TEC) derived from the global positioning system (GPS) at three different stations namely, Hyderabad (latitude 17.38°N, longitude 78.48°E), Singapore (latitude 1.37°N, longitude 103.84°E) and Port Blair (latitude 11.62°N, longitude 92.72°E) due to this EQ has also been presented. Significant perturbation in TEC data (enhancements and depletion) is noted before and after the main shock of the EQ. The possible mechanisms behind these perturbations due to EQ have also been discussed.     

Schumann resonance frequency and conductivity in the nighttime ionosphere of Mars: A source for lightning

We have solved the Maxwellian equations of electromagnetic waves which oscillate within the cavity formed in the lower ionosphere of Mars between 0 and 70?km. The electrical conductivity and Schumann Resonance (SR) frequencies are calculated in the lower ionosphere of Mars, in the presence of a major dust storm that occurred in Martian Year (MY) 25 at low latitude region (25°–35°S). It is found that the atmospheric conductivity reduced by one to two orders of magnitude in the presence of a dust storm. It represents a small dust layer at about 25–30?km altitudes where lightning can occur. We also found that the SR frequencies peak at?～18?km with values 19.9, 34.5 and 48.8?Hz for the modes l?=?1, 2 and 3, respectively, in the non-homogeneous medium. Our results indicate that practical or measurable values of SR are dependent on the altitudes.     

Balloon-borne,high-altitude gravimetry

Gravity measurements from a high-altitude balloon can verify global and upward-continued gravity models. A gravimeter suspended beneath a balloon is in a dynamic, and largely unpredictable, environment sensing accelerations due to gravity and balloon motions. Independent measurements of balloon motions using inertial navigation data combined with ground tracking data will allow for separation of balloon-induced accelerations from gravitational accelerations. Analysis of these data must estimate: 1) vertical gravimeter accelerations due to motion and gravity, 2) horizontal velocity to estimate the Eötvös effect, and 3) gravimeter position for comparison with gravity models. The first engineering test flight occurred on 11 October 1983, during the seasonal wind reversal and was very successful. Flight duration was approximately seven hours, with two hours of data collected at each of 30 km and 26 km altitudes. The results include gravity estimates, design criteria for future flights and feasibility analysis for vertical gravity profiles during ascent and descent.     

Ionospheric total electron content response to annular solar eclipse on June 21, 2020

The effects of physical events on the ionosphere structure is an important field of study, especially for navigation and radio communication. The paper presents the spatio-temporal ionospheric TEC response to the recent annular solar eclipse on June 21, 2020, which spans across two continents, Africa and Asia, and 14 countries. This eclipse took place on the same day as the June Solstice. The Global Navigation Satellite System (GNSS) based TEC data of the Global Ionosphere Maps (GIMs), 9 International GNSS Service (IGS) stations and FORMOSAT-7/COSMIC-2 (F7/C2) were utilized to analyze TEC response during the eclipse. The phases of the TEC time series were determined by taking the difference of the observed TEC values on eclipse day from the previous 5-day median TEC values. The results showed clear depletions in the TEC time series on June 21. These decreases were between 1 and 9 TECU (15–60%) depending on the location of IGS stations. The depletions are relatively higher at the stations close to the path of annular eclipse than those farther away. Furthermore, a reduction of about ?10 TECU in the form of an equatorial plasma bubble (EPB) was observed in GIMs at ~20° away from the equator towards northpole, between 08:00–11:00 UT where its maximum phase is located in southeast Japan. Additionally, an overall depletion of ~10% was observed in F7/C2 derived TEC at an altitude of 240 km (hmF2) in all regions affected by the solar eclipse, whereas, significant TEC fluctuations between the altitudes of 100 km ? 140 km were analyzed using the Savitzky-Golay smoothing filter. To prove TEC depletions are not caused by space weather, the variation of the sunspot number (SSN), solar wind (V_SW), disturbance storm-time (Dst), and Kp indices were investigated from 16th to 22nd June. The quiet space weather before and during the solar eclipse proved that the observed depletions in the TEC time series and profiles were caused by the annular solar eclipse.     

TS-InSAR analysis for monitoring ground deformation in Lanzhou New District,the loess Plateau of China,from 2017 to 2019

Large-scale land creation projects involving the cutting of mountains to infill gullies for construction have been carried out in Lanzhou New District (LZND). However, there is an urgent need for comprehensive and detailed research on the spatiotemporal evolution of ground deformation in LZND. Based on Sentinel-1A SAR data, combined with the urban geological background, the ground deformation in LZND from 2017 to 2019 was analysed. Two independent, multi-temporal techniques, persistent scatterers interferometry (PS-InSAR) and the small baseline subset (SBAS-InSAR), were used to calculate the deformation time series, and the results were cross-verified. The time series-monitoring results of the SBAS and PS calculations exhibited strong consistency in LZND and verified the high reliability of the experimental results. The results showed the whole surface of the LZND from March 2017 to October 2019 maintained stability, and the deformation rate was primarily in the range of ?10 to 10 mm/year. However, ground deformation in the Xicha area was evident. The maximum annual deformation rates monitored by SBAS-InSAR and PS-InSAR were ?52.48 mm/year and ?56.35 mm/year, respectively. The most typical deformation areas include the built-up area and the land creation area. The surface subsidence area was concentrated in the filling area. The ground deformation range of LZND kept expanding and accelerating from 2017 to 2019. Land creation, urban construction, geology and precipitation were the primary factors contributing to local severe ground deformation. The results of this study provide reference for the regional urban planning in LZND.     

Analysis of coastal altimetry in the Mexican Caribbean

A reprocessing of sea-level anomalies (SLA) resulting from X-TRACK coastal altimetry was carried out for the ENVISAT (2002–2010) and TOPEX/POSEIDON-Jason (1992–2019) satellite missions in the coastal area of the Mexican Caribbean. This consisted of applying a tidal correction to coastal altimetry sea level observations. Harmonic analysis of five coastal tide gauge records was performed to estimate the most important tidal components of the area, resulting on M2, N2, O1, S2, K1, MF, and MM. The tidal signal was reconstructed with the seven tidal components using the TPXO9 model. The SLA signals corrected with the seven tidal components were validated with in situ data from coastal tide gauges. The validation showed that the TPXO9 tidal barotropic model (1/30° grid) used to reconstruct the tidal signal with the seven representative tidal components performed better than the FES2012 global model (1/16° grid) that uses 33 tidal components. The reprocessed SLAs showed clear seasonality with significant signals at 4, 6, and 12 months, with the annual signal being the dominant one. In the Mexican Caribbean coastal zone, oceanographic processes with different scales (from coastal to mesoscale) converge, showing their complexity in the different SLA signals observed. The aim of this work is to contribute to the analysis of coastal altimetry data and understanding the sea level variations in the Mexican Caribbean. This work is the first step in the implementation of methodologies that take advantage of coastal satellite altimetry in the Caribbean Sea.