OMIP contribution to CMIP6: experimental and diagnostic protocol for the physical component of the Ocean Model Intercomparison Project

The Ocean Model Intercomparison Project (OMIP) is an endorsed project in the Coupled Model Intercomparison Project Phase 6 (CMIP6). OMIP addresses CMIP6 science questions, investigating the origins and consequences of systematic model biases. It does so by providing a framework for evaluating (including assessment of systematic biases), understanding, and improving ocean, sea-ice, tracer, and biogeochemical components of climate and earth system models contributing to CMIP6. Among the WCRP Grand Challenges in climate science (GCs), OMIP primarily contributes to the regional sea level change and near-term (climate/decadal) prediction GCs. OMIP provides (a) an experimental protocol for global ocean/sea-ice models run with a prescribed atmospheric forcing; and (b) a protocol for ocean diagnostics to be saved as part of CMIP6. We focus here on the physical component of OMIP, with a companion paper (Orr et al., 2016) detailing methods for the inert chemistry and interactive biogeochemistry. The physical portion of the OMIP experimental protocol follows the interannual Coordinated Ocean-ice Reference Experiments (CORE-II). Since 2009, CORE-I (Normal Year Forcing) and CORE-II (Interannual Forcing) have become the standard methods to evaluate global ocean/sea-ice simulations and to examine mechanisms for forced ocean climate variability. The OMIP diagnostic protocol is relevant for any ocean model component of CMIP6, including the DECK (Diagnostic, Evaluation and Characterization of Klima experiments), historical simulations, FAFMIP (Flux Anomaly Forced MIP), C4MIP (Coupled Carbon Cycle Climate MIP), DAMIP (Detection and Attribution MIP), DCPP (Decadal Climate Prediction Project), ScenarioMIP, HighResMIP (High Resolution MIP), as well as the ocean/sea-ice OMIP simulations.

The role of the Indian Ocean sector and sea surface salinity for prediction of the coupled Indo-Pacific system

The purpose of this dissertation is to evaluate the potential downstream influence of the Indian Ocean (IO) on El Niño/Southern Oscillation (ENSO) forecasts through the oceanic pathway of the Indonesian Throughflow (ITF), atmospheric teleconnections between the IO and Pacific, and assimilation of IO observations. Also the impact of sea surface salinity (SSS) in the Indo-Pacific region is assessed to try to address known problems with operational coupled model precipitation forecasts. The ITF normally drains warm fresh water from the Pacific reducing the mixed layer depths (MLD). A shallower MLD amplifies large-scale oceanic Kelvin/Rossby waves thus giving ~10% larger response and more realistic ENSO sea surface temperature (SST) variability compared to observed when the ITF is open. In order to isolate the impact of the IO sector atmospheric teleconnections to ENSO, experiments are contrasted that selectively couple/decouple the interannual forcing in the IO. The interannual variability of IO SST forcing is responsible for 3 month lagged widespread downwelling in the Pacific, assisted by off-equatorial curl, leading to warmer NINO3 SST anomaly and improved ENSO validation (significant from 3-9 months). Isolating the impact of observations in the IO sector using regional assimilation identifies large-scale warming in the IO that acts to intensify the easterlies of the Walker circulation and increases pervasive upwelling across the Pacific, cooling the eastern Pacific, and improving ENSO validation (r ~ 0.05, RMS~0.08C). Lastly, the positive impact of more accurate fresh water forcing is demonstrated to address inadequate precipitation forecasts in operational coupled models. Aquarius SSS assimilation improves the mixed layer density and enhances mixing, setting off upwelling that eventually cools the eastern Pacific after 6 months, counteracting the pervasive warming of most coupled models and significantly improving ENSO validation from 5-11 months. In summary, the ITF oceanic pathway, the atmospheric teleconnection, the impact of observations in the IO, and improved Indo-Pacific SSS are all responsible for ENSO forecast improvements, and so each aspect of this study contributes to a better overall understanding of ENSO. Therefore, the upstream influence of the IO should be thought of as integral to the functioning of ENSO phenomenon.

