Antarctic ocean and sea ice response to ozone depletion: a two timescale problem

The response of the Southern Ocean to a repeating seasonal cycle of ozone loss is studied in two coupled climate models and found to comprise both fast and slow processes. The fast response is similar to the inter-annual signature of the Southern Annular Mode (SAM) on Sea Surface Temperature (SST), on to which the ozone-hole forcing projects in the summer. It comprises enhanced northward Ekman drift inducing negative summertime SST anomalies around Antarctica, earlier sea ice freeze-up the following winter, and northward expansion of the sea ice edge year-round. The enhanced northward Ekman drift, however, results in upwelling of warm waters from below the mixed layer in the region of seasonal sea ice. With sustained bursts of westerly winds induced by ozone-hole depletion, this warming from below eventually dominates over the cooling from anomalous Ekman drift. The resulting slow-timescale response (years to decades) leads to warming of SSTs around Antarctica and ultimately a reduction in sea-ice cover year-round. This two-timescale behavior - rapid cooling followed by slow but persistent warming - is found in the two coupled models analysed, one with an idealized geometry, the other a complex global climate model with realistic geometry. Processes that control the timescale of the transition from cooling to warming, and their uncertainties are described. Finally we discuss the implications of our results for rationalizing previous studies of the effect of the ozone-hole on SST and sea-ice extent. %Interannual variability in the Southern Annular Mode (SAM) and sea ice covary such that an increase and southward shift in the surface westerlies (a positive phase of the SAM) coincides with a cooling of Sea Surface Temperature (SST) around 70-50$^\circ$S and an expansion of the sea ice cover, as seen in observations and models alike. Yet, in modeling studies, the Southern Ocean warms and sea ice extent decreases in response to sustained, multi-decadal positive SAM-like wind anomalies driven by 20th century ozone depletion. Why does the Southern Ocean appear to have disparate responses to SAM-like variability on interannual and multidecadal timescales? Here it is demonstrated that the response of the Southern Ocean to ozone depletion has a fast and a slow response. The fast response is similar to the interannual variability signature of the SAM. It is dominated by an enhanced northward Ekman drift, which transports heat northward and causes negative SST anomalies in summertime, earlier sea ice freeze-up the following winter, and northward expansion of the sea ice edge year round. The enhanced northward Ekman drift causes a region of Ekman divergence around 70-50$^\circ$S, which results in upwelling of warmer waters from below the mixed layer. With sustained westerly wind enhancement in that latitudinal band, the warming due to the anomalous upwelling of warm waters eventually dominates over the cooling from the anomalous Ekman drift. Hence, the slow response ultimately results in a positive SST anomaly and a reduction in the sea ice cover year round. We demonstrate this behavior in two models: one with an idealized geometry and another, more detailed, global climate model. However, the models disagree on the timescale of transition from the fast (cooling) to the slow (warming) response. Processes that controls this transition and their uncertainties are discussed.

The dependence of radiative forcing and feedback on evolving patterns of surface temperature change in climate models

Experiments with CO2 instantaneously quadrupled and then held constant are used to show that the relationship between the global-mean net heat input to the climate system and the global-mean surface-air-temperature change is nonlinear in Coupled Model Intercomparison Project phase 5 (CMIP5) Atmosphere-Ocean General Circulation Models (AOGCMs). The nonlinearity is shown to arise from a change in strength of climate feedbacks driven by an evolving pattern of surface warming. In 23 out of the 27 AOGCMs examined the climate feedback parameter becomes significantly (95% confidence) less negative – i.e. the effective climate sensitivity increases – as time passes. Cloud feedback parameters show the largest changes. In the AOGCM-mean approximately 60% of the change in feedback parameter comes from the topics (30N-30S). An important region involved is the tropical Pacific where the surface warming intensifies in the east after a few decades. The dependence of climate feedbacks on an evolving pattern of surface warming is confirmed using the HadGEM2 and HadCM3 atmosphere GCMs (AGCMs). With monthly evolving sea-surface-temperatures and sea-ice prescribed from its AOGCM counterpart each AGCM reproduces the time-varying feedbacks, but when a fixed pattern of warming is prescribed the radiative response is linear with global temperature change or nearly so. We also demonstrate that the regression and fixed-SST methods for evaluating effective radiative forcing are in principle different, because rapid SST adjustment when CO2 is changed can produce a pattern of surface temperature change with zero global mean but non-zero change in net radiation at the top of the atmosphere (~ -0.5 Wm-2 in HadCM3).

