ALBERT

Description: Observations and numerical simulations have shown that the meridional overturning circulation (MOC) exhibits substantial variability on sub- to interannual timescales. This variability is not fully understood. In particular it is not known what fraction of the MOC variability is caused by processes such as mesoscale ocean eddies and waves which are ubiquitous in the ocean. Here we analyse twin experiments performed with a global ocean model at eddying (1/4°) and non-eddying (1°) resolutions. The twin experiments are forced with the same surface fluxes for the 1958 to 2001 period but start from different initial conditions. Our results show that on subannual to interannual timescales a large fraction of MOC variability directly reflects variability in the surface forcing. Nevertheless, in the eddy-permitting case there is an initial-condition-dependent MOC variability (hereinafter referred to as "chaotic" variability) of several Sv (1Sv = 106 m3 s−1) in the Atlantic and the Indo-Pacific. In the Atlantic the chaotic MOC variability represents up to 30% of the total variability at the depths where the maximum MOC occurs. In comparison the chaotic MOC variability is only 5–10% in the non-eddying case. The surface forcing being almost identical in the twin experiments suggests that mesoscale ocean eddies are the most likely cause for the increased chaotic MOC variability in the eddying case. The exact formation time of eddies is determined by the initial conditions which are different in the two model passes, and as a consequence the mesoscale eddy field is decorrelated in the twin experiments. In regions where eddy activity is high in the eddy-permitting model, the correlation of sea surface height variability in the twin runs is close to zero. In the non-eddying case in contrast, we find high correlations (0.9 or higher) over most regions. Looking at the sub- and interannual MOC components separately reveals that most of the chaotic MOC variability is found on subannual timescales for the eddy-permitting model. On interannual timescales the amplitude of the chaotic MOC variability is much smaller and the amplitudes are comparable for both the eddy-permitting and non-eddy-permitting model resolutions. Whereas the chaotic MOC variability on interannual timescales only accounts for a small fraction of the total chaotic MOC variability in the eddy-permitting case, it is the main contributor to the chaotic variability in the non-eddying case away from the Equator.

Description: In a sensitivity experiment, an eddy-permitting ocean general circulation model is forced with freshwater fluxes from the Greenland Ice Sheet, averaged for the period 1991–2000. The fluxes are obtained with a mass balance model for the ice sheet, forced with the ERA-40 reanalysis dataset. The freshwater flux is distributed around Greenland as an additional term in prescribed runoff, representing seasonal melting of the ice sheet and a fixed year-round iceberg calving flux, for 8.5 model years. The impacts on regional hydrography and circulation are investigated by comparing the sensitivity experiment to a control experiment, without Greenland fluxes. By the end of the sensitivity experiment, the majority of additional fresh water has accumulated in Baffin Bay, and only a small fraction has reached the interior of the Labrador Sea, where winter mixed layer depth is sensitive to small changes in salinity. As a consequence, the impact on large-scale circulation is very slight. An indirect impact of strong freshening off the west coast of Greenland is a small anti-cyclonic circulation around Greenland which opposes the wind-driven cyclonic circulation and reduces net southward flow through the Canadian Archipelago by ~10%. Implications for the post-2000 acceleration of Greenland mass loss are discussed.

Description: Model studies of the Southern Ocean, reported here, show that the Antarctic Circumpolar Current responds within two days to changes in the zonal wind stress at the latitudes of Drake Passage. Further investigation shows that the response is primarily barotropic and that, as one might expect, it is controlled by topography. Analysis of the results show that the changes in the barotropic flow are sufficient to transfer the changed surface wind stress to the underlying topography and that during this initial phase baroclinic processes are not involved. The model results also show that the Deacon Cell responds to changes in the wind stress on the same rapid time scale. It is shown that the changes in the Deacon Cell can also be explained by the change in the barotropic velocity field, an increase in the zonal wind stress producing an increased northward flow in shallow regions and southward flow where the ocean is deep. This new explanation is unexpected as previously the Deacon Cell has been thought of as a baroclinic feature of the ocean. The results imply that where baroclinic processes do appear to be involved in either the zonal momentum balance of the Southern Ocean or the formation of the Deacon Cell, they are part of the long term baroclinic response of the ocean's density field to the changes in the barotropic flow.

