scholarly journals A near-global eddy-resolving OGCM for climate studies

Author(s):  
X. Zhang ◽  
P. R. Oke ◽  
M. Feng ◽  
M. A. Chamberlain ◽  
J. A. Church ◽  
...  

Abstract. Eddy-resolving global ocean models are highly desired for spatially-improved climate studies, but this is challenging because they require careful configuration and substantial computational resources. Model drift, partially related to insufficient model spin-up, imperfect model physics or bias in surface forcing, can be problematic, leading to contamination of climate change signals. In this study, we adapt a near-global eddy-resolving ocean general circulation model, originally developed for short-range ocean forecasting, for climate studies. The Ocean Forecasting Australia Model version 3 (OFAM3) is spun up for 20 years, with repeated year 1979 forcing and adaptive relaxation (Newtonian nudging) of temperature and salinity in the deep ocean to an observation-based climatology. In addition, surface heat fluxes from the JRA-55 atmospheric reanalysis are adjusted during the spin-up experiment to minimise excessive net heat uptake in the ocean. In the historical experiment, spanning 1979–2014, a non-adaptive relaxation is applied by repeating the same relaxation rates derived from the last five years of the spin-up experiment, and the surface heat flux adjustment diagnosed during the spinup experiment is also maintained. We demonstrate that the historical experiment driven by the JRA-55 reanalysis does not have significant drifts (e.g., as shown by simulated global ocean heat content), and also provides an eddy-resolving simulation of the global ocean circulation over the period 1979–2014. Decadal changes, such as the strengthening of the subtropical gyre circulation, are also reasonably simulated. A biogeochemical model is coupled with OFAM3 to produce patterns of primary productivity and carbon fluxes that are consistent with observations. Experiences gained from our numerical experiments will be helpful to other modelling groups who are interested in running global eddy-resolving OGCMs for climate studies.

1997 ◽  
Vol 25 ◽  
pp. 111-115 ◽  
Author(s):  
Achim Stössel

This paper investigates the long-term impact of sea ice on global climate using a global sea-ice–ocean general circulation model (OGCM). The sea-ice component involves state-of-the-art dynamics; the ocean component consists of a 3.5° × 3.5° × 11 layer primitive-equation model. Depending on the physical description of sea ice, significant changes are detected in the convective activity, in the hydrographic properties and in the thermohaline circulation of the ocean model. Most of these changes originate in the Southern Ocean, emphasizing the crucial role of sea ice in this marginally stably stratified region of the world's oceans. Specifically, if the effect of brine release is neglected, the deep layers of the Southern Ocean warm up considerably; this is associated with a weakening of the Southern Hemisphere overturning cell. The removal of the commonly used “salinity enhancement” leads to a similar effect. The deep-ocean salinity is almost unaffected in both experiments. Introducing explicit new-ice thickness growth in partially ice-covered gridcells leads to a substantial increase in convective activity, especially in the Southern Ocean, with a concomitant significant cooling and salinification of the deep ocean. Possible mechanisms for the resulting interactions between sea-ice processes and deep-ocean characteristics are suggested.


2012 ◽  
Vol 8 (5) ◽  
pp. 1581-1598 ◽  
Author(s):  
V. Mariotti ◽  
L. Bopp ◽  
A. Tagliabue ◽  
M. Kageyama ◽  
D. Swingedouw

Abstract. Marine sediments records suggest large changes in marine productivity during glacial periods, with abrupt variations especially during the Heinrich events. Here, we study the response of marine biogeochemistry to such an event by using a biogeochemical model of the global ocean (PISCES) coupled to an ocean-atmosphere general circulation model (IPSL-CM4). We conduct a 400-yr-long transient simulation under glacial climate conditions with a freshwater forcing of 0.1 Sv applied to the North Atlantic to mimic a Heinrich event, alongside a glacial control simulation. To evaluate our numerical results, we have compiled the available marine productivity records covering Heinrich events. We find that simulated primary productivity and organic carbon export decrease globally (by 16% for both) during a Heinrich event, albeit with large regional variations. In our experiments, the North Atlantic displays a significant decrease, whereas the Southern Ocean shows an increase, in agreement with paleo-productivity reconstructions. In the Equatorial Pacific, the model simulates an increase in organic matter export production but decreased biogenic silica export. This antagonistic behaviour results from changes in relative uptake of carbon and silicic acid by diatoms. Reasonable agreement between model and data for the large-scale response to Heinrich events gives confidence in models used to predict future centennial changes in marine production. In addition, our model allows us to investigate the mechanisms behind the observed changes in the response to Heinrich events.


1997 ◽  
Vol 25 ◽  
pp. 116-120 ◽  
Author(s):  
S. Legutke ◽  
E. Maier-Reimkr ◽  
A. Stössel ◽  
A. Hellbach

A global ocean general circulation model has been coupled with a dynamic thermodynamic sea-ice model. This model has been spun-up in a 1000 year integration using daily atmosphere model data. Main water masses and currents are reproduced as well as the seasonal characteristics of the ice cover of the Northern and Southern Hemispheres. Model results for the Southern Ocean, however, show the ice cover as too thin, and there are large permanent polynyas in the Weddell and Ross Seas. These polynyas are due to a large upward oceanic heat flux caused by haline rejection during the freezing of sea ice. Sensitivity studies were performed to test several ways of treating the sea-surface salinity and the rejected brine. The impact on the ice cover, water-mass characteristics, and ocean circulation are described.


