Resolution-dependent variations of sinking particle trajectories in general circulation models: Implications for data-model comparison in past climate

Author(s):  
Peter Nooteboom ◽  
Philippe Delandmeter ◽  
Peter Bijl ◽  
Erik van Sebille ◽  
Henk Dijkstra ◽  
...  

<p>Any type of non-buoyant material in the ocean is transported by currents during its sinking journey. This transport can be far from negligible for typical (plankton) particles with a low sinking velocity. To estimate the lateral transport, the material can be modelled as a set of Lagrangian particles advected by currents that are obtained from Ocean General Circulation Models (OGCMs). State-of-the-art OGCMs are often strongly eddying, providing flow fields with a horizontal resolution of  10km on a daily basis. However, many long term climate modelling studies (e.g. in palaeoclimate) rely on low resolution models that cannot capture mesoscale features. The lower model resolution could influence data-model comparisons using Lagrangian techniques, but this is not properly evaluated yet through a direct comparison.</p><p>In this study, we simulate the transport of sinking Lagrangian particles using low (1°; non-eddying)  and high (0.1°; eddying) horizontal resolution OGCMs of the present-day ocean, and evaluate the effect of the two resolutions on particle transport. We find major differences between the transport in the non-eddying versus the eddying OGCM (in terms of the divergence of particle trajectories and their mean trajectory). Addition of stochastic noise to the particle trajectory parameterizes the effect of eddies well in some regions (e.g. in the North Pacific gyre).</p><p>We recommend to apply sinking Lagrangian particles only in velocity fields with eddying OGCMs, which basically excludes all paleo-simulations. We are currently simulating the equilibrium Eocene (38Ma) climate using an eddying OGCM, to be able to apply these Lagrangian techniques in an eddying ocean of the past. We expect this to lead towards a better agreement between the OGCM and sedimentary fossil microplankton.</p>

2007 ◽  
Vol 3 (3) ◽  
pp. 499-512 ◽  
Author(s):  
S. Brewer ◽  
J. Guiot ◽  
F. Torre

Abstract. We present here a comparison between the outputs of 25 General Circulation Models run for the mid-Holocene period (6 ka BP) with a set of palaeoclimate reconstructions based on over 400 fossil pollen sequences distributed across the European continent. Three climate parameters were available (moisture availability, temperature of the coldest month and growing degree days), which were grouped together using cluster analysis to provide regions of homogenous climate change. Each model was then investigated to see if it reproduced 1) similar patterns of change and 2) the correct location of these regions. A fuzzy logic distance was used to compare the output of the model with the data, which allowed uncertainties from both the model and data to be taken into account. The models were compared by the magnitude and direction of climate change within the region as well as the spatial pattern of these changes. The majority of the models are grouped together, suggesting that they are becoming more consistent. A test against a set of zero anomalies (no climate change) shows that, although the models are unable to reproduce the exact patterns of change, they all produce the correct signs of change observed for the mid-Holocene.


2021 ◽  
Author(s):  
Enrico Scoccimarro ◽  
Daniele Peano ◽  
Silvio Gualdi ◽  
Alessio Bellucci ◽  
Tomas Lovato ◽  
...  

Abstract. The recent advancements in climate modelling partially build on the improvement of horizontal resolution in different components of the simulating system. A higher resolution is expected to provide a better representation of the climate variability, and in this work we are particularly interested in the potential improvements in representing extreme events of high temperature and precipitation. The two versions of the CMCC-CM2 model used here, adopt the highest horizontal resolutions available within the last family of the global coupled climate models de¬veloped at CMCC to participate in the CMIP6 effort. The main aim of this study is to document the ability of the CMCC-CM2 models in representing the spatial distribution of extreme events of temperature and precipitation, under the historical period, comparing model results to observations (ERA5 Reanalysis and CHIRPS observations). For a more detailed evaluation we investigate both 6 hourly and daily time series for the definition of the extreme conditions. In terms of mean climate, the two models are able to realistically reproduce the main patterns of temperature and precipitation. The very-high resolution version (¼ degree horizontal resolution) of the atmospheric model provides better results than the high resolution one (one degree), not only in terms of means but also in terms of extreme events of temperature defined at daily and 6-hourly frequency. This is also the case of average precipitation. On the other hand the extreme precipitation is not improved by the adoption of a higher horizontal resolution.