Evaluating the ability of process based models to project sea-level change

We evaluate the ability of process based models to reproduce observed global mean sea-level change. When the models are forced by changes in natural and anthropogenic radiative forcing of the climate system and anthropogenic changes in land-water storage, the average of the modelled sea-level change for the periods 1900–2010, 1961–2010 and 1990–2010 is about 80%, 85% and 90% of the observed rise. The modelled rate of rise is over 1 mm yr−1 prior to 1950, decreases to less than 0.5 mm yr−1 in the 1960s, and increases to 3 mm yr−1 by 2000. When observed regional climate changes are used to drive a glacier model and an allowance is included for an ongoing adjustment of the ice sheets, the modelled sea-level rise is about 2 mm yr−1 prior to 1950, similar to the observations. The model results encompass the observed rise and the model average is within 20% of the observations, about 10% when the observed ice sheet contributions since 1993 are added, increasing confidence in future projections for the 21st century. The increased rate of rise since 1990 is not part of a natural cycle but a direct response to increased radiative forcing (both anthropogenic and natural), which will continue to grow with ongoing greenhouse gas emissions

Intercomparison of general circulation models for hot extrasolar planets

We compare five general circulation models (GCMs) which have been recently used to study hot extrasolar planet atmospheres (BOB, CAM, IGCM, MITgcm, and PEQMOD), under three test cases useful for assessing model convergence and accuracy. Such a broad, detailed intercomparison has not been performed thus far for extrasolar planets study. The models considered all solve the traditional primitive equations, but employ di↵erent numerical algorithms or grids (e.g., pseudospectral and finite volume, with the latter separately in longitude-latitude and ‘cubed-sphere’ grids). The test cases are chosen to cleanly address specific aspects of the behaviors typically reported in hot extrasolar planet simulations: 1) steady-state, 2) nonlinearly evolving baroclinic wave, and 3) response to fast timescale thermal relaxation. When initialized with a steady jet, all models maintain the steadiness, as they should—except MITgcm in cubed-sphere grid. A very good agreement is obtained for a baroclinic wave evolving from an initial instability in pseudospectral models (only). However, exact numerical convergence is still not achieved across the pseudospectral models: amplitudes and phases are observably di↵erent. When subject to a typical ‘hot-Jupiter’-like forcing, all five models show quantitatively di↵erent behavior—although qualitatively similar, time-variable, quadrupole-dominated flows are produced. Hence, as have been advocated in several past studies, specific quantitative predictions (such as the location of large vortices and hot regions) by GCMs should be viewed with caution. Overall, in the tests considered here, pseudospectral models in pressure coordinate (PEBOB and PEQMOD) perform the best and MITgcm in cubed-sphere grid performs the worst.

Centennial changes in solar activity and the response of galactic cosmic rays

There is a growing consensus that the eleven year modulation of galactic cosmic rays (GCRs) resulting from solar activity is related to interplanetary propagating diffusive barriers (PDBs). The source of these PDBs is not well understood and numerical models describing GCR modulation simulate their effect by scaling the diffusion tensor to the interplanetary magnetic field strength (IMF). The implications of a century-scale change in solar wind speed and open solar flux, for numerical modelling of GCR modulation and the reconstruction of GCR variations over the last hundred years are discussed. The dominant role of the solar non-axisymmetric magnetic field in both forcing longitudinal solar wind speed fluctuations at solar maximum and in increasing the IMF is discussed in the context of a long-term rise in the open solar magnetic flux.

A numerical model of the ionospheric signatures of time-varying magnetic reconnection: II. Measuring expansions in the ionospheric flow response

