利用海表盐度等遥感观测产品重构三维盐度场的性能评估

Several remotely sensed sea surface salinity(SSS) retrievals with various resolutions from the soil moisture and ocean salinity(SMOS) and Aquarius/SAC-D missions are applied as inputs for retrieving salinity profiles(S) using multilinear regressions. The performance is evaluated using a total root mean square(RMS) error, different error sources, and the feature resolutions of the retrieved S fields. In the mixed layer of the salinity, the SSS-S regression coefficients are uniformly large. The SSS inputs yield smaller RMS errors in the retrieved S with respect to Argo profiles as their spatial or temporal resolution decreases. The projected SSS errors are dominant, and the retrieved S values are more accurate than those of climatology in the tropics except for the tropical Atlantic, where the regression errors are abnormally large. Below that level, because of the influence of a sea level anomaly, the areas of high-accuracy S values shift to higher latitudes except in the high-latitude southern oceans, where the projected SSS errors are abnormally large. A spectral analysis suggests that the CATDS-0.25° results are much noisier and that the BEC-L4-0.25° results are much smoother than those of the other retrievals. Aquarius-CAP-1° generates the smallest RMS errors, and Aquarius-V2-1° performs well in depicting large-scale phenomena. BEC-L3-0.25°,which has small RMS errors and remarkable mesoscale energy, is the best fit for portraying mesoscale features in the SSS and retrieved S fields. The current priority for retrieving S is to improve the reliability of satellite SSS especially at middle and high latitudes, by developing advanced algorithms, combining both sensors, or weighing between accuracy and resolutions.     

降雨导致的海面粗糙度对Aquarius卫星盐度反演的影响研究

Rainfall has two significant effects on the sea surface, including salinity decreasing and surface becoming rougher,which have further influence on L-band sea surface emissivity. Investigations using the Aquarius and TRMM 3B42 matchup dataset indicate that the retrieved sea surface salinity(SSS) is underestimated by the present Aquarius algorithm compared to numerical model outputs, especially in cases of a high rain rate. For example, the bias between satellite-observed SSS and numerical model SSS is approximately 2 when the rain rate is 25 mm/h. The bias can be eliminated by accounting for rain-induced roughness, which is usually modeled by rain-generated ring-wave spectrum. The rain spectrum will be input into the Small Slope Approximation(SSA) model for the simulation of sea surface emissivity influenced by rain. The comparison with theoretical model indicated that the empirical model of rain spectrumis more suitable to be used in the simulation. Further, the coefficients of the rain spectrum are modified by fitting the simulations with the observations of the 2–year Aquarius and TRMM matchup dataset. The calculations confirm that the sea surface emissivity increases with the wind speed and rain rate. The increase induced by the rain rate is rapid in the case of low rain rate and low wind speed. Finally, a modified model of sea surface emissivity including the rain spectrum is proposed and validated by using the matchup dataset in May 2014. Compared with observations, the bias of the rain-induced sea surface emissivity simulated by the modified modelis approximately 1e–4, and the RMSE is slightly larger than 1e–3. With using more matchup data, thebias between model retrieved sea surface salinities and observationsmay be further corrected,and the RMSE may be reduced to less than 1 in the cases of low rain rate and low wind speed.     

一种利用SMOS卫星观测的海表面盐度反演海水盐度廓线的方法

This paper proposes a new method to retrieve salinity profiles from the sea surface salinity(SSS) observed by the Soil Moisture and Ocean Salinity(SMOS) satellite. The main vertical patterns of the salinity profiles are firstly extracted from the salinity profiles measured by Argo using the empirical orthogonal function. To determine the time coefficients for each vertical pattern, two statistical models are developed. In the linear model, a transfer function is proposed to relate the SSS observed by SMOS(SMOS_SSS) with that measured by Argo, and then a linear relationship between the SMOS_SSS and the time coefficient is established. In the nonlinear model, the neural network is utilized to estimate the time coefficients from SMOS_SSS, months and positions of the salinity profiles. The two models are validated by comparing the salinity profiles retrieved from SMOS with those measured by Argo and the climatological salinities. The root-mean-square error(RMSE) of the linear and nonlinear model are 0.08–0.16 and 0.08–0.14 for the upper 400 m, which are 0.01–0.07 and 0.01–0.09 smaller than the RMSE of climatology. The error sources of the method are also discussed.     