2015 ◽  
Vol 28 (11) ◽  
pp. 4513-4524 ◽  
Author(s):  
Jeff Willison ◽  
Walter A. Robinson ◽  
Gary M. Lackmann

Abstract Mesoscale condensational heating can increase the sensitivity of modeled extratropical cyclogenesis to horizontal resolution. Here a pseudo global warming experiment is presented to investigate how this heating-enhanced sensitivity to resolution changes in a warmer and thus moister atmosphere. The Weather Research and Forecasting (WRF) Model with 120- and 20-km grid spacing is used to simulate current and future climates. It is found that the North Atlantic storm-track response to global warming is amplified at the higher model resolution. The most dramatic changes occur over the northeastern Atlantic, where resolution typical of current general circulation models (GCMs) results in a smaller global warming response in comparison with that in the 20-km simulations. These results suggest that caution is warranted when interpreting projections from coarse-resolution GCMs of future cyclone activity over the northeastern Atlantic.


1995 ◽  
Vol 43 (2) ◽  
pp. 147-158 ◽  
Author(s):  
Anatoly V. Lozhkin ◽  
Patricia M. Anderson

AbstractAlluvial, fluvial, and organic deposits of the last interglaciation are exposed along numerous river terraces in northeast Siberia. Although chronological control is often poor, the paleobotanical data suggest range extensions of up to 1000 km for the primary tree species. These data also indicate that boreal communities of the last interglaciation were similar to modern ones in composition, but their distributions were displaced significantly to the north-northwest. Inferences about climate of this period suggest that mean July temperatures were warmer by 4 to 8°C, and seasonal precipitation was slightly greater. Mean January temperatures may have been severely cooler than today (up to 12°C) along the Arctic coast, but similar or slightly warmer than present in other areas. The direction and magnitude of change in July temperatures agree with Atmospheric General Circulation Models, but the 126,000-year-B.P. model results also suggest trends opposite to the paleobotanical data, with simulated cooler winter temperatures and drier conditions than present during the climatic optimum.


2018 ◽  
Vol 31 (14) ◽  
pp. 5437-5459 ◽  
Author(s):  
Hui Ding ◽  
Matthew Newman ◽  
Michael A. Alexander ◽  
Andrew T. Wittenberg

Seasonal forecasts made by coupled atmosphere–ocean general circulation models (CGCMs) undergo strong climate drift and initialization shock, driving the model state away from its long-term attractor. Here we explore initializing directly on a model’s own attractor, using an analog approach in which model states close to the observed initial state are drawn from a “library” obtained from prior uninitialized CGCM simulations. The subsequent evolution of those “model-analogs” yields a forecast ensemble, without additional model integration. This technique is applied to four of the eight CGCMs comprising the North American Multimodel Ensemble (NMME) by selecting from prior long control runs those model states whose monthly tropical Indo-Pacific SST and SSH anomalies best resemble the observations at initialization time. Hindcasts are then made for leads of 1–12 months during 1982–2015. Deterministic and probabilistic skill measures of these model-analog hindcast ensembles are comparable to those of the initialized NMME hindcast ensembles, for both the individual models and the multimodel ensemble. In the eastern equatorial Pacific, model-analog hindcast skill exceeds that of the NMME. Despite initializing with a relatively large ensemble spread, model-analogs also reproduce each CGCM’s perfect-model skill, consistent with a coarse-grained view of tropical Indo-Pacific predictability. This study suggests that with little additional effort, sufficiently realistic and long CGCM simulations provide the basis for skillful seasonal forecasts of tropical Indo-Pacific SST anomalies, even without sophisticated data assimilation or additional ensemble forecast integrations. The model-analog method could provide a baseline for forecast skill when developing future models and forecast systems.


2021 ◽  
Author(s):  
Luca Famooss Paolini ◽  
Alessio Bellucci ◽  
Paolo Ruggieri ◽  
Panos Athanasiadis ◽  
Silvio Gualdi

