海洋模式中垂直混合参数化方案介绍

介绍了海洋垂直混合过程参数化方案的发展,以及不同参数化方案在海洋模式中的应用情况。首先,介绍不同垂直混合参数化方案的物理问题、理论依据、数学表达和特征,并对不同参数化方案进行了比较。其次,针对中尺度涡、亚中尺度涡以及波浪、潮流混合参数化的最新研究进展进行了总结并对垂直混合参数化的未来发展提出了一些建议。     

The mixing mechanism in the formation of ocean shear waves

Ｔｈｅｍｉｘｉｎｇｍｅｃｈａｎｉｓｍｉｎｔｈｅｆｏｒｍａｔｉｏｎｏｆｏｃｅａｎｓｈｅａｒｗａｖｅｓ￥ＱｉａｏＦａｎｇｌｉ（ＦｉｒｓｔＩｎｓｔｉｔｕｔｅｏｆＯｃｅａｎｏｇｒａｐｈｙ，ＳｔａｔｅｏｃｅａｎｉｃＡｄｍｉｎｉｓｔｒａｔｉｏｎ，Ｑｉｎｇｄａｏ２...     

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.     

Effects of horizontal mixing on the upper ocean temperature in the equatorial Pacific Ocean

The influence of horizontal mixing on the thermal structure of the equatorial Pacific Ocean is examined based on a sigma coordinate model.In general,the distributions of the temperature and currents si...     

Gabriel T. Csanady: Understanding the physics of the ocean

Gabriel T. Csanady turned 80 in December 2005 and we celebrate it with this special Progress in Oceanography issue. It comprises 20 papers covering some of the many areas that Gabe contributed significantly throughout his professional career. In this introductory paper we briefly review Gabe’s career as an engineer, meteorologist and oceanographer, and highlight some of his major contributions to oceanography, both as a scientist as well as an educator. But we also use this opportunity to remember and thank Gabe, and his wife Joyce, for being such good friends and mentors to several generations of oceanographers. The authors of the collection of papers in this volume deserve special thanks for their efforts. We also are pleased to acknowledge the support of Progress in Oceanography’s editor, Detlef Quadfasel, and the many anonymous reviewers who generously contributed their time and expertise.     

南极西北威德尔海上层海洋湍流混合研究

Turbulent mixing in the upper ocean(30-200 m) of the northwestern Weddell Sea is investigated based on profiles of temperature,salinity and microstructure data obtained during February 2014.Vertical thermohaline structures are distinct due to geographic features and sea ice distribution,resulting in that turbulent dissipation rates(ε) and turbulent diffusivity(K) are vertically and spatially non-uniform.On the shelf north of Antarctic Peninsula and Philip Ridge,with a relatively homogeneous vertical structure of temperature and salinity through the entire water column in the upper 200 m,both ε and K show significantly enhanced values in the order of O(10~(-7))-O(10~(-6)) W/kg and O(10~(-3))-O(10~(-2)) m~2/s respectively,about two or three orders of magnitude higher than those in the open ocean.Mixing intensities tend to be mild due to strong stratification in the Powell Basin and South Orkney Plateau,where s decreases with depth from O(10~(-8)) to O(10~(-9)) W/kg,while K changes vertically in an inverse direction relative to s from O(10~(-6)) to O(10~(-5)) m~2/s.In the marginal ice zone,K is vertically stable with the order of10~(-4) m~2/s although both intense dissipation and strong stratification occur at depth of 50-100 m below a cold freshened mixed layer.Though previous studies indentify wind work and tides as the primary energy sources for turbulent mixing in coastal regions,our results indicate weak relationship between K and wind stress or tidal kinetic energy.Instead,intensified mixing occurs with large bottom roughness,demonstrating that only when internal waves generated by wind and tide impinge on steep topography can the energy dissipate to support mixing.In addition,geostrophic current flowing out of the Weddell Sea through the gap west of Philip Passage is another energy source contributing to the local intense mixing.     

A generalized coordinate ocean model and a comparison of the bottom boundary layer dynamics in terrain-following and in z-level grids

Sensitivity studies with a new generalized coordinate ocean model are performed in order to compare the behavior of bottom boundary layers (BBLs) when terrain-following (sigma or combined sigma and z-level) or z-level vertical grids are used, but most other numerical aspects remain unchanged. The model uses a second-order turbulence closure scheme that provides surface and BBL mixing and results in a quite realistic climatology and deep water masses after 100 year simulations with a coarse resolution (1° × 1°) basin-scale terrain-following grid. However, with the same turbulence scheme but using a z-level grid, the model was unable to produce dense water masses in the deep ocean. The latter is a known problem for coarse resolution z-level models, unless they include highly empirical BBL schemes.A set of dense water overflow experiments with high-resolution grids (10 and 2.5 km) are used to investigate the influence of model parameters such as horizontal diffusivity, vertical mixing, horizontal resolution, and vertical resolution on the simulation of bottom layers for the different coordinate systems. Increasing horizontal diffusivity causes a thinner BBL and a bottom plume that extends further downslope in a sigma grid, but causes a thicker BBL and limited downslope plume extension in a z-level grid. A major difference in the behavior of the BBL in the two grids is due to the larger vertical mixing generated by the turbulence scheme over the step-like topography in the z-level grid, compared to a smaller vertical mixing and a more stably stratified BBL in the sigma grid. Therefore, the dense plume is able to maintain its water mass better and penetrates farther downslope in the sigma grid than in the z-level grid. Increasing horizontal and vertical resolution in the z-level grid converges the results toward those obtained by a much coarser resolution sigma coordinate grid, but some differences remain due to the basic differences in the mixing process in the BBL.     

Flow of winter-transformed Pacific water into the Western Arctic

The dynamics of the flow of dense water through Barrow Canyon is investigated using data from a hydrographic survey in summer 2002. The focus is on the winter-transformed Bering water—the highest volumetric mode of winter water in the Chukchi Sea—which drains northward through the canyon in spring and summer. The transport of this water mass during the time of the survey was 0.2–0.3 Sv. As the layer flowed from the head of the canyon to the mouth, it sank, decelerated, and stretched. Strong cyclonic relative vorticity was generated on the seaward side of the jet, which compensated for the stretching. This adjustment was incomplete, however, in that it did not extend across the entire current, possibly because of internal mixing due to shear instabilities. The resulting vorticity structure of the flow at the canyon mouth was conducive for baroclinic instability and eddy formation. Multiple eddies of winter-transformed Bering water were observed along the Chukchi–Beaufort shelfbreak. Those to the west of Barrow Canyon were in the process of being spawned by the eastward-flowing shelfbreak current emanating from Herald Canyon, while the single eddy observed to the east originated from the Barrow Canyon outflow. It is argued that such an eddy formation is a major source of the ubiquitous cold-core anti-cyclones observed historically throughout the Canada Basin. Implications for the ventilation of the upper halocline of the Western Arctic are discussed.     

Surface layer mixing during the SAGE ocean fertilization experiment

Vessel-based observations of the oceanic surface layer during the 14-day 2004 SAGE ocean fertilization experiment were conducted using ADCP, CTD and temperature microstructure in a frame of reference moving with a patch of injected SF₆ tracer. During the experiment the mixed layer depth z_mld ranged between 50 and 80 m, with several re-stratifying events that brought z_mld up to less than 40 m. These re-stratifying events were not directly attributable to local surface-down development of stratification and were more likely associated with horizontal variation in density structure. Comparison between the CTD and a one-dimensional model confirmed that the SAGE experiment was governed by 3-d processes. A new method for estimating z_mld was developed that incorporates a component that is proportional to density gradient. This highlighted the need for well-conditioned near-surface data which are not always available from vessel-based survey CTD profiles. A centred-displacement scale, L_c, equivalent to the Thorpe lengthscale, reached a maximum of 20 m, with the eddy-centroid located at around 40 m depth. Temperature gradient microstructure-derived estimates of the vertical turbulent eddy diffusivity of scalar (temperature) material yielded bin-averaged values around 10⁻³ m² s⁻¹ in the pycnocline rising to over 10⁻² m² s⁻¹ higher in the surface layer. This suggests transport rates of nitrate and silicate at the base of the surface layer generate mixed layer increases of the order of 38 and 13 mmol/m²/day, respectively, during SAGE. However, the variability in measured vertical transport processes highlights the importance of transient events like wind mixing and horizontal intrusions.     

Three-dimensional characteristics of mesoscale eddies simulated by a regional model in the northwestern Pacific Ocean during 2000–2008