Ocean Science ◽  
2009 ◽  
Vol 5 (4) ◽  
pp. 461-473 ◽  
Author(s):  
P. Oddo ◽  
M. Adani ◽  
N. Pinardi ◽  
C. Fratianni ◽  
M. Tonani ◽  
...  

Abstract. A new numerical general circulation ocean model for the Mediterranean Sea has been implemented nested within an Atlantic general circulation model within the framework of the Marine Environment and Security for the European Area project (MERSEA, Desaubies, 2006). A 4-year twin experiment was carried out from January 2004 to December 2007 with two different models to evaluate the impact on the Mediterranean Sea circulation of open lateral boundary conditions in the Atlantic Ocean. One model considers a closed lateral boundary in a large Atlantic box and the other is nested in the same box in a global ocean circulation model. Impact was observed comparing the two simulations with independent observations: ARGO for temperature and salinity profiles and tide gauges and along-track satellite observations for the sea surface height. The improvement in the nested Atlantic-Mediterranean model with respect to the closed one is particularly evident in the salinity characteristics of the Modified Atlantic Water and in the Mediterranean sea level seasonal variability.


2005 ◽  
Vol 35 (6) ◽  
pp. 933-948 ◽  
Author(s):  
Henk A. Dijkstra ◽  
Wilbert Weijer

Abstract A study of the stability of the global ocean circulation is performed within a coarse-resolution general circulation model. Using techniques of numerical bifurcation theory, steady states of the global ocean circulation are explicitly calculated as parameters are varied. Under a freshwater flux forcing that is diagnosed from a reference circulation with Levitus surface salinity fields, the global ocean circulation has no multiple equilibria. It is shown how this unique-state regime transforms into a regime with multiple equilibria as the pattern of the freshwater flux is changed in the northern North Atlantic Ocean. In the multiple-equilibria regime, there are two branches of stable steady solutions: one with a strong northern overturning in the Atlantic and one with hardly any northern overturning. Along the unstable branch that connects both stable solution branches (here for the first time computed for a global ocean model), the strength of the southern sinking in the South Atlantic changes substantially. The existence of the multiple-equilibria regime critically depends on the spatial pattern of the freshwater flux field and explains the hysteresis behavior as found in many previous modeling studies.


2016 ◽  
Vol 13 (15) ◽  
pp. 4359-4377 ◽  
Author(s):  
Zuo Xue ◽  
Ruoying He ◽  
Katja Fennel ◽  
Wei-Jun Cai ◽  
Steven Lohrenz ◽  
...  

Abstract. A three-dimensional coupled physical–biogeochemical model was used to simulate and examine temporal and spatial variability of sea surface pCO2 in the Gulf of Mexico (GoM). The model was driven by realistic atmospheric forcing, open boundary conditions from a data-assimilative global ocean circulation model, and observed freshwater and terrestrial nutrient and carbon input from major rivers. A 7-year model hindcast (2004–2010) was performed and validated against ship measurements. Model results revealed clear seasonality in surface pCO2 and were used to estimate carbon budgets in the Gulf. Based on the average of model simulations, the GoM was a net CO2 sink with a flux of 1.11 ± 0.84  ×  1012 mol C yr−1, which, together with the enormous fluvial inorganic carbon input, was comparable to the inorganic carbon export through the Loop Current. Two model sensitivity experiments were performed: one without biological sources and sinks and the other using river input from the 1904–1910 period as simulated by the Dynamic Land Ecosystem Model (DLEM). It was found that biological uptake was the primary driver making GoM an overall CO2 sink and that the carbon flux in the northern GoM was very susceptible to changes in river forcing. Large uncertainties in model simulations warrant further process-based investigations.


2008 ◽  
Vol 21 (22) ◽  
pp. 6015-6035 ◽  
Author(s):  
James A. Carton ◽  
Anthony Santorelli