Description: A global general circulation model coupled to a simple six-compartment ecosystem model is used to study the extent to which global variability in primary and export production can be realistically predicted on the basis of advanced parameterizations of upper mixed layer physics, without recourse to introducing extra complexity in model biology. The ''K profile parameterization'' (KPP) scheme employed, combined with 6-hourly external forcing, is able to capture short-term periodic and episodic events such as diurnal cycling and storm-induced deepening. The model realistically reproduces various features of global ecosystem dynamics that have been problematic in previous global modelling studies, using a single generic parameter set. The realistic simulation of deep convection in the North Atlantic, and lack of it in the North Pacific and Southern Oceans, leads to good predictions of chlorophyll and primary production in these contrasting areas. Realistic levels of primary production are predicted in the oligotrophic gyres due to high frequency external forcing of the upper mixed layer (accompanying paper Popova et al., 2006) and novel parameterizations of zooplankton excretion. Good agreement is shown between model and observations at various JFOFS time series sites: BATS, KERFIX, Papa and station India. One exception is that the high zooplankton grazing rates required to maintain low chlorophyll in high-nutrient low-chlorophyll and oligotrophic systems lessened agreement between model and data in the northern North Atlantic, where mesozooplankton with lower grazing rates may be dominant. The model is therefore not globally robust in the sense that additional parameterizations were needed to realistically simulate ecosystem dynamics in the North Atlantic. Nevertheless, the work emphasises the need to pay particular attention to the parameterization of mixed layer physics in global ocean ecosystem modelling as a prerequisite to increasing the complexity of ecosystem models.

Description: Model studies of the Southern Ocean, reported here, show that the Antarctic Circumpolar Current responds within two days to changes in the zonal wind stress at the latitudes of Drake Passage. Further investigation shows that the response is primarily barotropic and that, as one might expect, it is controlled by topography. Analysis of the results show that the changes in the barotropic flow are sufficient to transfer the changed surface wind stress to the underlying topography and that during this initial phase baroclinic processes are not involved. The model results also show that the Deacon Cell responds to changes in the wind stress on the same rapid time scale. It is shown that the changes in the Deacon Cell can also be explained by the change in the barotropic velocity field, an increase in the zonal wind stress producing an increased northward flow in shallow regions and southward flow where the ocean is deep. This new explanation is unexpected as previously the Deacon Cell has been thought of as a baroclinic feature of the ocean. The results imply that where baroclinic processes do appear to be involved in either the zonal momentum balance of the Southern Ocean or the formation of the Deacon Cell, they are part of the long term baroclinic response of the ocean's density field to the changes in the barotropic flow.

Description: The use of 6 h, daily, weekly and monthly atmospheric forcing resulted in dramatically different predictions of plankton productivity in a global 3-D coupled physical-biogeochemical model. Resolving the diurnal cycle of atmospheric variability by use of 6 h forcing, and hence also diurnal variability in UML depth, produced the largest difference, reducing predicted global primary and new production by 25% and 10% respectively relative to that predicted with daily and weekly forcing. This decrease varied regionally, being a 30% reduction in equatorial areas primarily because of increased light limitation resulting from deepening of the mixed layer overnight as well as enhanced storm activity, and 25% at moderate and high latitudes primarily due to increased grazing pressure resulting from late winter stratification events. Mini-blooms of phytoplankton and zooplankton occur in the model during these events, leading to zooplankton populations being sufficiently well developed to suppress the progress of phytoplankton blooms. A 10% increase in primary production was predicted in the peripheries of the oligotrophic gyres due to increased storm-induced nutrient supply end enhanced winter production during the short term stratification events that are resolved in the run forced by 6 h meteorological fields. By resolving the diurnal cycle, model performance was significantly improved with respect to several common problems: underestimated primary production in the oligotrophic gyres; overestimated primary production in the Southern Ocean; overestimated magnitude of the spring bloom in the subarctic Pacific Ocean, and overestimated primary production in equatorial areas. The result of using 6 h forcing on predicted ecosystem dynamics was profound, the effects persisting far beyond the hourly timescale, and having major consequences for predicted global and new production on an annual basis.

Description: A global general circulation model coupled to a simple six-compartment ecosystem model is used to study the extent to which global variability in primary and export production can be realistically predicted on the basis of advanced parameterizations of upper mixed layer physics, without recourse to introducing extra complexity in model biology. The "K profile parameterization" (KPP) scheme employed, combined with 6-hourly external forcing, is able to capture short-term periodic and episodic events such as diurnal cycling and storm-induced deepening. The model realistically reproduces various features of global ecosystem dynamics that have been problematic in previous global modelling studies, using a single generic parameter set. The realistic simulation of deep convection in the North Atlantic, and lack of it in the North Pacific and Southern Oceans, leads to good predictions of chlorophyll and primary production in these contrasting areas. Realistic levels of primary production are predicted in the oligotrophic gyres due to high frequency external forcing of the upper mixed layer (accompanying paper Popova et al., 2006) and novel parameterizations of zooplankton excretion. Good agreement is shown between model and observations at various JGOFS time series sites: BATS, KERFIX, Papa and HOT. One exception is the northern North Atlantic where lower grazing rates are needed, perhaps related to the dominance of mesozooplankton there. The model is therefore not globally robust in the sense that additional parameterizations are needed to realistically simulate ecosystem dynamics in the North Atlantic. Nevertheless, the work emphasises the need to pay particular attention to the parameterization of mixed layer physics in global ocean ecosystem modelling as a prerequisite to increasing the complexity of ecosystem models.

