Coupling winds to ocean surface currents over the global ocean

A Wind stress–Current Coupled System (WCCS) consisting of the HYbrid Coordinate Ocean Model (HYCOM) and an improved wind stress algorithm based on Donelan et al. [Donelan, W.M., Drennan, Katsaros, K.B., 1997. The air–sea momentum flux in mixed wind sea and swell conditions. J. Phys. Oceanogr. 27, 2087–2099] is developed by using the Earth System Modeling Framework (ESMF). The WCCS is applied to the global ocean to study the interactions between the wind stress and the ocean surface currents. In this study, the ocean surface current velocity is taken into consideration in the wind stress calculation and air–sea heat flux calculation. The wind stress that contains the effect of ocean surface current velocity will be used to force the HYCOM. The results indicate that the ocean surface velocity exerts an important influence on the wind stress, which, in turn, significantly affects the global ocean surface currents, air–sea heat fluxes, and the thickness of ocean surface boundary layer. Comparison with the TOGA TAO buoy data, the sea surface temperature from the wind–current coupled simulation showed noticeable improvement over the stand-alone HYCOM simulation.     

Measurement of ocean surface currents using a long-range,high-frequency ground wave radar

High-frequency (HF) ground wave radar (GWR) is emerging as a significant tool for monitoring ocean surface conditions at ranges well beyond the line-of-sight horizon that limits conventional systems. An experimental GWR system at Cape Race, Newfoundland, Canada that has been operational since 1991, has the ability to performing routine surveillance of oceanic surface parameters and surface target detection. Operating in the frequency range between 5 and 8 MHz, the frequency modulated interrupted continuous wave (FMICW) radar has a nominal range capability of 200 km. An experiment was performed during the period of October 20-November 21, 1992 to test the surface current measuring capability of the Cape Race system. Here, near real-time radial surface current information is extracted from the Doppler spectra of the radar time series data and a comparison is performed to the Lagrangian velocities derived from the position-time tracks of Accurate Surface Tracker (AST) drifters. A wide range of oceanic conditions were experienced during the experimental period, and favorable results were obtained from the comparison regardless of the sea state conditions. The analysis shows the standard deviation in the radar radial velocity component to be approximately 5 cm/s     

基于现场观测数据的HY-2A卫星有效波高校正研究

Significant wave height(SWH) can be computed from the returning waveform of radar altimeter, this parameter is only raw estimates if it does not calibrate. But accurate calibration is important for all applications, especially for climate studies. HY-2a altimeter has been operational since April 2012 and its products are available to the scientific community. In this work, SWH data from HY-2A altimeters are calibrated against in situ buoy data from the National Data Buoy Center(NDBC), Distinguished from previous calibration studies which generally regarded buoy data as "truth", the work of calibration for HY-2A altimeter wave data against in situ buoys was applied a more sophisticated statistical technique—the total least squares(TLS) method which can take into account errors in both variables. We present calibration results for HY-2A radar altimeter measurement of wave height against NDBC buoys. In addition, cross-calibration for HY-2A and Jason-2 wave data are talked over and the result is given.     

Mapping the sea surface using a GPS buoy

On May 22 and 24, 1995, a buoy, designed to float with the water surface and equipped with a GPS antenna, was deployed off the California coast at 16 locations near the Texaco oil platform, Harvest. The purpose of this deployment was threefold:.(1) to demonstrate the ability of this style of buoy to calibrate the TOPEXIPOSEIDON (TIP) altimeter range measurement as it overflew the platform: (2) to demonstrate the ability of the buoy to map the ocean's surface over a 10‐km‐diameter circle surrounding platform Harvest; and (3) to demonstrate the ability of the buoy to measure the sea state accurately. During the 1.6‐h period surrounding the time of the TIP overflight, the buoy‐measured sea level never differed by more than 1.5 cm from the sea level measured by the National Oceanic and Atmospheric Administration (NOAA) acoustic tide gauge on the platform. The good agreement demonstrated the capability of this style of buoy to calibrate altimetric satellites. A paraboloid was fitted to sea level from 16 buoy locations surrounding the platform with a 2.5‐cm rms residual. On a 10‐km‐diameter circle centered on the platform, the paraboloid was within 2.4‐cm rms of the Ohio State University Mean Sea Surface (OSUMSS95). H _u3 values calculated around the overflight times from the GPS buoy vertical positions had a mean difference of 2 cm and a standard deviation of 18 cm from values calculated from the University of Colorado (CU) pressure gauge system. At the time of the overflight, H _u3 was near 2 m, while 3‐m seas were observed by the CU pressure system during measurements later in the day. This experiment demonstrates that a simple wave‐rider buoy design can give comparable accuracies to that of more complex GPS platforms such as the University of Colorado's spar buoy, but is much easier to deploy and capable of being used in more severe weather conditions. Thus, such a buoy and derivative designs have great potential for calibrating altimetric experiments, and for oceanographic and geodetic mapping experiments.     