A numerical model embodying the concepts of the Cowley-Lockwood (Cowley and Lockwood, 1992, 1997) paradigm has been used to produce a simple Cowley– Lockwood type expanding flow pattern and to calculate the resulting change in ion temperature. Cross-correlation, fixed threshold analysis and threshold relative to peak are used to determine the phase speed of the change in convection pattern, in response to a change in applied reconnection. Each of these methods fails to fully recover the expansion of the onset of the convection response that is inherent in the simulations. The results of this study indicate that any expansion of the convection pattern will be best observed in time-series data using a threshold which is a fixed fraction of the peak response. We show that these methods used to determine the expansion velocity can be used to discriminate between the two main models for the convection response to a change in reconnection.

Non-replication of the association between 5HTTLPR and response to psychological therapy for child anxiety disorders

Background We previously reported an association between 5HTTLPR genotype and outcome following cognitive–behavioural therapy (CBT) in child anxiety (Cohort 1). Children homozygous for the low-expression short-allele showed more positive outcomes. Other similar studies have produced mixed results, with most reporting no association between genotype and CBT outcome. Aims To replicate the association between 5HTTLPR and CBT outcome in child anxiety from the Genes for Treatment study (GxT Cohort 2, n = 829). Method Logistic and linear mixed effects models were used to examine the relationship between 5HTTLPR and CBT outcomes. Mega-analyses using both cohorts were performed. Results There was no significant effect of 5HTTLPR on CBT outcomes in Cohort 2. Mega-analyses identified a significant association between 5HTTLPR and remission from all anxiety disorders at follow-up (odds ratio 0.45, P = 0.014), but not primary anxiety disorder outcomes. Conclusions The association between 5HTTLPR genotype and CBT outcome did not replicate. Short-allele homozygotes showed more positive treatment outcomes, but with small, non-significant effects. Future studies would benefit from utilising whole genome approaches and large, homogenous samples.

Deconstructing the climate change response of the Northern Hemisphere wintertime storm tracks

There are large uncertainties in the circulation response of the atmosphere to climate change. One manifestation of this is the substantial spread in projections for the extratropical storm tracks made by different state-of-the-art climate models. In this study we perform a series of sensitivity experiments, with the atmosphere component of a single climate model, in order to identify the causes of the differences between storm track responses in different models. In particular, the Northern Hemisphere wintertime storm tracks in the CMIP3 multi-model ensemble are considered. A number of potential physical drivers of storm track change are identified and their influence on the storm tracks is assessed. The experimental design aims to perturb the different physical drivers independently, by magnitudes representative of the range of values present in the CMIP3 model runs, and this is achieved via perturbations to the sea surface temperature and the sea-ice concentration forcing fields. We ask the question: can the spread of projections for the extratropical storm tracks present in the CMIP3 models be accounted for in a simple way by any of the identified drivers? The results suggest that, whilst the changes in the upper-tropospheric equator-to-pole temperature difference have an influence on the storm track response to climate change, the large spread of projections for the extratropical storm track present in the northern North Atlantic in particular is more strongly associated with changes in the lower-tropospheric equator-to-pole temperature difference.

Bond graph models of DC-DC converters operating in both CCM and DCM

In this paper, Bond Graphs are employed to develop a novel mathematical model of conventional switched-mode DC-DC converters valid for both continuous and discontinuous conduction modes. A unique causality bond graph model of hybrid models is suggested with the operation of the switch and the diode to be represented by a Modulated Transformer with a binary input and a resistor with fixed conductance causality. The operation of the diode is controlled using an if-then function within the model. The extracted hybrid model is implemented on a Boost and Buck converter with their operations to change from CCM to DCM and to return to CCM. The vector fields of the models show validity in a wide operation area and comparison with the simulation of the converters using PSPICE reveals high accuracy of the proposed model, with the Normalised Root Means Square Error and the Maximum Absolute Error remaining adequately low. The model is also experimentally tested on a Buck topology.

Vegetation and soils feedbacks on the response of the African monsoon response to orbital forcing in early to middle Holocene

