首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 842 毫秒
1.
《Ocean Modelling》1999,1(1):1-15
In regions of sloping isopycnals, isopycnal mixing acting in conjunction with biological cycling can produce patterns in the nutrient field which have negative values of tracer in light water and unrealistically large values of tracer in dense water. Under certain circumstances, these patterns can start to grow unstably. This paper discusses why such behavior occurs. Using a simple four-box model, it demonstrates that the instability appears when the isopycnal slopes exceed the grid aspect ratio (Δz/Δx). In contrast to other well known instabilities of the CFL type, this instability does not depend on the time step or time-stepping scheme. Instead it arises from a fundamental incompatibility between two requirements for isopycnal mixing schemes, namely that they should produce no net flux of passive tracer across an isopycnal and everywhere reduce tracer extrema. In order to guarantee no net flux of tracer across an isopycnal, some upgradient fluxes across certain parts of an isopycnal are required to balance downgradient fluxes across other parts of the isopycnal. However, these upgradient fluxes can cause local maxima in the nutrient field to become self-reinforcing. Although this is less of a problem in larger domains, there is still a strong tendency for isopycnal mixing to overconcentrate tracer in the dense water. The introduction of eddy-induced advection is shown to be capable of counteracting the upgradient fluxes of nutrient which cause problems, stabilizing the solution. The issue is not simply a numerical curiosity. When used in a GCM, different parameterizations of eddy mixing result in noticeably different distributions of nutrient and large differences in biological production. While much of this is attributable to differences in convection and circulation, the numerical errors described here may also play an important role in runs with isopycnal mixing alone.  相似文献   

2.
Ventilation of the deep basins of the North Aegean Sea takes place during relatively scarce events of massive dense water formation in that region. In the time intervals between such events, the bottom waters of each sub-basin are excluded from interaction with other water masses through advection or isopycnal mixing and the only process that changes their properties is diapycnal mixing with overlying waters. In this work we utilize a simple one-dimensional model in order to estimate the vertical eddy diffusion coefficient Kρ based on the observed rate of change of density and stratification. Vertical diffusivity is estimated for each of three sub-basins of the North Aegean, one of convex shape of the seabed and the other two of concave topography. It is noteworthy that the convex sub-basin exhibited much higher vertical diffusivity than the two concave sub-basins, a fact consistent with theoretical predictions that internal-wave-induced mixing is higher over the former shape of seabed. Furthermore, the estimates of Kρ are exploited in computing the vertical transport of dissolved oxygen through diffusion and the rate of oxygen consumption by decaying organic matter. The different levels of the estimated diffusion and oxygen consumption rates testify to the dynamical and biogeochemical characteristics of each basin.  相似文献   

3.
《Ocean Modelling》2002,4(1):55-88
A newly developed hybrid-coordinate ocean circulation model is documented and tested. Coordinate surfaces in this model adhere to isopycnals wherever this does not violate minimum layer thickness requirements; elsewhere, coordinate surfaces are geometrically constrained. The intent of this approach, some of whose features are reminiscent of the Arbitrary Lagrangian–Eulerian (ALE) technique, is to combine the best features of isopycnic-coordinate and fixed-grid circulation models within a single framework. The hybrid model is an offshoot of the Miami Isopycnic Coordinate Ocean Model whose solutions, obtained under identical geographic and forcing conditions, serve as reference. Century-scale simulations on a coarse-mesh near-global domain show considerable similarities in the modeled thermohaline-forced circulation. Certain architectural details, such as the choice of prognostic thermodynamic variables (ρ,S versus T,S) and the algorithm for moving coordinate surfaces toward their reference isopycnals, are found to only have a minor impact on the solution. Emphasis in this article is on the numerical resiliency of the hybrid coordinate approach. Exploitation of the model's flexible coordinate layout in areas of ocean physics where pure isopycnic coordinate models only have limited options, such as mixed-layer turbulence parameterization, will be the subject of forthcoming articles.  相似文献   