Abstract This paper examines nine analyses of global ocean 0-/700-m temperature and heat content during the 43-yr period of warming, 1960–2002. Among the analyses are two that are independent of any numerical model, six that rely on sequential data assimilation, including an ocean general circulation model, and one that uses four-dimensional variational data assimilation (4DVAR), including an ocean general circulation model and its adjoint. Most analyses show gradual warming of the global ocean with an ensemble trend of 0.77 × 108 J m−2 (10 yr)−1 (=0.24 W m−2) as the result of rapid warming in the early 1970s and again beginning around 1990. One proposed explanation for these variations is the effect of volcanic eruptions in 1963 and 1982. Examination of this hypothesis suggests that while there is an oceanic signal, it is insufficient to explain the observed heat content variations. A second potential cause of decadal variations in global heat content is the uncorrelated contribution of heat content variations in individual ocean basins. The subtropical North Atlantic is warming at twice the global average, with accelerated warming in the 1960s and again beginning in the late 1980s and extending through the end of the record. The Barents Sea region of the Arctic Ocean and the Gulf of Mexico have also warmed, while the western subpolar North Atlantic has cooled. Heat content variability in the North Pacific differs significantly from the North Atlantic. There the spatial and temporal patterns are consistent with the decadal variability previously identified through observational and modeling studies examining SST and surface winds. In the Southern Hemisphere large heat content anomalies are evident, and while there is substantial disagreement among analyses on average the band of latitudes at 30°–60°S contribute significantly to the global warming trend. Thus, the uncorrelated contributions of heat content variations in the individual basins are a major contributor to global heat content variations. A third potential contributor to global heat content variations is the effect of time-dependent bias in the set of historical observations. This last possibility is examined by comparing the analyses to the unbiased salinity–temperature–depth dataset and finding a very substantial warm bias in all analyses in the 1970s relative to the latter decades. This warm bias may well explain the rapid increase in analysis heat content in the early 1970s, but not the more recent increase, which began in the early 1990s. Finally, this study provides information about the similarities and differences between analyses that are independent of a model and those that use sequential assimilation and 4DVAR. The comparisons provide considerable encouragement for the use of the sequential analyses for climate research despite the presence of erroneous variability (also present in the no-model analyses) resulting from instrument bias. The strengths and weaknesses of each analysis need to be considered for a given application.


2013 ◽  
Vol 10 (11) ◽  
pp. 7219-7234 ◽  
Author(s):  
Z. Xue ◽  
R. He ◽  
K. Fennel ◽  
W.-J. Cai ◽  
S. Lohrenz ◽  
...  

Abstract. A three-dimensional coupled physical-biogeochemical model is applied to simulate and examine temporal and spatial variability of circulation and biogeochemical cycling in the Gulf of Mexico (GoM). The model is driven by realistic atmospheric forcing, open boundary conditions from a data assimilative global ocean circulation model, and observed freshwater and terrestrial nitrogen input from major rivers. A 7 yr model hindcast (2004–2010) was performed, and validated against satellite observed sea surface height, surface chlorophyll, and in situ observations including coastal sea level, ocean temperature, salinity, and dissolved inorganic nitrogen (DIN) concentration. The model hindcast revealed clear seasonality in DIN, phytoplankton and zooplankton distributions in the GoM. An empirical orthogonal function analysis indicated a phase-locked pattern among DIN, phytoplankton and zooplankton concentrations. The GoM shelf nitrogen budget was also quantified, revealing that on an annual basis the DIN input is largely balanced by the removal through denitrification (an equivalent of ~ 80% of DIN input) and offshore exports to the deep ocean (an equivalent of ~ 17% of DIN input).


2013 ◽  
Vol 10 (5) ◽  
pp. 7785-7830 ◽  
Author(s):  
Z. Xue ◽  
R. He ◽  
K. Fennel ◽  
W.-J. Cai ◽  
S. Lohrenz ◽  
...  

Abstract. A three-dimensional coupled physical-biogeochemical model is applied to simulate and examine temporal and spatial variability of circulation and biogeochemical cycling in the Gulf of Mexico (GoM). The model is driven by realistic atmospheric forcing, open boundary conditions from a data assimilative global ocean circulation model, and observed freshwater and terrestrial nutrient input from major rivers. A 7 yr model hindcast (2004–2010) was performed, and validated against satellite observed sea surface height, surface chlorophyll, and in-situ observations including coastal sea-level, ocean temperature, salinity, and nutrient concentration. The model hindcast revealed clear seasonality in nutrient, phytoplankton and zooplankton distributions in the GoM. An Empirical Orthogonal Function analysis indicated a phase-locked pattern among nutrient, phytoplankton and zooplankton concentrations. The GoM shelf nutrient budget was also quantified, revealing that on an annual basis ~80% of nutrient input was denitrified on the shelf and ~17% was exported to the deep ocean.


2012 ◽  
Vol 9 (1) ◽  
pp. 25-61
Author(s):  
A. M. Huerta-Casas ◽  
D. J. Webb

Abstract. The transport and storage of heat by the ocean is of crucial importance because of its effect on ocean dynamics and its impact on the atmosphere, climate and climate change. Unfortunately, limits to the amount of data that can be collected and stored means that many experimental and modelling studies of the heat budget have to make use of mean datasets where the effects of short term fluctuations are lost. In this paper we investigate the magnitude of the resulting errors making use of data from OCCAM, a high resolution global ocean model. The model carries out a proper heat balance every timestep so any imbalances that are found in the analysis must result from the use of mean fields. The study concentrates on two areas of the ocean affecting the El Nino. The first is the region of tropical instability waves north of the equator. The second is in the upwelling region along the equator. It is shown that in both cases, processes with a period of less than five days can have a significant impact on the heat budget. Thus analyses using data averaged over five days or more are likely to have significant errors. It is also shown that if a series of instantaneous values is available, reasonable estimates can be made of the size of the errors. In model studies such values are available in the form of the datasets used to restart the model. In experimental studies they may be in the form of individual unaveraged observations.


Sign in / Sign up

Export Citation Format

Share Document