A coupled cloud physics–radiation parameterization of the bulk optical properties of cirrus and its impact on the Met Office unified model global atmosphere 5.0 configuration

A new coupled cloud physics–radiation parameterization of the bulk optical properties of ice clouds is presented. The parameterization is consistent with assumptions in the cloud physics scheme regarding particle size distributions (PSDs) and mass–dimensional relationships. The parameterization is based on a weighted ice crystal habit mixture model, and its bulk optical properties are parameterized as simple functions of wavelength and ice water content (IWC). This approach directly couples IWC to the bulk optical properties, negating the need for diagnosed variables, such as the ice crystal effective dimension. The parameterization is implemented into the Met Office Unified Model Global Atmosphere 5.0 (GA5) configuration. The GA5 configuration is used to simulate the annual 20-yr shortwave (SW) and longwave (LW) fluxes at the top of the atmosphere (TOA), as well as the temperature structure of the atmosphere, under various microphysical assumptions. The coupled parameterization is directly compared against the current operational radiation parameterization, while maintaining the same cloud physics assumptions. In this experiment, the impacts of the two parameterizations on the SW and LW radiative effects at TOA are also investigated and compared against observations. The 20-yr simulations are compared against the latest observations of the atmospheric temperature and radiative fluxes at TOA. The comparisons demonstrate that the choice of PSD and the assumed ice crystal shape distribution are as important as each other. Moreover, the consistent radiation parameterization removes a long-standing tropical troposphere cold temperature bias but slightly warms the southern midlatitudes by about 0.5 K.

The Met Office Unified Model Global Atmosphere 3.0/3.1 and JULES Global Land 3.0/3.1 configurations

We describe Global Atmosphere 3.0 (GA3.0): a configuration of the Met Office Unified Model (MetUM) developed for use across climate research and weather prediction activities. GA3.0 has been formulated by converging the development paths of the Met Office's weather and climate global atmospheric model components such that wherever possible, atmospheric processes are modelled or parametrized seamlessly across spatial resolutions and timescales. This unified development process will provide the Met Office and its collaborators with regular releases of a configuration that has been evaluated, and can hence be applied, over a variety of modelling régimes. We also describe Global Land 3.0 (GL3.0): a configuration of the JULES community land surface model developed for use with GA3.0.

Further considerations of cosmic ray modulation of infra-red radiation in the atmosphere

Understanding effects of ionisation in the lower atmosphere is a new interdisciplinary area, crossing the traditionally distinct scientific boundaries between astro-particle and atmospheric physics and also requiring understanding of both heliospheric and magnetospheric influences on cosmic rays. Following the paper of Erlykin et al. (2014) we develop further the interpretation of our observed changes in long-wave (LW) radiation, Aplin and Lockwood (2013) by taking account of both cosmic ray ionisation yields and atmospheric radiative transfer. To demonstrate this, we show that the thermal structure of the whole atmosphere needs to be considered along with the vertical profile of ionisation. Allowing for, in particular, ionisation by all components of a cosmic ray shower and not just by the muons, reveals that the effect we have detected is certainly not inconsistent with laboratory observations of the LW absorption cross section. The analysis presented here, although very different from that of Erlykin et al., does come to the same conclusion that the events detected by AL were not caused by individual cosmic ray primaries – not because it is impossible on energetic grounds, but because events of the required energy are too infrequent for the 12 h_1 rate at which they were seen by the AL experiment. The present paper numerically models the effect of three different scenario changes to the primary GCR spectrum which all reproduce the required magnitude of the effect observed by AL. However, they cannot solely explain the observed delay in the peak effect which, if confirmed, would appear to open up a whole new and interesting area in the study of water oligomers and their effects on LW radiation. We argue that a technical artefact in the AL experiment is highly unlikely and that our initial observations merit both a wide-ranging follow-up experiment and more rigorous, self-consistent, three-dimensional radiative transfer modelling

Exploring strategies for coupled 4D-Var data assimilation using an idealised atmosphere-ocean model

