ALBERT

Description: The total alkalinity–pH equation, which relates total alkalinity and pH for a given set of total concentrations of the acid–base systems that contribute to total alkalinity in a given water sample, is reviewed and its mathematical properties established. We prove that the equation function is strictly monotone and always has exactly one positive root. Different commonly used approximations are discussed and compared. An original method to derive appropriate initial values for the iterative solution of the cubic polynomial equation based upon carbonate-borate-alkalinity is presented. We then review different methods that have been used to solve the total alkalinity–pH equation, with a main focus on biogeochemical models. The shortcomings and limitations of these methods are made out and discussed. We then present two variants of a new, robust and universally convergent algorithm to solve the total alkalinity–pH equation. This algorithm does not require any a priori knowledge of the solution. SolveSAPHE (Solver Suite for Alkalinity-PH Equations) provides reference implementations of several variants of the new algorithm in Fortran 90, together with new implementations of other, previously published solvers. The new iterative procedure is shown to converge from any starting value to the physical solution. The extra computational cost for the convergence security is only 10–15% compared to the fastest algorithm in our test series.

Description: During glacial-interglacial cycles, atmospheric CO2 concentration varied by about 100 ppmv in amplitude. While testing mechanisms that have led to the low glacial CO2 level could be done in equilibrium model experiments, an ultimate goal is to explain CO2 changes in transient simulations through the complete glacial-interglacial cycle. The computationally efficient Earth System model of intermediate complexity CLIMBER-2 is used to simulate global biogeochemistry over the last glacial cycle (126 kyr). The physical core of the model (atmosphere, ocean, land and ice sheets) is driven by orbital changes and reconstructed radiative forcing from greenhouses gases, ice, and aeolian dust. The carbon cycle model is able to reproduce the main features of the CO2 changes: a 50 ppmv CO2 drop during glacial inception, a minimum concentration at the last glacial maximum 80 ppmv lower than the Holocene value, and an abrupt 60 ppmv CO2 rise during the deglaciation. The model deep ocean δ13C also resembles reconstructions from deep-sea cores. The main drivers of atmospheric CO2 evolve in time: changes in sea surface temperatures and in the volume of bottom water of southern origin control atmospheric CO2 during the glacial inception and deglaciation; changes in carbonate chemistry and marine biology are dominant during the first and second parts of the glacial cycle, respectively. These feedback mechanisms could also significantly impact the ultimate climate response to the anthropogenic perturbation.

Description: The Middle Miocene was one of the last warm periods of the Neogene, culminating with the Middle Miocene Climatic Optimum (MMCO, approximatively 17–15 Ma). Several proxy-based reconstructions support warmer and more humid climate during the MMCO. The mechanisms responsible for the warmer climate at the MMCO and particularly the role of the atmospheric carbon dioxide are still highly debated. Here we carried out a series of sensitivity experiments with the model of intermediate complexity Planet Simulator, investigating the contributions of the absence of ice on the continents, the opening of the Central American and Eastern Tethys Seaways, the lowering of the topography on land, the effect of various atmospheric CO2 concentrations and the vegetation feedback. Our results show that a higher than present-day CO2 concentration is necessary to generate a warmer climate at all latitudes at the Middle Miocene, in agreement with the terrestrial proxy reconstructions which suggest high atmospheric CO2 concentrations at the MMCO. Nevertheless, the changes in sea-surface conditions, the lowering of the topography on land and the vegetation feedback also produce significant local warming that may, locally, even be stronger than the CO2 induced temperature increases. The lowering of the topography leads to a more zonal atmospheric circulation and allows the westerly flow to continue over the lowered Plateaus at mid-latitudes. The reduced height of the Tibetan Plateau notably prevents the development of a monsoon-like circulation, whereas the reduction of elevations of the North American and European reliefs strongly increases precipitation from northwestern to eastern Europe. The changes in vegetation cover contribute to maintain and even to intensify the warm and humid conditions produced by the other factors, suggesting that the vegetation-climate interactions could help to improve the model-data comparison.