Validation and application of soil moisture active passive sea surface salinity observation over the Changjiang River Estuary

Using sea surface salinity(SSS)observation from the soil moisture active passive(SMAP)mission,we analyzed the spatial distribution and seasonal variation of SSS around Changjiang River(Yangtze River)Estuary for the period of September 2015 to August 2018.First,we found that the SSS from SMAP is more accurate than soil moisture and ocean salinity(SMOS)mission observation when comparing with the in situ observations.Then,the SSS signature of the Changjiang River freshwater was analyzed using SMAP data and the river discharge data from the Datong hydrological station.The results show that the SSS around the Changjiang River Estuary is significantly lower than that of the open ocean,and shows significant seasonal variation.The minimum value of SSS appears in July and maximum SSS in December.The root mean square difference of daily SSS between SMAP observation and in situ observation is around 3 in both summer and winter,which is much lower than the annual range of SSS variation.In summer,the diffusion direction of the Changjiang River freshwater depicted by SSS from SMAP is consistent with the path of freshwater from in situ observation,suggesting that SMAP observation may be used in coastal seas in monitoring the diffusion and advection of freshwater discharge.     

Variation of Indo-Pacific upper ocean heat content during 2001–2012 revealed by Argo

Understanding of the temporal variation of oceanic heat content(OHC) is of fundamental importance to the prediction of climate change and associated global meteorological phenomena. However, OHC characteristics in the Pacific and Indian oceans are not well understood. Based on in situ ocean temperature and salinity profiles mainly from the Argo program, we estimated the upper layer(0–750 m) OHC in the Indo-Pacific Ocean(40°S–40°N, 30°E–80°W). Spatial and temporal variability of OHC and its likely physical mechanisms are also analyzed. Climatic distributions of upper-layer OHC in the Indian and Pacific oceans have a similar saddle pattern in the subtropics, and the highest OHC value was in the northern Arabian Sea. However, OHC variabilities in the two oceans were different. OHC in the Pacific has an east-west see-saw pattern, which does not appear in the Indian Ocean. In the Indian Ocean, the largest change was around 10°S. The most interesting phenomenon is that, there was a long-term shift of OHC in the Indo-Pacific Ocean during 2001–2012. Such variation coincided with modulation of subsurface temperature/salinity. During 2001–2007, there was subsurface cooling(freshening)nearly the entire upper 400 m layer in the western Pacific and warming(salting) in the eastern Pacific. During2008–2012, the thermocline deepened in the western Pacific but shoaled in the east. In the Indian Ocean, there was only cooling(upper 150 m only) and freshening(almost the entire upper 400 m) during 2001–2007. The thermocline deepened during 2008–2012 in the Indian Ocean. Such change appeared from the equator to off the equator and even to the subtropics(about 20°N/S) in the two oceans. This long-term change of subsurface temperature/salinity may have been caused by change of the wind field over the two oceans during 2001–2012, in turn modifying OHC.     

环南澳大利亚太平洋-印度洋海洋间南极中层水环流

On the basis of the salinity distribution of isopycnal(σ_0=27.2 kg/m~3) surface and in salinity minimum, the Antarctic Intermediate Water(AAIW) around South Australia can be classified into five types corresponding to five regions by using in situ CTD observations. Type 1 is the Tasman AAIW, which has consistent hydrographic properties in the South Coral Sea and the North Tasman Sea. Type 2 is the Southern Ocean(SO) AAIW, parallel to and extending from the Subantarctic Front with the freshest and coldest AAIW in the study area. Type 3 is a transition between Type 1 and Type 2. The AAIW transforms from fresh to saline with the latitude declining(equatorward). Type 4, the South Australia AAIW, has relatively uniform AAIW properties due to the semienclosed South Australia Basin. Type 5, the Southeast Indian AAIW, progressively becomes more saline through mixing with the subtropical Indian intermediate water from south to north. In addition to the above hydrographic analysis of AAIW, the newest trajectories of Argo(Array for real-time Geostrophic Oceanography) floats were used to constructed the intermediate(1 000 m water depth) current field, which show the major interocean circulation of AAIW in the study area. Finally, a refined schematic of intermediate circulation shows that several currents get together to complete the connection between the Pacific Ocean and the Indian Ocean. They include the South Equatorial Current and the East Australia Current in the Southwest Pacific Ocean, the Tasman Leakage and the Flinders Current in the South Australia Basin, and the extension of Flinders Current in the southeast Indian Ocean.     