Operational forecasting centres are currently developing data assimilation systems for coupled atmosphere-ocean models. Strongly coupled assimilation, in which a single assimilation system is applied to a coupled model, presents significant technical and scientific challenges. Hence weakly coupled assimilation systems are being developed as a first step, in which the coupled model is used to compare the current state estimate with observations, but corrections to the atmosphere and ocean initial conditions are then calculated independently. In this paper we provide a comprehensive description of the different coupled assimilation methodologies in the context of four dimensional variational assimilation (4D-Var) and use an idealised framework to assess the expected benefits of moving towards coupled data assimilation. We implement an incremental 4D-Var system within an idealised single column atmosphere-ocean model. The system has the capability to run both strongly and weakly coupled assimilations as well as uncoupled atmosphere or ocean only assimilations, thus allowing a systematic comparison of the different strategies for treating the coupled data assimilation problem. We present results from a series of identical twin experiments devised to investigate the behaviour and sensitivities of the different approaches. Overall, our study demonstrates the potential benefits that may be expected from coupled data assimilation. When compared to uncoupled initialisation, coupled assimilation is able to produce more balanced initial analysis fields, thus reducing initialisation shock and its impact on the subsequent forecast. Single observation experiments demonstrate how coupled assimilation systems are able to pass information between the atmosphere and ocean and therefore use near-surface data to greater effect. We show that much of this benefit may also be gained from a weakly coupled assimilation system, but that this can be sensitive to the parameters used in the assimilation.

Origin and impact of initialisation shocks in coupled atmosphere-ocean forecasts

Current methods for initialising coupled atmosphere-ocean forecasts often rely on the use of separate atmosphere and ocean analyses, the combination of which can leave the coupled system imbalanced at the beginning of the forecast, potentially accelerating the development of errors. Using a series of experiments with the European Centre for Medium-range Weather Forecasts coupled system, the magnitude and extent of these so-called initialisation shocks is quantified, and their impact on forecast skill measured. It is found that forecasts initialised by separate ocean and atmospheric analyses do exhibit initialisation shocks in lower atmospheric temperature, when compared to forecasts initialised using a coupled data assimilation method. These shocks result in as much as a doubling of root-mean-square error on the first day of the forecast in some regions, and in increases that are sustained for the duration of the 10-day forecasts performed here. However, the impacts of this choice of initialisation on forecast skill, assessed using independent datasets, were found to be negligible, at least over the limited period studied. Larger initialisation shocks are found to follow a change in either the atmospheric or ocean model component between the analysis and forecast phases: changes in the ocean component can lead to sea surface temperature shocks of more than 0.5K in some equatorial regions during the first day of the forecast. Implications for the development of coupled forecast systems, particularly with respect to coupled data assimilation methods, are discussed.

Energetic particle influence on the Earth's atmosphere

This manuscript gives an up-to-date and comprehensive overview of the effects of energetic particle precipitation (EPP) onto the whole atmosphere, from the lower thermosphere/mesosphere through the stratosphere and troposphere, to the surface. The paper summarizes the different sources and energies of particles, principally galactic cosmic rays (GCRs), solar energetic particles (SEPs) and energetic electron precipitation (EEP). All the proposed mechanisms by which EPP can affect the atmosphere are discussed, including chemical changes in the upper atmosphere and lower thermosphere, chemistry-dynamics feedbacks, the global electric circuit and cloud formation. The role of energetic particles in Earth’s atmosphere is a multi-disciplinary problem that requires expertise from a range of scientific backgrounds. To assist with this synergy, summary tables are provided, which are intended to evaluate the level of current knowledge of the effects of energetic particles on processes in the entire atmosphere.

Atmosphere-ocean coupled processes in the Madden-Julian oscillation

The Madden-Julian oscillation (MJO) is a convectively coupled 30-70 day (intraseasonal) tropical atmospheric mode that drives variations in global weather, but which is poorly simulated in most atmospheric general circulation models. Over the past two decades, field campaigns and modeling experiments have suggested that tropical atmosphere-ocean interactions may sustain or amplify the pattern of enhanced and suppressed atmospheric convection that defines the MJO, and encourage its eastward propagation through the Indian and Pacific Oceans. New observations collected during the past decade have advanced our understand of the ocean response to atmospheric MJO forcing and the resulting intraseasonal sea surface temperature (SST) fluctuations. Numerous modeling studies have revealed a considerable impact of the mean state on MJO ocean-atmosphere coupled processes, as well as the importance of resolving the diurnal cycle of atmosphere--upper-ocean interactions. New diagnostic methods provide insight to atmospheric variability and physical processes associated with the MJO, but offer limited insight on the role of ocean feedbacks. Consequently, uncertainty remains concerning the role of the ocean in MJO theory. Our understanding of how atmosphere-ocean coupled processes affect the MJO can be improved by collecting observations in poorly sampled regions of MJO activity, assessing oceanic and atmospheric drivers of surface fluxes, improving the representation of upper-ocean mixing in coupled-model simulations, designing model experiments that minimize mean-state differences, and developing diagnostic tools to evaluate the nature and role of coupled ocean-atmosphere processes over the MJO cycle.

Regionally coupled atmosphere-ocean-sea ice-marine biogeochemistry model ROM: 1. Description and validation

The general circulation models used to simulate global climate typically feature resolution too coarse to reproduce many smaller-scale processes, which are crucial to determining the regional responses to climate change. A novel approach to downscale climate change scenarios is presented which includes the interactions between the North Atlantic Ocean and the European shelves as well as their impact on the North Atlantic and European climate. The goal of this paper is to introduce the global ocean-regional atmosphere coupling concept and to show the potential benefits of this model system to simulate present-day climate. A global ocean-sea ice-marine biogeochemistry model (MPIOM/HAMOCC) with regionally high horizontal resolution is coupled to an atmospheric regional model (REMO) and global terrestrial hydrology model (HD) via the OASIS coupler. Moreover, results obtained with ROM using NCEP/NCAR reanalysis and ECHAM5/MPIOM CMIP3 historical simulations as boundary conditions are presented and discussed for the North Atlantic and North European region. The validation of all the model components, i.e., ocean, atmosphere, terrestrial hydrology, and ocean biogeochemistry is performed and discussed. The careful and detailed validation of ROM provides evidence that the proposed model system improves the simulation of many aspects of the regional climate, remarkably the ocean, even though some biases persist in other model components, thus leaving potential for future improvement. We conclude that ROM is a powerful tool to estimate possible impacts of climate change on the regional scale.

Current model capabilities for simulating black carbon and sulfate concentrations in the Arctic atmosphere: a multi-model evaluation using a comprehensive measurement data set

The concentrations of sulfate, black carbon (BC) and other aerosols in the Arctic are characterized by high values in late winter and spring (so-called Arctic Haze) and low values in summer. Models have long been struggling to capture this seasonality and especially the high concentrations associated with Arctic Haze. In this study, we evaluate sulfate and BC concentrations from eleven different models driven with the same emission inventory against a comprehensive pan-Arctic measurement data set over a time period of 2 years (2008–2009). The set of models consisted of one Lagrangian particle dispersion model, four chemistry transport models (CTMs), one atmospheric chemistry-weather forecast model and five chemistry climate models (CCMs), of which two were nudged to meteorological analyses and three were running freely. The measurement data set consisted of surface measurements of equivalent BC (eBC) from five stations (Alert, Barrow, Pallas, Tiksi and Zeppelin), elemental carbon (EC) from Station Nord and Alert and aircraft measurements of refractory BC (rBC) from six different campaigns. We find that the models generally captured the measured eBC or rBC and sulfate concentrations quite well, compared to previous comparisons. However, the aerosol seasonality at the surface is still too weak in most models. Concentrations of eBC and sulfate averaged over three surface sites are underestimated in winter/spring in all but one model (model means for January–March underestimated by 59 and 37 % for BC and sulfate, respectively), whereas concentrations in summer are overestimated in the model mean (by 88 and 44 % for July–September), but with overestimates as well as underestimates present in individual models. The most pronounced eBC underestimates, not included in the above multi-site average, are found for the station Tiksi in Siberia where the measured annual mean eBC concentration is 3 times higher than the average annual mean for all other stations. This suggests an underestimate of BC sources in Russia in the emission inventory used. Based on the campaign data, biomass burning was identified as another cause of the modeling problems. For sulfate, very large differences were found in the model ensemble, with an apparent anti-correlation between modeled surface concentrations and total atmospheric columns. There is a strong correlation between observed sulfate and eBC concentrations with consistent sulfate/eBC slopes found for all Arctic stations, indicating that the sources contributing to sulfate and BC are similar throughout the Arctic and that the aerosols are internally mixed and undergo similar removal. However, only three models reproduced this finding, whereas sulfate and BC are weakly correlated in the other models. Overall, no class of models (e.g., CTMs, CCMs) performed better than the others and differences are independent of model resolution.

Stimulation of "stress-regulated" mitogen-activated protein kinases (stress-activated protein kinases/c-Jun N-terminal kinases and p38-mitogen-activated protein kinases) in perfused rat hearts by oxidative and other stresses.

"Stress-regulated" mitogen-activated protein kinases (SR-MAPKs) comprise the stress-activated protein kinases (SAPKs)/c-Jun N-terminal kinases (JNKs) and the p38-MAPKs. In the perfused heart, ischemia/reperfusion activates SR-MAPKs. Although the agent(s) directly responsible is unclear, reactive oxygen species are generated during ischemia/reperfusion. We have assessed the ability of oxidative stress (as exemplified by H2O2) to activate SR-MAPKs in the perfused heart and compared it with the effect of ischemia/reperfusion. H2O2 activated both SAPKs/JNKs and p38-MAPK. Maximal activation by H2O2 in both cases was observed at 0.5 mM. Whereas activation of p38-MAPK by H2O2 was comparable to that of ischemia and ischemia/reperfusion, activation of the SAPKs/JNKs was less than that of ischemia/reperfusion. As with ischemia/reperfusion, there was minimal activation of the ERK MAPK subfamily by H2O2. MAPK-activated protein kinase 2 (MAPKAPK2), a downstream substrate of p38-MAPKs, was activated by H2O2 to a similar extent as with ischemia or ischemia/reperfusion. In all instances, activation of MAPKAPK2 in perfused hearts was inhibited by SB203580, an inhibitor of p38-MAPKs. Perfusion of hearts at high aortic pressure (20 kilopascals) also activated the SR-MAPKs and MAPKAPK2. Free radical trapping agents (dimethyl sulfoxide and N-t-butyl-alpha-phenyl nitrone) inhibited the activation of SR-MAPKs and MAPKAPK2 by ischemia/reperfusion. These data are consistent with a role for reactive oxygen species in the activation of SR-MAPKs during ischemia/reperfusion.

Stimulation of multiple mitogen-activated protein kinase sub-families by oxidative stress and phosphorylation of the small heat shock protein, HSP25/27, in neonatal ventricular myocytes.

We investigated the activation of three subfamilies of mitogen-activated protein kinases (MAPKs), namely the stress-activated protein kinases/c-Jun N-terminal kinases (SAPKs/JNKs), the extracellularly responsive kinases (ERKs) and p38-MAPK, by oxidative stress as exemplified by H2O2 in primary cultures of neonatal rat ventricular myocytes. The 46 and 54 kDa species of SAPKs/JNKs were activated 5- and 10-fold, respectively, by 0.1 mM H2O2 (the maximally effective concentration). Maximal activation occurred at 15-30 min, but was still detectable after 2 h. Both ERK1 and ERK2 were activated 16-fold by 0.1 mM H2O2 with a similar time course to the SAPKs/JNKs, and this was comparable with their activation by 1 microM PMA, the most powerful activator of ERKs that we have so far identified in these cells. The activation of ERKs by H2O2 was inhibited by PD98059, which inhibits the activation of MAPK (or ERK) kinases, and by the protein kinase C (PKC) inhibitor, GF109203X. ERK activation was also inhibited by down-regulation of PMA-sensitive PKC isoforms. p38-MAPK was activated by 0.1 mM H2O2 as shown by an increase in its phosphorylation. However, maximal phosphorylation (activation) was more rapid (<5 min) than for the SAPKs/JNKs or the ERKs. We studied the downstream consequences of p38-MAPK activation by examining activation of MAPK-activated protein kinase 2 (MAPKAPK2) and phosphorylation of the MAPKAPK2 substrate, the small heat shock protein HSP25/27. As with p38-MAPK, MAPKAPK2 was rapidly activated (maximal within 5 min) by 0.1 mM H2O2. This activation was abolished by 10 microM SB203580, a selective inhibitor of certain p38-MAPK isoforms. The phosphorylation of HSP25/27 rapidly followed activation of MAPKAPK2 and was also inhibited by SB203580. Phosphorylation of HSP25/27 was associated with a decrease in its aggregation state. These data indicate that oxidative stress is a powerful activator of all three MAPK subfamilies in neonatal rat ventricular myocytes. Activation of all three MAPKs has been associated with the development of the hypertrophic phenotype. However, stimulation of p38-MAPK and the consequent phosphorylation of HSP25/27 may also be important in cardioprotection.