Description: During glacial-interglacial cycles, atmospheric CO2 concentration varied by about 100 ppmv in amplitude. While testing mechanisms that had led to the low glacial CO2 level could be done in equilibrium model experiments, an ultimate goal is to explain CO2 changes in transient simulations through the complete glacial-interglacial cycle. A computationally efficient Earth System model of intermediate complexity CLIMBER-2 is used to simulate global biogeochemistry over the last glacial cycle (126 kyr). The physical core of the model (atmosphere, ocean, land and ice sheets) is driven by orbital changes and reconstructed radiative forcing from greenhouses gases, ice, and aeolian dust. The carbon cycle model is able to reproduce the main features of the CO2 changes: a 50 ppmv CO2 drop during glacial inception, a minimum concentration at the last glacial maximum by 80 ppmv lower than the Holocene value, and an abrupt 60 ppmv CO2 rise during the deglaciation. The model deep ocean δ13C also resembles reconstructions from deep-sea cores. The main drivers of atmospheric CO2 evolve with time: changes in sea surface temperatures and in the volume of bottom water of southern origin controls atmospheric CO2 during the glacial inception and deglaciation, while changes in carbonate chemistry and marine biology are dominant during the first and second parts of the glacial cycle, respectively. These feedback mechanisms could also significantly impact the ultimate climate response to the anthropogenic perturbation.

Description: Many sensitivity studies have been carried out, using climate models of different degrees of complexity to test the climate response to Last Glacial Maximum boundary conditions. Here, instead of adding the forcings successively as in most previous studies, we applied the separation method of U. Stein et P. Alpert 1993, in order to determine rigorously the different contributions of the boundary condition modifications, and isolate the pure contributions from the interactions among the forcings. We carried out a series of sensitivity experiments with the model of intermediate complexity Planet Simulator, investigating the contributions of the ice sheet expansion and elevation, the lowering of the atmospheric CO2 and of the vegetation cover change on the LGM climate. The separation of the ice cover and orographic contributions shows that the ice albedo effect is the main contributor to the cooling of the Northern Hemisphere, whereas orography has only a local cooling impact over the ice sheets. The expansion of ice cover in the Northern Hemisphere causes a disruption of the tropical precipitation, and a southward shift of the ITCZ. The orographic forcing mainly contributes to the disruption of the atmospheric circulation in the Northern Hemisphere, leading to a redistribution of the precipitation, but weakly impacts the tropics. The isolated vegetation contribution also induces strong cooling over the continents of the Northern Hemisphere that further affects the tropical precipitation and reinforce the southward shift of the ITCZ, when combined with the ice forcing. The combinations of the forcings generate many non-linear interactions that reinforce or weaken the pure contributions, depending on the climatic mechanism involved, but they are generally weaker than the pure contributions. Finally, the comparison between the LGM simulated climate and climatic reconstructions over Eurasia suggests that our results reproduce well the south-west to north-east temperature gradients over Eurasia.

Description: Many sensitivity studies have been carried out, using simplified GCMs to test the climate response to Last Glacial Maximum boundary conditions. Here, instead of adding the forcings successively as in previous studies, we applied the separation method of Stein and Alpert (1993), in order to determine rigourously the different contributions of the boundary condition modifications, and isolate the pure contributions from the interactions among the forcings. We carried out a series of sensitivity experiments with the model of intermediate complexity Planet Simulator, investigating the contributions of the ice sheet expansion and elevation, the lowering of the atmospheric CO2 and of the vegetation cover change on the LGM climate. The results clearly identify the ice cover forcing as the main contributor to the cooling of the Northern Hemisphere, and also to the tropical precipitation disruption, leading to the shouthward shift of the ITCZ, while the orographic forcing mainly contributes to the disruption of the atmospheric circulation in the Northern Hemisphere. The isolated vegetation contribution also induces strong cooling over the continents of the Northern Hemisphere, that is further sufficient to affect the tropical precipitation and reinforce the southwards shift of the ITCZ, when combined with the ice forcing. The combinations of the forcings generate many non linear interactions, that reinforce or weaken the pure contributions, depending on the climatic mechanism involved, but they are generally weaker than the pure contributions. Finally, the comparison between the LGM simulated climate and climatic reconstructions over Eurasia suggests that our results reproduce well the south-west to north-east temperature gradients over Eurasia.

