中国沿海验潮站及其邻近地区陆地垂直运动分析

利用中国沿岸验潮站GNSS和邻近地区陆态网络GNSS基准站观测数据,结合卫星高度计和验潮站海平面观测数据分析了中国沿海验潮站及其邻近地区陆地垂直运动特征。中国沿海海平面观测以及验潮站和陆态网GNSS基准站观测结果显示,中国沿海省区市及沿海验潮站陆地垂直运动总体表现为：辽宁至江苏沿海上升、上海至福建泉州沿海沉降、福建厦门至广西沿海升降交替的格局,局部滨海平原地区如华北平原天津南部、河北平原的沧县则表现出显著的沉降特征。验潮站陆地的抬升与沉降是沿海相对海平面变化的重要组成部分,准确掌握验潮站及其邻近区域的陆地垂直运动特征,可为沿海相对海平面变化分析、海平面变化影响评估以及未来海平面上升预测提供依据。     

Human impacts overwhelm the effects of sea-level rise on Basque coastal habitats (N Spain) between 1954 and 2004

According to coastal measurements, global mean sea-level has risen at a rate of 1.8 mm yr⁻¹ between 1950 and 2000, with large spatial variability at regional scales. Within the Bay of Biscay, trends computed from coastal tide gauges records have revealed that sea-level rise is accelerating over this period of time; this is in agreement with rates obtained from satellite imagery in the open ocean since 1993. The objectives of the present study are: (1) to assess the evidence of the relative sea-level rise on coastal morphology and habitats in the Gipuzkoan littoral zone (Basque coast, northern Spain) for the period 1954–2004, and (2) to evaluate the relative contribution of local anthropogenic versus sea-level rise impacts for explaining inter-supratidal habitat changes. A high-resolution airborne laser altimetry data (LIDAR) has been used to derive a Digital Terrain Model (DTM) of 15-cm vertical resolution. Coastal habitats were mapped for two periods, using historic airborne photography (1954) and high-resolution imagery (2004). Analysis of tide gauge records from Santander (northern Spain) has revealed that relative mean sea-level has been rising at a rate of 2.08 ± 0.33 mm yr⁻¹ from 1943 to 2004; this is consistent with sea-level trends from other measurements within the area (St. Jean de Luz and Bilbao), obtained over shorter periods of time, and with previous results obtained in the Bay of Biscay. Based upon this sea-level trend and by means of a LIDAR-based DTM, the results have indicated that the predicted change along the Gipuzkoan coast due to sea-level rise was of 11.1 ha within the 50-yr period. In contrast, comparison of historical and recent orthophotography has detected only 2.95 ha of change, originated possibly from sea-level rise, and 98 ha transformed by anthropogenic impacts. Hence, coastal changes due to sea-level rise might be overwhelmed by excessive human impacts, at the spatial and temporal scales of the analysis. This work highlights that local anthropogenic impact is the major threat to Basque coastal and estuarine habitats, compared with natural erosive processes and global climate change driving forces over recent times.     

Absolute sea level variability of Arctic Ocean in 1993–2018 from satellite altimetry and tide gauge observations

Arctic absolute sea level variations were analyzed based on multi-mission satellite altimetry data and tide gauge observations for the period of 1993–2018. The range of linear absolute sea level trends were found ?2.00 mm/a to 6.88 mm/a excluding the central Arctic, positive trend rates were predominantly located in shallow water and coastal areas, and negative rates were located in high-latitude areas and Baffin Bay. Satellite-derived results show that the average secular absolute sea level trend was (2.53±0.42) mm/a in the Arctic region. Large differences were presented between satellite-derived and tide gauge results, which are mainly due to low satellite data coverage, uncertainties in tidal height processing and vertical land movement (VLM). The VLM rates at 11 global navigation satellite system stations around the Arctic Ocean were analyzed, among which 6 stations were tide gauge co-located, the results indicate that the absolute sea level trends after VLM corrected were of the same magnitude as satellite altimetry results. Accurately calculating VLM is the primary uncertainty in interpreting tide gauge measurements such that differences between tide gauge and satellite altimetry data are attributable generally to VLM.     