Global air-sea surface carbon-dioxide transfer velocity and flux estimated using ERS-2 data and a new parametric formula

Using data from the European remote sensing scatterometer（ERS-2） from July 1997 to August 1998,global distributions of the air-sea CO₂ transfer velocity and flux are retrieved.A new model of the air-sea CO₂ transfer velocity with surface wind speed and wave steepness is proposed.The wave steepness（5） is retrieved using a neural network（NN） model from ERS-2 scatterometer data,while the wind speed is directly derived by the ERS-2 scatterometer.The new model agrees well with the formulations based on the wind speed and the variation in the wind speed dependent relationships presented in many previous studies can be explained by this proposed relation with variation in wave steepness effect.Seasonally global maps of gas transfer velocity and llux are shown on the basis of the new model and the seasonal variations of the transfer velocity and llux during the 1 a period.The global mean gas transfer velocity is 30 cm/h after area-weighting and Schmidt number correction and its accuracy remains calculation with in situ data.The highest transfer velocity occurs around 60°N and 60°S,while the lowest on the equator.The total air to sea CO₂ llux（calculated by carbon） in that year is 1.77 Pg.The strongest source of CO₂ is in the equatorial east Pacific Ocean, while the strongest sink is in the 68°N.Full exploration of the uncertainty of this estimate awaits further data.An effectual method is provided to calculate the effect of waves on the determination of air-sea CO₂ transfer velocity and fluxes with ERS-2 scatterometer data.     

Argo时代东南印度洋上层盐度变化对海平面年代际演变的影响

In the past nearly two decades, the Argo Program has created an unprecedented global observing array with continuous in situ salinity observations, providing opportunities to extend our knowledge on the variability and effects of ocean salinity. In this study, we utilize the Argo data during 2004–2017, together with the satellite observations and a newly released version of ECCO ocean reanalysis, to explore the decadal salinity variability in the Southeast Indian Ocean(SEIO) and its impacts on the regional sea level changes. Both the observations and ECCO reanalysis show that during the Argo era, sea level in the SEIO and the tropical western Pacific experienced a rapid rise in 2005–2013 and a subsequent decline in 2013–2017. Such a decadal phase reversal in sea level could be explained, to a large extent, by the steric sea level variability in the upper 300 m. Argo data further show that, in the SEIO, both the temperature and salinity changes have significant positive contributions to the decadal sea level variations. This is different from much of the Indo-Pacific region, where the halosteric component often has minor or negative contributions to the regional sea level pattern on decadal timescale. The salinity budget analyses based on the ECCO reanalysis indicate that the decadal salinity change in the upper 300 m of SEIO is mainly caused by the horizontal ocean advection. More detailed decomposition reveals that in the SEIO, there exists a strong meridional salinity front between the tropical low-salinity and subtropical high salinity waters. The meridional component of decadal circulation changes will induce strong cross-front salinity exchange and thus the significant regional salinity variations.     

Upper ocean high resolution regional modeling of the Arabian Sea and Bay of Bengal

In this paper, effort is made to demonstrate the quality of high-resolution regional ocean circulation model in realistically simulating the circulation and variability properties of the northern Indian Ocean(10°S–25°N,45°–100°E) covering the Arabian Sea(AS) and Bay of Bengal(BoB). The model run using the open boundary conditions is carried out at 10 km horizontal resolution and highest vertical resolution of 2 m in the upper ocean.The surface and sub-surface structure of hydrographic variables(temperature and salinity) and currents is compared against the observations during 1998–2014(17 years). In particular, the seasonal variability of the sea surface temperature, sea surface salinity, and surface currents over the model domain is studied. The highresolution model's ability in correct estimation of the spatio-temporal mixed layer depth(MLD) variability of the AS and BoB is also shown. The lowest MLD values are observed during spring(March-April-May) and highest during winter(December-January-February) seasons. The maximum MLD in the AS(BoB) during December to February reaches 150 m (67 m). On the other hand, the minimum MLD in these regions during March-April-May becomes as low as 11–12 m. The influence of wind stress, net heat flux and freshwater flux on the seasonal variability of the MLD is discussed. The physical processes controlling the seasonal cycle of sea surface temperature are investigated by carrying out mixed layer heat budget analysis. It is found that air-sea fluxes play a dominant role in the seasonal evolution of sea surface temperature of the northern Indian Ocean and the contribution of horizontal advection, vertical entrainment and diffusion processes is small. The upper ocean zonal and meridional volume transport across different sections in the AS and BoB is also computed. The seasonal variability of the transports is studied in the context of monsoonal currents.     

