全球海洋环流模式中上层海洋对表面强迫的响应和调整 I. 年际变率

利用一个较高分辨率的全球海洋环流模式在COADS1945～1993年逐月平均资料的强迫下对海温和环流场进行了模拟试验,研究了全球热带海洋(主要是热带太平洋)海温和环流场的年际变化特征及模式ENSO冷暖事件演变的控制机理.结果表明,模式成功地再现了和观测一致的海温和环流的年际变化以及ENSO演变特征.其中热带印度洋年际SST变率的主要模态表现为与ENSO相联系的海盆尺度的一致性增暖或变冷现象,次级模态为热带印度洋偶极子模态;热带大西洋的SST年际变率表现为类ENSO的年际振荡现象.在热带太平洋,SST年际变化主要表现为ENSO型,环流的年际变率表现为与ENSO相对应的热带海洋质量循环圈的年际异常.对应于暖(冷)事件,前期赤道海洋垂直环流圈显示出减弱(增强)的特征.其中南赤道流异常的位相较Nino3区海温总体要超前5个月左右的时间;赤道上翻流异常的位相在表层要超前4个月,并随时间由上至下扩展;赤道潜流的异常则显示出东传特征,其中最早的较为显著的异常发生ENSO成熟前3个月180°附近.在模式ENSO冷暖事件的演变过程中,次表层海温异常沿赤道的东传起了关键作用,模式的ENSO模态主要表现为"时滞振子"模态.     

Local oceanic response to atmospheric forcing in the Gulf Stream region

The dominance of shifts in the location of the Gulf Stream (GS) in the local heat balance was observed in an hourly 15-month record of unprecedented surface mooring measurements at a site in the western North Atlantic occupied from November 2005 to January 2007. Instrumentation on the buoy provided a high quality record of air-sea exchanges of momentum, heat, and freshwater flux; and oceanographic sensors recorded ocean variability in the upper 640 m. The mooring was at times in the GS and at other times north of the GS. Our intent was to isolate the local oceanic response to the atmosphere from the influence of the GS shifts. A one-dimensional heat budget analysis indicated that the advective contribution from the GS shifts dwarfed the heat contribution by atmospheric forcing and therefore played the dominant role for upper oceanic thermal variability during the whole time record. A GS case study (i.e., when the surface mooring was in the GS), isolated the upper oceanic response to the atmospheric forcing in the GS and supported the critical role of GS shifts in total oceanic heat content. Through both an Empirical Orthogonal Function (EOF) analysis and by referencing temperatures to that observed at 200 m, the impact of GS shifts and atmospheric forcing were decomposed, allowing the local oceanic thermal response to be isolated. This local oceanic response was particularly prominent during the period of sustained heating during summer. A case study of summer conditions revealed a near surface flow consistent with Ekman dynamics within a shallow, warm ocean mixed layer.     

西北太平洋上层海洋对台风响应的个例研究

基于数字台风网、欧洲中心ERA-Interim、美国国家海洋与大气局以及中国Argo实时资料中心的资料研究了西北太平洋上层海洋对台风"奥鹿"的响应。研究结果表明,当"奥鹿"移动速度在2 m/s以下时,强风应力产生的Ekman泵是上层海洋响应的主要机制,移动速度越慢,Ekman抽吸速率(EPV)越大,海表温度(SST)降温持续时间短,冷尾迹出现在台风中心位置处。当"奥鹿"移动速度达到6 m/s以上时,持续风应力驱动的惯性泵是主导机制,SST降温持续时间长,冷尾迹出现在台风路径的右侧。惯性泵比Ekman泵持续的时间长,但Ekman泵影响深度比惯性泵大得多。在"奥鹿"经过西北太平洋时,混合层深度(MLD)变浅并伴随着"冷抽吸"作用的出现。上层海洋中"冷抽吸"现象较"热泵"现象影响深度深,持续时间长,在"奥鹿"过境后可持续20天以上。     

The upper ocean response to a moving typhoon