Model-based clustering of Baltic sea-level

Long (>30 years) monthly records of relative sea-level heights from tide gauges in the Baltic sea are analyzed. Time series clustering based on forecast densities is applied in order to describe regional sea-level variability in the Baltic Sea in terms of future relative heights. The tide gauge records are clustered on the basis of forecasts at 3-month and 6-month horizons. For the 3-month horizon, the results of the cluster analysis show a fairly spatial coherency in terms of grouping together locations from the same sub-basin, with the northern records in the Bothnian Sea and Gulf of Finland clustering together, followed by the tide gauges in the Baltic Proper and lastly the southernmost stations in the western Baltic. For the 6-month horizon, the results show a higher degree of homogeneity between different locations, but a clear separation between the stations at the Baltic entrance and the tide gauges inside the Baltic basin. Moreover, when considering detrended records, reflecting mainly the seasonal cycle, the clustering results are more homogeneous and suggest a distinct response of coastal sea-level in spring and in summer.     

Ku– and Ka-band Altimeter Data in the Northwestern Mediterranean Sea: Impact on the Observation of the Coastal Ocean Variability

The strong increase in altimeter measurement errors near land surfaces is a limiting factor for coastal applications. We analyze the performance of the new Ka-band SARAL/AltiKa (SRL) mission in the northwestern Mediterranean Sea. SRL sea surface height (SSH) measurements are compared with those from the Jason-2 Ku-band satellite mission. The results show a significant increase in both quantity and quality of SSH data available near coastlines when using SRL data. Available edited data are 95.1% of SRL compared with 88.6% for Jason-2. Closer than 10 km to the coastline, available SRL data are still about 60% and only about 31% for Jason-2. Comparisons of the altimeter sea level variations are made with available coastal tide gauge data. The differences obtained between altimeter and tide gauge SLA time series are reduced for SRL (3.3 cm in average) compared with Jason-2 (4.2 cm in average), especially closer than 30 km to the land. It results in higher correlations (by 30%) obtained with SRL data. The coastal circulation derived from altimetry using SRL data shows an offshore meandering, which is more stable in time and with larger velocities close to the coast than that derived from Jason-2 observations.     

Absolute Sea Level Surface Modeling for the Mediterranean from Satellite Altimeter and Tide Gauge Measurements

Changes in the height of the ocean can be described through the relative and absolute sea level changes depending on the geodetic reference the sea level records are related to. Satellite altimetry provides absolute sea level (ASL) measurements related to the global geodetic reference, whereas tide gauges provide relative sea level (RSL) measurements related to the adjacent land. This study aims at computing the ASL surfaces for different time epochs from combined satellite altimeter and tide gauge records. A method of sea level data fusion is proposed to enable modeling of the impact of present and future sea level changes on the coast. Sea surface modeling was investigated for ten different gridding methods commonly used for the interpolation of altimeter data over the open ocean and extrapolation over the coastal zones. The performance of gridding methods was assessed based on the comparison of the gridded altimeter data and corrected tide gauge measurements. Finally, the sea level surfaces related to the GRS80 global reference ellipsoid were computed for the Mediterranean Sea over the altimeter period. In addition, the current sea level trends were estimated from both sea level measurements.     

Inferring the global mean sea level from a global tide gauge network