<p>Western boundary currents transport a large amount of heat from the Tropics toward higher latitudes; furthermore they are characterized by a strong sea surface temperature (SST) gradient, which anchors zones of intense upward motion extending up to the upper-troposphere and shapes zones of intense baroclinic eddy activity (storm tracks). For such reasons they have been shown to be fundamental in influencing the climate of the Northern Hemisphere and its variability, and a potentially relevant source of atmospheric predictability. </p><p> </p><p>General circulation models show deficiencies in simulating the observed atmospheric response to SST front variability. The atmospheric horizontal resolution has been recently proposed as a key element in understanding such differences. However, the number of studies on this subject is still limited. Furthermore, a multi-model analysis to systematically investigate differences between low-resolution and high-resolution atmospheric response to oceanic forcing is still lacking. </p><p> </p><p>The present work has the objective to fill this gap, analysing the atmospheric response to Gulf Stream SST front shifting using data from recent High Resolution Model Intercomparison Project (HighResMIP). This project was designed with the specific objective of investigating the impact of increased model horizontal resolution on the representation of the observed climate. Ensembles of historical simulations performed with three atmospheric general circulation models (AGCMs) have been analysed, each conducted with a low-resolution (LR, about 1°) and a high-resolution (HR, about 0.25°) configuration. AGCMs have been forced with observed SSTs (HadISST2 dataset), available at daily frequency on a 0.25° grid, during 1950–2014. </p><p><br>Results show atmospheric responses to the SST-induced diabatic heating anomalies that are strongly resolution dependent. In LR simulations a low-pressure anomaly is present downstream of the SST anomaly, while the diabatic heating anomaly is mainly balanced by meridional advection of air coming from higher latitudes, as expected for an extra-tropical shallow heat source. In contrast, HR simulations generate a high-pressure anomaly downstream of the SST anomaly, thus driving positive temperature advection from lower latitudes (not balancing diabatic heating). Along the vertical direction, both in LR and HR simulation, the diabatic heating in the interior of the atmosphere is balanced by upward motion south of GS SST front and downward motion north and further south of the Gulf Stream. Finally, LR simulations show a reduction in storm-track activity over the North Atlantic, whereas HR simulations show a meridional displacement of the storm-track considerably larger (yet in the same direction) than that of the SST front. HR simulations reproduce the atmospheric response obtained from observations, albeit weaker. This is a hint for the existence of a positive feedback between ocean and atmosphere, as proposed in previous studies. These findings are qualitatively consistent with previous results in literature and, leveraging on recent coordinated modelling efforts, shed light on the effective role of atmospheric horizontal resolution in modelling the atmospheric response to extra-tropical oceanic forcing.</p>


2021 ◽  
Author(s):  
Peter Nooteboom ◽  
Michiel Baatsen ◽  
Peter Bijl ◽  
Erik van Sebille ◽  
Appy Sluijs ◽  
...  

<p>Simulations of the geological past using General Circulation Models (GCMs) are computationally expensive. Mainly because of the long equilibration time scales, most of these GCMs have ocean components with a horizontal resolution of 1° or coarser. Such models are non-eddying and the effects of mesoscale ocean eddies on the transport of heat and salt are parameterized. However, from present-day ocean modeling studies, it is known that eddying ocean models better represent regional and time-mean ocean flows compared to non-eddying models. At the same time, proxy data from sediment sample sites represent climate at specific locations. Hence, the coarse ocean resolution of typical palaeo-GCMs lead to a challenge for model-data comparison in past climates.</p><p>Here we present the first simulations of a global eddying Eocene ocean with a 0.1° (horizontal) resolution model, which are initialized and forced with data from a coarser resolution (1° horizontally) equilibrated coupled ocean-atmosphere GCM. We investigate the response of the model equilibrium state to the change in ocean resolution and the consequences this has for model-data comparison in the middle-late Eocene (38Ma). We find that, compared to the non-eddying model, the eddying ocean resolution of palaeomodels reduce the biases in both sea surface temperatures and biogeographic patterns which are derived from proxy data.</p>


