Improving Ocean State by Assimilating SARAL/AltiKa Derived Sea Level and Other Satellite-Derived Data in MITGCM

Assimilation of satellite-derived surface datasets has been explored in the study. Three types of surface data, namely sea level anomaly, sea surface temperature and sea surface salinity, have been used in various data assimilation experiments. The emphasis has been on the extra benefit arising out of the additional sea level assimilation and hence there are two parallel runs, in one of which sea level assimilation has been withheld. The model used is a state-of-the art ocean general circulation model (OGCM) and the assimilation method is the widely used singular evolutive extended Kalman filter (SEEK). Evaluation of the assimilation skill has been carried out by comparing the simulated depth of the 20°C isotherm with the same quantity measured by buoys and Argo floats. Simulated subsurface temperature and salinity profiles have also been compared with the same profiles measured by Argo floats. Finally, surface currents in the assimilation runs have been compared with currents measured by several off-equatorial buoys. Addition of sea level has been found to substantially improve the quality of simulation. An important feature that has been effectively simulated by the addition of sea level in the assimilation scheme is the near-surface temperature inversion (2-3°C) in the northern Bay of Bengal.     

Verification of an operational ocean circulation-surface wave coupled forecasting system for the China's seas

An operational ocean circulation-surface wave coupled forecasting system for the seas off China and adjacent areas(OCFS-C) is developed based on parallelized circulation and wave models. It has been in operation since November 1, 2007. In this paper we comprehensively present the simulation and verification of the system, whose distinguishing feature is that the wave-induced mixing is coupled in the circulation model. In particular, with nested technique the resolution in the China's seas has been updated to(1/24)° from the global model with(1/2)°resolution. Besides, daily remote sensing sea surface temperature(SST) data have been assimilated into the model to generate a hot restart field for OCFS-C. Moreover, inter-comparisons between forecasting and independent observational data are performed to evaluate the effectiveness of OCFS-C in upper-ocean quantities predictions, including SST, mixed layer depth(MLD) and subsurface temperature. Except in conventional statistical metrics, non-dimensional skill scores(SS) is also used to evaluate forecast skill. Observations from buoys and Argo profiles are used for lead time and real time validations, which give a large SS value(more than 0.90). Besides, prediction skill for the seasonal variation of SST is confirmed. Comparisons of subsurface temperatures with Argo profiles data indicate that OCFS-C has low skill in predicting subsurface temperatures between 100 m and 150 m. Nevertheless, inter-comparisons of MLD reveal that the MLD from model is shallower than that from Argo profiles by about 12 m, i.e., OCFS-C is successful and steady in MLD predictions. Validation of 1-d, 2-d and 3-d forecasting SST shows that our operational ocean circulation-surface wave coupled forecasting model has reasonable accuracy in the upper ocean.     

Seasonal variability of sonic layer depth in the Central Arabian Sea

The seasonal variability of sonic layer depth (SLD) in the central Arabian Sea (CAS) (0 to 25°N and 62-66°E) was studied using the temperature and salinity (T/S) profiles from Argo floats for the years 2002–2006. The atmospheric forcing responsible for the observed changes was explored using the meteorological data from NCEP/NCAR and Quickscat winds. SLD was obtained from sound velocity profiles computed from T/S data. Net heat flux and wind forcing regulated SLD in the CAS. Up-welling and down-welling (Ekman dynamics) associated with the Findlater Jet controlled SLD during the summer monsoon. While in winter monsoon, cooling and convective mixing regulated SLD in the study region. Weak winds, high insolation and positive net heat flux lead to the formation of thin, warm and stratified sonic layer during pre and post summer monsoon periods, respectively.     

Mixed layer depth climatology of the North Pacific based on Argo observations

