Ocean tides for satellite geodesy

Abstract

Spherical harmonic tidal solutions have been obtained at the frequencies of the 32 largest luni‐solar tides using prior theory of the author. That theory was developed for turbulent, nonglobal, self‐gravitating, and loading oceans possessing realistic bathymetry and linearized bottom friction; the oceans satisfy no‐flow boundary conditions at coastlines. In this theory the eddy viscosity and bottom drag coefficients are treated as spatially uniform. Comparison of the predicted degree‐2 components of the Mf, PI, and M2 tides with those from numerical and satellite‐based tide models allows the ocean friction parameters to be estimated at long and short periods. Using the 32 tide solutions, the frequency dependence of tidal admittance is investigated, and the validity of sideband tide models used in satellite orbit analysis is examined. The implications of admittance variability for oceanic resonances are also explored. By extending the theory to include a second constraint derived from tide observations or data‐constrained tide models, it is possible to assess those models from a fluid dynamic perspective. One general conclusion from such exercises is that the large higher‐degree admittances of current short‐period tide models are dynamically incompatible with their degree‐2 admittances. Eventually it may prove possible to produce dynamically sound, observationally consistent tide models by combining the author's tide theory with satellite orbit determination.     

Non-stationary extreme value models to account for trends and shifts in the extreme wave climate due to climate change

The extreme values of wave climate data are of great interest in a number of different ocean engineering applications, including the design and operation of ships and offshore structures, marine energy generation, aquaculture and coastal installations. Typically, the return values of certain met-ocean parameters such as significant wave height are of particular importance. There exist many methods for estimating such return values, including the initial distribution approach, the block maxima approach and the peaks-over threshold approach. In a climate change perspective, projections of such return values to a future climate are of great importance for risk management and adaptation purposes. However, many approaches to extreme value modelling assume stationary conditions and it is not straightforward how to include non-stationarity of the extremes due to for example climate change. In this paper, various non-stationary GEV-models for significant wave height are developed that account for trends and shifts in the extreme wave climate due to climate change. These models are fitted to block maxima in a particular set of wave data obtained for a historical control period and two future projections for a future period corresponding to different emission scenarios. These models are used to investigate whether there are trends in the data within each period that influence the extreme value analysis and need to be taken into account. Moreover, it will be investigated whether there are significant inter-period shifts or trends in the extreme wave climate from the historical period to the future periods. The results from this study suggest that the intra-period trends are not statistically significant and that it might be reasonable to ignore these in extreme value analyses within each period. However, when it comes to comparing the different data sets, i.e. the historical period and the future projections, statistical significant inter-period changes are detected. Hence, the accumulated effect of a climatic trend may not be negligible over longer time periods. Interestingly enough, such statistically significant shifts are not detected if stationary extreme value models are fitted to each period separately. Therefore, the non-stationary extreme value models with inter-period shifts in the parameters are proposed as an alternative for extreme value modelling in a climate change perspective, in situations where historical data and future projections are available.     

A reduced-dynamics variational approach for the assimilation of altimeter data into eddy-resolving ocean models

A new method of assimilating sea surface height (SSH) data into ocean models is introduced and tested. Many features observable by satellite altimetry are approximated by the first baroclinic mode over much of the ocean, especially in the lower (but non-equatorial) and mid latitude regions. Based on this dynamical trait, a reduced-dynamics adjoint technique is developed and implemented with a three-dimensional model using vertical normal mode decomposition. To reduce the complexity of the variational data assimilation problem, the adjoint equations are based on a one-active-layer reduced-gravity model, which approximates the first baroclinic mode, as opposed to the full three-dimensional model equations. The reduced dimensionality of the adjoint model leads to lower computational cost than a traditional variational data assimilation algorithm. The technique is applicable to regions of the ocean where the SSH variability is dominated by the first baroclinic mode. The adjustment of the first baroclinic mode model fields dynamically transfers the SSH information to the deep ocean layers. The technique is developed in a modular fashion that can be readily implemented with many three-dimensional ocean models. For this study, the method is tested with the Navy Coastal Ocean Model (NCOM) configured to simulate the Gulf of Mexico.     

Turbulent dispersion in the ocean

