scholarly journals Assessment of the sea-ice carbon pump: Insights from a three-dimensional ocean-sea-ice biogeochemical model (NEMO-LIM-PISCES)

Elem Sci Anth ◽  
2016 ◽  
Vol 4 ◽  
Author(s):  
Sébastien Moreau ◽  
Martin Vancoppenolle ◽  
Laurent Bopp ◽  
Oliver Aumont ◽  
Gurvan Madec ◽  
...  

Abstract The role of sea ice in the carbon cycle is minimally represented in current Earth System Models (ESMs). Among potentially important flaws, mentioned by several authors and generally overlooked during ESM design, is the link between sea-ice growth and melt and oceanic dissolved inorganic carbon (DIC) and total alkalinity (TA). Here we investigate whether this link is indeed an important feature of the marine carbon cycle misrepresented in ESMs. We use an ocean general circulation model (NEMO-LIM-PISCES) with sea-ice and marine carbon cycle components, forced by atmospheric reanalyses, adding a first-order representation of DIC and TA storage and release in/from sea ice. Our results suggest that DIC rejection during sea-ice growth releases several hundred Tg C yr−1 to the surface ocean, of which < 2% is exported to depth, leading to a notable but weak redistribution of DIC towards deep polar basins. Active carbon processes (mainly CaCO3 precipitation but also ice-atmosphere CO2 fluxes and net community production) increasing the TA/DIC ratio in sea-ice modified ocean-atmosphere CO2 fluxes by a few Tg C yr−1 in the sea-ice zone, with specific hemispheric effects: DIC content of the Arctic basin decreased but DIC content of the Southern Ocean increased. For the global ocean, DIC content increased by 4 Tg C yr−1 or 2 Pg C after 500 years of model run. The simulated numbers are generally small compared to the present-day global ocean annual CO2 sink (2.6 ± 0.5 Pg C yr−1). However, sea-ice carbon processes seem important at regional scales as they act significantly on DIC redistribution within and outside polar basins. The efficiency of carbon export to depth depends on the representation of surface-subsurface exchanges and their relationship with sea ice, and could differ substantially if a higher resolution or different ocean model were used.

Elem Sci Anth ◽  
2016 ◽  
Vol 4 ◽  
Author(s):  
R. Grimm ◽  
D. Notz ◽  
R.N. Glud ◽  
S. Rysgaard ◽  
K.D. Six

Abstract It has been suggested that geochemical processes related to sea-ice growth and melt might be important for the polar carbon cycle via the so called sea-ice carbon pump (SICP). The SICP affects the air-sea CO2 exchange by influencing the composition of dissolved inorganic carbon (DIC) and total alkalinity (TA) in the surface ocean. Here we quantify the strength of the SICP-induced air-sea CO2 flux using the global three-dimensional ocean-sea-ice-biogeochemical model MPIOM/HAMOCC. Simulations prescribing the range of observed DIC and TA concentrations in the sea ice were performed under two idealized climate scenarios for the present-day and the future oceanic and sea-ice state, both forced with a fixed atmospheric CO2 concentration. Model results indicate that the SICP-induced air-sea CO2 uptake increases with higher ratios of TA:DIC prescribed in the sea ice relative to the basic oceanic TA:DIC ratios. Independent of the modeled scenario, the simulated strength of the SICP is larger in the Antarctic than in the Arctic, because of more efficient export of brine-associated DIC from the Antarctic mixed layer. On an annual basis, we generally find an enhanced SICP-induced oceanic CO2 uptake in regions with net sea-ice melt, and enhanced SICP-induced oceanic CO2 out-gassing in regions with net sea-ice growth. These general regional patterns are modified further by the blockage of air-sea gas exchange through sea-ice coverage. Integrated over the sea-ice zones of both hemispheres, the SICP-induced oceanic CO2 uptake ranges from 2 to 14 Tg C yr−1, which is up to 7% of the simulated net CO2 uptake in polar regions, but far less than 1% of the current global oceanic CO2 uptake. Hence, while we find that the SICP plays a minor role in the modern global carbon cycle, it is of importance for the regional carbon cycle at high latitudes.


2019 ◽  
Vol 12 (8) ◽  
pp. 3745-3758 ◽  
Author(s):  
François Massonnet ◽  
Antoine Barthélemy ◽  
Koffi Worou ◽  
Thierry Fichefet ◽  
Martin Vancoppenolle ◽  
...  