Mesoscale eddies play vital roles in ocean processes. Although previous studies focused on eddy surface features and individual three-dimensional (3D) eddy cases in the northwestern Pacific Ocean, the analysis of unique eddy 3D regional characteristics is still lacking. A 3D eddy detection scheme is applied to 9 years (2000–2008) of eddy-resolving Regional Ocean Modeling System (ROMS) output to obtain a 3D eddy dataset from the surface to a depth of 1 000 m in the northwestern Pacific Ocean (15°–35°N, 120°–145°E). The 3D characteristics of mesoscale eddies are analyzed in two regions, namely, Box1 (Subtropical Countercurrent, 15°–25°N, 120°–145°E) and Box2 (Southern Kuroshio Extension, 25°–35°N, 120°–145°E). In Box1, the current is characterized by strong vertical shear and weak horizontal shear. In Box2, the current is characterized by the strong Kuroshio, topographic effect, and the westward propagation of Rossby waves. The results indicate the importance of baroclinic instability in Box1, whereas in Box2, both the barotropic and baroclinic instability are important. Moreover, the mesoscale eddies’ properties in Box1 and Box2 are distinct. The eddies in Box1 have larger number and radius but a shorter lifetime. By contrast, Box2 has fewer eddies, which have smaller radius but longer lifetime. Vertically, more eddies are detected at the subsurface than at the surface in both regions; the depth of 650 m is the turning point in Box1. Above this depth, the number of cyclonic eddies (CEs) is larger than that of anticyclonic eddies (AEs). In Box2, the number of CEs is dominant vertically. Eddy kinetic energy (EKE) and mean normalized relative vorticity in Box2 are significantly higher than those in Box1. With increasing depth, the attenuation trend of EKE and relative vorticity of Box1 become greater than those of Box2. Furthermore, the upper ocean (about 300 m in depth) contains 68.6% of the eddies (instantaneous eddy). Only 16.6% of the eddies extend to 1 000 m. In addition, about 87% of the eddies are bowl-shaped eddies in the two regions. Only about 3% are cone-shaped eddies. With increasing depth of the eddies, the proportion of bowl-shaped eddies gradually decreases. Conversely, the cone- and lens-shaped eddies are equal in number at 700–1 000 m, accounting for about 30% each. Studying the 3D characteristics of eddies in two different regions of the northwestern Pacific Ocean is an important stepping stone for discussing the different eddy generation mechanisms.     

On the growth of ocean waves

The availability of 10 h of continuous, uninterrupted field measurements of wind waves recorded in the western Pacific and containing a complete wave growth episode, has provided a distinct opportunity for us to make a novel, unprecedented examination of detailed wave growth processes. We found that the significance of the size of data used in the measurement, which can only be addressed with continuous and uninterrupted measurements, reflected the ineptness of the conventional approach toward further detailed understanding of realistic wave growth processes, as the conventional 20 min data size essentially stamped out any dynamics with time scale below 20 min. While our conventional understanding and modeling were generally operative and useful, they left no real vestige on time localized mechanisms such as wave grouping or wave breaking processes all with time scales much less than 20 min.     

The mesoscale variability in the Caribbean Sea. Part II: Energy sources

The processes which drive the production and the growth of the strong mesoscale eddy field in the Caribbean Sea are examined using a general circulation model. Diagnostics of the simulations suggest that:(1) The mean currents in the Caribbean Sea are intrinsically unstable. The nature of the instability and its strength vary spatially due to strong differences of current structure among basins.(2) The greatest and most energetic eddies of the Caribbean Sea originate in the Venezuela Basin by mixed barotropic-baroclinic instability of an intense jet, formed with waters mostly from the surface return flow of the Meridional Overturning Circulation and the North Equatorial Current which converge and accelerate through the Grenada Passage. The vertical shear of this inflow is enhanced by an eastward undercurrent, which flows along the south American Coast between 100 and 250 m depth. The shallow eddies (less than 200 m depth) formed in the vicinity of the Grenada Passage get rapidly deeper (down to 1000 m depth) and stronger by their interaction with the deep interior flow of the Subtropical Gyre, which enters through passages north of St. Lucia. These main eastern Caribbean inflows merge and form the southern Caribbean Current, whose baroclinic instability is responsible for the westward growth and strengthening of these eddies from the Venezuela to the Colombia Basin.(3) Eddies of lesser strength are produced in other regions of the Caribbean Sea. Their generation and growth is also linked with instability of the local currents. First, cyclones are formed in the cyclonic shear of the northern Caribbean Current, but appear to be rapidly dissipated or absorbed by the large anticyclones coming from the southern Caribbean. Second, eddies in the Cayman Sea, which impact the Yucatan region, are locally produced and enhanced by barotropic instability of the deep Cayman Current.(4) The role of the North Brazil Current (NBC) rings is mostly to act as a finite perturbation for the instability of the mean flow. Their presence near the Lesser Antilles is ubiquitous and they appear to be linked with most of the Caribbean eddies. There are some evidences that the frequency at which they form near the Grenada Passage is influenced by the frequency at which the NBC rings impinge the Lesser Antilles. But large Caribbean eddies also form without a close influence of any ring, and comparison between simulations shows that mean eddy kinetic energy and eddy population in the Caribbean Sea are not substantially different in absence or presence of NBC rings: their presence is not a necessary condition for the generation and growth of the Caribbean eddies.     