Description: The Middle Miocene was one of the last warm periods of the Neogene, culminating with the Middle Miocene Climatic Optimum (MMCO, approximatively 17–15 Ma). Several proxy-based reconstructions support warmer and more humid climate during the MMCO. The mechanisms responsible for the warming at MMCO and particulary the role of the atmospheric carbon dioxide CO2 are still highly debated. Here we carried out a series of sensitivity experiments with the model of intermediate complexity Planet Simulator, investigating the contributions of the absence of ice on the continents, the opening of the Central American and Eastern Tethys Seaways, the lowering of the topography on land, the effect of various atmospheric CO2 concentrations and the vegetation retroaction. Our results show that a higher than present-day CO2 concentration is necessary to generate a warmer climate at all latitudes at the Middle Miocene, in agreement with the terrestrial proxy reconstructions which suggest high atmospheric CO2 concentrations at MMCO. Nevertheless, the changes in sea-surface conditions and the lowering of the topography on land also produce significant local warming that may, locally, even be stronger than the CO2 induced temperature increases. The lowering of the topography leads to a more zonal atmospheric circulation and allows the westerly flow to continue over the lowered Plateaus at mid-latitudes. The reduced height of the Tibetan Plateau notably prevents the development of a monsoon-like circulation, whereas the reduction of elevations of the North American and European reliefs strongly increases precipitation from northwestern to eastern Europe. The changes in vegetation cover contributes to maintain and even to intensify the the warm and humid conditions produced by the other factors, suggesting that the vegetation-climate interactions could help to improve the model-data comparison.

Description: Trends in the atmospheric concentration of CO2 during three recent interglacials, the Holocene, the Eemian and Marine Isotope Stage (MIS) 11, are investigated using an Earth system Model of Intermediate Complexity, which we extended with modules to dynamically determine two slow carbon cycle processes – peat accumulation and shallow-water CaCO3 sedimentation (coral reef formation). For all three interglacials, model simulations considering peat accumulation and shallow water CaCO3 sedimentation substantially improve the agreement between model results and ice core CO2 reconstructions in comparison to a carbon cycle setup neglecting these processes. This enables us to model the trends in atmospheric CO2, with modelled trends similar to the ice core data, forcing the model only with orbital and sea level changes. During the Holocene, anthropogenic CO2 emissions are required to match the observed rise in atmospheric CO2 after 3 ka BP, but are not relevant before this time. Therefore our model experiments show for the first time how the CO2 evolution during the Holocene and two recent interglacials can be explained consistently using an identical model setup.

Description: The total alkalinity–pH equation, which relates total alkalinity and pH for a given set of total concentrations of the acid-base systems that contribute to total alkalinity in a given water sample, is reviewed and its mathematical properties established. We prove that the equation function is strictly monotone and always has exactly one positive root. Different commonly used approximations are discussed and compared. An original method to derive appropriate initial values for the iterative solution of the cubic polynomial equation based upon carbonate-borate-alkalinity is presented. We then review different methods that have been used to solve the total alkalinity–pH equation, with a main focus on biogeochemical models. The shortcomings and limitations of these methods are made out and discussed. We then present two variants of a new, robust and universally convergent algorithm to solve the total alkalinity–pH equation. This algorithm does not require any a priori knowledge of the solution. The iterative procedure is shown to converge from any starting value to the physical solution. The extra computational cost for the convergence security is only 10–15% compared to the fastest algorithm in our test series.

Description: The main characteristics of the new version 1.2 of the three-dimensional Earth system model of intermediate complexity LOVECLIM are briefly described. LOVECLIM 1.2 includes representations of the atmosphere, the ocean and sea ice, the land surface (including vegetation), the ice sheets, the icebergs and the carbon cycle. The atmospheric component is ECBilt2, a T21, 3-level quasi-geostrophic model. The oceanic component is CLIO3, which is made up of an ocean general circulation model coupled to a comprehensive thermodynamic-dynamic sea-ice model. Its horizontal resolution is 3° by 3°, and there are 20 levels in the ocean. ECBilt-CLIO is coupled to VECODE, a vegetation model that simulates the dynamics of two main terrestrial plant functional types, trees and grasses, as well as desert. VECODE also simulates the evolution of the carbon cycle over land while the oceanic carbon cycle is represented in LOCH, a comprehensive model that takes into account both the solubility and biological pumps. The ice sheet component AGISM is made up of a three-dimensional thermomechanical model of the ice sheet flow, a visco-elastic bedrock model and a model of the mass balance at the ice-atmosphere and ice ocean interfaces. For both the Greenland and Antarctic ice sheets, calculations are made on a 10 km by 10 km resolution grid with 31 sigma levels. LOVECLIM 1.2 reproduces well the major characteristics of the observed climate both for present-day conditions and for key past periods such as the last millennium, the mid-Holocene and the Last Glacial Maximum. However, despite some improvements compared to earlier versions, some biases are still present in the model. The most serious ones are mainly located at low latitudes with an overestimation of the temperature there, a too symmetric distribution of precipitation between the two hemispheres, an overestimation of precipitation and vegetation cover in the subtropics. In addition, the atmospheric circulation is too weak. The model also tends to underestimate the surface temperature changes (mainly at low latitudes) and to overestimate the ocean heat uptake observed over the last decades.