Ocean Science ◽  
2011 ◽  
Vol 7 (3) ◽  
pp. 389-404 ◽  
Author(s):  
I. Medhaug ◽  
T. Furevik

Abstract. Output from a total of 24 state-of-the-art Atmosphere-Ocean General Circulation Models is analyzed. The models were integrated with observed forcing for the period 1850–2000 as part of the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report. All models show enhanced variability at multi-decadal time scales in the North Atlantic sector similar to the observations, but with a large intermodel spread in amplitudes and frequencies for both the Atlantic Multidecadal Oscillation (AMO) and the Atlantic Meridional Overturning Circulation (AMOC). The models, in general, are able to reproduce the observed geographical patterns of warm and cold episodes, but not the phasing such as the early warming (1930s–1950s) and the following colder period (1960s–1980s). This indicates that the observed 20th century extreme in temperatures are due to primarily a fortuitous phasing of intrinsic climate variability and not dominated by external forcing. Most models show a realistic structure in the overturning circulation, where more than half of the available models have a mean overturning transport within the observed estimated range of 13–24 Sverdrup. Associated with a stronger than normal AMOC, the surface temperature is increased and the sea ice extent slightly reduced in the North Atlantic. Individual models show potential for decadal prediction based on the relationship between the AMO and AMOC, but the models strongly disagree both in phasing and strength of the covariability. This makes it difficult to identify common mechanisms and to assess the applicability for predictions.


2019 ◽  
Author(s):  
Allison C. Michaelis ◽  
Gary M. Lackmann ◽  
Walter A. Robinson

Abstract. We present multi-seasonal simulations representative of present-day and future thermodynamic environments using the global Model for Prediction Across Scales-Atmosphere (MPAS) version 5.1 with high resolution (15 km) throughout the Northern Hemisphere. We select ten simulation years with varying phases of El Niño-Southern Oscillation (ENSO) and integrate each for 14.5 months. We use analysed sea surface temperature (SST) patterns for present-day simulations. For the future climate simulations, we alter present-day SSTs by applying monthly-averaged temperature changes derived from a 20-member ensemble of Coupled Model Intercomparison Project Phase 5 (CMIP5) general circulation models (GCMs) following the Representative Concentration Pathway (RCP) 8.5 emissions scenario. Daily sea ice fields, obtained from the monthly-averaged CMIP5 ensemble mean sea ice, are used for present-day and future simulations. The present-day simulations provide a reasonable reproduction of large-scale atmospheric features in the Northern Hemisphere such as the wintertime midlatitude storm tracks, upper-tropospheric jets, and maritime sea-level pressure features as well as annual precipitation patterns across the tropics. The simulations also adequately represent tropical cyclone (TC) characteristics such as strength, spatial distribution, and seasonal cycles for most of Northern Hemispheric basins. These results demonstrate the applicability of these model simulations for future studies examining climate change effects on various Northern Hemispheric phenomena, and, more generally, the utility of MPAS for studying climate change at spatial scales generally unachievable in GCMs.


Author(s):  
J.D Annan ◽  
J.C Hargreaves

In this paper, we review progress towards efficiently estimating parameters in climate models. Since the general problem is inherently intractable, a range of approximations and heuristic methods have been proposed. Simple Monte Carlo sampling methods, although easy to implement and very flexible, are rather inefficient, making implementation possible only in the very simplest models. More sophisticated methods based on random walks and gradient-descent methods can provide more efficient solutions, but it is often unclear how to extract probabilistic information from such methods and the computational costs are still generally too high for their application to state-of-the-art general circulation models (GCMs). The ensemble Kalman filter is an efficient Monte Carlo approximation which is optimal for linear problems, but we show here how its accuracy can degrade in nonlinear applications. Methods based on particle filtering may provide a solution to this problem but have yet to be studied in any detail in the realm of climate models. Statistical emulators show great promise for future research and their computational speed would eliminate much of the need for efficient sampling techniques. However, emulation of a full GCM has yet to be achieved and the construction of such represents a substantial computational task in itself.


Sign in / Sign up

Export Citation Format

Share Document