春季黄海海-气湍流通量交换的特征

本文利用１９９６年４月中国和韩国第一次同时在黄海海域的调查资料,选用两种整体参数化公式分析春季黄海海气湍流通量交换特征。分析结果表明：动量通量和感热通量的高值区皆位于山东半岛和朝鲜西岸的近岸海域,而低值区则出现在长江口东北近海;在平稳天气条件下,表层各通量具有日变化特征,而强风条件下的感热、潜热和动量通量分别是平稳天气条件下的５．７倍、６．６倍和７倍。     

A fast inversion method for ocean parameters based on dispersion curves with a single hydrophone

The dispersion characteristics of shallow water can be described by the dispersion curves, which contain substantial ocean parameter information. A fast ocean parameter inversion method based on dispersion curves with a single hydrophone is presented in this paper. The method is achieved through Bayesian theory. Several sets of dispersion curves extracted from measured data are used as the input function. The inversion is performed by matching a replica calculated with a dispersion formula. The bottom characteristics can be described by the bottom reflection phase shift parameter P. The propagation range and the depth can be inverted quickly when the seabed parameters are represented by on parameter P. The inversion results improve the inversion efficiency of the seabed parameters. Consequently, the inversion efficiency and accuracy are improved while the number of inversion parameters is decreased and the computational speed of replica is increased. The inversion results have lower error than the reference values, and the dispersion curves calculated with inversion parameters are also in good agreement with extracted curves from measured data; thus, the effectiveness of the inversion method is demonstrated.     

The release of dissolved actinium to the ocean: A global comparison of different end-members

The measurement of short-lived ²²³Ra often involves a second measurement for supported activities, which represents ²²⁷Ac in the sample. Here we exploit this fact, presenting a set of 284 values on the oceanic distribution of ²²⁷Ac, which was collected when analyzing water samples for short-lived radium isotopes by the radium delayed coincidence counting system. The present work compiles ²²⁷Ac data from coastal regions all over the northern hemisphere, including values from ground water, from estuaries and lagoons, and from marine end-members. Deep-sea samples from a continental slope off Puerto Rico and from an active vent site near Hawaii complete the overview of ²²⁷Ac near its potential sources.The average ²²⁷Ac activities of nearshore marine end-members range from 0.4 dpm m^− 3 at the Gulf of Mexico to 3.0 dpm m^− 3 in the coastal waters of the Korean Strait. In analogy to ²²⁸Ra, we find the extension of adjacent shelf regions to play a substantial role for ²²⁷Ac activities, although less pronounced than for radium, due to its weaker shelf source. Based on previously published values, we calculate an open ocean ²²⁷Ac inventory of 1.35 * 10¹⁸ dpm ²²⁷Ac_ex in the ocean, which corresponds to 37 moles, or 8.4 kg. This implies a flux of 127 dpm m⁻² y^− 1 from the deep-sea floor. For the shelf regions, we obtain a global inventory of ²²⁷Ac of 4.5 * 10¹⁵ dpm, which cannot be converted directly into a flux value, as the regional loss term of ²²⁷Ac to the open ocean would have to be included.Ac has so far been considered to behave similarly to Ra in the marine environment, with the exception of a strong Ac source in the deep-sea due to ²³¹Pa_ex. Here, we present evidence of geochemical differences between Ac, which is retained in a warm vent system, and Ra, which is readily released [Moore, W.S., Ussler, W. and Paull, C.K., 2008-this issue. Short-lived radium isotopes in the Hawaiian margin: Evidence for large fluid fluxes through the Puna Ridge. Marine Chemistry]. Another potential mechanism of producing deviations in ²²⁷Ac/²²⁸Ra and daughter isotope ratios from the expected production value of lithogenic material is observed at reducing environments, where enrichment in uranium may occur. The presented data here may serve as a reference for including ²²⁷Ac in circulation models, and the overview provides values for some end-members that contribute to the global Ac distribution.     