Description: Water mass transformation in the North Atlantic is examined in an eddy-permitting simulation with the OCCAM ocean general circulation model, forced by realistic surface fluxes over the period 1985–2002. Three regions are considered: the Subtropics, the Mid-latitudes and the Northeast Atlantic. The oceanic boundaries of each region coincide with hydrographic sections occupied in recent years. These regions broadly represent the formation sites of Subtropical Mode Water (STMW) and Subpolar Mode Water (SPMW). A water mass budget is obtained for each region and year. Terms in the budget comprise surface-forced transformation rates, boundary exchanges and unsteadiness. Unsteadiness is relatively small, so that regional net water mass transformation is largely balanced by net boundary exchanges. Transformation rates due to mixing are then obtained as the difference between net and surface transformation rates. Transports at the boundaries are compared with recent observations, and reasonable agreement is obtained. For the period 1985-1993, model surface transformation rates are broadly in agreement with equivalent rates computed using the globally-balanced SOC fluxes of heat and freshwater, derived from ship observations. In each Atlantic region, surface transformation rates reach 10-15 Sv, based on both model and SOC fluxes. Higher spatial and temporal resolution in the OCCAM surface fluxes may improve the realism of surface transformation rates in some regions, notably the Labrador Sea. However, the unrealistic location of the North Atlantic Current too far south leads to spurious surface heating east of the Grand Banks in OCCAM. Period-mean transformation rates due to mixing reveal the formation of intermediate waters in each region due to the "consumption" of lighter and denser waters formed by surface fluxes, and are comparable with recent inverse estimates. There is, however, strong interannual-to-decadal variability in the consumption rates. In the subtropics, STMW consumption rates co-vary with STMW formation rates. In the subpolar gyre, anomalies in SPMW consumption rate appear to lag anomalies in the surface formation rate by up to 4 years.

Description: The use of 6 h, daily, weekly and monthly atmospheric forcing resulted in dramatically different predictions of plankton productivity in a global 3-D coupled physical-biogeochemical model. Resolving the diurnal cycle of atmospheric variability by use of 6 h forcing, and hence also diurnal variability in UML depth, produced the largest difference, reducing predicted global primary and new production by 25% and 10% respectively relative to that predicted with daily and weekly forcing. This decrease varied regionally, being a 30% reduction in equatorial areas and 25% at moderate and high latitudes. A 10% increase in the primary production was seen in the peripheries of the oligotrophic gyres. By resolving the diurnal cycle, model performance was significantly improved with respect to several common problems: underestimated primary production in the oligotrophic gyres; overestimated primary production in the Southern Ocean; overestimated magnitude of the spring bloom in the subarctic Pacific Ocean, and overestimated primary production in equatorial areas. The result of using 6 h forcing on predicted ecosystem dynamics was profound, the effects persisting far beyond the hourly timescale, and having major consequences for predicted global and new production on an annual basis.

Description: It is well-established that, at periods shorter than a year, variations in Antarctic Circumpolar Transport are reflected in a barotropic mode, known as the Southern Mode, in which sea level and bottom pressure varies coherently around Antarctica. Here, we use two multidecadal ocean model runs to investigate the behaviour of the Southern Mode at time scales on which density changes become important, leading to a baroclinic component to the adjustment. We find that the concept of a Southern Mode in bottom pressure remains valid, and remains a direct measure of the circumpolar transport, with changes at the northern boundary playing only a small role even on decadal time scales. However, at periods longer than about 5 yr, density changes start to play a role, leading to a surface intensification of the vertical profile of the transport. We also find that barotropic currents on the continental slope account for a significant fraction of the variability, and produce surface intensification in the meridional-integral flow. The role of density variations results in a sea level signal which, although reflecting transport changes at all time scales, has a ratio of sea level to transport which becomes larger at longer time scales. This means that any long-term transport monitoring strategy based on present measurement systems must involve multiplying the observed quantity by a factor which depends on frequency.