4.
Study of oceanic circulation and climate requires models which can simulate tracer eddy diffusion and ad vection accurately. It is shown that the traditional Eulerian coordinates can introduce large artificial hori zontal diffusivity/viscosity due to the incorrect alignment of the axis. Therefore, such models can smear sharp fronts and introduce other numerical artifacts. For simulation with relatively low resolution, large lateral diffusion was explicitly used in models; therefore, such numerical diffusion may not be a problem. However, with the increase of horizontal resolution, the artificial diffusivity/viscosity associated with hori zontal advection in the commonly used Eulerian coordinates may become one of the most challenging ob stacles for modeling the ocean circulation accurately. Isopycnal eddy diffusion （mixing） has been widely used in numerical models. The common wisdom is that mixing along isopycnal is energy free. However, a careful examination reveals that this is not the case. In fact, eddy diffusion can be conceptually separated into two steps： stirring and subscale diffusion. Due to the thermobaric effect, stirring, or exchanging water masses, along isopycnal surface is associated with the change of GPE in the mean state. This is a new type of instability, called the thermobaric instability. In addition, due to cabbeling subscale diffusion of water parcels always leads to the release of GPE. The release of GPE due to isopycnal stirring and subscale diffusion may lead to the thermobaric instability.  相似文献   

5.
An air-sea buoyancy flux out of the ocean between the surface outcroppings of different isopycnals must be balanced by a convergence of advective and diffusive fluxes of buoyancy across those isopycnals (Walin, 1982; Tziperman, 1986; Garrett et al., 1995). For steady conditions, the diapycnal diffusive flux due to vertical mixing in the surface mixed layer is very small, so that the advective buoyancy flux dominates (Speer, 1993; Garrett et al., 1995). The associated advective buoyancy flux can then be used to estimate the volume flux of water out of the base of the surface mixed layer. The resulting thermodynamic algorithm provides a valuable estimate of water mass formation in the ocean.In contrast, for the time-dependent real ocean with horizontal and vertical gradients of the horizontal buoyancy gradient, diurnal and seasonal mixed layer deepening and entrainment in the presence of a buoyancy jump at the base of the mixed layer contributes to the annual volume flux out of the base of the deepest (wintertime) mixed layer. The mismatch between the predictions of the ideal algorithm and measured rates of water mass formation (Speer, 1997) may thus be partly due to mixed layer processes rather than diapycnal mixing in the thermocline.  相似文献   

6.
Few basins in the world exhibit such a wide range of water properties as those of the Nordic Seas with cold freshwaters from the Arctic in the western basins and warm saline waters from the Atlantic in the eastern basins. In this study we present a 50-year hydrographic climatology of the Nordic Seas in terms of depth and temperature patterns on four upper ocean specific volume anomaly surfaces. This approach allows us to better distinguish between change due to variations along such surfaces and change due to depth variations of the stratified water column. Depth variations indicate changes in the mass field while property variations along isopycnals give insight into isopycnal advection and mixing, as well as diapycnal processes. We find that the warmest waters on each surface are found in the north, close to where the isopycnal outcrops, a clear indication of downward mixing of the warmer, more saline waters on shallower isopycnals due to convective cooling at the surface. These saline waters come from the Norwegian Atlantic Slope Current by means of a very high level of eddy activity in the Lofoten Basin.The isopycnal analyses further show that the principal water mass boundary between the waters of Arctic origin in the west and Atlantic waters in the east aligns quite tightly with the Jan Mayen, Mohn, Knipovich Ridge system suggesting little cross-ridge exchange. Instead, the main routes of exchange between the eastern and western basins appear to be limited to the northern and southern ends of ridge system: Atlantic waters into the Greenland Sea in the Fram St and Artic waters into the southern Norwegian Sea just north of the Iceland-Faroe Ridge.Analysis of a representative isopycnal in the main pycnocline shows it to be stable over time with only small variations with season (except where it outcrops in winter in the Greenland and Iceland Seas). However, two very cold winters, 1968–1969, led to greater than average heat losses across the entire Lofoten Basin that eroded away much of the Lofoten eddy and induced the greatest temperature anomaly in the entire 50-year record. Interannual variations in isopycnal layer temperature correlate with the NAO index such that waters in the Iceland Sea become warmer than average with warming air temperatures and conversely in the Lofoten Basin.  相似文献   