Abstract. The ice thickness distribution (ITD) is one of the core constituents of modern sea ice models. The ITD accounts for the unresolved spatial variability of sea ice thickness within each model grid cell. While there is a general consensus on the added physical realism brought by the ITD, how to discretize it remains an open question. Here, we use the ocean–sea ice general circulation model, Nucleus for European Modelling of the Ocean (NEMO) version 3.6 and Louvain-la-Neuve sea Ice Model (LIM) version 3 (NEMO3.6-LIM3), forced by atmospheric reanalyses to test how the ITD discretization (number of ice thickness categories, positions of the category boundaries) impacts the simulated mean Arctic and Antarctic sea ice states. We find that winter ice volumes in both hemispheres increase with the number of categories and attribute that increase to a net enhancement of basal ice growth rates. The range of simulated mean winter volumes in the various experiments amounts to ∼30 % and ∼10 % of the reference values (run with five categories) in the Arctic and Antarctic, respectively. This suggests that the way the ITD is discretized has a significant influence on the model mean state, all other things being equal. We also find that the existence of a thick category with lower bounds at ∼4 and ∼2 m for the Arctic and Antarctic, respectively, is a prerequisite for allowing the storage of deformed ice and therefore for fostering thermodynamic growth in thinner categories. Our analysis finally suggests that increasing the resolution of the ITD without changing the lower limit of the upper category results in small but not negligible variations of ice volume and extent. Our study proposes for the first time a bi-polar process-based explanation of the origin of mean sea ice state changes when the ITD discretization is modified. The sensitivity experiments conducted in this study, based on one model, emphasize that the choice of category positions, especially of thickest categories, has a primary influence on the simulated mean sea ice states while the number of categories and resolution have only a secondary influence. It is also found that the current default discretization of the NEMO3.6-LIM3 model is sufficient for large-scale present-day climate applications. In all cases, the role of the ITD discretization on the simulated mean sea ice state has to be appreciated relative to other influences (parameter uncertainty, forcing uncertainty, internal climate variability).


2019 ◽  
Author(s):  
François Massonnet ◽  
Antoine Barthélemy ◽  
Koffi Worou ◽  
Thierry Fichefet ◽  
Martin Vancoppenolle ◽  
...  

Abstract. The Ice Thickness Distribution (ITD) is one of the core constituents of modern sea ice models. The ITD accounts for the unresolved spatial variability of sea ice thickness within each model grid cell. While there is a general consensus on the added physical realism brought by the ITD, how to implement it remains an open question. Here, we use the ocean--sea ice general circulation model NEMO3.6-LIM3 forced by atmospheric reanalyses to test how the ITD discretization (number of ice thickness categories, positions of the category boundaries) impacts the simulated mean Arctic and Antarctic sea ice states. We find that winter ice volumes in both hemispheres increase with the number of categories, and attribute that increase to a net enhancement of basal ice growth rates. The range of simulated mean winter volumes in the various experiments amounts to ~ 30 % and ~ 10 % of the reference values (run with 5 categories) in the Arctic and Antarctic, respectively. This suggests that the way the ITD is discretized has a significant influence on the model mean state, all other things being equal. We also find that the existence of a thick category with lower bounds at ~ 4 m and ~ 2 m for the Arctic and Antarctic, respectively, is a prerequisite for allowing the storage of deformed ice, and therefore for fostering thermodynamic growth in thinner categories. Our analysis finally suggests that increasing the resolution of the ITD without changing the lower limit of the upper category results in small but not negligible variations of ice volume and extent. Our study proposes for the first time a bi-polar process-based explanation of the origin of mean state changes when the ITD discretization is modified. The sensitivity experiments conducted in this study, based on one model, emphasize that the choice of category positions, especially of thickest categories, has a primary influence on the simulated mean sea ice states while the number of categories and resolution have only a secondary influence.


2020 ◽  
Author(s):  
Miguel Angel Morales Maqueda ◽  
Alethea Sara Mountford

<p>The presence of microplastics in the Arctic sea ice cover and water column, as well as on land, has raised the already high concerns about the dispersion of litter in the global environment. We present a 50-year simulation carried out with the NEMO ocean general circulation model of the dispersion of buoyant and neutrally buoyant microplastics in the global ocean that includes a simple formulation of microplastic accumulation in, and advection by, sea ice. Microplastics enter the Arctic predominantly through the Barents Sea, with a smaller input through the Bering Strait, although the simulation also takes into account small plastic sources along the Arctic coastline. Microplastics become trapped in newly formed sea ice chiefly on the Eurasian shelves and the Chukchi Sea, but a still significant amount is transferred from the mixed layer to the ice base through congelation in the central Arctic, where microplastics congregate nearer to the surface than elsewhere in the global ocean due to the strong stratification and the relatively small levels of vertical turbulence underneath multiyear sea ice. In the model, the maximum average residence time of sea ice in the Arctic is about six years, and this is also, therefore, the typical timescale for maximum microplastic accumulation in the ice cover. Plastics trapped in sea ice undergo a seasonal cycle of accumulation and release in consonance with the freeze and melt sea ice cycle but ultimately are release back into the ocean in the Greenland and Labrador seas, from where they will be subsequently transported into the North Atlantic.</p>


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.


Ocean Science ◽  
2017 ◽  
Vol 13 (1) ◽  
pp. 61-75 ◽  
Author(s):  
Arash Bigdeli ◽  
Brice Loose ◽  
An T. Nguyen ◽  
Sylvia T. Cole