SMAP卫星海表盐度产品的检验与校正

In this study, sea surface salinity(SSS) Level 3(L3) daily product derived from soil moisture active passive(SMAP)during the year 2016, was validated and compared with SSS daily products derived from soil Moisture and ocean salinity(SMOS) and in-situ measurements. Generally, the root mean square error(RMSE) of the daily SSS products is larger along the coastal areas and at high latitudes and is smaller in the tropical regions and open oceans. Comparisons between the two types of daily satellite SSS product revealed that the RMSE was higher in the daily SMOS product than in the SMAP, whereas the bias of the daily SMOS was observed to be less than that of the SMAP when compared with Argo floats data. In addition, the latitude-dependent bias and RMSE of the SMAP SSS were found to be primarily influenced by the precipitation and the sea surface temperature(SST). Then, a regression analysis method which has adopted the precipitation and SST data was used to correct the larger bias of the daily SMAP product. It was confirmed that the corrected daily SMAP product could be used for assimilation in high-resolution forecast models, due to the fact that it was demonstrated to be unbiased and much closer to the in-situ measurements than the original uncorrected SMAP product.     

热带印度洋降水、蒸发的时空特征及其对海表盐度的影响

本文利用降水、蒸发等资料分析热带印度洋年降水量、蒸发量、净淡水通量的分布特征，并选取4个典型海域来分析降水量、蒸发量、净淡水通量的季节变化和年际变化。结果表明:东印度洋的苏门答腊岛西部海域年降水量最大，季节变化较小，属全年降雨型；孟加拉湾的东北部和安达曼海的北部海域年降水量较大，其年际变化以4.2 mm/a的速率增长，强降水出现在5-9月；阿拉伯海的西部海域年降水量较小；南印度洋东部（20°~30°S，80°~110°E）海域年降水量较小，年蒸发量较大，年蒸发量在2000年之前以5.1 mm/a的速率增长，之后以4.5 mm/a的速率减小。本文还采用Argo盐度等资料探讨降水、蒸发对海表盐度的影响，研究结果表明:降水量远大于蒸发量的海域，海表盐度较低；降水量远小于蒸发量的海域，海表盐度较高。表层水平环流是导致高净淡水通量中心与低盐中心并不重合的主要原因，也是导致强蒸发中心与高盐中心并不重合的主要原因。选取的4个典型海域海表盐度的季节变化与净淡水通量关系不大，而是与表层水平环流有关。孟加拉湾强降水对表层盐度的影响显著，强降水发生后表层盐度降低0.2~0.8，其影响深度为30~50 m。     

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.     

印度洋海表盐度预报的线性马尔可夫模型及其改进办法

海洋盐度在水循环、海洋环流、海洋生态系统、全球天气和气候变化等方面起着至关重要的作用。然而, 受观测的限制, 以往对海洋盐度的研究相对匮乏, 对其进行预报的工作更为少见。本文采用线性马尔可夫模型对印度洋海表面盐度(sea surface salinity, SSS)开展初步的预报工作。根据混合层盐度收支方程, 选择海表面高度(sea surface height, SSH)、海表面温度 (sea surface temperature, SST)、SSS等物理量的异常值作为模型的组成部分, 对印度洋SSS开展预报工作。结果表明, 马尔可夫模型可提前9个月对印度洋SSS进行较好的预报。此外, 南太平洋海表面温度异常(sea surface temperature anomaly, SSTA), 海表面高度异常(sea surface height anomaly, SSHA)和印度洋偶极子(Indian Ocean dipole, IOD)系数等遥相关因素的加入可将线性马尔可夫预报对印度洋SSS的预报效果(相关系数)平均提高10%。利用改进的模型对印度洋SSS进行提前1~11个月的“实时”预测, 得出预报的SSS时空变化特征与观测场相吻合。综上所述, 改进的线性马尔可夫模型对印度洋SSS具有一定的预测能力, 未来可进一步完善。     