A monthly mean climatology of the mixed layer depth (MLD) in the North Pacific has been produced by using Argo observations. The optimum method and parameter for evaluating the MLD from the Argo data are statistically determined. The MLD and its properties from each density profile were calculated with the method and parameter. The monthly mean climatology of the MLD is computed on a 2° × 2° grid with more than 30 profiles for each grid. Two bands of deep mixed layer with more than 200 m depth are found to the north and south of the Kuroshio Extension in the winter climatology, which cannot be reproduced in some previous climatologies. Early shoaling of the winter mixed layer between 20–30°N, which has been pointed out by previous studies, is also well recognized. A notable feature suggested by our climatology is that the deepest mixed layer tends to occur about one month before the mixed layer density peaks in the middle latitudes, especially in the western region, while they tend to coincide with each other in higher latitudes.     

Information-and-measuring drifter system for the monitoring of currents in the coastal zone of the sea

The autonomous drifting buoys equipped with satellite link turn into one of the most important components of the global system of operative observations of the ocean and the surface layer of the atmosphere. However, on the regional level, the problem of analysis of the surface circulation of waters in the coastal zone and sea straits remains quite urgent because the available drifters cannot be used in this case due to their sizes and long intervals between measurements. We present the results of development and testing of a new drifter system aimed at measuring currents. The system is based on the use of buoys operating at depths less than 1 m. To improve the space-and-time resolution of measurements, the buoys are equipped with receivers of the global positioning system (GPS) and GSM modems for the data transfer via cellular communication networks. The drifter system guarantees the possibility of determination of the coordinates of buoys with a resolution of 3 min in time and 14 m in space. We describe the specific features of the design of the proposed information-and-measuring drifter system and present the first results of application of new buoys called “minidrifters” for the pilot monitoring of currents in the Kerch Strait.     

基于观测的南海越南沿岸次表层涡旋

In this study, subsurface eddies near the Vietnam coast of the South China Sea were observed with in situ observations, including Argo, CTD, XBT and some processed and quality controlled data. Based on temperature profiles from four Argo floats near the coast of Vietnam, a subsurface warm eddy was identified in spring and summer. The multi-year Argo and Global Temperature and Salinity Profile Programme(GTSPP) data were merged on a seasonal basis based on the data interpolating variational analysis(DIVA) method to reconstruct the three-dimensional temperature structure. There is a warm eddy in the central subsurface at 12.5°N, 111°E below300 m depth in spring, which does not exist in autumn and is weak in winter and summer. From CSIRO Atlas of Regional Seas(CARS) and Generalized Digital Environment Model(GDEM) reanalysis data, this subsurface warm eddy is also verified in spring.     

Evidence for the Trapping and Amplification of Near-Inertial Motions in a Large Anticyclonic Ring in the Oyashio

Three ARGOS drift buoys were deployed in the Oyashio Current off the Kuril Islands near 45°N in fall, 1990, during a joint Russia/Canada study of western boundary current dynamics in the Subarctic Pacific Ocean. We here report on one buoy deployed within an anticyclonic warm core ring (WCR86B) which shows evidence of large amplitude inertial motions of near-diurnal frequency. During its first week within the ring the buoy drifted with a mean azimuthal current speed of 0.40–0.45 m s⁻¹ and a radius of rotation of 15–20 km. However, superimposed on the mean rotation of the ring at this time were “loops” of near-diurnal period, radius 7–8 km and speeds exceeding 1 m s⁻¹. During successive rotations the buoy spiraled outward, its mean period of rotation increased and the amplitude of the near-diurnal motions decreased. The large motions are explained by inertial wave trapping and amplification within the extremely large and weakly stratified eddy, wherein the negative vorticity of the eddy reduces the local inertial frequency to near-diurnal frequency. We here suggest that either tidal or wind forcing may generate these high-amplitude “loop” motions. This revised version was published online in August 2006 with corrections to the Cover Date.     

The relationship between the behavioral activity and the underwater vocalization of the beluga whale (<Emphasis Type="Italic">Delphinapterus leucas</Emphasis>)