Fossil pollen, ancient lake sediments and archaeological evidence from Africa indicate that the Sahel and Sahara regions were considerably wetter than today during the early to middle Holocene period, about 12,000 to 5,000 years ago1–4. Vegetation associated with the modern Sahara/Sahel boundary was about 5° farther north, and there were more and larger lakes between 15 and 30° N. Simulations with climate models have shown that these wetter conditions were probably caused by changes in Earth's orbital parameters that increased the amplitude of the seasonal cycle of solar radiation in the Northern Hemisphere, enhanced the land-ocean temperature contrast, and thereby strengthened the African summer monsoon5–7. However, these simulations underestimated the consequent monsoon enhancement as inferred from palaeorecords4. Here we use a climate model to show that changes in vegetation and soil may have increased the climate response to orbital forcing. We find that replacing today's orbital forcing with that of the mid-Holocene increases summer precipitation by 12% between 15 and 22° N. Replacing desert with grassland, and desert soil with more loamy soil, further enhances the summer precipitation (by 6 and 10% respectively), giving a total precipitation increase of 28%. When the simulated climate changes are applied to a biome model, vegetation becomes established north of the current Sahara/Sahel boundary, thereby shrinking the area of the Sahara by 11% owing to orbital forcing alone, and by 20% owing to the combined influence of orbital forcing and the prescribed vegetation and soil changes. The inclusion of the vegetation and soil feedbacks thus brings the model simulations and palaeovegetation observations into closer agreement.

Clinical predictors of response to cognitive-behavioural therapy in pediatric anxiety disorders: the Genes for Treatment (GXT) Study

Objective The Genes for Treatment study is an international, multisite collaboration exploring the role of genetic, demographic, and clinical predictors in response to cognitive-behavioral therapy (CBT) in pediatric anxiety disorders. The current article, the first from the study, examined demographic and clinical predictors of response to CBT. We hypothesized that the child’s gender, type of anxiety disorder, initial severity and comorbidity, and parents’ psychopathology would significantly predict outcome. Method A sample of 1,519 children 5 to 18 years of age with a primary anxiety diagnosis received CBT across 11 sites. Outcome was defined as response (change in diagnostic severity) and remission (absence of the primary diagnosis) at each time point (posttreatment, 3-, 6-, and/or 12-month follow-up) and analyzed using linear and logistic mixed models. Separate analyses were conducted using data from posttreatment and follow-up assessments to explore the relative importance of predictors at these time points. Results Individuals with social anxiety disorder (SoAD) had significantly poorer outcomes (poorer response and lower rates of remission) than those with generalized anxiety disorder (GAD). Although individuals with specific phobia (SP) also had poorer outcomes than those with GAD at posttreatment, these differences were not maintained at follow-up. Both comorbid mood and externalizing disorders significantly predicted poorer outcomes at posttreatment and follow-up, whereas self-reported parental psychopathology had little effect on posttreatment outcomes but significantly predicted response (although not remission) at follow-up. Conclusion SoAD, nonanxiety comorbidity, and parental psychopathology were associated with poorer outcomes after CBT. The results highlight the need for enhanced treatments for children at risk for poorer outcomes.

The dependence of wintertime Mediterranean precipitation on the atmospheric circulation response to climate change

Climate models indicate a future wintertime precipitation reduction in the Mediterranean region but there is large uncertainty in the amplitude of the projected change. We analyse CMIP5 climate model output to quantify the role of atmospheric circulation in the Mediterranean precipitation change. It is found that a simple circulation index, i.e. the 850 hPa zonal wind (U850) in North Africa, well describes the year to year fluctuations in the area-averaged Mediterranean precipitation, with positive (i.e. westerly) U850 anomalies in North Africa being associated with positive precipitation anomalies. Under climate change, U850 in North Africa and the Mediterranean precipitation are both projected to decrease consistently with the relationship found in the inter-annual variability. This enables us to estimate that about 85% of the CMIP5 mean precipitation response and 80% of the variance in the inter-model spread are related to changes in the atmospheric circulation. In contrast, there is no significant correlation between the mean precipitation response and the global-mean surface warming across the models. It follows that the uncertainty in cold-season Mediterranean precipitation projection will not be narrowed unless the uncertainty in the atmospheric circulation response is reduced.

HPA axis related genes and response to psychological therapies: genetics and epigenetics