利用卫星观测海面信息反演三维温度场

基于历史观测的温盐剖面资料,采用回归分析方法统计出海面温度异常、海面动力高度异常与温度剖面异常之间的相关关系;然后利用高分辨率的卫星遥感海表面温度(SST)和卫星观测海面高度(SSH)信息重构了三维海洋温度场。在台湾岛周边海域建立了时间分辨率为天、空间分辨率为0.25°×0.25°的三维温度分析场。通过与实测资料的比较分析,文章所构建的分析场能够较好地描述海洋三维温度场的结构特征,能够较为真实地反映海洋的中尺度变化过程。该分析场可以作为海洋数值模式的初始场,也可以作为伪观测同化到海洋数值再分析和预报系统中,进而改善三维温、盐、流的数值再分析和预报。     

利用卫星测高、GRACE和GOCE资料估计全球海洋表面地转流

重力恢复和气候试验GRACE(gravity recovery and climate experiment)卫星极大地提高了地球重力场的精度和分辨率,特别是中长波分量,联合卫星测高数据可获得全球海洋表面大尺度洋流循环。另外,新一代地球重力和海洋环流探测卫星GOCE(gravity field and steady-state ocean circulation explorer)于2009年3月成功发射,采用卫星重力梯度测量原理,对重力场的高频部分非常敏感,使其高分辨率监测全球海洋循环成为可能。本文利用1~7年GRACE观测数据确定的重力场模型和18个月GOCE观测数据确定的地球重力场模型GO_CONS_GCF_2_TIM_R3,联合卫星测高确定的平均海面高模型MSS_CNES_CLS_11,分别估计全球海洋表面地转流,并且与实测浮标数据结果进行比较。分析表明GOCE重力卫星确定的重力场模型具有更高的空间分辨率,能够确定高精度和高空间分辨率的全球海洋地转流,如墨西哥湾暖流的细节和特征,并且与实测浮标结果基本一致。而基于1~4年GRACE观测资料的模型不能很好估计全球地转流特征,基于7年GRACE观测资料的重力场模型ITG-Grace2010s确定的全球地转流的精度仍低于18个月GOCE观测数据确定的地球重力场模型GO_CONS_GCF_2_TIM_R3的结果,估计的全球地转流仍含有较大的噪声,不能很好地反应中小尺度地转流细节特征。并计算ITG_Grace2010s和GOCE_TIM3的稳态海面地形和全球平均地转流的内符合精度,结果显示,在全球范围内,GOCE_TIM3的稳态海面地形和全球平均地转流的精度都比ITG_Grace2010s结果的精度有着很大的改善,其中ITG_Grace2010s的稳态海面地形的精度为21.6cm,而GOCE_TIM3的结果则为7.45cm,ITG_Grace2010s的全球平均地转流的精度为40.7cm/s,而GOCE_TIM3的结果则为19.6cm/s。     

High-frequency ocean radar derived characteristics of sea surface currents in the Ariake Sea,Japan