The underwater vocalizations of the beluga whale summering in Onega Bay (64°24′N, 35°49′E) were recorded in June–July of 2008. The vocalizations were classified into five major whistle types, four types of pulsed tones, click series, and noise vocalizations. To determine the relationship between the behavioral activity and the underwater vocalizations, a total of fifty-one 2 minute-long samples of the audio records were analyzed in the next six behavioral contexts: directional movements, quiet swimming, resting, social interactions, individual hunting behavior, and the exploration of hydrophones by beluga whales. The overall vocalization rate and the percentage of the main types of signals depend on the behavior of the belugas. We suggest that one of the whistle types (the “stereotype whistle”) is used by belugas for long-distance communications, while other whistle types (with the exception of “squeaks”) and three types of pulsed tones (with the exception of “vowels”) are used for short distance communication. The percentage of “squeaks” and “vowels” was equally high in all the behavioral situations. Thus, we assume that “squeaks” are contact signals. “Vowels” have a specific physical structure and probably play a role in identification signals. A high rate of the click series was observed in the process of social interactions.     

The calibration of wind and wave model data in the Mediterranean Sea

Three sources of long-term wind and wave data are available in the Mediterranean Sea: numerical models, satellites and buoys. We make use of the overall information to obtain calibrated decadal time series at a large number of points, distributed at 0.5° intervals. We discuss the accuracy of the three sources and point out the errors that affect the accuracy of the final results.     

Estimation of eddy heat transport in the global ocean from Argo data

The Argo data are used to calculate eddy(turbulence)heat transport(EHT)in the global ocean and analyze its horizontal distribution and vertical structure.We calculate the EHT by averaging all the v′,T′profiles within each 2×2 bin.The velocity and temperature anomalies are obtained by removing their climatological values from the Argo"instantaneous"values respectively.Through the Student’s t-test and an error evaluation,we obtained a total of 87%Argo bins with significant depth-integrated EHTs(D-EHTs).The results reveal a positive-and-negative alternating D-EHT pattern along the western boundary currents(WBC)and Antarctic Circumpolar Current(ACC).The zonally-integrated D-EHT(ZI-EHT)of the global ocean reaches 0.12 PW in the northern WBC band and–0.38 PW in the ACC band respectively.The strong ZI-EHT across the ACC in the global ocean is mainly caused by the southern Indian Ocean.The ZI-EHT in the above two bands accounts for a large portion of the total oceanic heat transport,which may play an important role in regulating the climate.The analysis of vertical structures of the EHT along the 35 N and45 S section reveals that the oscillating EHT pattern can reach deep in the northern WBC regions and the Agulhas Return Current(ARC)region.It also shows that the strong EHT could reach 600 m in the WBC regions and 1 000 m in the ARC region,with the maximum mainly located between 100 and 400 m depth.The results would provide useful information for improving the parameterization scheme in models.     

Impact of wind/wind-stress field in the North Pacific constructed by ADEOS/NSCAT data

Data sets of surface wind and wind-stress fields in the North Pacific from September 1996 to June 1997 have been constructed using NASA Scatterometer (NSCAT) data on-board ADEOS to investigate their variability and implications for the wind-driven oceanic circulation. Using a weighting function decreasing with the distance between each grid and data points, and of Gaussian type for time, daily, 10-day and monthly averages are calculated for each 1°×1° grid. Products are validated by comparison with those calculated from in-situ measurement data at oceanic buoys around Japan (JMA) and in the equatorial area (TAO). The RMS differences for wind direction and speed never exceed 20° and 2 ms⁻¹, respectively, for the TAO buyos. This does not hold for data taken by JMA buoys, suggesting that the reliability in the mid-latitudes is not good for time averages shorter than several days. Zonal integration of the Sverdrup transport in a zone of 28°–30°N calculated from the monthly-mean products ranges between 25 and 60×10⁶ m³s⁻¹ (Sv) around its mean of 38 Sv. These are not so different from the Kuroshio transport values calculated from oceanic measurements.     

Impact of SARAL/AltiKa-Derived Sea Level Anomaly in a Data Assimilative Ocean Prediction System for the Indian Ocean