Abstract. In ice-covered regions it is challenging to determine constituent budgets – for heat and momentum, but also for biologically and climatically active gases like carbon dioxide and methane. The harsh environment and relative data scarcity make it difficult to characterize even the physical properties of the ocean surface. Here, we sought to evaluate if numerical model output helps us to better estimate the physical forcing that drives the air–sea gas exchange rate (k) in sea ice zones. We used the budget of radioactive 222Rn in the mixed layer to illustrate the effect that sea ice forcing has on gas budgets and air–sea gas exchange. Appropriate constraint of the 222Rn budget requires estimates of sea ice velocity, concentration, mixed-layer depth, and water velocities, as well as their evolution in time and space along the Lagrangian drift track of a mixed-layer water parcel. We used 36, 9 and 2 km horizontal resolution of regional Massachusetts Institute of Technology general circulation model (MITgcm) configuration with fine vertical spacing to evaluate the capability of the model to reproduce these parameters. We then compared the model results to existing field data including satellite, moorings and ice-tethered profilers. We found that mode sea ice coverage agrees with satellite-derived observation 88 to 98 % of the time when averaged over the Beaufort Gyre, and model sea ice speeds have 82 % correlation with observations. The model demonstrated the capacity to capture the broad trends in the mixed layer, although with a significant bias. Model water velocities showed only 29 % correlation with point-wise in situ data. This correlation remained low in all three model resolution simulations and we argued that is largely due to the quality of the input atmospheric forcing. Overall, we found that even the coarse-resolution model can make a modest contribution to gas exchange parameterization, by resolving the time variation of parameters that drive the 222Rn budget, including rate of mixed-layer change and sea ice forcings.


2017 ◽  
Vol 30 (11) ◽  
pp. 3945-3962 ◽  
Author(s):  
James A. Screen

Abstract The loss of Arctic sea ice is already having profound environmental, societal, and ecological impacts locally. A highly uncertain area of scientific research, however, is whether such Arctic change has a tangible effect on weather and climate at lower latitudes. There is emerging evidence that the geographical location of sea ice loss is critically important in determining the large-scale atmospheric circulation response and associated midlatitude impacts. However, such regional dependencies have not been explored in a thorough and systematic manner. To make progress on this issue, this study analyzes ensemble simulations with an atmospheric general circulation model prescribed with sea ice loss separately in nine regions of the Arctic, to elucidate the distinct responses to regional sea ice loss. The results suggest that in some regions, sea ice loss triggers large-scale dynamical responses, whereas in other regions sea ice loss induces only local thermodynamical changes. Sea ice loss in the Barents–Kara Seas is unique in driving a weakening of the stratospheric polar vortex, followed in time by a tropospheric circulation response that resembles the North Atlantic Oscillation. For October–March, the largest spatial-scale responses are driven by sea ice loss in the Barents–Kara Seas and the Sea of Okhotsk; however, different regions assume greater importance in other seasons. The atmosphere responds very differently to regional sea ice losses than to pan-Arctic sea ice loss, and the response to pan-Arctic sea ice loss cannot be obtained by the linear addition of the responses to regional sea ice losses. The results imply that diversity in past studies of the simulated response to Arctic sea ice loss can be partly explained by the different spatial patterns of sea ice loss imposed.


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.


2016 ◽  
Vol 12 (9) ◽  
pp. 1907-1918 ◽  
Author(s):  
Rasmus A. Pedersen ◽  
Peter L. Langen ◽  
Bo M. Vinther

Abstract. Insolation changes during the Eemian (the last interglacial period, 129 000–116 000 years before present) resulted in warmer than present conditions in the Arctic region. The NEEM ice core record suggests warming of 8 ± 4 K in northwestern Greenland based on stable water isotopes. Here we use general circulation model experiments to investigate the causes of the Eemian warming in Greenland. Simulations of the atmospheric response to combinations of Eemian insolation and preindustrial oceanic conditions and vice versa are used to disentangle the impacts of the insolation change and the related changes in sea surface temperatures and sea ice conditions. The changed oceanic conditions cause warming throughout the year, prolonging the impact of the summertime insolation increase. Consequently, the oceanic conditions cause an annual mean warming of 2 K at the NEEM site, whereas the insolation alone causes an insignificant change. Taking the precipitation changes into account, however, the insolation and oceanic changes cause more comparable increases in the precipitation-weighted temperature, implying that both contributions are important for the ice core record at the NEEM site. The simulated Eemian precipitation-weighted warming of 2.4 K at the NEEM site is low compared to the ice core reconstruction, partially due to missing feedbacks related to ice sheet changes and an extensive sea ice cover. Surface mass balance calculations with an energy balance model further indicate that the combination of temperature and precipitation anomalies leads to potential mass loss in the north and southwestern parts of the ice sheet. The oceanic conditions favor increased accumulation in the southeast, while the insolation appears to be the dominant cause of the expected ice sheet reduction. Consequently, the Eemian is not a suitable analogue for future ice sheet changes.


Sign in / Sign up

Export Citation Format

Share Document