7.
Gravitational Potential Energy （GPE） change due to horizontal/isopycnal eddy diffusion and advection is examined. Horizontal/isopycnal eddy diffusion is conceptually separated into two steps： stirring and sub scale diffusion. GPE changes associated with these two steps are analyzed. In addition, GPE changes due to stirring and subscale diffusion associated with horizontal/isopycnal advection in the Eulerian coordinates are analyzed. These formulae are applied to the SODA data for the world oceans. Our analysis indicates that horizontal/isopycnal advection in Eulerian coordinates can introduce large artificial diffusion in the model. It is shown that GPE source/sink in isopycnal coordinates is closely linked to physical property distribution, such as temperature, salinity and velocity. In comparison with z-coordinates, GPE source/sink due to stir ring/cabbeling associated with isopycnal diffusion/advection is much smaller. Although isopycnal coordi nates may be a better choice in terms of handling lateral diffusion, advection terms in the traditional Eule rian coordinates can produce artificial source of GPE due to cabbeling associated with advection. Reducing such numerical errors remains a grand challenge.  相似文献   

8.
Fifty years of measurements at Ocean Station Papa (OSP, 50°N, 145°W) show trends in the interior waters of the subarctic Pacific that are both impacted by short term (few years to bi-decadal) atmospheric or ocean circulation oscillations and by persistent climate trends. Between 1956 and 2006, waters below the ocean mixed layer to a depth of at least 1000 m have been warming and losing oxygen. On density surfaces found in the depth range 100-400 m (σθ = 26.3-27.0), the ocean is warming at 0.005-0.012 °C y−1, whereas oxygen is declining at 0.39-0.70 μmol kg−1 y−1 or at an integrated rate of 123 mmol m−2 y−1 (decrease of 22% over 50 years). During this time, the hypoxic boundary (defined as 60 μmol O2 kg−1) has shoaled from ∼400 to 300 m. In the Alaska Gyre, the 26.2 isopycnal occasionally ventilates, whereas at OSP 26.0σθ has not been seen at the ocean surface since 1971 as the upper ocean continues to stratify. To interpret the 50 year record at OSP, the isopycnal transport of oxygenated waters within the interior of the subarctic Pacific is assessed by using a slightly modified “NO” parameter [Broecker, W., 1974. “NO” a conservative water-mass tracer. Earth and Planetary Science Letters 23, 100-107]. The highest nitrate-oxygen signature in interior waters of the North Pacific is found in the Bering Sea Gyre, Western Subarctic Gyre and East Kamchatka Current region as a consequence of winter mixing to the ∼26.6 isopycnal. By mixing with low NO waters found in the subtropics and Okhotsk Sea, this signature is diluted as waters flow eastward across the Pacific. Evidence of low NO waters flowing north from California is seen along the coasts of British Columbia and SE Alaska. Oxygen in the subsurface waters of the Alaskan Gyre was supplied ∼60% by subarctic and 40% by subtropical waters during WOCE surveys, whereas such estimates are shown to periodically vary by 20% at OSP. Other features discernable in the OSP data include periods of increased ventilation of deeper isopycnals on an ∼18 year cycle and strong, short term (few month) variability caused by passing mesoscale eddies. The potential impacts of declining oxygen on coastal ecosystems are discussed.  相似文献   

9.
Distributions and characteristics of water mass and chlorofluorocarbons (CFCs) in the North Pacific are investigated by using a General Circulation Model (GCM). The anthropogenic CO2 uptake by the ocean is estimated with velocity fields derived from the GCM experiments. The sensitivity of the uptake to different diffusion parameterizations and different surface forcing used in the GCM is investigated by conducting the three GCM experiments; the diffusive processes are parameterized by horizontal and vertical eddy diffusion which is used in many previous models (RUN1), parameterized by isopycnal diffusion (RUN2), and isopycnal diffusion and perpetual winter forcing for surface temperature and salinity (RUN3). Realistic features for water masses and CFCs can be simulated by the isopycnal diffusion models. The horizontal and vertical diffusion model fails to simulate the salinity minimum and realistic penetration of CFCs into the ocean. The depth of the salinity minimum layer is better simulated under the winter forcing. The results suggest that both isopycnal parameterization and winter forcing are crucial for the model water masses and CFCs simulations. The oceanic uptake of anthropogenic CO2 in RUN3 is about 19.8 GtC in 1990, which is larger by about 10% than that in RUN1 with horizontal and vertical diffusive parameterization. RUN3 well simulates the realistic water mass structure of the intermediate layer considered as a candidate of oceanic sink for anthropogenic CO2. The results suggest that the previous models with horizontal and vertical diffusive parameterization may give the oceanic uptake of anthropogenic CO2 underestimated. This revised version was published online in August 2006 with corrections to the Cover Date.  相似文献   