利用BP神经网络标准化东南太平洋茎柔鱼(Dosidicus gigas)资源丰度

东南太平洋茎柔鱼(Dosidicus gigas)是短生命周期大洋性经济鱼类,其资源量受环境因素变化的影响较大。根据我国鱿钓船队2013~2017年在东南太平洋的生产统计数据,结合海洋环境数据包括海表面温度(SST)、海表面盐度(SSS)、叶绿素a浓度(chl a),运用BP神经网络(back propagation network)模型来标准化单位捕捞努力量渔获量(catch per unit effort, CPUE,也称名义CPUE)。以均方误差(mean square errors, MSE)和平均相对变动值(average relative variances, ARV)为最优模型判断依据,比较隐含层节点3-10的神经网络模型,发现6-9-1结构为最优模型。用Garson算法解释模型结果,发现各输入层因子对东南太平洋茎柔鱼资源丰度影响重要度排序为chl a、SST、经度(Lon)、SSS、纬度(Lat)、月份(Month)。并作名义CPUE和标准化CPUE资源丰度对比分布图,结果显示CPUE与标准化CPUE总体分布状况基本一致,但局部区域存在明显差异, 80°~85°W及10°~20°S海域适宜鱿钓生产,表明BP神经网络模型可以适用于东南太平洋茎柔鱼的CPUE标准化,从而为鱿钓渔业生产提供一定参考依据。     

基于现场观测资料的卫星海表盐度网格化产品误差分析与质量评估

自欧洲土壤湿度和盐度卫星SMOS和美国宝瓶座盐度卫星Aquarius相继发射之后,多个数据中心发布了两颗卫星的海表盐度网格化产品,其中包括法国海洋研究院SMOS卫星数据小组发布SMOS Locean L3盐度产品、西班牙巴塞罗那专家中心发布SMOS BEC L4盐度产品和美国宇航局喷气动力实验室发布AquariusV3.0 CAP L3盐度产品。本文利用精确盐度现场观测资料从产品精度和模拟海洋现象能力两个方面对以上3种产品质量进行了评估。研究表明:(1) 在精度方面,与盐度现场资料相比,Aquarius CAP 产品质量最高,产品盐度偏差和均方根误差全年稳定且偏差较小,部分海域达到了设计精度;SMOS两种卫星产品在全球海域偏差较不稳定,个别月份出现异常偏差值;SMOS产品在低纬和开阔海域的数据质量相对较高,但在高纬海域仍存在较大误差,需要进一步提升;(2) 在刻画海洋现象方面,Aquarius产品在热带太平洋较好刻画了淡池东缘盐度锋,SMOS BEC产品的刻画能力次之,SMOS Locean产品在热带太平洋充满了小尺度噪音,描述物理现象方面表现偏差。     

基于SMOS卫星数据的海表面盐度模型

海表面盐度是研究海洋对全球气候影响以及大洋环流的重要参量之一,而卫星遥感技术是获取海表面盐度数据的最有效方法.目前,L波段的SMOS和Aquarius/SAC-D遥感卫星正在用于探测海表面盐度,并根据卫星观测数据和物理机制反演出海表面盐度的产品.但在某些近陆地区域,由于淡水流入及陆地射频(RFI)等因素影响,卫星反演盐度的产品精度较低.文中利用“东方红2号”科学考察船的实测数据、SMOS卫星数据,首次针对中国南海海域提出了用贝叶斯网络模型计算海表面盐度,并用验证数据集(实测Argo盐度)对模型进行适应性评估.经过计算,模型误差和验证误差分别为0.47 psu和0.45 psu,而相应的SMOS Level 2产品的精度分别为1.90 psu和1.82 psu.此模型为海表面盐度的计算提供了一个新方法.