The upper ocean response to the translation speed of typhoons is studied using a three-dimensional primitive equation model. Similar models studied previously have applied stability criteria rather than the diffusion term to simulate the vertical mixing process. This study retains the diffusion term and uses the level-2 turbulence closure scheme to estimate the vertical eddy viscosity. The model results indicate that in the forced period, the mixed-layer temperature decrease is greater for a slow-moving storm due to stronger upwelling caused by the longer residence time. A fast-moving storm can attain a similar cooling intensity in the wake period if its residence time allows the wind to resonate with the current. The significant downward momentum diffusion and advection in the first few inertial periods of these events leads to strong, persistent inertial pumping throughout the upper ocean in the wake period. The mixed layer is further cooled by turbulent mixing supported by vertical current shears. Meanwhile, the upper thermocline exhibits a compensating temperature increase. The vertical transfer magnitude and penetration scale are smaller in the slow-moving case, when the inertial motion decays rapidly. The model results also indicate that the dominant cooling process can be inferred from the non-dimensional storm speed. However, this value may be misleading for rapidly moving storms in which the current response is so distant from the storm that little wind work is performed on the ocean.     

The contribution of atmospheric acid deposition to ocean acidification in the subtropical North Atlantic Ocean

The absorption of anthropogenic CO₂ and atmospheric deposition of acidity can both contribute to the acidification of the global ocean. Rainfall pH measurements and chemical compositions monitored on the island of Bermuda since 1980, and a long-term seawater CO₂ time-series (1983–2005) in the subtropical North Atlantic Ocean near Bermuda were used to evaluate the influence of acidic deposition on the acidification of oligotrophic waters of the North Atlantic Ocean and coastal waters of the coral reef ecosystem of Bermuda. Since the early 1980's, the average annual wet deposition of acidity at Bermuda was 15 ± 14 mmol m^− 2 year^− 1, while surface seawater pH decreased by 0.0017 ± 0.0001 pH units each year. The gradual acidification of subtropical gyre waters was primarily due to uptake of anthropogenic CO₂. We estimate that direct atmospheric acid deposition contributed 2% to the acidification of surface waters in the subtropical North Atlantic Ocean, although this value likely represents an upper limit. Acidifying deposition had negligible influence on seawater CO₂ chemistry of the Bermuda coral reef, with no evident impact on hard coral calcification.     

西北太平洋上层海洋对台风响应的合成分析研究

基于台风最佳路径资料和20a的Argo浮标资料,采用一种合成分析方法,得到了西北太平洋上层海洋的温度、盐度、溶解氧等三维要素场对台风的平均响应特征,结果表明:台风移动路径右侧海表面温度呈明显的负异常,即具有"右偏性",次表层在"热泵"和"冷抽吸"共同影响下,温度异常呈正负相间的分布,次表层以下大部区域以温度负异常为主,但"右偏性"已不明显;表层至1 000 m深度,各层占主导的盐度异常值呈"正-负-正"的垂向分布,次表层以下,各层的盐度异常在平行和垂直于台风路径方向上的分布较为均匀,不具有明显的各向性差异;表层至1 000 m深度层,溶解氧浓度呈"小-大-小"的垂向分布,各层的平面分布特征较为相似,均在平行于台风路径-2～2 R50范围存在一带状高值区域;温度、盐度、溶解氧浓度等要素对台风的平均响应深度可达到1 000 m以上.     

阿拉伯海上层海温的双峰特征

利用全球海洋Argo网格数据集、SODA月平均海洋数据集和CCMP风场数据,通过EOF分析,揭示了阿拉伯海5、50、100、200 m层海温全年2次增温、2次降温的双峰变化特征.结果表明,5 m层温度变化双峰信号出现在第一模态,其方差贡献率为75.79%,该信号主要受风场、太阳辐射及风生环流影响;50 m层温度变化双峰信号出现在第三模态,其方差贡献率为11.95%,该信号主要受风生环流影响;100 m层温度变化双峰信号出现在第一模态和第三模态,其中第一模态方差贡献率为52.03%,第三模态方差贡献率为9.55%.由100 m层第一模态可知,100 m层温度变化幅度最大、变化范围最广,是由于100 m层处于海洋温度变化最为剧烈的温跃层中.100 m层海温变化主要受风应力旋度(方向:向上为正)影响,风应力旋度为负时,大气对海洋的强迫导致局地海水辐合,温跃层加深,100 m层部分海域温度升高;风应力旋度为正时,大气对海洋的抽吸导致局地海水辐散,海洋深处的冷水上升,100 m层部分海域温度降低.     

Oceanic response to atmospheric forcing derived from different climatic data sets. Intercomparison study for the Black Sea