10.
Data archives from four global coupled ocean–atmosphere models are used to construct maps of diapycnal mass flux through selected isopycnal surfaces in the model oceans. The maps illustrate location and strength of the up and downwelling limbs of thermohaline-forced overturning loops whose stability in the face of rising atmospheric carbon dioxide (CO2) concentrations is of major concern in century-scale climate prediction. The up and downwelling limbs simulated by the four models for present-day greenhouse gas concentrations are compared with observational estimates. Predicted changes in the overturning brought about by gradually rising atmospheric CO2 content are compared model-to-model. While all four models predict some decline in the rate of Atlantic overturning during CO2-induced global warming, the geographic layout of the overturning circulations in each model is found to be insensitive to the changing climate.  相似文献   

11.
《Ocean Modelling》2008,20(3):223-239
A turbulence closure for the effect of mesoscale eddies in non-eddy-resolving ocean models is proposed. The closure consists of a prognostic equation for the eddy kinetic energy (EKE) that is integrated as an additional model equation, and a diagnostic relation for an eddy length scale (L), which is given by the minimum of Rhines scale and Rossby radius. Combining EKE and L using a standard mixing length assumption gives a diffusivity (K), corresponding to the thickness diffusivity in the [Gent, P.R., McWilliams, J.C. 1990. Isopycnal mixing in ocean circulation models. J. Phys. Oceanogr. 20, 150–155] parameterisation. Assuming downgradient mixing of potential vorticity with identical diffusivity shows how K is related to horizontal and vertical mixing processes in the horizontal momentum equation, and also enables us to parameterise the source of EKE related to eddy momentum fluxes.The mesoscale eddy closure is evaluated using synthetic data from two different eddy-resolving models covering the North Atlantic Ocean and the Southern Ocean, respectively. The diagnosis shows that the mixing length assumption together with the definition of eddy length scales is valid within certain limitations. Furthermore, implementation of the closure in non-eddy-resolving models of the North Atlantic and the Southern Ocean shows consistently that the closure has skill at reproducing the results of the eddy-resolving model versions in terms of EKE and K.  相似文献   

12.
Vertical profiles of tritium in seawater were determined for samples collected during the period from 1988 to 1990 at fourteen stations in the northwestern North Pacific (the Oyashio region) including the Okhotsk Sea and the Bering Sea. The profiles usually had a maximum in the surface layer and decreased gradually with depth down to 1,000 m. The water column inventory of tritium averaged 63% of the total atmospheric input in this region.The horizontal distribution of tritium showed a maximum in the region facing the Okhotsk Sea near 45°N for every isopycnal surface of 0 ranging from 26.60 to 27.40. The ages of the intermediate water were calculated for the respective isopycnal surfaces in the maximum region. This calculation assumed that the intermediate water was formed by the isopycnal mixing of two water masses—the Okhotsk Sea and the Bering Sea Component Waters, which had been produced in wintertime by the diapycnal mixing of the surface and the deep waters in the respective marginal seas. The results show that the intermediate water in this region was formed in the late 1980's for the water which has 0 of 26.60 to 26.80 and about 1970 for the water which has 0 of 27.00 to 27.40. Although we have estimated the mean ages of the intermediate water, the horizontal profile of dissolved oxygen suggests that the Okhotsk Sea Component Water is younger than the mean age.  相似文献   