Regulation of bcl-2 family proteins during development and in response to oxidative stress in cardiac myocytes: association with changes in mitochondrial membrane potential.

Cardiac myocyte apoptosis is potentially important in many cardiac disorders. In other cells, Bcl-2 family proteins and mitochondrial dysfunction are probably key regulators of the apoptotic response. In the present study, we characterized the regulation of antiapoptotic (Bcl-2, Bcl-xL) and proapoptotic (Bad, Bax) Bcl-2 family proteins in the rat heart during development and in oxidative stress-induced apoptosis. Bcl-2 and Bcl-xL were expressed at high levels in the neonate, and their expression was sustained during development. In contrast, although Bad and Bax were present at high levels in neonatal hearts, they were barely detectable in adult hearts. We confirmed that H(2)O(2) induced cardiac myocyte cell death, stimulating poly(ADP-ribose) polymerase proteolysis (from 2 hours), caspase-3 proteolysis (from 2 hours), and DNA fragmentation (from 8 hours). In unstimulated neonatal cardiac myocytes, Bcl-2 and Bcl-xL were associated with the mitochondria, but Bad and Bax were predominantly present in a crude cytosolic fraction. Exposure of myocytes to H(2)O(2) stimulated rapid translocation of Bad (<5 minutes) to the mitochondria. This was followed by the subsequent degradation of Bad and Bcl-2 (from approximately 30 minutes). The levels of the mitochondrial membrane marker cytochrome oxidase remained unchanged. H(2)O(2) also induced translocation of cytochrome c from the mitochondria to the cytosol within 15 to 30 minutes, which was indicative of mitochondrial dysfunction. Myocytes exposed to H(2)O(2) showed an early loss of mitochondrial membrane potential (assessed by fluorescence-activated cell sorter analysis) from 15 to 30 minutes, which was partially restored by approximately 1 hour. However, a subsequent irreversible loss of mitochondrial membrane potential occurred that correlated with cell death. These data suggest that the regulation of Bcl-2 and mitochondrial function are important factors in oxidative stress-induced cardiac myocyte apoptosis.

Regulation of protein kinase B and 4E-BP1 by oxidative stress in cardiac myocytes.

Stimulation of phosphatidylinositol 3'-kinase (PI3K) and protein kinase B (PKB) is implicated in the regulation of protein synthesis in various cells. One mechanism involves PI3K/PKB-dependent phosphorylation of 4E-BP1, which dissociates from eIF4E, allowing initiation of translation from the 7-methylGTP cap of mRNAs. We examined the effects of insulin and H(2)O(2) on this pathway in neonatal cardiac myocytes. Cardiac myocyte protein synthesis was increased by insulin, but was inhibited by H(2)O(2). PI3K inhibitors attenuated basal levels of protein synthesis and inhibited the insulin-induced increase in protein synthesis. Insulin or H(2)O(2) increased the phosphorylation (activation) of PKB through PI3K, but, whereas insulin induced a sustained response, the response to H(2)O(2) was transient. 4E-BP1 was phosphorylated in unstimulated cells, and 4E-BP1 phosphorylation was increased by insulin. H(2)O(2) stimulated dephosphorylation of 4E-BP1 by increasing protein phosphatase (PP1/PP2A) activity. This increased the association of 4E-BP1 with eIF4E, consistent with H(2)O(2) inhibition of protein synthesis. The effects of H(2)O(2) were sufficient to override the stimulation of protein synthesis and 4E-BP1 phosphorylation induced by insulin. These results indicate that PI3K and PKB are important regulators of protein synthesis in cardiac myocytes, but other factors, including phosphatase activity, modulate the overall response.