The mathematical framework for turbulent transport in the ocean is reasonably well established. It may be applied to large-scale fields of scalars in the ocean and to the instantaneous or continuous discharge from a point. The theory and its physical basis can also provide an interpretation of passive scalar spectra. Spatial variations in the rate of turbulent transfer can be related to the movement of the center of mass of a scalar and to a formulation in terms of entrainment. The relative dispersion of a scalar with respect to its center of mass and the streakiness of the concentration field within the relative dispersion domain need to be considered. In many of these problems it is valuable to think in terms of simple models for individual streaks, as well as overall statistical properties.     

Vertical transport in the ocean due to sub-mesoscale structures: Impacts in the Kerguelen region

The summertime phytoplankton bloom near the Kerguelen Plateau is in marked contrast to the low-chlorophyll conditions typical of the Southern Ocean and is thought to arise from natural iron fertilisation. The mechanisms of iron supply to the euphotic zone in this region are poorly understood, and numerical studies of iron transport have until now omitted fine-scale (sub-mesoscale) dynamics which have been shown to significantly increase vertical transport in other parts of the ocean.We present the first sub-mesoscale-resolving study of the flow and vertical transport in this region. The modelled transport and flow structure agree well with observations. We find that an increase in horizontal resolution from mesoscale-resolving (1/20°) to 1/80° resolves sub-mesoscale filamentary frontal structures in which vertical velocities are dramatically higher and are consistent with available observations. Lagrangian tracking shows that water is advected to the surface from much greater depth in the sub-mesoscale-resolving experiment, and that vertical exchange is far more rapid and frequent. This study of sub-mesoscale vertical velocities sets the foundation for subsequent investigation of iron transport in this environment.     

A singular evolutive extended Kalman filter to assimilate ocean color data in a coupled physical–biochemical model of the North Atlantic ocean

Within the European DIADEM project, a data assimilation system for coupled ocean circulation and marine ecosystem models has been implemented for the North Atlantic and the Nordic Seas. One objective of this project is to demonstrate the relevance of sophisticated methods to assimilate satellite data such as altimetry, surface temperature and ocean color, into realistic ocean models. In this paper, the singular evolutive extended Kalman (SEEK) filter, which is an advanced assimilation scheme where three-dimensional, multivariate error statistics are taken into account, is used to assimilate ocean color data into the biological component of the coupled system. The marine ecosystem model, derived from the FDM model [J. Mar. Res. 48 (1990) 591], includes 11 nitrogen and carbon compartments and describes the synthesis of organic matter in the euphotic zone, its consumption by animals of upper trophic levels, and the recycling of detritic material in the deep ocean. The circulation model coupled to the ecosystem is the Miami isopycnic coordinate ocean model (MICOM), which covers the Atlantic and the Arctic Oceans with an enhanced resolution in the North Atlantic basin. The model is forced with realistic ECMWF ocean/atmosphere fluxes, which permits to resolve the seasonal variability of the circulation and mixed layer properties. In the twin assimimation experiments reported here, the predictions of the coupled model are corrected every 10 days using pseudo-measurements of surface phytoplankton as a substitute to chlorophyll concentrations measured from space. The diagnostics of these experiments indicate that the assimilation is feasible with a reduced-order Kalman filter of small rank (of order 10) as long as a sufficiently good identification of the error structure is available. In addition, the control of non-observed quantities such as zooplankton and nitrate concentrations is made possible, owing to the multivariate nature of the analysis scheme. However, a too severe truncation of the error sub-space downgrades the propagation of surface information below the mixed layer. The reduction of the actual state vector to the surface layers is therefore investigated to improve the estimation process in the perspective of sea-viewing wide field-of-view sensor (SeaWiFS) data assimilation experiments.     

Impact of Aquarius Sea-Surface Salinity Assimilation in Improving the Ocean Analysis Over Indian Ocean

The distribution of ocean salinity controls the density field and thereby plays a major role in influencing the ocean dynamics. It has been a challenging task to understand the variability of salinity structure in the regions of large fresh water discharge and high precipitation such as Bay of Bengal (BoB). Recent advancement in satellite technology has made possible the measurement of sea surface salinity (SSS). Aquarius is the satellite which measured the global SSS for the period 2011 to 2015. In the present study, we assimilated Aquarius SSS in the Global Ocean Data Assimilation System based on 3DVAR technique. The assimilation of Aquarius SSS resulted in reduced biases in salinity not only at the surface, but also in the vertical distribution of salinity and better captured the temporal variations of salinity structure in sensitive regions, such as the Bay of Bengal. In addition, the assimilation of SSS showed marginal improvement in ocean thermal structure over data sparse regions of Indian Ocean. It is also shown that the assimilation of Aquarius SSS has improved the stratification in the upper Ocean which is the key factor in the observed improvement in ocean analysis.     