The impact of SARAL/AltiKa derived sea level anomaly (SLA) has been studied by assimilating it along with Jason-2 and Cryosat-2 SLA in the Princeton Ocean model (POM) using ensemble optimal interpolation (EnOI) technique. For isolating the extra benefit brought by SARAKL/Altika, a parallel run with assimilation of only Jason-2 and Cryosat-2 SLA has also been conducted. The importance of SARAL SLA in a data assimilative ocean prediction system has been evaluated with special emphasis on the improvement in thermocline depth, depth of the 20° isotherm, subsurface temperature and currents. Comparison with RAMA buoy has shown a positive impact of up to 13% for 20°C isotherm and up to 17% for thermocline depth after assimilating SARAL SLA. An overall improvement in temperature profile is also observed when compared with analogous profiles from RAMA buoys and Argo floats. Improvement in zonal currents away from the equator has also been noticed.     

Using Argo data to investigate the Meridional Overturning Circulation in the North Atlantic

Using a variety of oceanographic data, including direct volume transports in the Florida Strait, and Argo float profiles and drift velocities at 24°N and 36°N in the North Atlantic, inverse calculations are presented in which the net meridional transport, down to a depth of approximately 1600 m, is estimated at both latitudes for a 5-year period 2003–2007. The upper ocean is divided into seven layers using neutral density, and mass conservation constraints have been applied to a closed box bounded by these latitudes, including the Florida Strait. Ekman layer transports have been included in the top-most layer, and the inverse calculation has solved for changes from the initial reference velocities, Ekman and Florida Strait transports, given a priori estimates on the accuracy of each of these quantities. Solutions with and without transformations due to Mediterranean Water (MW) formation are made. Our results indicate that (1) time-averaged transport estimates derived from Argo have significant less eddy noise than individual hydrographic sections, (2) Argo drift velocities provide information to the inverse solution for the ocean interior, and (3) comparison of the total integrated interior mass transports in the thermocline waters for the period 2003–2007 with the previous estimates based on trans-ocean hydrographic sections shows that, within the errors of our estimation, the upper limb of the Atlantic Meridional Overturning Circulation has not significantly changed since 1957.     

Outbreak of warm water from the Kuroshio South of Japan

During the concentrated observation (April–May 1988) conducted as a part of the Ocean Mixed Layer Experiment (OMLET) in the sea area south of Japan, a conspicuous outbreak of warm water occurred from the large-meander region of the Kuroshio toward the southwest in the direction of the former Ocean Weather Station “T”. A series of NOAA-AVHRR infrared images clearly showed the process of this event. A surface buoy-mooring system deployed in this experiment recorded the arrival of this outbreak of water, in terms of the rise of sea-surface temperature (SST) of 1.5°C and the flow of warm water of 1.5kt toward the northwest at “T”. We studied this phenomenon by combining time series of infrared SST images with the oceanographic data obtained by two research vessels. The warm water was about 100 m deep in the section at 137°E along the edge of the Off-Shikoku Warm Water. It was estimated that about twenty outbreaks of this kind in a year can compensate a large heat loss to the atmosphere above this ocean region.     

Anticyclonic eddies in the northeastern South China Sea during winter 2003/2004

The origins and evolutions of two anticyclonic eddies in the northeastern South China Sea (SCS) were examined using multi-satellite remote sensing data, trajectory data of surface drifting buoys, and in-situ hydrographic data during winter 2003/2004. The results showed that buoy 22918 tracked an anti-cyclonic warm-core eddy (AE1) for about 20 days (December 4–23, 2003) in the northeastern SCS, and then escaped from AE1 eventually. Subsequently to that, buoy 22517 remained within a different anti-cyclonic warm-core eddy (AE2) for about 78 days (from January 28 to April 14, 2004) in the same area. It drifted southwestward for about 540 km, and finally entered into the so-called “Luzon Gyre”. Using inference from sea level anomaly (SLA), sea surface temperature (SST), geostrophic currents and the buoys’ trajectories, it is shown that both eddies propagated southwestward along the continental slope of the northern SCS. The mean speeds of AE1 and AE2 movements were 9.7 cm/s and 10.5 cm/s, respectively, which are similar to the phase speed of Rossby waves in the northern SCS. The variation of instantaneous speeds of the eddy movement and intensity of anticyclonic eddy may suggest complex interactions between an anticyclonic eddy and its ambient fluids in the northern SCS, where the eddy propagated southwestward with Rossby waves. Furthermore, SLA and SST images in combination with the temperature and salinity profiles obtained during a cruise suggested that AE1 was generated in the interior SCS and AE2 was shed from the “Kuroshio meander”.     