Composite analysis was conducted using high-frequency radar data obtained during 2006–2015 in order to gain a better understanding of the current field in the Ariake Sea. The seasonally averaged surface current in the Ariake Sea was directed southward in all seasons, except around river mouths during summer. Heavy rainfall enhanced the outflow along the eastern coast of the Shimabara Peninsula from Isahaya Bay to the southern area 2–5 days after heavy rainfall. Spring–neap differences were clearly seen in the southward current along the Shimabara Peninsula. Interannual variation in the M₂ tidal current amplitude was synchronized with the lunar nodal cycle.     

The long-term trend of the sea surface wind speed and the wave height (wind wave, swell, mixed wave)in global ocean during the last 44 a

Utilizing the 45 a European Centre for Medium-Range Weather Forecasts(ECMWF)reanalysis wave data(ERA-40),the long-term trend of the sea surface wind speed and(wind wave,swell,mixed wave)wave height in the global ocean at grid point 1.5×1.5 during the last 44 a is analyzed.It is discovered that a majority of global ocean swell wave height exhibits a significant linear increasing trend(2–8 cm/decade),the distribution of annual linear trend of the significant wave height(SWH)has good consistency with that of the swell wave height.The sea surface wind speed shows an annually linear increasing trend mainly concentrated in the most waters of Southern Hemisphere westerlies,high latitude of the North Pacific,Indian Ocean north of 30 S,the waters near the western equatorial Pacific and low latitudes of the Atlantic waters,and the annually linear decreasing mainly in central and eastern equator of the Pacific,Juan.Fernandez Archipelago,the waters near South Georgia Island in the Atlantic waters.The linear variational distribution characteristic of the wind wave height is similar to that of the sea surface wind speed.Another find is that the swell is dominant in the mixed wave,the swell index in the central ocean is generally greater than that in the offshore,and the swell index in the eastern ocean coast is greater than that in the western ocean inshore,and in year-round hemisphere westerlies the swell index is relatively low.     

基于浮标数据的卫星雷达高度计海浪波高数据评价与校正

卫星雷达高度计是海浪有效波高（significant wave height,SWH）观测的重要手段之一,本文利用时空匹配方法对T/P、Jason-1、Envisat、Jason-2、Cryosat-2和HY-2A共6颗卫星雷达高度计SWH数据与NDBC（National Data Buoy Center,NDBC）浮标SWH数据进行对比验证,并对雷达高度计SWH数据进行校正。全部卫星雷达高度计SWH数据时间跨度为1992年9月25日到2015年9月1日,对比验证NDBC浮标共53个,包括7个大洋浮标。精度评价发现除T/P外,各卫星雷达高度计SWH的RMSE都在0.4~0.5 m之间,经过校正后,RMSE都有显著下降,下降程度最大为13.82%;对于大洋浮标,评价结果RMSE在0.20~0.28 m之间,结果明显优于全部NDBC浮标的精度评价结果;HY-2A卫星雷达高度计SWH在经过校正后数据质量与国外其他5颗卫星雷达高度计SWH数据质量差异较小。     

卫星高度计海面风速的校准与验证

为了改善不同卫星高度计海面风速数据之间的一致性,以浮标数据为基准,对国内的HY-2A和国外的T/P、GFO、Jason-1、Envisat、Jason-2、CryoSat-2共7颗卫星高度计的海面风速数据进行了分析,给出了各个卫星高度计的海面风速校准公式,并对其校准效果进行了验证。验证结果表明:各个卫星高度计的海面风速在经过校准后,与浮标海面风速差异的均值和均方根都有所降低,其中HY-2A最为显著。经过校准后所有卫星高度计的海面风速与浮标海面风速差异的均值都在±0.2m/s以内。除了HY-2A、GFO和Jason-1,其余4颗卫星高度计校准后的海面风速与浮标海面风速差异的均方根都在1.6m/s以下。由此可以得出结论,利用本文的校准公式对各个卫星高度计(特别是HY-2A卫星高度计)的海面风速进行校准,可以有效减少其与浮标海面风速之间的差异。     

Maximum entropy principle and statistical distribution of ocean wave heights