Comparison of the nonlocal transport characteristics of a series of one-dimensional oceanic boundary layer models

Seven one-dimensional oceanic boundary layer models are investigated to assess the possible nonlocal transport characteristics of mass and heat in the upper ocean. The dynamical models have been chosen from the diffusion and bulk types currently in use plus two modifications of the transilient type that have been used extensively for atmospheric work by Stull. The models are forced using wind speed and insolation conditions recorded during a 9-day oceanographic cruise near Bermuda in March 1993 during the decline of the spring bloom. The attenuation of sunlight in the upper ocean is calculated using a full spectral model for downwelling irradiance. The vertical heat transport characteristics are reported and compared. A series of spectral diagnostic tests (Green's function analysis, process spectra, and overall mixing lengths) reveal significant differences in the vertical transport characteristics of the models that are not observed in commonly used diagnostics such as sea surface temperature or mixed layer depth. Age spectra (or modal time since last surface contact) are calculated from Green's functions and reveal how a water mass can be cut off from the surface within a short time period. The large differences in the vertical mixing characteristics of the upper ocean boundary layer models have potential implications on the vertical distributions of short-lived chemical tracers and phytoplankton.     

卫星海洋水色遥感的辐射模式研究

海洋水体向上的辐射由水中多种成分的浓度所决定，同时辐射的强度十分微弱并随着时间和空间在变化。因此，海洋水色遥感需要从水体和大气辐射传递机理出发，研究到达卫星遥感器的辐射模式，以达到模拟仿真卫星水色遥感图像和卫星图像质量预测的目的。本研究通过叶绿素、悬浮泥沙和黄色物质等海洋水色的主要因子，分别从可见光波段和荧光波段对辐射的贡献出发，发展了海洋水体离水辐射率模式。在卫星水色遥感机理的基础上，同时运用Iqbal，Gordon和Sturm等的大气程辐射模式和太阳耀光模式，使卫星海洋水色遥感的辐射模式系统化，模拟仿真了我国1990年9月3日发射的FYIB卫星两个海洋水色通道的全轨道辐射图像，同时利用从德国柏林大学气象中心接收到的1990年9月23－25日卫星图像进行验证。结果认为，本文所发展的模式可应用于卫星图像质量的预测和水色信息的提取。     

Evaluation of GODAS Using RAMA Mooring Observations from the Indian Ocean

We present a comparison of the Global Ocean Data Assimilation System (GODAS) five-day ocean analyses against in situ daily data from Research Moored Array for African-Asian-Australian Monsoon Analysis and Prediction (RAMA) moorings at locations 90°E, 12°N; 90°E, 8°N; 90°E, 0°N and 90°E, 1.5°S in the equatorial Indian Ocean and the Bay of Bengal during 2002–2008. We find that the GODAS temperature analysis does not adequately capture a prominent signal of Indian Ocean dipole mode of 2006 seen in the mooring data, particularly at 90°E 0°N and 90°E 1.5°S in the eastern India Ocean. The analysis, using simple statistics such as bias and root-mean-square deviation, indicates that standard GODAS temperature has definite biases and significant differences with observations on both subseasonal and seasonal scales. Subsurface salinity has serious deficiencies as well, but this may not be surprising considering the poorly constrained fresh water forcing, and possible model deficiencies in subsurface vertical mixing. GODAS reanalysis needs improvement to make it more useful for study of climate variability and for creating ocean initial conditions for prediction.     

论厦门海洋经济的持续发展

本文回顾了厦门海洋经济发展的历史,评述厦门海洋经济发展的现状和趋势,认为厦门已经形成了沿海城市经济、已经具有相当完备的海洋经济体系,但传统产业依然占据较大的比重,新兴产业还处于起步阶段,未来产业基本空白,城市经济的发展与海洋资源的利用结合不够紧密。海洋资源的开发利用尚为粗放式,科技进步的贡献率在海洋增长中不够大。海洋经济未走上开发和保护相结合的可持续发展道路。在发展重点海洋产业中,必须树立“科技兴     

Sea ice–ocean coupling using a rescaled vertical coordinate z

Realistic representation of sea ice in ocean models involves the use of a non-linear free-surface, a real freshwater flux and observance of requisite conservation laws. We show here that these properties can be achieved in practice through use of a rescaled vertical coordinate “z^*” in z-coordinate models that allows one to follow undulations in the free-surface under sea ice loading. In particular, the adoption of “z^*” avoids the difficult issue of vanishing levels under thick ice.Details of the implementation within MITgcm are provided. A high resolution global ocean sea ice simulation illustrates the robustness of the z^* formulation and reveals a source of oceanic variability associated with sea ice dynamics and ice-loading effects. The use of the z^* coordinate allows one to achieve perfect conservation of fresh water, heat and salt, as shown in extended integration of coupled ocean sea ice atmospheric model.     

海浪破碎对海洋上混合层中湍能量收支的影响

海浪破碎产生一向下输入的湍动能通量,在近海表处形成一湍流生成明显增加的次层,加强了海洋上混合层中的湍流垂向混合。为了研究海浪破碎对混合层中湍能量收支的影响,文中分析了海浪破碎对海洋上混合层中湍流生成的影响机制,采用垂向一维湍封闭混合模式,通过改变湍动能方程的上边界条件,引入了海浪破碎产生的湍动能通量,并分别对不同风速下海浪破碎的影响进行了数值研究,分析了混合层中湍能量收支的变化。当考虑海浪破碎影响时,近海表次层中的垂直扩散项和耗散项都有显著的增加,该次层中被耗散的湍动能占整个混合层中耗散的总的湍能量的92.0%,比无海浪破碎影响的结果增加了近1倍;由于平均流场切变减小,混合层中的湍流剪切生成减小了3.5%,形成一种存在于湍动能的耗散和垂直扩散之间的局部平衡关系。在该次层以下,局部平衡关系与壁层定律的结论一致,即湍动能的剪切生成与耗散相平衡。研究结果表明,海浪破碎在海表产生的湍动能通量影响了海洋上混合层中的各项湍能量收支间的局部平衡关系。     

印度洋大西洋上升流的长期变化及其与爱尔尼诺南方涛动的关系

用Ｘｉｅ和Ｈｓｉｅｈ从三十九年（１９５０－１９８８）Ｃｏａｄｓ月平均风应力资料算出的全球海洋Ｅｋｍａｎ层底垂直速度资料，分别对印度洋和大西洋的垂直速度场作ＥＯＦ分析。在３８４个模中集中讨论了前四个模，印度洋的四个模被发现均与Ｅｌｎｉｎｏ事件有关，模一与模二表现为Ｅｌｎｉｎｏ前后的剧降与剧升，（只是模二的位相较模一滞后半年）而模三与模四的位相与模一、二反相、时间振幅的曲线在Ｅｌｎｉｎｏ前达到正峰值，     

大洋科学钻探特点与发展趋势——基于国际大洋发现计划科学框架的对比分析

大洋科学钻探计划是地球科学领域迄今为止历时最长、成效最大的国际科学合作计划。自1968年以来,先后经历了深海钻探计划、国际大洋钻探计划、综合大洋钻探计划和国际大洋发现计划4个阶段。进入国际大洋发现计划,IODP工作组两次发布科学框架,文章对两个科学框架进行了对比分析。大洋科学钻探的特点和发展趋势是坚持解决重大科学问题,由地球拓展到太阳系,重视技术和大数据,以及扩大社会影响和宣传。这给我国发起国际大科学计划提供了启示,一要注重科技创新,运用大数据等先进技术;二要加强科学管理,促进开放和共享;三要重视提升影响力,扩大宣传和影响。     

2011年3月日本福岛核电站核泄漏在海洋中的传输

使用全球版本的迈阿密等密度海洋环流模式对2011年3月日本福岛核电站泄漏在海洋中的传输以及扩散进行了数值模拟。数值模式中核废料(示踪物)排放情景采取等通量连续排放,排放时间从3月25日开始,分别持续20 d以及1 a,两种情形分别积分20 a。为了减少海洋环流年际变化带来的数值模拟的的不确定性,20 a的模式积分分别用2010年、1991-2011年、1971-1991年以及1951-1971年4个不同时段的NCEP/NCAR逐日再分析资料作为大气强迫场,因此每种排放情形包含4个数值试验。模拟结果的分析表明,不同核废料排放情景及其在不同时段大气资料对海洋模式的驱动下,模拟的示踪物总体的传输扩散路径(包括表层以及次表层)、传输速率以及垂直扩展的范围没有显著的差异。集合平均数值模拟的结果显示:在两种排放情景下,日本福岛核泄漏在海洋的传输路径受北太平洋副热带涡旋洋流系统主导,其传输路径首先主要向东,到达东太平洋后,再向南向西扩散至西太平洋,可能在10~15 a左右影响到我国东部沿海海域,且海洋次表层的传输信号比表层信号早5 a左右。通过进一步分析模式积分过程中最大示踪物浓度随时间变化发现,在积分第20 a(2031年3月),海洋表层和次表层浓度的最高值分别只有模式积分第一年浓度的0.1%和1%。在积分的20 a里,排放的核废料主要滞留在北太平洋海域(超过86%±1.5%的核废料在积分结束时,滞留在北太平洋),而在积分的前10 a(2021年之前),几乎所有的核废料滞留在北太平洋;在核废料的垂直分布上,主要集中在海洋表层至600 m的深度,在积分的20 a时间里,没有核废料信号扩散至1 000 m以下的深度。数值模拟的结果也表明核废料浓度减弱的强度以及演变的时间特征主要受洋流系统的影响,与排放源的排放时间长短关系不大。值得指出的是,更加准确地评估一个真实的核泄漏事故对海洋环境所造成的可能影响,还需要考虑大气中的放射性物质的沉降以及海洋生态对核物质的响应。     

2011年3月日本福岛核电站核泄漏在海洋中的传输

使用全球版本的迈阿密等密度海洋环流模式对2011年3月日本福岛核电站泄漏在海洋中的传输以及扩散进行了数值模拟。数值模式中核废料(示踪物)排放情景采取等通量连续排放,排放时间从3月25日开始,分别持续20 d以及1 a,两种情形分别积分20 a。为了减少海洋环流年际变化带来的数值模拟的的不确定性,20 a的模式积分分别用2010年、1991-2011年、1971-1991年以及1951-1971年4个不同时段的NCEP/NCAR逐日再分析资料作为大气强迫场,因此每种排放情形包含4个数值试验。模拟结果的分析表明,不同核废料排放情景及其在不同时段大气资料对海洋模式的驱动下,模拟的示踪物总体的传输扩散路径(包括表层以及次表层)、传输速率以及垂直扩展的范围没有显著的差异。集合平均数值模拟的结果显示:在两种排放情景下,日本福岛核泄漏在海洋的传输路径受北太平洋副热带涡旋洋流系统主导,其传输路径首先主要向东,到达东太平洋后,再向南向西扩散至西太平洋,可能在10~15 a左右影响到我国东部沿海海域,且海洋次表层的传输信号比表层信号早5 a左右。通过进一步分析模式积分过程中最大示踪物浓度随时间变化发现,在积分第20 a(2031年3月),海洋表层和次表层浓度的最高值分别只有模式积分第一年浓度的0.1%和1%。在积分的20 a里,排放的核废料主要滞留在北太平洋海域(超过86%±1.5%的核废料在积分结束时,滞留在北太平洋),而在积分的前10 a(2021年之前),几乎所有的核废料滞留在北太平洋;在核废料的垂直分布上,主要集中在海洋表层至600 m的深度,在积分的20 a时间里,没有核废料信号扩散至1 000 m以下的深度。数值模拟的结果也表明核废料浓度减弱的强度以及演变的时间特征主要受洋流系统的影响,与排放源的排放时间长短关系不大。值得指出的是,更加准确地评估一个真实的核泄漏事故对海洋环境所造成的可能影响,还需要考虑大气中的放射性物质的沉降以及海洋生态对核物质的响应。     

Determination of Ocean Tide Loading Displacement by GPS PPP with Priori Information Constraint of NAO99b Global Ocean Tide Model

The published global ocean tide models show good agreement in deep oceans and exhibit differences in complex coastal areas, along with subsequent Ocean Tide Loading Displacement (OTLD) modeling differences. Meanwhile, OTLD parameters (amplitudes and phase lags) derived by Global Positioning System (GPS) Precise Point Positioning (PPP) approach need long time to converge to a stable state and show poor precision of S2, K1, and K2 constituents. Based on the fact that no constraint is imposed in the current kinematic solution, a new method is put forward, in which global ocean tide model predictions are taken as the priori information constraints to speed up the convergence rate and improve the accuracy of the GPS-derived OTLD parameters. First, the data of tide gauge from 01 January 2014 to 31 December 2016 are used to generate the harmonic parameters to evaluate the accuracy of six global ocean tide models and a regional ocean tide model (osu.chinesea.2010). Osu.chinesea.2010 model shows good agreement with the tide gauge results, while NAO99b model presents relatively large difference. The predictions from osu.chinesea.2010 and NAO99b model are employed as reference and the prior information, respectively. Second, continuous observations of 12 GPS sites during 2006–2013 in Hong Kong are collected to generate three dimensional OTLD amplitudes and phase lags of eight constituents using PPP with prior information constraints approach and harmonic analysis. Third, comparing the convergence time of eight constituents from PPP without and with priori information constraints approaches, the results show that the new method can significantly improve the convergence rate of OTLD amplitude estimates which obtain a certain level of stability seven years earlier than that derived by the PPP without priori information constraints. Precision of OTLD parameters derived by the new method is about 1 mm. By comparing with the precision of single PPP approach, the accuracy of eight constituents has been improved, especially for S2, K1, and K2 constituents. Finally, through comparing the different correction effects of OTLD estimates on the coordinates and their time series of the ground GPS stations, the results show that OTLD estimates derived by the new approach have similar influence as the osu.chinasea.2010 ocean tide model. The new method provides an effective means to improve the convergence and precision of the GPS-derived OTLD parameters, and achieve a similar correction as the high precision ocean tide model.     

Simulated Heat Content Variability in the Upper Layers of the Tropical Indian Ocean

Ocean General Circulation Model (OGCM) simulations from 1970–2007 are used to study the upper ocean heat content variability in the Tropical Indian Ocean (TIO). Model computed heat contents up to 50 m (denoted by HC50 m hereafter) representing upper ocean heat content and 300 m (HC300 m) representing heat content up to thermocline depth are first compared with heat contents computed from observations of two buoys in the TIO. It is found that there is good agreement between the model and observations. Fourier analysis of heat content is carried out in different regions of TIO. The amplitudes of semi-annual variability for HC50 m and HC300 m are observed to be greater than those for the annual variability in the Bay of Bengal, while in the Arabian Sea there is a mixed result. Heat content tendency is known to be governed by net surface heat flux and horizontal as well as vertical heat transports. For understanding the relative importance of these processes, a detailed analysis of these terms in the tendency equation is carried out. Rossby wave is observed in the annual mode of heat transport while equatorial jet and Kelvin waves are observed in the semi-annual mode of heart transport. Finally, the correlation between heat content and sea surface temperature (SST) and sea level anomaly (SLA), taken one at a time, is computed. It is found that the correlation improves significantly when both these quantities are together taken into account.     

Accurate polynomial expressions for the density and specific volume of seawater using the TEOS-10 standard

A new set of approximations to the standard TEOS-10 equation of state are presented. These follow a polynomial form, making it computationally efficient for use in numerical ocean models. Two versions are provided, the first being a fit of density for Boussinesq ocean models, and the second fitting specific volume which is more suitable for compressible models. Both versions are given as the sum of a vertical reference profile (6th-order polynomial) and an anomaly (52-term polynomial, cubic in pressure), with relative errors of ∼0.1% on the thermal expansion coefficients. A 75-term polynomial expression is also presented for computing specific volume, with a better accuracy than the existing TEOS-10 48-term rational approximation, especially regarding the sound speed, and it is suggested that this expression represents a valuable approximation of the TEOS-10 equation of state for hydrographic data analysis. In the last section, practical aspects about the implementation of TEOS-10 in ocean models are discussed.