Available climatic and atmospheric analysis data have been used to prepare forcing functions for the Black Sea numerical model, based on the Bryan-Semtner-Cox Modular Ocean Model and including parameterizations for the atmosphere-ocean exchange, inflow through the strait of Bosphorus and the Mediterranean plume. Atmospheric data from different sources are compared and the drawbacks of some of them illustrated. A new wind stress data set, based on ship observations, is prepared. Compared to the existing wind stress estimates, the present ones use additional data and more accurate parameterization of the boundary layer physics. The intercomparison between forcing data sets is focused on the heat flux and freshwater flux at the sea surface.The model simulates adequately vertical stratification, seasonal variability and horizontal patterns. Five data sets for heat flux, freshwater flux and wind stress are used in different combinations to study the model response. The large differences between the simulations, forced by different wind stress and identical thermohaline forcing, justify the computation of the new wind stress. Though the forcing data used are perhaps close to the best available at the moment for the Black Sea, the model simulations range in large intervals and some of them are very poor. The model responses to forcing functions of different origin give rough estimates on the possible errors in present-day simulations. Some inconsistencies give clear indications that further verifications, improvements of the forcing functions, and intercomparisons between the responses simulated by the ocean circulation models are needed.     

Inferring upwelling rates in the equatorial Atlantic using Be measurements in the upper ocean

Ocean upwelling rates are difficult to measure because of the relatively small velocities involved, and therefore are typically inferred from indirect methods such as heat budget estimates or tracer observations. Here we present the first results using a novel technique, based on the isotope ⁷Be, to infer rates of upwelling along the equator. Beryllium-7 (half-life=53.3 d) is a cosmic-ray produced radioactive nuclide that is deposited by rainfall upon the ocean surface and subsequently enriched and homogenized within the mixed layer. Previous investigations have utilized the penetration of characteristically high mixed layer concentrations into the upper thermocline to trace ocean ventilation and subduction over seasonal timescales. Here, the tracer is used in a reverse sense; that is, the ⁷Be concentration in the usually ⁷Be-rich surface mixed layer will be diluted from penetration of ⁷Be “dead” water upwelled from below. This dilution provides a means to infer upwelling rates. Furthermore, with knowledge of upwelling rates, ⁷Be profiles can be used to constrain vertical diffusivity within the upper thermocline. These ideas were tested with samples collected during the Tropical Atlantic Climate Experiment (TACE) cruise (May 22-June 27, 2009). The observations indicated a nearly linear relationship between ⁷Be inventory and mixed layer temperature, as with increased upwelling, lower mixed layer temperatures correspond to greater ⁷Be dilution from depth. With this data, upwelling rates were estimated at a number of stations near the equator between 0°E and 30°W within and adjacent to the equatorial cold tongue. The derived upwelling rates ranged from 0 to 2.2 m/d, with maximum values found between the equator and 2°S. The corresponding K_z values derived for the upper thermocline were in the range 1-4×10⁻⁴ m²/s.     

A description of short period temperature fluctuations in the upper ocean

Dynamical properties of short-period temperature fluctuations are studied. Water temperature was measured continuously at several depths at the following stations: at 38°29.5′N, 141°35.8′E (100 m depth) on the continental shelf off Miyagi Prefecture in the summer of 1967, at 35°01.8′N, 139°0.8.5′E (100 m depth) in Sagami Bay in the summer of 1968, and at 32°32.2′N, 129°53.7′E (74 m depth) in Tachibana Bay in the summer of 1970. These measurements were made with a thermistor array laid down from the R. V.Tanseimaru (Ocean Research Institute, University of Tokyo) which was fixed with bow and stern anchors. Significant temperature fluctuations found at the first and the third stations are thought to be due to first mode internal waves having amplitude 3 to 5 m and period 5 to 20 minutes. The wave length of the waves is estimated to be 25 m to 400 m from the observed density structure. At the second station, we found second-mode internal waves. The period, amplitude and wave length of the waves are about 30 minutes, 1.3 m and 600 m, respectively. In all cases, the spectral density of the temperature fluctuations decreases with increase in frequency. However, the decrease obey neither the ?3 power law nor the ?5/3 power law. Coherences in the temperature fluctuations between two depths of measurement in the seasonal thermocline are significantly high in the range of frequencies lower than the local Brunt-Väisälä frequency, but are low in the higher frequency range. At the first and the third stations, the difference in the level of coherences between the lower frequencies and the higher frequencies are large. Phase differences between two depths in the thermocline are small in the lower frequency range. This suggests that the first-mode internal waves are predominant over higher-mode internal waves and over other disturbances.     

南海西北部上层海洋对台风“杜苏芮”的响应分析