Ｍａｘｉｍｕｍｅｎｔｒｏｐｙｐｒｉｎｃｉｐｌｅａｎｄｓｔａｔｉｓｔｉｃａｌｄｉｓｔｒｉｂｕｔｉｏｎｏｆｏｃｅａｎｗａｖｅｈｅｉｇｈｔｓ￥ＷｕＫｅｊｉａｎａｎｄＳｕｎＦｕ（ＲｅｃｅｉｖｅｄＭａｙ１５，１９９５；ａｃｃｅｐｔｅｄＪｕｎｅ２１，１９９５）Ａ...     

利用GPS拖体在万山区域测量平均海面技术研究

Wanshan area has been chosen to be the specified field to calibrate and validate(Cal/Val) the HY-2 altimeter and its follow-on satellites. In March 2018, an experiment has been conducted to determine the sea surface height(SSH) under the HY-2 A ground track(Pass No. 203). A GPS towing-body(GPS-TB) was designed to measure the SSH covering an area of about 6 km×28 km wide centered on the HY-2 A altimeter satellite ground track. Three GPS reference stations, one tide gauge and a GPS buoy were placed in the research area, in order to process and resolve the kinematic solution and check the precision of the GPS-TB respectively. All the GPS data were calculated by the GAMIT/GLOBK software and TRACK module. The sea surface was determined by the GPS-TB solution and the tide gauge placed on Zhiwan Island. Then the sea surface of this area was interpolated by Arc GIS10.2 with ordinary Kriging method. The results showed that the precision of the GPS-TB is about 1.10 cm compared with the tide gauge placed nearby, which has an equivalent precision with the GPS buoy. The interpolated sea surface has a bias of –1.5–4.0 cm with standard deviation of 0.2–2.4 cm compared with the checking line. The gradient of the measured sea surface is about 1.62 cm/km along the HY-2 orbit which shows a good agreement compared with the CLS11 mean sea surface(MSS). In the Cal/Val of satellites, the sea surface between the tide gauge/GPS buoy and the footprint of altimeter can be improved by this work.     

Absolute calibration of HY-2A and Jason-2 altimeters for sea surface height using GPS buoy in Qinglan,China

GPS buoy methodology is one of the main calibration methodologies for altimeter sea surface height calibration. This study introduces the results of the Qinglan calibration campaign for the HY-2A and Jason-2 altimeters. It took place in two time slices;one was from August to September 2014, and the other was in July 2015. One GPS buoy and two GPS reference stations were used in this campaign. The GPS data were processed using the real-time kinematic (RTK) technique. The fi nal error budget estimate when measuring the sea surface height (SSH) with a GPS buoy was better than 3.5 cm. Using the GPS buoy, the altimeter bias estimate was about -2.3 cm for the Jason-2 Geophysical Data Record (GDR) Version ‘D' and from -53.5 cm to -75.6 cm for the HY-2A Interim Geophysical Data Record (IGDR). The bias estimates for Jason-2 GDR-D are similar to the estimates from dedicated calibration sites such as the Harvest Platform, the Crete Site and the Bass Strait site. The bias estimates for HY-2A IGDR agree well with the results from the Crete calibration site. The results for the HY-2A altimeter bias estimated by the GPS buoy were verifi ed by cross-calibration, and they agreed well with the results from the global analysis method.     

Analysis of the global relationship of biennial variation of sea surface temperature and air-sea heat flux using satellite data