A coupled ocean-atmosphere system for short term forecasting in the southern european seas

This coupled model combines two state-of-the-art numerical models, NEMO for the oceanic component and WRF for the atmospheric component and implements them at an appropriate resolution. The oceanic model has been implemented starting from the Mediterranean Forecasting System with a resolution of 1/24° and the domain was extended to exactly match the grid of a newly implemented atmospheric model for the same area. The uncoupled ocean model has been validated against SST observed data, both in the simulation of an extreme event and in the short-term forecast of two seasonal periods. A new setup of the model was successfully tested in which the downward radiative fluxes were prescribed from atmospheric forecasts. Various physical schemes, domain, boundary, and initial conditions were tested with the atmospheric model to obtain the best representation of medicane Ianos. The heat fluxes calculated by the uncoupled models were compared to determine which setup gave the best energy balance between the components of the coupled model. The coupling strategy used is the traditional one, where the ocean is driven by the surface stress, heat fluxes, and radiative fluxes computed in the atmospheric component, which in turn receives the SST and surface currents. As expected, the overall skills of the coupled model are slightly degraded compared to the uncoupled models, even though the positioning and timing of the cyclone at the time of the landfall is enhanced. The mean heat fluxes do not change compared to the uncoupled model, whereas the pattern of the shortwave radiation and latent heat is changed. Moreover, the two energy fluxes are larger in absolute values than those calculated with the MFS formulas. The fact that they have opposite signs give raise to a compensation error that limits the overall degradation of the coupled simulation.

A state-space modeling approach for active structural acoustic control