An attempt is made to infer the global mean sea level(GMSL) from a global tide gauge network and frame the problem in terms of the limitations of the network. The network,owing to its limited number of gauges and poor geographical distribution complicated further by unknown vertical land movements,is ill suited for measuring the GMSL. Yet it remains the only available source for deciphering the sea level rise over the last 100 a. The poor sampling characteristics of the tide gauge network have necessitated the usage of statistical inference. A linear optimal estimator based on the Gauss-Markov theorem seems well suited for the job. This still leaves a great deal of freedom in choosing the estimator. GMSL is poorly correlated with tide gauge measurements because the small uniform rise and fall of sea level are masked by the far larger regional signals. On the other hand,a regional mean sea level(RMSL) is much better correlated with the corresponding regional tide gauge measurements. Since the GMSL is simply the sum of RMSLs,the problem is transformed to one of estimating the RMSLs from regional tide gauge measurements. Specifically for the annual heating and cooling cycle,we separate the global ocean into 10-latitude bands and compute for each 10-latitude band the estimator that predicts its RMSL from tide gauges within. In the future,the statistical correlations are to be computed using satellite altimetry. However,as a first attempt,we have used numerical model outputs instead to isolate the problem so as not to get distracted by altimetry or tide gauge errors. That is,model outputs for sea level at tide gauge locations of the GLOSS network are taken as tide gauge measurements,and the RMSLs are computed from the model outputs. The results show an estimation error of approximately 2 mm versus an error of 2.7 cm if we simply average the tide gauge measurements to estimate the GMSL,caused by the much larger regional seasonal cycle and mesoscale variation plaguing the individual tide gauges. The numerical model,Los Alamos POP model Run 11 lasting 3 1/4 a,is one of the best eddy-resolving models and does a good job simulating the annual heating and cooling cycle,but it has no global or regional trend. Thus it has basically succeeded in estimating the seasonal cycle of the GMSL. This is still going to be the case even if we use the altimetry data because the RMSLs are dominated by the seasonal cycle in relatively short periods. For estimating the GMSL trend,longer records and low-pass filtering to isolate the statistical relations that are of interest. Here we have managed to avoid the much larger regional seasonal cycle plaguing individual tide gauges to get a fairly accurate estimate of the much smaller seasonal cycle in the GMSL so as to enhance the prospect of an accurate estimate of GMSL trend in short periods. One should reasonably expect to be able to do the same for longer periods during which tide gauges are plagued by much larger regional interannual(e. g.,ENSO events) and decadal sea level variations. In the future,with the availability of the satellite altimeter data,we could use the same approach adopted here to estimate the seasonal variations of GMSL and RMSL accurately and remove these seasonal variations accordingly so as to get a more accurate statistical inference between the tide gauge data and the RMSLs(therefore the GMSL) at periods longer than 1 a,i. e.,the long-term trend.     

Improving the Coastal Mean Dynamic Topography by Geodetic Combination of Tide Gauge and Satellite Altimetry

Abstract

The ocean mean dynamic topography (MDT) is the surface representation of the ocean circulation. The MDT may be determined by the ocean approach, which involves temporal averaging of numerical ocean circulation model information, or by the geodetic approach, wherein the MDT is derived using the ellipsoidal height of the mean sea surface (MSS), or mean sea level (MSL) minus the geoid as the geoid. The ellipsoidal height of the MSS might be estimated either by satellite or coastal tide gauges by connecting the tide gauge datum to the Earth-centred reference frame. In this article we present a novel approach to improve the coastal MDT, where the solution is based on both satellite altimetry and tide gauge data using new set of 302 tide gauges with ellipsoidal heights through the SONEL network. The approach was evaluated for the Northeast Atlantic coast where a dense network of GNSS-surveyed tide gauges is available. The typical misfit between tide gauge and satellite or oceanographic MDT was found to be around 9?cm. This misfit was found to be mainly due to small scale geoid errors. Similarly, we found, that a single tide gauge places only weak constraints on the coastal dynamic topography.     

Influence of the Kuroshio Fluctuations on Sea Level Variations along the South Coast of Japan

The correlation between the Kuroshio and coastal sea level south of Japan has been examined using the altimetry and tide gauge data during the period 1992–2000. The sea level varies uniformly in a region bounded by the coast and the mean Kuroshio axis, which stretches for several hundred kilometers along the coast. These variations are related with the Kuroshio velocity, as coastal sea level decreases (or increases) when the Kuroshio is faster (or slower). To the east of the Kii Peninsula, where sea level variations are different from these to the west, movement of the Kuroshio axis additionally affects coastal sea level variations.     

A Century of Sea Level Measurements at Newlyn,Southwest England

The Newlyn Tidal Observatory is the most important sea level station in the United Kingdom. It commenced operations in 1915 as part of the Second Geodetic Levelling of England and Wales, and the mean sea level determined from the tide gauge during the first six years (May 1915–April 1921) defined Ordnance Datum Newlyn (ODN) which became the national height datum for the whole of Great Britain. The 100 years of sea level data now available have contributed significantly to many studies in oceanography, geology, and climate change. This paper marks the centenary of this important station by reviewing the sea level (and, more recently, detailed land level) measurements and Newlyn's contributions to UK cartography, geodesy and sea-level science in general. Recommendations are made on how sea and land level measurements at Newlyn might be enhanced.     