We investigated the phase difference and the cross correlation coefficient between the band-pass filtered biennial variations of sea surface temperature (SST) and air-sea heat flux estimated by the monthly mean 2°×2° satellite data of Advanced Very High Resolution Radiometer (AVHRR) and Special Sensor Microwave/Imager (SSM/I) from July 1987 to June 1991. Judging from the phase difference, it can be determined whether the biennial variation of SST is controlled by local thermal air-sea interaction or oceanic processes of horizontal transport. When the local air-sea heat flux controls the biennial variation of SST, the phase of SST advances /2 (6 months) against that of the air-sea heat flux. In contrast, when the biennial variation of SST is controlled by the oceanic process, the phase difference between the SST and the air-sea heat flux becomes 0 or (12 months). In this case, two types of the phase differences are determined, depending on which variability of SST and air-sea heat flux is larger. The close thermal air-sea interaction is noticeable in the tropics and in the western boundary current region. The phase difference of /2 appears mainly in the north Pacific, the southeast Indian Ocean, and the western tropical Pacific; zero in the eastern tropical Pacific and the northeast and equatorial Atlantic; and that of in the central equatorial Pacific and north of the intertropical convergence zone (ITCZ) of the Atlantic. Phase differences of 0, , or /2 are possible in the western boundary current regions. This fact indicates that each current plays a different role to the biennial variation of SST. It is inferred that SST anomalies in the tropics are mutually correlated, and the process in which marked SST anomalies in the tropics are transferred to the remote area was probed. In the equatorial Pacific, the SST anomaly is transferred by the long planetary wave. On the other hand, it is found from the phase relationship and the horizontal correlation of SST that the SST anomaly in the central and western equatorial Pacific is connected through atmospheric mediation. It is suggested that the biennial variation of SST in the eastern Indian Ocean is affected by heat transport due to the Indonesian throughflow from the western tropical Pacific. It is found that the mentioned pattern of the interannual variation of SST in the tropical Atlantic as a dipole is not tenable.     

Advanced sea surface temperature retrieval using the Japanese geostationary satellite,Himawari-6

Producing high-quality match-ups coupling the Japanese geostationary satellite, Himawari-6 (H6), and buoy SST observations, we have developed the new SST retrieval method. Kawamura et al. (2010) developed the previous version of SST product called MTSAT SST, which left several scientific/technical questions. For solving them, 6,711 algorithm tuning match-ups with precise navigation and 240,476 validation match-ups are generated for covering all seasons and wide ocean coverage. For discriminating the previous MTSAT SST, we call the new version of SST H6 SST. It is found that the SZA dependences of MTSAT SST algorithm are different from area to area of SZA > 40–50° N/S. The regionally different SZA dependences are treated by dividing the H6 disk coverage into five areas by the latitude lines of 40° N/S first and the longitude lines of 100° K and 180° K. Using the algorithm tuning match-ups, Nonlinear SST (NLSST) equations are derived for all of the five areas. Though the sun zenith angle dependency correction term is also examined, there is no significant regional difference. Therefore, this term is used in the H6 SST algorithm again. The new H6 SST equation is formed by the areal NLSST and the sun zenith angle dependency term for each area. The statistical evaluation of H6 SST using the validation match-ups show the small negative biases and the RMS errors of about 0.74° K for each area. For the full H6 disk, the bias is −0.1° K and the RMS error 0.74° K. The histogram of H6 SST minus the in situ SST for each area has a similar Gaussian shape with small negative skewness, and the monthly validation of H6 SST for each area is consistent with those for the whole period and the histograms     

卫星遥感海表温度资料和高度计资料的变分同化

利用国家气候中心正在发展的第二代全球海洋资料同化系统(BCC_GODAS2.0),针对多变量同化的协调性问题,发展了一种基于三维变分框架(3DVAR)下的高度计和海表温度(SST)相互约束的同化方法。该方法使海面高度和SST资料在同一个动力约束关系下进行同化。在一般方法中,海面高度和SST观测项是代价函数中2个独立的观测项,海面高度项引入动力高度计算公式,海表温度项用统计关系进行垂向投影。在代价函数的实际求解的计算过程中,虽然其总体积分效应受海面高度观测的约束,但整个水柱中各层温盐分析变量的调整是无序的。针对这个问题,文章提出一种新的同化方案。该方案将SST的观测项并入海面高度观测项中,海面高度的一部分,确切说是上层海洋部分,由SST决定,因此至少在SST的统计关系能影响到深度的上层海洋,在代价函数的求解过程中,温盐的调整是受较强的统计关系约束的,而这种统计关系的有效性已经在很多SST的同化试验中被其他学者广泛应用并证明。利用该方法,对1993—1997年的AVHRR卫星遥感海表温度资料进行变分同化试验,用TAO、OISST和SODA数据集进行检验证明,通过对卫星遥感资料的同化能够有效改进对海洋温度和盐度的估...     

