Formulation of a new explicit tidal scheme in ocean general circulation model

Tides play an important role in ocean energy transfer and mixing, and provide major energy for maintaining thermohaline circulation. This study proposes a new explicit tidal scheme and assesses its performance in a global ocean model. Instead of using empirical specifications of tidal amplitudes and frequencies, the new scheme directly uses the positions of the Moon and Sun in a global ocean model to incorporate tides. Compared with the traditional method that has specified tidal constituents, the new scheme can better simulate the diurnal and spatial characteristics of the tidal potential of spring and neap tides as well as the spatial patterns and magnitudes of major tidal constituents (K1 and M2). It significantly reduces the total errors of eight tidal constituents (with the exception of N2 and Q1) in the traditional explicit tidal scheme. Relative to the control simulation without tides, both the new and traditional tidal schemes can lead to better dynamic sea level (DSL) simulation in the North Atlantic, reducing significant negative biases in this region. The new tidal scheme also shows smaller positive bias than the traditional scheme in the Southern Ocean. The new scheme is suited to calculate regional distributions of sea level height in addition to tidal mixing.

Regional aspects of global climate change simulations : Validation and assessment of climate response over Indian monsoon region to transient increase of greenhouse gases and sulfate aerosols

The regional climatic impacts associated with global climatic change and their assessment are very important since agriculture, water resources, ecology etc., are all vulnerable to climatic changes on regional scale. Coupled Atmosphere-Ocean general circulation model (AOGCM) simulations provide a range of scenarios, which can be used, for the assessment of impacts and development of adaptive or mitigative strategies. Validation of the models against the observations and establishing the sensitivity to climate change forcing are essential before the model projections are used for assessment of possible impacts. Moreover model simulated climate projections are often of coarse resolution while the models used for impact assessment, (e.g. crop simulation models, or river runoff models etc.) operate on a higher spatial resolution. This spatial mismatch can be overcome by adopting an appropriate strategy of downscaling the GCM output. This study examines two AOGCM (ECHAM4/OPYC3 and HadCM2) climate change simulations for their performance in the simulation of monsoon climate over India and the sensitivity of the simulated monsoon climate to transient changes in the atmospheric concentrations of greenhouse gases and sulfate aerosols. The results show that the two models simulate the gross features of climate over India reasonably well. However the inter-model differences in simulation of mean characteristics, sensitivity to forcing and in the simulation of climate change suggest need for caution. Further an empirical downscaling approach in used to assess the possibility of using GCM projections for preparation of regional climate change scenario for India.

Climate and Ocean Circulation in the Aftermath of a Marinoan Snowball Earth

When a snowball Earth deglaciates through a very high atmospheric CO2 concentration, the resulting inflow of freshwater leads to a stably stratified ocean, and the strong greenhouse conditions drive the climate into a very warm state. Here, we use a coupled atmosphere-ocean general circulation model, applying different scenarios for the evolution of atmospheric CO2, to conduct the first simulation of the climate and the three-dimensional ocean circulation in the aftermath of the Marinoan snowball Earth. The simulations show that the strong freshwater stratification breaks up on a timescale in the order of 103 years, mostly independent of the applied CO2 scenario. This is driven by the upwelling of salty waters in high latitudes, mainly the northern hemisphere, where a strong circumpolar current dominates the circulation. In the warmest CO2 scenario, the simulated Marinoan supergreenhouse climate reaches a global mean surface temperature of about 30 °C under an atmospheric CO2 concentration of 15 × 103 parts per million by volume, which is a moderate temperature compared to previous estimates. Consequently, the thermal expansion of seawater causes a sea-level rise of only 8 m, with most of it occurring during the first 3000 years. Our results imply that the surface temperatures of that time were potentially not as threatening for early metazoa as previously assumed. Furthermore, the short destratification timescale found in this study implies a very rapid accumulation of Marinoan cap dolostones, given that they were deposited in a freshwater environment.

Next generation of Bluelink ocean reanalysis with multiscale data assimilation: BRAN2020