13.
Analysis of wintertime CLIMODE data for 2007 indicates that a substantial portion of new Eighteen Degree Water (EDW) is likely ventilated within the eastward flowing Gulf Stream (GS) between 67°W and 52°W longitudes, possibly exceeding that formed elsewhere in the northern Sargasso Sea. Use of some global air–sea interaction data sets applied to the study region for Feb/Mar of 2007 indicate that this winter may have been anomalously energetic in air–sea exchange compared to the mean of the prior 19 yr. The largest heat and freshwater fluxes found directly over the meandering warm core of the Gulf Stream are capable of removing most of the subtropical heat anomaly of the GS, but cross-frontal fluxes of salinity are required to account for the observed regional salinity structure. An isopycnal diffusivity of ∼100 m2 s−1 is inferred from the salinity balance. This mixing would also account for the observation that EDW formed in the GS is slightly fresher than that formed in northern Sargasso Sea. The lateral flux of heat across the GS north wall also acts to cool the resulting EDW water, but the heat balance for EDW production is largely determined from GS advection and air–sea fluxes, in contrast to salinity. Based on oxygen saturation data, we estimate that 1.8–3.0 Sv-yr of new EDW is formed in the GS for the winter of 2007. EDW originating from the GS is generated in a separate location from where it is accumulated in the northern Sargasso Sea. This manner of EDW formation will produce unique characteristics of EDW found in the northern Sargasso Sea: ones that differ in T/S properties from that formed south of the GS under the more traditional 1D, cooling-driven convection process.  相似文献   

14.
Gravitational potential energy （GPE） source and sink due to stirring and cabbeling associated with sigma dif fusion/ advection is analyzed. It is shown that GPE source and sink is too big, and they are not closely linked to physical property distribution, such as temperature, salinity and velocity. Although the most frequently quoted advantage of sigma coordinate models are their capability of dealing with topography; the exces sive amount of GPE source and sink due to stirring and cabbeling associated with sigma diffusion/advec tion diagnosed from our analysis raises a very serious question whether the way lateral diffusion/advection simulated in the sigma coordinates model is physically acceptable. GPE source and sink in three coordinates is dramatically different in their magnitude and patterns. Overall, in terms of simulating lateral eddy diffu sion and advection isopycnal coordinates is the best choice and sigma coordinates is the worst. The physical reason of the excessive GPE source and sink in sigma coordinates is further explored in details. However, even in the isopycnal coordinates, simulation based on the Eulerian coordinates can be contaminated by the numerical errors associated with the advection terms.  相似文献   

15.
Newly formed North Pacific Tropical Water (NPTW) is carried to the Philippine Sea (PS) by the North Equatorial Current (NEC) as a subsurface salinity maximum. In this study its spreading and salinity change processes are explored using existing hydrographic data of the World Ocean Database 2009 and Argo floats. Spreading of NPTW is closely associated with the transports of the NEC, Mindanao Current (MC), and Kuroshio. Estimated for subsurface water with salinity S greater than 34.8?psu, the southward (northward) geostrophic transport of NPTW by the MC (Kuroshio) at 8°N (18°N) is about 4.4 (5.7)?Sv (1?Sv?=?106?m3?s?1), which is not sensitive to reference level choice. Fields of salinity maximum, geostrophic current, sea level variation, and potential vorticity suggest that the equatorward spreading of NPTW to the tropics is primarily afforded by the MC, whereas its poleward spreading is achieved by both the Kuroshio transport along the coast and open-ocean mesoscale eddy fluxes in the northern PS. The NPTW also undergoes a prominent freshening in the PS. Lying beneath fresh surface water, salinity decreases quicker in the upper part of the NPTW, which gradually lowers the salinity maximum of NPTW to denser isopycnals. Salinity decrease is especially fast in the MC, with along-path decreasing rate reaching O (10?7?psu?s?1). Both diapycnal and isopycnal mixing effects are shown to be elevated in the MC owing to enhanced salinity gradient near the Mindanao Eddy. These results suggest intensive dispersion of thermal anomalies along the subtropical-to-tropical thermocline water pathway near the western boundary.  相似文献   