Satellite monitoring of the sea surface temperature using the NOAA satellite data

A method of obtaining the operative sea surface temperature (SST)t using satellite scanner observations in the spectral ranges 3.53–3.94 m and 10.3–11.3 m is realized. The method represents a combination of McClainet al.'s formula (1983) and expressions suggested by the authors which describe the universal angular structures of the radiation temperature fields. The RMS error of reconstructingt at scanning angles of 0–55^o is equal to 0.2–0.3°C for atmospheric states corresponding to the SST variation within the limit 6–28°C. An atlas of temperature maps on the grid 0.5×0.5^o with temporal averaging from 5 days to 1 month is compiled using the data obtained on board the RVAkademik Vernadsky in the Atlantic Ocean in 1987–1989.Translated by Mikhail M. Trufanov.     

Effects of the seasonal variation in chlorophyll concentration on sea surface temperature in the global ocean

The effects of biological heating on the upper-ocean temperature of the global ocean are investigated using two ocean-only experiments forced by prescribed atmospheric fields during 1990–2007, on with fixed constant chlorophyll concentration, and the other with seasonally varying chlorophyll concentration. Although the existence of high chlorophyll concentrations can trap solar radiation in the upper layer and warm the surface, cooling sea surface temperature (SST) can be seen in some regions and seasons. Seventeen regions are selected and classified according to their dynamic processes, and the cooling mechanisms are investigated through heat budget analysis. The chlorophyll-induced SST variation is dependent on the variation in chlorophyll concentration and net surface heat flux and on such dynamic ocean processes as mixing, upwelling and advection. The mixed layer depth is also an important factor determining the effect. The chlorophyll-induced SST warming appears in most regions during the local spring to autumn when the mixed layer is shallow, e.g., low latitudes without upwelling and the mid-latitudes. Chlorophyll-induced SST cooling appears in regions experiencing strong upwelling, e.g., the western Arabian Sea, west coast of North Africa, South Africa and South America, the eastern tropical Pacific Ocean and the Atlantic Ocean, and strong mixing (with deep mixed layer depth), e.g., the mid-latitudes in winter.     

Surveillance of waste disposal activity at sea using satellite ocean color imagers: GOCI and MODIS

Korean Geostationary Ocean Color Imager (GOCI) and Moderate Resolution Imaging Spectroradiometer (MODIS) Aqua observations of the variation in ocean color at the sea surface were utilized to monitor the impact of nutrient-rich sewage sludge disposal in the oligotrophic area of the Yellow Sea. MODIS revealed that algal blooms persisted in the spring annually at the dump site in the Yellow Sea since year 2000 to the present. A number of implications of using products of the satellite ocean color imagers were exploited here based on the measurements in the Yellow Sea. GOCI observes almost every hour during the daylight period, every day since June 2011. Therefore, GOCI provides a powerful tool to monitor waste disposal at sea in real time. Tracking of disposal activity from a large tanker was possible hour by hour from the GOCI timeseries images compared to MODIS. Smaller changes in the color of the ocean surface can be easily observed, as GOCI resolves images at smaller scales in space and time in comparison to polar orbiting satellites, e.g., MODIS. GOCI may be widely used to monitor various marine activities in the sea, including waste disposal activity from ships.     

Studying the sea surface slopes using an array of wave gauge sensors

The deviations of the marine surface slope spectra (measured using an array of wave gauge sensors) from the theoretical estimates obtained using the linear spectral model of the wave field are analyzed. It has been indicated that the average measured full slope spectra (the sum of the slope component spectra in the orthogonal directions) is higher than the theoretical estimates by 6% at frequencies from the surface wave spectral peak (f _m ) to 4.5 f _m . The difference between the measured and theoretical estimates of the full slope spectrum rapidly increases at frequencies of f < f _m . At f _m ≈ 0.75 f _m , the average measured full slope spectrum is higher than the theoretical estimate by a factor of more than 5.