Assessment of SARAL-ALTIKA Tidal Corrections in the Coastal Oceans Around India

The present tidal correction of sea level records of Satellite with ARgoes and ALtimeter (SARAL) is based on the finite element solution (FES) of global tide model FES2012 tidal solution. In this study, we examined the validity of the tidal corrections in the coastal oceans around India using tide gauge measurements and a regional tidal model. Our regional model is based on the barotropic version of the Princeton Ocean Model that is forced by the time-varying tidal levels at the open ocean end based on the global FES99 tidal solution. Tide charts prepared from the simulated tidal levels are very similar to the FES tidal solutions. Comparison with the tide gauge measurement shows close agreement with the regional tidal solutions. On the other hand, the agreement with the FES tide models differ significantly in the Gulf of Khambhat and the Gulf of Kutch on the northwest, and in the Hooghly estuary on the northeast continental shelf. However, the agreement is exceptional in other parts of the study domain. These tidal solutions are used in the SARAL-ALTIKA X-track data to assess the FES tidal correction and to draw some inferences associated with the coastal processes. It is revealed that these corrections are reasonably accurate for the coastal oceans around India except the aforementioned converging channels.     

Interannual Sea Level Variations in the South China Sea Over 1950–2009

Spatial patterns of interannual sea level variations in the South China Sea (SCS) are investigated by analyzing an EOF-based 2-dimensional past sea level reconstruction from 1950 to 2009 and satellite altimetry data from 1993 to 2009. Long-term tide gauge records from 14 selected stations in this region are also used to assess the quality of reconstructed sea levels and determine the rate of sea level along the coastal area. We found that the rising rate of sea levels derived from merged satellite altimetry data during 1993–2009 and past sea level reconstruction over 1950–2009 is about 3.9 ± 0.6 mm/yr and 1.7 ± 0.1 mm/yr, respectively. For the longer period, this rate is not significantly different from the global mean rate (of 1.8 ± 0.3 mm/yr). The interannual mean sea level of the SCS region appears highly correlated with Niño 4 indices (a proxy of El Niño-Southern Oscillation/ENSO), suggesting that the interannual sea level variations over the SCS region is driven by ENSO events. Interpolation of the reconstructed sea level data for 1950–2009 at sites where tide gauge records are of poor quality (either short or gapped) show that sea level along the Chinese coastal area is rising faster than the global mean rate of 1.8 mm/yr. At some sites, the rate is up to 2.5 mm/yr.     

崖城油气田附近海域溢油漂移轨迹数值模拟预测

通过采用不规则的三角网格和有限体积法的FVCOM模式,建立三维潮流数值模型。利用大海域计算得到的调和常数值作为开边界的输入值,模拟出崖城附近海域的潮流和潮位变化情况。在潮流、潮位验证正确的前提下,利用欧拉—拉格朗日追踪方法,建立了溢油轨迹预测模型,进行崖城油气田附近海域溢油中心轨迹的预测,同时预测了溢油漂移的平均速率和油膜抵达敏感区的时间,为油气田实施应急措施提供技术支持。     

Sea-Level Variations Measured by the New Altimetry Mission SARAL/AltiKa and its Validation Based on Spatial Patterns and Temporal Curves Using Jason-2, Tide Gauge Data and an Overview of the Annual Sea Level Budget

High-precision satellite altimeters help in measuring the variations in sea level since the early 1990s. After a number of such successful altimetry missions such as Topex/Poseidon, Jason-1, Jason-2, and Envisat, SARAL/AltiKa, a high resolution altimetry mission based on the Ka frequency band that can also cover high latitudinal zones, was launched in February 2013. Even though the data set available from this recent mission is not yet suitable for climate research owing to its short duration, in this study we perform a preliminary validation of SARAL/AltiKa sea-level data. The first part of the validation is the comparison of SARAL/AltiKa and Jason-2 sea-level data between March 2013 and August 2014 in terms of temporal mean spatial pattern. Comparisons in terms of global mean sea-level time series and latitudinal band-based mean time series are also performed. The second part of the validation is the comparison of the SARAL/AltiKa sea-level based time series with several tide gauge records covering the period of our study. Finally, an analysis of the annual sea-level budget with SARAL/AltiKa data, steric sea level, and ocean mass is performed. Results of these preliminary comparisons show good agreement with other sea-level data.     