The mesoscale variability in the Caribbean Sea. Part I: Simulations and characteristics with an embedded model

The variability in the Caribbean Sea is investigated using high resolution (1/15°) general circulation model experiments. For the first time in this region, simulations were carried out with a 2-way nested configuration of the NEMO primitive equation model. A coarse North Atlantic grid (1/3°) reproduces the main features of the North Atlantic and Equatorial circulation capable of influencing ocean dynamics in the Caribbean Sea. This numerical study highlights strong dynamical differences among basins and modifies the view that dynamics are homogeneous over the whole Caribbean Basin. The Caribbean mean flow is shown to organize in two intense jets flowing westward along the northern and southern boundaries of the Venezuela Basin, which merge in the center of the Colombia Basin. Diagnostics of model outputs show that width, depth and strength of baroclinic eddies increase westward from the Lesser Antilles to the Colombia Basin. The widening and strengthening to the west is consistent with altimetry data and drifter observations. Although influenced by the circulation in the Colombia Basin, the variability in the Cayman Basin (which also presents a westward growth from the Chibcha Channel) is deeper and less energetic than the variability in the Colombia/Venezuela Basins. Main frequency peaks for the mesoscale variability present a westward shift, from roughly 50 days near the Lesser Antilles to 100 days in the Cayman Basin, which is associated with growth and merging of eddies.     

The role of surface waves in the ocean mixed layer

Previously, most ocean circulation models have overlooked the role of the surface waves. As a result, these models have produced insufficient vertical mixing, with an under - prediction of the ,nixing layer （ML） depth and an over - prediction of the sea surface temperature （SST）, particularly during the summer season. As the ocean surface layer determines the lower boundary conditions of the atmosphere, this deficiency has severely limited the performance of the coupled ocean - atmospheric models and hence the climate studies. To overcome this shortcoming, a new parameterization for the wave effects in the ML model that will correct this systematic error of insufficient mixing. The new scheme has enabled the mixing layer to deepen, the surface excessive heating to be corrected, and an excellent agreement with observed global climatologic data. The study indicates that the surface waves are essential for ML formation, and that they are the primer drivers of the upper ocean dynamics; therefore, they are critical for climate studies.     

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.     

National ocean policy—new opportunities for Malaysian ocean development

Being a littoral state, Malaysian maritime sectors provide the basis for the growth and development of marine-related activities. As the country enters the new millennium along with the ‘2020 Vision’ which serves as the basic template against which the economic performance of these sectors should be judged, it poses a new challenge for Malaysia in developing these sectors. This paper attempts to address three basic issues: Firstly, is there any national ocean policy in place? Secondly, are the maritime sectors adequately organized to play an important role in achieving the objectives? Thirdly, what are the opportunities available in developing the ocean?It can be seen that existing policies to develop comprehensive ocean governance have not received the full attention they deserve. Organizational structures governing the ocean for implementing national policies are well in place but in a fragmented and uncoordinated fashion. As a result, sectoral and intersectoral management problems were created such as multiple-use conflicts, overlapping of jurisdiction and duplication of efforts. Environmental problems have also not been properly addressed. A few sectors have been identified as offering opportunities to further develop the Malaysian maritime areas. Among these are marine education and human resources development, marine tourism and the seafood industry.     

Modelling ocean waves in ice-covered seas

Over the past decade there has been a rapid growth of interest in wave propagation through ice covers. This paper summarizes the author’s observation of the modeling efforts on this topic. Models can be theory-based, data-driven, or a combination of the two. A pure data-driven model relies on a large amount of observations and is only becoming available recently. Theory-based models on the other hand have a long history. They are always a simplified version of the reality. As our knowledge grows, theories become more complicated. A theory for waves-in-ice that captures all possible processes does not exist. However, when integrated with observation through calibration, these combined theory + data-based models may be used with some confidence. In this paper, different models, their basic concepts, their calibration and validation are discussed. The present theory-based models do not have the correct spectral attenuation trend as observed from field or laboratory experiments. Hence, through calibration they may fit different parts of the wave spectra but not all. Pure data-driven models can reproduce the correct trend, but its dependability outside the situation where the data are collected is uncertain. In addition to offering tools to forecast waves-in-ice, these model building and validating efforts point to missing mechanisms that should be carefully studied. Despite the many challenges towards building a satisfactory general waves-in-ice model, significant progress has been made for models that work reasonably well in the marginal ice zone. We anticipate much more data will become available in the coming years to help us improve the existing models.