16.
Ａｑｕａｓｉ－ｔｈｒｅｅ－ｄｉｍｅｎｓｉｏｎａｌｎｕｍｅｒｉｃａｌｐｒｅｄｉｃｔｉｏｎｍｏｄｅｌｏｆｓａｌｉｎｉｔｙｓｔｒｕｃｔｕｒｅｉｎＢｏｈａｉＳｅａａｎｄＨｕａｎｇｈａｉＳｅａ￥ＳｕｎＷｅｉｙａｎｇａｎｄＷａｎｇＺｏｎｇｓｈａｎ（Ｒｅｃｅｉｖｅ...  相似文献   

17.
We have combined the available total CO2, temperature, salinity and oxygen data from the TTO, SAVE and WOCE programs in the Atlantic Ocean to parameterize TCO2 below 500 m depth as a function of potential temperature, salinity and apparent oxygen utilization. We then use the Levitus data set of temperature, salinity and oxygen to compute the TCO2 profiles at the resolution of the Levitus data set on a 1° × 1° grid with a vertical resolution of 33 layers, more densely spaced in the upper 1500 m than below. Depending on the method used to interpolate the data (along isopycnals or vertically by station), the estimated random uncertainty of the computed TC02 values in the Atlantic Ocean throughout the water column below the wintertime mixed layer depth ranges from ± 7.1 μmol kg−1 to t 5.9 μmol kg−1.  相似文献   

18.
Predictions of maximum wave height Hmax are made at Cromer, Happisburgh and Lowestoft on the East Anglian Coast using the formula Hmax = CKRKS U2/g where C is a constant, KR, KS are refraction and shoaling cofficients, U is wind speed and g is the acceleration due to gravity. Comparisonsare made with the models of Darbyshire Draper (1963) and Bretschneider (1958) Using this wave prediction formula, an estimate of the wave climate in the southern North Sea is deduced for the gales of 2–3 January 1976.  相似文献   

19.
Both horizontal and vertical heat exchanges and feedbacks between air temperature and anthropogenic heat fluxes significantly affect the characteristics of the urban heat island (UHI). The UHI intensity depends, in particular, on the ratio between the scales LA (area of anthropogenic forcing) and Lγ (distance passed by an air particle of the oncoming stably stratified flow before its temperature approaches air temperature within the UHI). Both advection and feedback effects may be estimated based on the equation for the local heat balance of the underlying surface. In this case, heat advection is taken into account by calculating temperatures individually for the atmospheric boundary layer and the surface of the urban canopy layer. The estimates show that the asymptotics of strong advection is more characteristic of a typical city. However, under weak winds, with consideration for the feedback between air temperature and anthropogenic heat flux, some deviations from this asymptotics are probable.  相似文献   

20.
Multi-parameter mixing analysis is applied to 26 stations of Auftrieb 75, a joint expedition of R.R.S. Discovery and FS Meteor. North Atlantic Central Water (NACW) and South Atlantic Central Water (SACW) are represented by two water types each, and phosphate, silicate and nitrate are all used as an additional parameter with temperature and salinity, yielding three sets of independent solutions. Solutions from different parameters show similar vertical distributions but can differ considerably in absolute terms. This is believed to be caused by insufficient knowledge of nutrient values of pure NACW and SACW and by low range/precision ratios of nutrient determination.It is shown that isopycnal mixing analysis correctly describes the water mass distribution to lowest order, i.e. there is agreement between results from isopycnal and multi-parameter analysis at most stations within the error bounds which are inherent in the method and which are presently very high. Significant non-isopycnal mixing could be confirmed for only one station, situated on the continental slope in the area of the undercurrent of the Canary Current upwelling regime, and may be caused by active upwelling, increased lateral and vertical shear in the under-current, or internal waves interacting with the bottom. Five stations do not yield oceanographically acceptable solutions, and an attempt is made to include a coastal water mass in the analysis in order to improve the result for these stations, i.e. three-parameter analysis for NACW and SACW is replaced by four-parameter analysis for NACW, SACW and coastal water by combining temperature and salinity with phosphate and silicate or with nitrate and silicate simultaneously (the combination phosphate and nitrate is degenerate because both nutrients are linearly correlated). This indeed improves the results, indicating the presence of coastal water at those stations, but accumulation of errors inhibits quantitative conclusions.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号