Decadal seesaw of the Central and Subtropical Mode Water formation associated with the Kuroshio Extension variability

Available Argo profiling float data from 2002 to 2011 were analyzed to examine the effect of the Kuroshio Extension (KE) current system variability on the formation of the Central Mode Water. Just north of the upstream portion of the KE at 140–152°E, formation of a lighter variety of the Central Mode Water in winter was active during the unstable period of the upstream KE in 2006–2009 and was reduced when the upstream KE was in the stable period of 2002–2005 and 2010–2011. This decadal formation variability is out of phase with that of the Subtropical Mode Water just south of the KE.     

基于Argo浮标和历史资料的热带西太平洋次表层与中层水年代变化分析

利用20世纪80年代和90年代WOD01(World Ocean Database2001)中的CTD温盐剖面资料和2000年以后Argo资料,对比分析了热带西太平洋次表层和中层水团分布的年代变化特征。分析结果表明,在这两个时期,起源于南北太平洋中高纬度海域的各次表层水和中层水,在热带西太平洋分布特征和交织在一起的总体态势基本一致,水团性质的年代变化不大。这与上述两个时段全球海洋-大气耦合系统趋于正常状态相吻合。通过辨识和跟踪表征次表层水性质的盐度极大值,发现南太平洋热带水沿西边界向北扩散程度有所加大,由前一时期的5°N,进一步扩散到6°~7°N;北太平洋热带水在西边界附近的向南扩散程度有所削弱,在2002-2005年间只向南扩散到4°N,而前一个时期则可向南扩散到2°N。通过辨识表征中层水性质的盐度极小值,南极中层水在西边界附近向北扩散程度有所加大,在2002-2005年到达13°N附近,而前一个时期只到达11°N;同期,北太平洋中层水在西边界附近的向南扩散程度有所削弱。上述年代变化与全球水循环强度的变化之间有何关系有待进一步研究。     

基于中国北极考察布放浮标观测和数值模拟的海冰和积雪厚度变化分析

Sea ice and the snow pack on top of it were investigated using Chinese National Arctic Research Expedition(CHINARE) buoy data.Two polar hydrometeorological drifters,known as Zeno? ice stations,were deployed during CHINARE 2003.A new type of high-resolution Snow and Ice Mass Balance Arrays,known as SIMBA buoys,were deployed during CHINARE 2014.Data from those buoys were applied to investigate the thickness of sea ice and snow in the CHINARE domain.A simple approach was applied to estimate the average snow thickness on the basis of Zeno~ temperature data.Snow and ice thicknesses were also derived from vertical temperature profile data based on the SIMBA buoys.A one-dimensional snow and ice thermodynamic model(HIGHTSI) was applied to calculate the snow and ice thickness along the buoy drift trajectories.The model forcing was based on forecasts and analyses of the European Centre for Medium-Range Weather Forecasts(ECMWF).The Zeno~ buoys drifted in a confined area during 2003–2004.The snow thickness modelled applying HIGHTSI was consistent with results based on Zeno~ buoy data.The SIMBA buoys drifted from 81.1°N,157.4°W to 73.5°N,134.9°W in 15 months during2014–2015.The total ice thickness increased from an initial August 2014 value of 1.97 m to a maximum value of2.45 m before the onset of snow melt in May 2015;the last observation was approximately 1 m in late November2015.The ice thickness based on HIGHTSI agreed with SIMBA measurements,in particular when the seasonal variation of oceanic heat flux was taken into account,but the modelled snow thickness differed from the observed one.Sea ice thickness derived from SIMBA data was reasonably good in cold conditions,but challenges remain in both snow and ice thickness in summer.