基于锚碇观测资料,本文分析了南海西北部陆坡区上层海洋对台风“杜苏芮”的动力学和热力学响应特征。在动力学响应方面,台风“杜苏芮”期间上层流速显著增强,混合层纬向流速可达1.20 m/s;“杜苏芮”经过后上层海水运动以近惯性振荡为主(流向顺时针旋转周期在36～40 h之间)。近惯性能量在垂向分布上存在两个高值中心,分别位于混合层和温跃层深度上。近惯性能量耗散过程的e折时间尺度约为3.7 d,我们认为能量的向下传播在局地近惯性能量衰减过程中起主要作用。对能量谱的分析表明,“杜苏芮”作用期间近惯性频率能量相对于其作用前增大了约29.4倍,而全日和半日频率(K₁和M₂)能量有所减弱。此外,能量谱显示近惯性频率存在明显的“蓝移”现象,即对于纬向和经向流速分量在400 m以浅平均的近惯性振荡频率分别为1.167 f₀和1.170 f₀(f₀为局地惯性频率)。蓝移与近惯性内波的向下传播及正的相对涡度的输入有关。在热力学响应方面,上层海洋在台风的搅拌作用下,40～250 m深度均出现较小增温,最大增温幅度接近1°C;此外70 m以浅盐度的降低可能与台风过境时的降水相关,而Ekman抽吸引起的上升流则可能对70～100 m深度盐度的升高具有重要作用。     

上层海洋对热带气旋的响应与反馈研究进展

对60年来有关上层海洋与热带气旋(Tropical Cyclone,TC)的响应与反馈的研究进行了回顾,通过观测手段的完善和改进模式的应用,人们的认识不断提高:TC直接激发的近惯性流最大可达1 m/s,其导致的强烈的剪切造成混合层对下层冷水的夹卷是引起混合层降温的主要原因,并往往伴随着混合层深度的增加,这一影响在TC右侧更为最著,并可延续几天到几十天不等。TC导致的混合层降温会使得海洋输出的热通量减少,反过来削弱TC的强度,形成一个负反馈,而海洋特殊的热力和环流结构(如暖涡、洋流等)则对TC有正反馈。所以了解TC经过前的海洋初始场对研究TC与海洋之间的相互作用、对预测,TC的强度、路径变化等尤其重要;通过准确的初始场结合越来越完善的模式可以对TC进行更真实的模拟和预测,使得对TC准确的预报和预警成为可能。     

Persistent thermal fronts in the Black Sea: Existence,variability, and response to atmospheric forcing

High-resolution satellite sea surface temperature (SST) measurements (PATHFINDER dataset) have been used to confirm an existence of persistent thermal frontal zones in the Black Sea. Fronts have been revealed in the winter season between the intensively cooling shallow northwestern area and warmer waters of central part, and in spring-summer season between the actively warming waters of the southeastern part and cooler central area. Interannual and quasi-decadal variability of the northwestern front have been documented, as well as a high negative correlation of its intensity with the winter-mean SST.     

Interdecadal temperature variations in the North Pacific Central Mode Water simulated by an OGCM

The North Pacific Central Mode Water (CMW) is a water mass that forms in the Kuroshio-Oyashio Extension (KOE) region with characteristic low potential vorticity. Recent studies have suggested that the CMW, as low potential vorticity water, plays an important role in the adjustment of the subtropical gyre and subsurface variability on decadal to interdecadal timescales. We have forced a realistic ocean general circulation model (OGCM) with observed wind stress and sea surface temperature (SST) forcing to investigate the decadal variations of the CMW. Associated with the large atmospheric changes after the mid-1970s climate regime shift, the upper thermocline experiences a cooling as negative SST anomalies in the central North Pacific are subducted and advected southward. In addition to this thermodynamic response, the CMW’s path shifts anomalously eastward in response to anomalous Ekman pumping. This eastward shift of the core of the CMW produces a lowering of the isotherms, and a consequent warming, on the path of the CMW core. This warming partially counteracts the cooling associated with subducted surface anomalies, and it may be responsible for the reduced temperature variations at the climatological position of the CMW when both anomalous wind and heat fluxes are given. Lateral induction across the sloping bottom of the winter mixed layer in the KOE is critical to the formation of the low potential vorticity CMW. Coarse resolution models, which are widely used in climate modeling, underestimate the horizontal gradient of the mixed layer depth and form only a weak CMW or none at all. We have conducted a coarse resolution experiment with the same OGCM, showing that the subsurface response is much reduced. In particular, there is no dynamic warming in the CMW and the thermodynamic response to the SST cooling dominates. The resultant total response differs substantially from that in the finer resolution run where a strong CMW forms. This sensitivity to the model resolution corroborates the important dynamical role that the CMW may play with its distinctive low potential vorticity character and calls for its improved simulation.