The demands for improvement in sound quality and reduction of noise generated by vehicles are constantly increasing, as well as the penalties for space and weight of the control solutions. A promising approach to cope with this challenge is the use of active structural-acoustic control. Usually, the low frequency noise is transmitted into the vehicle`s cabin through structural paths, which raises the necessity of dealing with vibro-acoustic models. This kind of models should allow the inclusion of sensors and actuators models, if accurate performance indexes are to be accessed. The challenge thus resides in deriving reasonable sized models that integrate structural, acoustic, electrical components and the controller algorithm. The advantages of adequate active control simulation strategies relies on the cost and time reduction in the development phase. Therefore, the aim of this paper is to present a methodology for simulating vibro-acoustic systems including this coupled model in a closed loop control simulation framework that also takes into account the interaction between the system and the control sensors/actuators. It is shown that neglecting the sensor/actuator dynamics can lead to inaccurate performance predictions.

Rocket tracking and decoupling eigenstructure control law

A control law was designed for a satellite launcher ( rocket ) vehicle using eigenstructure assignment in order that the vehicle tracks a reference attitude and also to decouple the yaw response from roll and pitch manoeuvres and to decouple the pitch response from roll and yaw manoeuvres. The design was based on a complete linear coupled model obtained from the complete vehicle non linear model by linearization at each trajectory point. After all, the design was assessed with the vehicle time varying non-linear model showing a good performance and robustness. The used design method is explained and a case study for the Brazilian satellite launcher ( VLS Rocket ) is reported.

The frequency and dynamics of stratospheric sudden warmings in the 21st century

Changes to stratospheric sudden warmings (SSWs) over the coming century, as predicted by the Geophysical Fluid Dynamics Laboratory (GFDL) chemistry climate model [Atmospheric Model With Transport and Chemistry (AMTRAC)], are investigated in detail. Two sets of integrations, each a three-member ensemble, are analyzed. The first set is driven with observed climate forcings between 1960 and 2004; the second is driven with climate forcings from a coupled model run, including trace gas concentrations representing a midrange estimate of future anthropogenic emissions between 1990 and 2099. A small positive trend in the frequency of SSWs is found. This trend, amounting to 1 event/decade over a century, is statistically significant at the 90% confidence level and is consistent over the two sets of model integrations. Comparison of the model SSW climatology between the late 20th and 21st centuries shows that the increase is largest toward the end of the winter season. In contrast, the dynamical properties are not significantly altered in the coming century, despite the increase in SSW frequency. Owing to the intrinsic complexity of our model, the direct cause of the predicted trend in SSW frequency remains an open question.

Subseasonal extremes of precipitation and active-break cycles of the Indian summer monsoon in a climate-change scenario

Changes to the behaviour of subseasonal precipitation extremes and active-break cycles of the Indian summer monsoon are assessed in this study using pre-industrial and 2 × CO2 integrations of the Hadley Centre coupled model HadCM3, which is able to simulate the monsoon seasonal cycle reasonably. At 2 × CO2, mean summer rainfall increases slightly, especially over central and northern India. The mean intensity of daily precipitation during the monsoon is found to increase, consistent with fewer wet days, and there are increases to heavy rain events beyond changes in the mean alone. The chance of reaching particular thresholds of heavy rainfall is found to approximately double over northern India, increasing the likelihood of damaging floods on a seasonal basis. The local distribution of such projections is uncertain, however, given the large spread in mean monsoon rainfall change and associated extremes amongst even the most recent coupled climate models. The measured increase of the heaviest precipitation events over India is found to be broadly in line with the degree of atmospheric warming and associated increases in specific humidity, lending a degree of predictability to changes in rainfall extremes. Active-break cycles of the Indian summer monsoon, important particularly due to their effect on agricultural output, are shown to be reasonably represented in HadCM3, in particular with some degree of northward propagation. We note an intensification of both active and break events, particularly when measured against the annual cycle, although there is no suggestion of any change to the duration or likelihood of monsoon breaks. Copyright © 2009 Royal Meteorological Society

Decadal Prediction. Can It Be Skillful?

A new field of study, “decadal prediction,” is emerging in climate science. Decadal prediction lies between seasonal/interannual forecasting and longer-term climate change projections, and focuses on time-evolving regional climate conditions over the next 10–30 yr. Numerous assessments of climate information user needs have identified this time scale as being important to infrastructure planners, water resource managers, and many others. It is central to the information portfolio required to adapt effectively to and through climatic changes. At least three factors influence time-evolving regional climate at the decadal time scale: 1) climate change commitment (further warming as the coupled climate system comes into adjustment with increases of greenhouse gases that have already occurred), 2) external forcing, particularly from future increases of greenhouse gases and recovery of the ozone hole, and 3) internally generated variability. Some decadal prediction skill has been demonstrated to arise from the first two of these factors, and there is evidence that initialized coupled climate models can capture mechanisms of internally generated decadal climate variations, thus increasing predictive skill globally and particularly regionally. Several methods have been proposed for initializing global coupled climate models for decadal predictions, all of which involve global time-evolving three-dimensional ocean data, including temperature and salinity. An experimental framework to address decadal predictability/prediction is described in this paper and has been incorporated into the coordinated Coupled Model Intercomparison Model, phase 5 (CMIP5) experiments, some of which will be assessed for the IPCC Fifth Assessment Report (AR5). These experiments will likely guide work in this emerging field over the next 5 yr.

Using idealized snow forcing to test teleconnections with the Indian summer monsoon in the Hadley Centre GCM

Anomalous heavy snow during winter or spring has long been regarded as a possible precursor of deficient Indian monsoon rainfall during the subsequent summer. However previous work in this field is inconclusive, in terms of the mechanism that communicates snow anomalies to the monsoon summer, and even the region from which snow has the most impact. In this study we explore these issues in coupled and atmosphere-only versions of the Hadley Centre model. A 1050-year control integration of the HadCM3 coupled model, which well represents the seasonal cycle of snow cover over the Eurasian continent, is analysed and shows evidence for weakened monsoons being preceded by strong snow forcing (in the absence of ENSO) over either the Himalaya/Tibetan Plateau or north/west Eurasia regions. However, empirical orthogonal function (EOF) analysis of springtime interannual variability in snow depth shows the leading mode to have opposite signs between these two regions, suggesting that competing mechanisms may be possible. To determine the dominant region, ensemble integrations are carried out using HadAM3, the atmospheric component of HadCM3, and a variety of anomalous snow forcing initial conditions obtained from the control integration of the coupled model. Forcings are applied during spring in separate experiments over the Himalaya/Tibetan Plateau and north/west Eurasia regions, in conjunction with climatological SSTs in order to avoid the direct effects of ENSO. With the aid of idealized forcing conditions in sensitivity tests, we demonstrate that forcing from the Himalaya region is dominant in this model via a Blanford-type mechanism involving reduced surface sensible heat and longwave fluxes, reduced heating of the troposphere over the Tibetan Plateau and consequently a reduced meridional tropospheric temperature gradient which weakens the monsoon during early summer. Snow albedo is shown to be key to the mechanism, explaining around 50% of the perturbation in sensible heating over the Tibetan Plateau, and accounting for the majority of cooling through the troposphere.

The interaction of the Madden-Julian Oscillation with the Maritime Continent in a GCM

Previous studies using the Hadley Centre coupled model (HadCM3) have shown that the islands of the Maritime Continent act as an unrealistic block to the eastward propagation of the Madden-Julian Oscillation (MJO). This blocking effect is investigated using a simplified, aqua-planet version of this GCM, with various idealized configurations of the Maritime Continent islands placed on the equator, and an MJO-like convective signal forced by a propagating sea-surface temperature anomaly dipole. Results suggest that it is the orography of the islands, rather than the presence of the islands themselves, which results in the blocking of the MJO. Although the peak elevation of the orography in the GCM is very much lower than in reality, it appears to act as effective block to the eastward propagation of the low-level Kelvin wave signal which accompanies the MJO. In particular, the representation of Sumatra in the GCM, as a north-south oriented ridge straddling the equator, seems to be particularly effective at blocking the Kelvin wave signal, which in a full GCM would result in the weakening or complete extinction of the MJO signal to the east of the Maritime Continent.

Investigating the causes of the response of the thermohaline circulation to past and future climate changes

The Atlantic thermohaline circulation (THC) is an important part of the earth's climate system. Previous research has shown large uncertainties in simulating future changes in this critical system. The simulated THC response to idealized freshwater perturbations and the associated climate changes have been intercompared as an activity of World Climate Research Program (WCRP) Coupled Model Intercomparison Project/Paleo-Modeling Intercomparison Project (CMIP/PMIP) committees. This intercomparison among models ranging from the earth system models of intermediate complexity (EMICs) to the fully coupled atmosphere-ocean general circulation models (AOGCMs) seeks to document and improve understanding of the causes of the wide variations in the modeled THC response. The robustness of particular simulation features has been evaluated across the model results. In response to 0.1-Sv (1 Sv equivalent to 10(6) ms(3) s(-1)) freshwater input in the northern North Atlantic, the multimodel ensemble mean THC weakens by 30% after 100 yr. All models simulate sonic weakening of the THC, but no model simulates a complete shutdown of the THC. The multimodel ensemble indicates that the surface air temperature could present a complex anomaly pattern with cooling south of Greenland and warming over the Barents and Nordic Seas. The Atlantic ITCZ tends to shift southward. In response to 1.0-Sv freshwater input, the THC switches off rapidly in all model simulations. A large cooling occurs over the North Atlantic. The annual mean Atlantic ITCZ moves into the Southern Hemisphere. Models disagree in terms of the reversibility of the THC after its shutdown. In general, the EMICs and AOGCMs obtain similar THC responses and climate changes with more pronounced and sharper patterns in the AOGCMs.

Metrics for solar wind prediction models: Comparison of empirical, hybrid, and physics-based schemes with 8 years of L1 observations

Space weather effects on technological systems originate with energy carried from the Sun to the terrestrial environment by the solar wind. In this study, we present results of modeling of solar corona-heliosphere processes to predict solar wind conditions at the L1 Lagrangian point upstream of Earth. In particular we calculate performance metrics for (1) empirical, (2) hybrid empirical/physics-based, and (3) full physics-based coupled corona-heliosphere models over an 8-year period (1995–2002). L1 measurements of the radial solar wind speed are the primary basis for validation of the coronal and heliosphere models studied, though other solar wind parameters are also considered. The models are from the Center for Integrated Space-Weather Modeling (CISM) which has developed a coupled model of the whole Sun-to-Earth system, from the solar photosphere to the terrestrial thermosphere. Simple point-by-point analysis techniques, such as mean-square-error and correlation coefficients, indicate that the empirical coronal-heliosphere model currently gives the best forecast of solar wind speed at 1 AU. A more detailed analysis shows that errors in the physics-based models are predominately the result of small timing offsets to solar wind structures and that the large-scale features of the solar wind are actually well modeled. We suggest that additional “tuning” of the coupling between the coronal and heliosphere models could lead to a significant improvement of their accuracy. Furthermore, we note that the physics-based models accurately capture dynamic effects at solar wind stream interaction regions, such as magnetic field compression, flow deflection, and density buildup, which the empirical scheme cannot.

Ocean circulation and climate in an idealised Pangean OAGCM

An idealised Pangean configuration is integrated in a coupled ocean atmosphere general circulation model to investigate the form of the ocean circulation and its impacts on the large scale climate system. A vigorous, hemispherically symmetric overturning is found, driven by deep water formation at high latitudes. Whilst the peak mass transport is around 100Sv, a low vertical temperature gradient in the ocean means that the maximum heat transport is only 1.2PW. The geographical change in the coupled model is found to produce a global average warming of 2°C, despite an increase in global surface albedo. This occurs through changes in the atmospheric water vapour and cloud distributions. There is also reduction in the equator-pole temperature gradient, largely attributable to the same causes, avoiding the paradox of low meridional temperature gradients without increased polar heat transport.

The impact of global freshwater forcing on the thermohaline circulation: adjustment of North Atlantic convection sites in a CGCM

On the time scale of a century, the Atlantic thermohaline circulation (THC) is sensitive to the global surface salinity distribution. The advection of salinity toward the deep convection sites of the North Atlantic is one of the driving mechanisms for the THC. There is both a northward and a southward contributions. The northward salinity advection (Nsa) is related to the evaporation in the subtropics, and contributes to increased salinity in the convection sites. The southward salinity advection (Ssa) is related to the Arctic freshwater forcing and tends on the contrary to diminish salinity in the convection sites. The THC changes results from a delicate balance between these opposing mechanisms. In this study we evaluate these two effects using the IPSL-CM4 ocean-atmosphere-sea-ice coupled model (used for IPCC AR4). Perturbation experiments have been integrated for 100 years under modern insolation and trace gases. River runoff and evaporation minus precipitation are successively set to zero for the ocean during the coupling procedure. This allows the effect of processes Nsa and Ssa to be estimated with their specific time scales. It is shown that the convection sites in the North Atlantic exhibit various sensitivities to these processes. The Labrador Sea exhibits a dominant sensitivity to local forcing and Ssa with a typical time scale of 10 years, whereas the Irminger Sea is mostly sensitive to Nsa with a 15 year time scale. The GIN Seas respond to both effects with a time scale of 10 years for Ssa and 20 years for Nsa. It is concluded that, in the IPSL-CM4, the global freshwater forcing damps the THC on centennial time scales.

Forecast calibration and combination: a simple Bayesian approach for ENSO

This study presents a new simple approach for combining empirical with raw (i.e., not bias corrected) coupled model ensemble forecasts in order to make more skillful interval forecasts of ENSO. A Bayesian normal model has been used to combine empirical and raw coupled model December SST Niño-3.4 index forecasts started at the end of the preceding July (5-month lead time). The empirical forecasts were obtained by linear regression between December and the preceding July Niño-3.4 index values over the period 1950–2001. Coupled model ensemble forecasts for the period 1987–99 were provided by ECMWF, as part of the Development of a European Multimodel Ensemble System for Seasonal to Interannual Prediction (DEMETER) project. Empirical and raw coupled model ensemble forecasts alone have similar mean absolute error forecast skill score, compared to climatological forecasts, of around 50% over the period 1987–99. The combined forecast gives an increased skill score of 74% and provides a well-calibrated and reliable estimate of forecast uncertainty.