Abstract
An anisotropic generalization of the Gent–McWilliams (GM) parameterization is presented for eddy-induced tracer transport and diffusion in ocean models, and it is implemented in an ocean ...general circulation model using a functional formalism to derive the spatial discretization. This complements the anisotropic viscosity parameterization recently developed by Smith and McWilliams. The anisotropic GM operator is potentially useful in both coarse- and high-resolution ocean models, and in this study the focus is on its application in high-resolution eddying solutions, for which it provides an adiabatic alternative to the more commonly used biharmonic horizontal diffusion operators. It is shown that realistically high levels of eddy energy can be simulated using harmonic anisotropic diffusion and friction operators. Isotropic forms can also be used, but these tend either to overly damp the solution when a large diffusion coefficient is used or to introduce unacceptable levels of numerical noise when a small coefficient is used. A series of numerical simulations of the North Atlantic Ocean are conducted at 0.2° resolution using anisotropic viscosity, anisotropic GM, and biharmonic mixing operators to investigate the effects of the anisotropic forms and to isolate changes in the solutions specifically associated with anisotropic GM. A high-resolution 0.1° simulation is then conducted using both anisotropic forms, and the results are compared with a similar run using biharmonic mixing. Modest improvements are seen in the mean wind-driven circulation with the anisotropic forms, but the largest effects are due to the anisotropic GM parameterization, which eliminates the spurious diapycnal diffusion inherent in horizontal tracer diffusion. This leads to significant improvements in the model thermohaline circulation, including the meridional heat transport, meridional overturning circulation, and deep-water formation and convection in the Labrador Sea.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Ocean models routinely used in simulations of the Earth's climate do not resolve mesoscale eddies because of the immense computational cost. A new parameterization of the effects of these eddies has ...been implemented in a widely used model. A comparison of its solution with that of the conventional parameterization shows significant improvements in the global temperature distribution, the poleward and surface heat fluxes, and the locations of deep-water formation.
A nonlocal K-profile parameterization (KPP) of the upper-ocean boundary layer is tested for the equatorial regions. First, the short-term performance of a one-dimensional model with KPP is found to ...compare favorably to large eddy simulations (LES), including nonlocal countergradient heat flux. The comparison is clean because both the surface forcing and the large-scale flow are identical in the two models. The comparison is direct because the parameterized turbulent flux profiles are explicitly computed in LES. A similar comparison is less favorable when KPP is replaced by purely downgradient diffusion with Richardson-number-dependent viscosity and diffusivity because of the absence of intense convection after sunset. Sensitivity experiments are used to establish parameter values in the interior mixing of KPP. Second, the impact of the parameterization on annual means and the seasonal cycle in a general circulation model of the upper, equatorial Pacific Ocean is described. The results of GCM runs with and without KPP are compared to annual mean profiles of zonal velocity and temperature from the TOGA-TAO array. The two GCM solutions are closer to each other than to the observations, with biases in zonal velocity in the western Pacific and in subsurface temperature in the eastern Pacific. Such comparisons are never clean because neither the wind stress and the surface heat flux nor the forcing by the large-scale flow are known to sufficient accuracy. Finally, comparisons are made of the equatorial Pacific Ocean GCM results when different heat flux formulations are used. These include bulk forcing where prescribed air temperature and humidity are used, SST forcing where the use of such ocean-controlled parameters is avoided, and a fully coupled atmospheric general circulation model where there is no prescribed control over any surface fluxes. It is concluded, especially in the eastern Pacific, that the use of specified air temperature and humidity does not overly constrain the model sea surface temperature.
The North Atlantic Ocean thermohaline circulation is examined in three experiments using the Climate System Model. They are a control integration for 1870 conditions, and particular emission ...scenarios for the 20th and 21st centuries. It is found that the strength of the thermohaline circulation does not change significantly over the 21st century. This is in contrast to several other recent studies, which have projected a significant reduction over the 21st century. The reason for the difference is that the Northwest Atlantic becomes warmer and more saline in the Climate System Model. These changes combine to make little change to the surface ocean density in this region, and hence to the rate of deep water formation. Caveats about the Climate System Model and other coupled climate models are then discussed.
The Community Climate System Model version 3 is used to analyse changes in water mass subduction rates in the South Atlantic Ocean over the 21st century. The model results are first compared to ...observations over 1950–2000, and shown to be rather good. The subduction rates do not change significantly over the 21st century, but the densities at which water masses form become significantly lighter. The strong westerly winds in this region do not change much, which suggests small changes to the rate at which the Atlantic sector of the Southern Ocean takes up heat and carbon dioxide over the 21st century.
An ensemble of nine simulations for the climate of the twentieth century has been run using the Community Climate System Model version 3 (CCSM3). Three of these runs also simulate the uptake of ...chlorofluorocarbon-11 (CFC-11) into the ocean using the protocol from the Ocean Carbon Model Intercomparison Project (OCMIP). Comparison with ocean observations taken between 1980 and 2000 shows that the global CFC-11 uptake is simulated very well. However, there are regional biases, and these are used to identify where too much deep-water formation is occurring in the CCSM3. The differences between the three runs simulating CFC-11 uptake are also briefly documented.
The variability in ocean heat content in the 1870 control runs is shown to be only a little smaller than estimates using ocean observations. The ocean heat uptake between 1957 and 1996 in the ensemble is compared to the recent observational estimates of the secular trend. The trend in ocean heat uptake is considerably larger than the natural variability in the 1870 control runs. The heat uptake down to 300 m between 1957 and 1996 varies by a factor of 2 across the ensemble. Some possible reasons for this large spread are discussed. There is much less spread in the heat uptake down to 3 km. On average, the CCSM3 twentieth-century ensemble runs take up 25% more heat than the recent estimate from ocean observations. Possible explanations for this are that the model heat uptake is calculated over the whole ocean, and not just in the regions where there are many observations and that there is no parameterization of the indirect effects of aerosols in CCSM3.
Celotno besedilo
Dostopno za:
BFBNIB, DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Twelve experiments with two coarse resolutions of a global ocean model using a variety of surface forcings are analyzed to address the question of what sets the mean transport through Drake Passage. ...Seven of the experiments do not have an active sea‐ice model, but the remaining five do. Previous theories have suggested that the Drake Passage transport is governed by the Cape Horn Sverdrup transport or, alternatively, is proportional to the square root of the meridional Ekman transport at the latitude of Drake Passage. The results presented here do not support either of these theories. The Drake Passage transport depends quite strongly on the isopycnal diffusivity parameter in the model and less strongly on the background vertical diffusivity and horizontal viscosity parameters. However, when the magnitudes of these parameters are fixed, the results show a very strong correlation between Drake Passage transport and both the strength of the meridional Ekman transport at the latitude of Drake Passage and the thermohaline circulation off the Antarctic shelf. The relationships are monotonic, but not linear. The best estimate is that the meridional Ekman transport drives ∼100 Sv of Drake Passage transport, while the remaining 30 Sv are driven by the global thermohaline circulation.
Ocean heat uptake and the thermohaline circulation are analyzed in present-day control, 1% increasing CO₂, and doubled CO₂ runs of the Community Climate System Model, version 2 (CCSM2). It is ...concluded that the observed 40-yr trend in the global heat content to 300 m, found by Levitus et al., is somewhat larger than the natural variability in the CCSM2 control run. The observed 40-yr trend in the global heat content down to a depth of 3 km is much closer to trends found in the control run and is not so clearly separated from the natural model variability. It is estimated that, in a 0.7% increasing CO₂ scenario that approximates the effect of increasing greenhouse gases between 1958 and 1998, the CCSM2 40-yr trend in the global heat content to 300 m is about the same as the observed value. This gives support for the CCSM2 climate sensitivity, which is 2.2°C.
Both the maximum of the meridional overturning streamfunction and the vertical flow across 1-km depth between 60° and 65°N decrease monotonically during the 1% CO₂ run. However, the reductions are quite modest, being 3 and 2 Sv, respectively, when CO₂ has quadrupled. The reason for this is that the surface potential density in the northern North Atlantic decreases steadily throughout the 1% CO₂ run. In the latter part of the doubled CO₂ run, the meridional overturning streamfunction recovers in strength back toward its value in the control run, but the deep-water formation rate across 1-km depth between 60° and 65°N remains at 85% of the control run value. The maximum northward heat transport at 22°N is governed by the maximum of the overturning, but the transport poleward of 62°N appears to be independent of the deep-water formation rate.
Celotno besedilo
Dostopno za:
BFBNIB, DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK