Global climate and ocean circulation on an aquaplanet ocean-atmosphere general circulation model

A low resolution coupled ocean-atmosphere general circulation model OAGCM is used to study the characteristics of the large scale ocean circulation and its climatic impacts in a series of global coupled aquaplanet experiments. Three configurations, designed to produce fundamentally different ocean circulation regimes, are considered. The first has no obstruction to zonal flow, the second contains a low barrier that blocks zonal flow in the ocean at all latitudes, creating a single enclosed basin, whilst the third contains a gap in the barrier to allow circumglobal flow at high southern latitudes. Warm greenhouse climates with a global average air surface temperature of around 27C result in all cases. Equator to pole temperature gradients are shallower than that of a current climate simulation. Whilst changes in the land configuration cause regional changes in temperature, winds and rainfall, heat transports within the system are little affected. Inhibition of all ocean transport on the aquaplanet leads to a reduction in global mean surface temperature of 8C, along with a sharpening of the meridional temperature gradient. This results from a reduction in global atmospheric water vapour content and an increase in tropical albedo, both of which act to reduce global surface temperatures. Fitting a simple radiative model to the atmospheric characteristics of the OAGCM solutions suggests that a simpler atmosphere model, with radiative parameters chosen a priori based on the changing surface configuration, would have produced qualitatively different results. This implies that studies with reduced complexity atmospheres need to be guided by more complex OAGCM results on a case by case basis.

Constraining a North Atlantic Ocean general circulation model with chlorofluorocarbon observations

Measurements of anthropogenic tracers such as chlorofluorocarbons and tritium must be quantitatively combined with ocean general circulation models as a component of systematic model development. The authors have developed and tested an inverse method, using a Green's function, to constrain general circulation models with transient tracer data. Using this method chlorofluorocarbon-11 and -12 (CFC-11 and -12) observations are combined with a North Atlantic configuration of the Miami Isopycnic Coordinate Ocean Model with 4/3 degrees resolution. Systematic differences can be seen between the observed CFC concentrations and prior CFC fields simulated by the model. These differences are reduced by the inversion, which determines the optimal gas transfer across the air-sea interface, accounting for uncertainties in the tracer observations. After including the effects of unresolved variability in the CFC fields, the model is found to be inconsistent with the observations because the model/data misfit slightly exceeds the error estimates. By excluding observations in waters ventilated north of the Greenland-Scotland ridge (sigma (0) < 27.82 kg m(-3); shallower than about 2000 m), the fit is improved, indicating that the Nordic overflows are poorly represented in the model. Some systematic differences in the model/data residuals remain and are related, in part, to excessively deep model ventilation near Rockall and deficient ventilation in the main thermocline of the eastern subtropical gyre. Nevertheless, there do not appear to be gross errors in the basin-scale model circulation. Analysis of the CFC inventory using the constrained model suggests that the North Atlantic Ocean shallower than about 2000 m was near 20% saturated in the mid-1990s. Overall, this basin is a sink to 22% of the total atmosphere-to-ocean CFC-11 flux-twice the global average value. The average water mass formation rates over the CFC transient are 7.0 and 6.0 Sv (Sv = 10(6) m(3) s(-1)) for subtropical mode water and subpolar mode water, respectively.

Impact of the Atlantic Ocean on the multidecadal fluctuation of El Nino–Southern Oscillation–South Asian monsoon relationship in a coupled general circulation model

The multidecadal variability of El Niño–Southern Oscillation (ENSO)–South Asian monsoon relationship is elucidated in a 1000 year control simulation of a coupled general circulation model. The results indicate that the Atlantic Multidecadal Oscillation (AMO), resulting from the natural fluctuation of the Atlantic Meridional Overturning Circulation (AMOC), plays an important role in modulating the multidecadal variation of the ENSO-monsoon relationship. The sea surface temperature anomalies associated with the AMO induce not only significant climate impact in the Atlantic but also the coupled feedbacks in the tropical Pacific regions. The remote responses in the Pacific Ocean to a positive phase of the AMO which is resulted from enhanced AMOC in the model simulation and are characterized by statistically significant warming in the North Pacific and in the western tropical Pacific, a relaxation of tropical easterly trades in the central and eastern tropical Pacific, and a deeper thermocline in the eastern tropical Pacific. These changes in mean states lead to a reduction of ENSO variability and therefore a weakening of the ENSO-monsoon relationship. This study suggests a nonlocal mechanism for the low-frequency fluctuation of the ENSO-monsoon relationship, although the AMO explains only a fraction of the ENSO–South Asian monsoon variation on decadal-multidecadal timescale. Given the multidecadal variation of the AMOC and therefore of the AMO exhibit decadal predictability, this study highlights the possibility that a part of the change of climate variability in the Pacific Ocean and its teleconnection may be predictable.

A simple relationship between volcanic sulfate aerosol optical depth and surface temperature change simulated in an atmosphere-ocean general circulation model.

In this study we quantify the relationship between the aerosol optical depth increase from a volcanic eruption and the severity of the subsequent surface temperature decrease. This investigation is made by simulating 10 different sizes of eruption in a global circulation model (GCM) by changing stratospheric sulfate aerosol optical depth at each time step. The sizes of the simulated eruptions range from Pinatubo‐sized up to the magnitude of supervolcanic eruptions around 100 times the size of Pinatubo. From these simulations we find that there is a smooth monotonic relationship between the global mean maximum aerosol optical depth anomaly and the global mean temperature anomaly and we derive a simple mathematical expression which fits this relationship well. We also construct similar relationships between global mean aerosol optical depth and the temperature anomaly at every individual model grid box to produce global maps of best‐fit coefficients and fit residuals. These maps are used with caution to find the eruption size at which a local temperature anomaly is clearly distinct from the local natural variability and to approximate the temperature anomalies which the model may simulate following a Tambora‐sized eruption. To our knowledge, this is the first study which quantifies the relationship between aerosol optical depth and resulting temperature anomalies in a simple way, using the wealth of data that is available from GCM simulations.

Classifying general nonlinear force laws in cell-based models via the continuum limit

though discrete cell-based frameworks are now commonly used to simulate a whole range of biological phenomena, it is typically not obvious how the numerous different types of model are related to one another, nor which one is most appropriate in a given context. Here we demonstrate how individual cell movement on the discrete scale modeled using nonlinear force laws can be described by nonlinear diffusion coefficients on the continuum scale. A general relationship between nonlinear force laws and their respective diffusion coefficients is derived in one spatial dimension and, subsequently, a range of particular examples is considered. For each case excellent agreement is observed between numerical solutions of the discrete and corresponding continuum models. Three case studies are considered in which we demonstrate how the derived nonlinear diffusion coefficients can be used to (a) relate different discrete models of cell behavior; (b) derive discrete, intercell force laws from previously posed diffusion coefficients, and (c) describe aggregative behavior in discrete simulations.

Pungent general anesthetics activate transient receptor potential-A1 to produce hyperalgesia and neurogenic bronchoconstriction

BACKGROUND: Volatile anesthetics such as isoflurane and halothane have been in clinical use for many years and represent the group of drugs most commonly used to maintain general anesthesia. However, despite their widespread use, the molecular mechanisms by which these drugs exert their effects are not completely understood. Recently, a seemingly paradoxical effect of general anesthetics has been identified: the activation of peripheral nociceptors by irritant anesthetics. This mechanism may explain the hyperalgesic actions of inhaled anesthetics and their adverse effects in the airways. METHODS: To test the hypothesis that irritant inhaled anesthetics activate the excitatory ion-channel transient receptor potential (TRP)-A1 and thereby contribute to hyperalgesia and irritant airway effects, we used the measurement of intracellular calcium concentration in isolated cells in culture. For our functional experiments, we used models of isolated guinea pig bronchi to measure bronchoconstriction and withdrawal threshold to mechanical stimulation with von Frey filaments in mice. RESULTS: Irritant inhaled anesthetics activate TRPA1 expressed in human embryonic kidney cells and in nociceptive neurons. Isoflurane induces mechanical hyperalgesia in mice by a TRPA1-dependent mechanism. Isoflurane also induces TRPA1-dependent constriction of isolated bronchi. Nonirritant anesthetics do not activate TRPA1 and fail to produce hyperalgesia and bronchial constriction. CONCLUSIONS: General anesthetics induce a reversible loss of consciousness and render the patient unresponsive to painful stimuli. However, they also produce excitatory effects such as airway irritation and they contribute to postoperative pain. Activation of TRPA1 may contribute to these adverse effects, a hypothesis that remains to be tested in the clinical setting.

Influence of the quasi-biennial oscillation on the extratropical winter stratosphere in an atmospheric general circulation model and in reanalysis data

The interannual variability of the stratospheric polar vortex during winter in both hemispheres is observed to correlate strongly with the phase of the quasi-biennial oscillation (QBO) in tropical stratospheric winds. It follows that the lack of a spontaneously generated QBO in most atmospheric general circulation models (AGCMs) adversely affects the nature of polar variability in such models. This study examines QBO–vortex coupling in an AGCM in which a QBO is spontaneously induced by resolved and parameterized waves. The QBO–vortex coupling in the AGCM compares favorably to that seen in reanalysis data [from the 40-yr ECMWF Re-Analysis (ERA-40)], provided that careful attention is given to the definition of QBO phase. A phase angle representation of the QBO is employed that is based on the two leading empirical orthogonal functions of equatorial zonal wind vertical profiles. This yields a QBO phase that serves as a proxy for the vertical structure of equatorial winds over the whole depth of the stratosphere and thus provides a means of subsampling the data to select QBO phases with similar vertical profiles of equatorial zonal wind. Using this subsampling, it is found that the QBO phase that induces the strongest polar vortex response in early winter differs from that which induces the strongest late-winter vortex response. This is true in both hemispheres and for both the AGCM and ERA-40. It follows that the strength and timing of QBO influence on the vortex may be affected by the partial seasonal synchronization of QBO phase transitions that occurs both in observations and in the model. This provides a mechanism by which changes in the strength of QBO–vortex correlations may exhibit variability on decadal time scales. In the model, such behavior occurs in the absence of external forcings or interannual variations in sea surface temperatures.

Potential Predictability and AMIP Implications of Midlatitude Climate Variability in Two General Circulation Models

Ensembles of extended Atmospheric Model Intercomparison Project (AMIP) runs from the general circulation models of the National Centers for Environmental Prediction (formerly the National Meteorological Center) and the Max-Planck Institute (Hamburg, Germany) are used to estimate the potential predictability (PP) of an index of the Pacific–North America (PNA) mode of climate change. The PP of this pattern in “perfect” prediction experiments is 20%–25% of the index’s variance. The models, particularly that from MPI, capture virtually all of this variance in their hindcasts of the winter PNA for the period 1970–93. The high levels of internally generated model noise in the PNA simulations reconfirm the need for an ensemble averaging approach to climate prediction. This means that the forecasts ought to be expressed in a probabilistic manner. It is shown that the models’ skills are higher by about 50% during strong SST events in the tropical Pacific, so the probabilistic forecasts need to be conditional on the tropical SST. Taken together with earlier studies, the present results suggest that the original set of AMIP integrations (single 10-yr runs) is not adequate to reliably test the participating models’ simulations of interannual climate variability in the midlatitudes.

Toward monitoring the tropospheric temperature by means of a general circulation model

The recent global tropospheric temperature trend can be reproduced by climate models that are forced only by observed sea surface temperature (SST) anomalies. In this study, simulations with the Hamburg climate model (ECHAM) are compared to temperatures from microwave sounding units (MSU) and to reanalyses from the European Centre for Medium-Range Weather Forecasts. There is overall agreement of observed and simulated tropospheric temperature anomalies in many regions, in particular in the tropics and over the oceans, which lack conventional observing systems. This provides the opportunity to link physically different quantities, such as surface observations or analyses (SST) and satellite soundings (MSU) by means of a general circulation model. The proposed method can indicate inconsistencies between MSU temperatures and SSTs and has apparently done so. Differences between observed and simulated tropospheric temperature anomalies can partly be attributed to stratospheric aerosol variations due to major volcanic eruptions.

Correlations of long-lived chemical species in a middle atmosphere general circulation model

Correlations between various chemical species simulated by the Canadian Middle Atmosphere Model, a general circulation model with fully interactive chemistry, are considered in order to investigate the general conditions under which compact correlations can be expected to form. At the same time, the analysis serves to validate the model. The results are compared to previous work on this subject, both from theoretical studies and from atmospheric measurements made from space and from aircraft. The results highlight the importance of having a data set with good spatial coverage when working with correlations and provide a background against which the compactness of correlations obtained from atmospheric measurements can be confirmed. It is shown that for long-lived species, distinct correlations are found in the model in the tropics, the extratropics, and the Antarctic winter vortex. Under these conditions, sparse sampling such as arises from occultation instruments is nevertheless suitable to define a chemical correlation within each region even from a single day of measurements, provided a sufficient range of mixing ratio values is sampled. In practice, this means a large vertical extent, though the requirements are less stringent at more poleward latitudes.

Understanding the transition to seizure by modeling the epileptiform activity of general anesthetic agents

A central difficulty in modeling epileptogenesis using biologically plausible computational and mathematical models is not the production of activity characteristic of a seizure, but rather producing it in response to specific and quantifiable physiologic change or pathologic abnormality. This is particularly problematic when it is considered that the pathophysiological genesis of most epilepsies is largely unknown. However, several volatile general anesthetic agents, whose principle targets of action are quantifiably well characterized, are also known to be proconvulsant. The authors describe recent approaches to theoretically describing the electroencephalographic effects of volatile general anesthetic agents that may be able to provide important insights into the physiologic mechanisms that underpin seizure initiation.

The contributions of domain-general and numerical factors to third-grade arithmetic skills and mathematical learning disability

Explanations of the marked individual differences in elementary school mathematical achievement and mathematical learning disability (MLD or dyscalculia) have involved domain-general factors (working memory, reasoning, processing speed and oral language) and numerical factors that include single-digit processing efficiency and multi-digit skills such as number system knowledge and estimation. This study of third graders (N = 258) finds both domain-general and numerical factors contribute independently to explaining variation in three significant arithmetic skills: basic calculation fluency, written multi-digit computation, and arithmetic word problems. Estimation accuracy and number system knowledge show the strongest associations with every skill and their contributions are both independent of each other and other factors. Different domain-general factors independently account for variation in each skill. Numeral comparison, a single digit processing skill, uniquely accounts for variation in basic calculation. Subsamples of children with MLD (at or below 10th percentile, n = 29) are compared with low achievement (LA, 11th to 25th percentiles, n = 42) and typical achievement (above 25th percentile, n = 187). Examination of these and subsets with persistent difficulties supports a multiple deficits view of number difficulties: most children with number difficulties exhibit deficits in both domain-general and numerical factors. The only factor deficit common to all persistent MLD children is in multi-digit skills. These findings indicate that many factors matter but multi-digit skills matter most in third grade mathematical achievement.

Response of tropical cyclones to idealized climate change experiments in a global high resolution coupled general circulation model

We present an assessment of how tropical cyclone activity might change due to the influence of increased atmospheric carbon dioxide concentrations, using the UK’s High Resolution Global Environment Model (HiGEM) with N144 resolution (~90 km in the atmosphere and ~40 km in the ocean). Tropical cyclones are identified using a feature tracking algorithm applied to model output. Tropical cyclones from idealized 30-year 2×CO2 (2CO2) and 4×CO2 (4CO2) simulations are compared to those identified in a 150-year present-day simulation, which is separated into a 5-member ensemble of 30-year integrations. Tropical cyclones are shown to decrease in frequency globally by 9% in the 2CO2 and 26% in the 4CO2. Tropical cyclones only become more intese in the 4CO2, however uncoupled time slice experiments reveal an increase in intensity in the 2CO2. An investigation into the large-scale environmental conditions, known to influence tropical cyclone activity in the main development regions, is used to determine the response of tropical cyclone activity to increased atmospheric CO2. A weaker Walker circulation and a reduction in zonally averaged regions of updrafts lead to a shift in the location of tropical cyclones in the northern hemisphere. A decrease in mean ascent at 500 hPa contributes to the reduction of tropical cyclones in the 2CO2 in most basins. The larger reduction of tropical cyclones in the 4CO2 arises from further reduction of mean ascent at 500 hPa and a large enhancement of vertical wind shear, especially in the southern hemisphere, North Atlantic and North East Pacific.

General circulation model estimates of aerosol transport and radiative forcing during the Indian Ocean Experiment

Aerosol sources, transport, and sinks are simulated, and aerosol direct radiative effects are assessed over the Indian Ocean for the Indian Ocean Experiment (INDOEX) Intensive Field Phase during January to March 1999 using the Laboratoire de Me´te´orologie Dynamique (LMDZT) general circulation model. The model reproduces the latitudinal gradient in aerosol mass concentration and optical depth (AOD). The model-predicted aerosol concentrations and AODs agree reasonably well with measurements but are systematically underestimated during high-pollution episodes, especially in the month of March. The largest aerosol loads are found over southwestern China, the Bay of Bengal, and the Indian subcontinent. Aerosol emissions from the Indian subcontinent are transported into the Indian Ocean through either the west coast or the east coast of India. Over the INDOEX region, carbonaceous aerosols are the largest contributor to the estimated AOD, followed by sulfate, dust, sea salt, and fly ash. During the northeast winter monsoon, natural and anthropogenic aerosols reduce the solar flux reaching the surface by 25 W m�2, leading to 10–15% less insolation at the surface. A doubling of black carbon (BC) emissions from Asia results in an aerosol single-scattering albedo that is much smaller than in situ measurements, reflecting the fact that BC emissions are not underestimated in proportion to other (mostly scattering) aerosol types. South Asia is the dominant contributor to sulfate aerosols over the INDOEX region and accounts for 60–70% of the AOD by sulfate. It is also an important but not the dominant contributor to carbonaceous aerosols over the INDOEX region with a contribution of less than 40% to the AOD by this aerosol species. The presence of elevated plumes brings significant quantities of aerosols to the Indian Ocean that are generated over Africa and Southeast and east Asia.