Sea level changes at Ascension Island in the last half century

An exercise in ‘data archaeology’ at Ascension Island has provided an estimate of sea level change between 1955 and 2001.5 (the mid-point of a recent dataset spanning 1993–2009). That average trend of 0.93 mm y^?1 (SE 0.69) compares to a larger rate during 1993–2009 itself of 2.55 (SE 0.13) and 2.07 (SE 0.30) mm y^?1 from tide gauge and altimeter data respectively, suggesting a recent acceleration in sea level rise. An ocean model and steric height datasets have been used for comparison to the measurements, with the conclusion that the acceleration was probably at least partly due to a steric height increase. This exercise is based on only one month of historical tide gauge data and is admittedly on the useful limit for long-term sea level studies. In addition, it is unfortunate that the tide gauge benchmark installed in 1955 has disappeared, even if one can estimate its height relative to modern marks. Nevertheless, the study does provide information of interest to climate studies, enables limits to be inferred on the real changes, and provides background information for other coastal studies. Most importantly, it is intended as a demonstration of the value of similar exercises where short historical records exist.     

Fluctuation of the sea surface dynamic topography southeast of Japan as estimated from Seasat altimetry data

The sea surface dynamic topography (the sea surface height relative to the geoid; hereafter abbreviated SSDT) can be divided into the temporal mean SSDT and the fluctuation SSDT around the mean. We use the optimal interpolation method to reduce the satellite radial orbit error and estimate the fluctuation SSDT southeast of Japan from Seasat altimetry data during the 17-day near-repeat mission. The fluctuation SSDT is further combined with the mean geopotential anomalies estimated from hydrographic data during the Seasat mission in order to give the approximated total SSDT, called here the composite SSDT (the approximated mean plus fluctuation SSDT's). The fluctuation SSDT is in accord with the low-frequency sea-level fluctuation recorded at tide gauge stations in the Japanese islands. The composite SSDT describes thoroughly variations of the location of the Kuroshio axis south of Japan determined on the basis of the GEK (Geomagnetic Electro-Kinematograph) surface velocities and the horizontal temperature distribution. The composite SSDT also agrees with oceanic variations east of Japan found in the temperature distribution at the depth of 200 m. These results confirm that the SSDT derived from altimetry data can provide fairly precise synoptic views of low-frequency oceanic phenomena.     

Ocean Tide Models Developed by Assimilating TOPEX/POSEIDON Altimeter Data into Hydrodynamical Model: A Global Model and a Regional Model around Japan

A global ocean tide model (NAO.99b model) representing major 16 constituents with a spatial resolution of 0.5° has been estimated by assimilating about 5 years of TOPEX/POSEIDON altimeter data into barotropic hydrodynamical model. The new solution is characterized by reduced errors in shallow waters compared to the other two models recently developed; CSR4.0 model (improved version of Eanes and Bettadpur, 1994) and GOT99.2b model (Ray, 1999), which are demonstrated in comparison with tide gauge data and collinear residual reduction test. This property mainly benefits from fine-scale along-track tidal analysis of TOPEX/POSEIDON data. A high-resolution (1/12°) regional ocean tide model around Japan (NAO.99Jb model) by assimilating both TOPEX/POSEIDON data and 219 coastal tide gauge data is also developed. A comparison with 80 independent coastal tide gauge data shows the better performance of NAO.99Jb model in the coastal region compared with the other global models. Tidal dissipation around Japan has been investigated for M₂ and K₁ constituents by using NAO.99Jb model. The result suggests that the tidal energy is mainly dissipated by bottom friction in localized area in shallow seas; the M₂ ocean tidal energy is mainly dissipated in the Yellow Sea and the East China Sea at the mean rate of 155 GW, while the K₁ energy is mainly dissipated in the Sea of Okhotsk at the mean rate of 89 GW. TOPEX/POSEIDON data, however, detects broadly distributed surface manifestation of M₂ internal tide, which observationally suggests that the tidal energy is also dissipated by the energy conversion into baroclinic tide.     

基于海气耦合模型的渤海两种类型风暴潮数值模拟与对比分析

渤海一年四季都易受到由温带风暴和热带气旋所致风暴潮的影响。为了缓解风暴潮灾害对海岸地区人员生命财产的影响,十分有必要了解大型风暴潮的发生过程和机制。目前大部分研究主要局限于单一的温带风暴潮或台风风暴潮。本文利用所构建的海气耦合数值模型研究了发生于渤海的两种类型的风暴潮,对发生在渤海的2次典型强风暴潮过程进行了模拟。由WRF模型模拟得到的风场强度和最低海平面气压与实测数据吻合较好,由ROMS模型模拟得到的风暴潮期间水位变化过程与潮位站观测结果也吻合较好。对两种类型风暴潮期间的风场结钩、海面风应力、海洋表面平均流场以及水位分布进行了分析对比,并将耦合模型结果与非耦合模型结果进行了对比。研究表明,渤海两种类型风暴潮期间的风场结钩、海面风应力、海洋表面平均流场以及水位分布等均存在巨大差异。渤海风暴潮的强度主要由海洋表面的驱动力所决定,但同时也受海岸地形地貌的影响。     

An Improved Calibration of Satellite Altimetric Heights Using Tide Gauge Sea Levels with Adjustment for Land Motion

Several major improvements to an existing method for calibrating satellite altimeters using tide gauge data are described. The calibration is in the sense of monitoring and correcting temporal drift in the altimetric time series, which is essential in efforts to use the altimetric data for especially demanding applications. Examples include the determination of the rate of change of global mean sea level and the study of the relatively subtle, but climatically important, decadal variations in basin scale sea levels. The improvements are to the method described by Mitchum (1998a), and the modifications are of two basic types. First, since the method depends on the cancellation of true ocean signals by differencing the altimetric data from the tide gauge sea level time series, improvements are made that produce a more complete removal of the ocean signals that comprise the noise for the altimetric drift estimation problem. Second, a major error source in the tide gauge data, namely land motion, is explicitly addressed and corrections are developed that incorporate space-based geodetic data (continuous GPS and DORIS measurements). The long-term solution, having such geodetic measurements available at all the tide gauges, is not yet a reality, so an interim solution is developed. The improved method is applied to the TOPEX altimetric data. The Side A data (August 1992?February 1999) are found to have a linear drift component of 0.55 + / 0.39 mm/yr, but there is also a significant quadratic component to the drift that is presently unexplained. The TOPEX Side B altimeter is estimated to be biased by 7.0 + / 0.7 mm relative to the Side A altimeter based on an analysis of the first 350 days of Side B data.     

Book Reviews

The UK Hydrographic Office (UKHO)-sponsored Vertical Offshore Reference Frames (VORF) project aims to develop tidal level transformation models that are referenced to the GRS80 ellipsoid and thus compatible with GNSS positioning; in particular, heighting. Benefits include increasing the efficiency of hydrographic surveying, providing a stable consistent reference frame and enabling integration with land data in the coastal zone. Seven contemporary global ocean tide models are used to derive Lowest Astronomical Tide (LAT) surfaces which are each assessed by comparison with LAT values from the 7,389-strong UKHO tide gauge database, with the results correlated with distance from land. The proportion of truly offshore and pelagic gauges is relatively limited; however, the transition zone whereby the global ocean tide models commence to deteriorate in accuracy is evident at approximately 30km from the coast. The DTU10 model was selected as the strongest candidate overall. Subsequently, a thin plate spline method is used with the tide gauge dataset to enhance the DTU10 LAT surface in the coastal zone, creating a high resolution global LAT surface with respect to mean sea level. It is seen by cross-validation that the method may be used to predict LAT in near-shore locations with a standard error of 0.23 m.