Background Hypothalamic–pituitary–adrenal (HPA) axis functioning has been implicated in the development of stress-related psychiatric diagnoses and response to adverse life experiences. This study aimed to investigate the association between genetic and epigenetics in HPA axis and response to cognitive behavior therapy (CBT). Methods Children with anxiety disorders were recruited into the Genes for Treatment project (GxT, N = 1,152). Polymorphisms of FKBP5 and GR were analyzed for association with response to CBT. Percentage DNA methylation at the FKBP5 and GR promoter regions was measured before and after CBT in a subset (n = 98). Linear mixed effect models were used to investigate the relationship between genotype, DNA methylation, and change in primary anxiety disorder severity (treatment response). Results Treatment response was not associated with FKBP5 and GR polymorphisms, or pretreatment percentage DNA methylation. However, change in FKBP5 DNA methylation was nominally significantly associated with treatment response. Participants who demonstrated the greatest reduction in severity decreased in percentage DNA methylation during treatment, whereas those with little/no reduction in severity increased in percentage DNA methylation. This effect was driven by those with one or more FKBP5 risk alleles, with no association seen in those with no FKBP5 risk alleles. No significant association was found between GR methylation and response. Conclusions Allele-specific change in FKBP5 methylation was associated with treatment response. This is the largest study to date investigating the role of HPA axis related genes in response to a psychological therapy. Furthermore, this is the first study to demonstrate that DNA methylation changes may be associated with response to psychological therapies in a genotype-dependent manner.

Predicting how many animals will be where: how to build, calibrate and evaluate individual-based models

Individual-based models (IBMs) can simulate the actions of individual animals as they interact with one another and the landscape in which they live. When used in spatially-explicit landscapes IBMs can show how populations change over time in response to management actions. For instance, IBMs are being used to design strategies of conservation and of the exploitation of fisheries, and for assessing the effects on populations of major construction projects and of novel agricultural chemicals. In such real world contexts, it becomes especially important to build IBMs in a principled fashion, and to approach calibration and evaluation systematically. We argue that insights from physiological and behavioural ecology offer a recipe for building realistic models, and that Approximate Bayesian Computation (ABC) is a promising technique for the calibration and evaluation of IBMs. IBMs are constructed primarily from knowledge about individuals. In ecological applications the relevant knowledge is found in physiological and behavioural ecology, and we approach these from an evolutionary perspective by taking into account how physiological and behavioural processes contribute to life histories, and how those life histories evolve. Evolutionary life history theory shows that, other things being equal, organisms should grow to sexual maturity as fast as possible, and then reproduce as fast as possible, while minimising per capita death rate. Physiological and behavioural ecology are largely built on these principles together with the laws of conservation of matter and energy. To complete construction of an IBM information is also needed on the effects of competitors, conspecifics and food scarcity; the maximum rates of ingestion, growth and reproduction, and life-history parameters. Using this knowledge about physiological and behavioural processes provides a principled way to build IBMs, but model parameters vary between species and are often difficult to measure. A common solution is to manually compare model outputs with observations from real landscapes and so to obtain parameters which produce acceptable fits of model to data. However, this procedure can be convoluted and lead to over-calibrated and thus inflexible models. Many formal statistical techniques are unsuitable for use with IBMs, but we argue that ABC offers a potential way forward. It can be used to calibrate and compare complex stochastic models and to assess the uncertainty in their predictions. We describe methods used to implement ABC in an accessible way and illustrate them with examples and discussion of recent studies. Although much progress has been made, theoretical issues remain, and some of these are outlined and discussed.

Trends in the CERES dataset, 2000–13: the effects of sea ice and jet shifts and comparison to climate models

We review the effects of dynamical variability on clouds and radiation in observations and models and discuss their implications for cloud feedbacks. Jet shifts produce robust meridional dipoles in upper-level clouds and longwave cloud-radiative effect (CRE), but low-level clouds, which do not simply shift with the jet, dominate the shortwave CRE. Because the effect of jet variability on CRE is relatively small, future poleward jet shifts with global warming are only a second-order contribution to the total CRE changes around the midlatitudes, suggesting a dominant role for thermodynamic effects. This implies that constraining the dynamical response is unlikely to reduce the uncertainty in extratropical cloud feedback. However, we argue that uncertainty in the cloud-radiative response does affect the atmospheric circulation response to global warming, by modulating patterns of diabatic forcing. How cloud feedbacks can affect the dynamical response to global warming is an important topic of future research.