BRAN2020 (2020 version of the Bluelink ReANalysis) is an ocean reanalysis that combines observations with an eddy-resolving, near-global ocean general circulation model to produce a four-dimensional estimate of the ocean state. The data assimilation system employed is ensemble optimal interpolation, implemented with a new multiscale approach that constrains the broad-scale ocean properties and the mesoscale circulation in two steps. There is a separation in the scales that are corrected in the two steps: the high-resolution step corrects the mesoscale dynamics in the same way as previous versions of BRAN, while the extra coarse step is effective at correcting biases that develop at large scales. The reanalysis currently spans January 1993 to December 2019 and assimilates observations of in situ temperature and salinity, as well as of satellite sea-level anomaly and sea surface temperature. BRAN2020 is planned to be updated to within months of real time after this initial release, until an updated version of BRAN is available. Reanalysed fields from BRAN2020 generally show much closer agreement to observations than all previous versions with misfits between reanalysed and observed fields reduced by over 30 % for some variables, for subsurface temperature and salinity in particular. The BRAN2020 dataset is comprised of daily averaged fields of temperature, salinity, velocity, mixed-layer depth and sea level. Reanalysed fields realistically represent all of the major current systems within 75∘ S and 75∘ N, excluding processes relating to sea ice but including boundary currents, equatorial circulation, Southern Ocean variability and mesoscale eddies. BRAN2020 is publicly available at https://doi.org/10.25914/6009627c7af03 (Chamberlain et al., 2021b) and is intended for use by the research community.

A salty deep ocean as a prerequisite for glacial termination

Deglacial transitions of the middle to late Pleistocene (terminations) are linked to gradual changes in insolation accompanied by abrupt shifts in ocean circulation. However, the reason these deglacial abrupt events are so special compared with their sub-glacial-maximum analogues, in particular with respect to the exaggerated warming observed across Antarctica, remains unclear. Here we show that an increase in the relative importance of salt versus temperature stratification in the glacial deep South Atlantic decreases the potential cooling effect of waters that may be upwelled in response to abrupt perturbations in ocean circulation, as compared with sub-glacial-maximum conditions. Using a comprehensive coupled atmosphere–ocean general circulation model, we then demonstrate that an increase in deep-ocean salinity stratification stabilizes relatively warm waters in the glacial deep ocean, which amplifies the high southern latitude surface ocean temperature response to an abrupt weakening of the Atlantic meridional overturning circulation during deglaciation. The mechanism can produce a doubling in the net rate of warming across Antarctica on a multicentennial timescale when starting from full glacial conditions (as compared with interglacial or subglacial conditions) and therefore helps to explain the large magnitude and rapidity of glacial terminations during the late Quaternary.

A discrete interaction numerical model for coagulation and fragmentation of marine detritic particulate matter (Coagfrag v.1)

A simplified model, representing the dynamics of marine organic particles in a given size range experiencing coagulation and fragmentation reactions, is developed. The framework is based on a discrete size spectrum on which reactions act to exchange properties between different particle sizes. The reactions are prescribed according to triplet interactions. Coagulation combines two particle sizes to yield a third one, while fragmentation breaks a given particle size into two (i.e. the inverse of the coagulation reaction). The complete set of reactions is given by all the permutations of two particle sizes associated with a third one. Since, by design, some reactions yield particle sizes that are outside the resolved size range of the spectrum, a closure is developed to take into account this unresolved range and satisfy global constraints such as mass conservation. In order to minimize the number of tracers required to apply this model to an ocean general circulation model, focus is placed on the robustness of the model to the particle size resolution. Thus, numerical experiments were designed to study the dependence of the results on (i) the number of particle size bins used to discretize a given size range (i.e. the resolution) and (ii) the type of discretization (i.e. linear vs. nonlinear). The results demonstrate that in a linearly size-discretized configuration, the model is independent of the resolution. However, important biases are observed in a nonlinear discretization. A first attempt to mitigate the effect of nonlinearity of the size spectrum is then presented and shows significant improvement in reducing the observed biases.

Next generation of Bluelink ocean reanalysis with multiscale data assimilation: BRAN2020

BRAN2020 is an ocean reanalysis that combines ocean observations with an eddy-resolving, near-global ocean general circulation model, to produce four-dimensional estimates of the ocean state. The data assimilation system employed is ensemble optimal interpolation, implemented with a new multiscale approach that constrains the broad-scale ocean properties and the mesoscale circulation in two steps. The reanalysis spans January 1993 to December 2019, and assimilates observations of in situ temperature and salinity, as well as satellite sea-level anomaly and sea surface temperature. Reanalysed fields from BRAN2020 generally show much closer agreement to observations than all previous versions with mis-fits between reanalysed and observed fields reduced by over 30 % for some variables. The BRAN2020 dataset is comprised of daily-averaged fields of temperature, salinity, velocity, mixed-layer depth, and sea-level. Reanalysed fields realistically represent all of the major current systems within 75° S and 75° N, excluding processes relating to sea ice, but including boundary currents, equatorial circulation, Southern Ocean variability, and mesoscale eddies. BRAN2020 is publicly-available at https://doi.org/10.25914/6009627c7af03 (Chamberlain et al., 2021b) and is intended for use by the research community.