Dynamical and observational constraints on tropical Pacific sea surface temperatures at the last glacial maximum

Asynchronously coupled atmosphere and ocean general circulation model simulations are used to examine the consequences of changes in the west/east sea-surface temperature (SST) gradient across the equatorial Pacific at the last glacial maximum (LGM). Simulations forced by the CLIMAP SST for the LGM, where the west/east SST gradient across the Pacific is reduced compared to present, produce a reduction in the strength of the trade winds and a decrease in the west/east slope of the equatorial thermocline that is incompatible with thermocline depths newly inferred from foraminiferal assemblages. Stronger-than-present trade winds, and a more realistic simulation of the thermocline slope, are produced when eastern Pacific SSTs are 2°C cooler than western Pacific SSTs. Our study highlights the importance of spatial heterogeneity in tropical SSTs in determining key features of the glacial climate.

Tropical palaeoclimates at the Last Glacial Maximum: comparison of Palaeoclimate Modeling Intercomparison Project (PMIP)simulations and paleodata.

Seventeen simulations of the Last Glacial Maximum (LGM) climate have been performed using atmospheric general circulation models (AGCM) in the framework of the Paleoclimate Modeling Intercomparison Project (PMIP). These simulations use the boundary conditions for CO2, insolation and ice-sheets; surface temperatures (SSTs) are either (a) prescribed using CLIMAP data set (eight models) or (b) computed by coupling the AGCM with a slab ocean (nine models). The present-day (PD) tropical climate is correctly depicted by all the models, except the coarser resolution models, and the simulated geographical distribution of annual mean temperature is in good agreement with climatology. Tropical cooling at the LGM is less than at middle and high latitudes, but greatly exceeds the PD temperature variability. The LGM simulations with prescribed SSTs underestimate the observed temperature changes except over equatorial Africa where the models produce a temperature decrease consistent with the data. Our results confirm previous analyses showing that CLIMAP (1981) SSTs only produce a weak terrestrial cooling. When SSTs are computed, the models depict a cooling over the Pacific and Indian oceans in contrast with CLIMAP and most models produce cooler temperatures over land. Moreover four of the nine simulations, produce a cooling in good agreement with terrestrial data. Two of these model results over ocean are consistent with new SST reconstructions whereas two models simulate a homogeneous cooling. Finally, the LGM aridity inferred for most of the tropics from the data, is globally reproduced by the models with a strong underestimation for models using computed SSTs.

Tropical climates at the Last Glacial Maximum: a new synthesis of terrestrial palaeoclimate data. I. Vegetation, lake-levels and geochemistry.

Palaeodata in synthesis form are needed as benchmarks for the Palaeoclimate Modelling Intercomparison Project (PMIP). Advances since the last synthesis of terrestrial palaeodata from the last glacial maximum (LGM) call for a new evaluation, especially of data from the tropics. Here pollen, plant-macrofossil, lake-level, noble gas (from groundwater) and δ18O (from speleothems) data are compiled for 18±2 ka (14C), 32 °N–33 °S. The reliability of the data was evaluated using explicit criteria and some types of data were re-analysed using consistent methods in order to derive a set of mutually consistent palaeoclimate estimates of mean temperature of the coldest month (MTCO), mean annual temperature (MAT), plant available moisture (PAM) and runoff (P-E). Cold-month temperature (MAT) anomalies from plant data range from −1 to −2 K near sea level in Indonesia and the S Pacific, through −6 to −8 K at many high-elevation sites to −8 to −15 K in S China and the SE USA. MAT anomalies from groundwater or speleothems seem more uniform (−4 to −6 K), but the data are as yet sparse; a clear divergence between MAT and cold-month estimates from the same region is seen only in the SE USA, where cold-air advection is expected to have enhanced cooling in winter. Regression of all cold-month anomalies against site elevation yielded an estimated average cooling of −2.5 to −3 K at modern sea level, increasing to ≈−6 K by 3000 m. However, Neotropical sites showed larger than the average sea-level cooling (−5 to −6 K) and a non-significant elevation effect, whereas W and S Pacific sites showed much less sea-level cooling (−1 K) and a stronger elevation effect. These findings support the inference that tropical sea-surface temperatures (SSTs) were lower than the CLIMAP estimates, but they limit the plausible average tropical sea-surface cooling, and they support the existence of CLIMAP-like geographic patterns in SST anomalies. Trends of PAM and lake levels indicate wet LGM conditions in the W USA, and at the highest elevations, with generally dry conditions elsewhere. These results suggest a colder-than-present ocean surface producing a weaker hydrological cycle, more arid continents, and arguably steeper-than-present terrestrial lapse rates. Such linkages are supported by recent observations on freezing-level height and tropical SSTs; moreover, simulations of “greenhouse” and LGM climates point to several possible feedback processes by which low-level temperature anomalies might be amplified aloft.

Dust sources and deposition during the last glacial maximum and current climate: a comparison of model results with paleodata from ice cores and marine sediments

Mineral dust aerosols in the atmosphere have the potential to affect the global climate by influencing the radiative balance of the atmosphere and the supply of micronutrients to the ocean. Ice and marine sediment cores indicate that dust deposition from the atmosphere was at some locations 2–20 times greater during glacial periods, raising the possibility that mineral aerosols might have contributed to climate change on glacial-interglacial time scales. To address this question, we have used linked terrestrial biosphere, dust source, and atmospheric transport models to simulate the dust cycle in the atmosphere for current and last glacial maximum (LGM) climates. We obtain a 2.5-fold higher dust loading in the entire atmosphere and a twenty-fold higher loading in high latitudes, in LGM relative to present. Comparisons to a compilation of atmospheric dust deposition flux estimates for LGM and present in marine sediment and ice cores show that the simulated flux ratios are broadly in agreement with observations; differences suggest where further improvements in the simple dust model could be made. The simulated increase in high-latitude dustiness depends on the expansion of unvegetated areas, especially in the high latitudes and in central Asia, caused by a combination of increased aridity and low atmospheric [CO2]. The existence of these dust source areas at the LGM is supported by pollen data and loess distribution in the northern continents. These results point to a role for vegetation feedbacks, including climate effects and physiological effects of low [CO2], in modulating the atmospheric distribution of dust.

Adaptive nonlinear equalizer using a mixture of gaussians based on-line density estimator

This paper introduces a new adaptive nonlinear equalizer relying on a radial basis function (RBF) model, which is designed based on the minimum bit error rate (MBER) criterion, in the system setting of the intersymbol interference channel plus a co-channel interference. Our proposed algorithm is referred to as the on-line mixture of Gaussians estimator aided MBER (OMG-MBER) equalizer. Specifically, a mixture of Gaussians based probability density function (PDF) estimator is used to model the PDF of the decision variable, for which a novel on-line PDF update algorithm is derived to track the incoming data. With the aid of this novel on-line mixture of Gaussians based sample-by-sample updated PDF estimator, our adaptive nonlinear equalizer is capable of updating its equalizer’s parameters sample by sample to aim directly at minimizing the RBF nonlinear equalizer’s achievable bit error rate (BER). The proposed OMG-MBER equalizer significantly outperforms the existing on-line nonlinear MBER equalizer, known as the least bit error rate equalizer, in terms of both the convergence speed and the achievable BER, as is confirmed in our simulation study

Simulated last glacial maximum D14Catm and the deep glacial ocean carbon reservoir

∆14Catm has been estimated as 420 ± 80‰ (IntCal09) during the Last Glacial Maximum (LGM) compared to preindustrial times (0‰), but mechanisms explaining this difference are not yet resolved. ∆14Catm is a function of both cosmogenic production in the high atmosphere and of carbon cycling and partitioning in the Earth system. 10Be-based reconstructions show a contribution of the cosmogenic production term of only 200 ± 200‰ in the LGM. The remaining 220‰ have thus to be explained by changes in the carbon cycle. Recently, Bouttes et al. (2010, 2011) proposed to explain most of the difference in pCO2atm and δ13C between glacial and interglacial times as a result of brine-induced ocean stratification in the Southern Ocean. This mechanism involves the formation of very saline water masses that contribute to high carbon storage in the deep ocean. During glacial times, the sinking of brines is enhanced and more carbon is stored in the deep ocean, lowering pCO2atm. Moreover, the sinking of brines induces increased stratification in the Southern Ocean, which keeps the deep ocean well isolated from the surface. Such an isolated ocean reservoir would be characterized by a low ∆14C signature. Evidence of such 14C-depleted deep waters during the LGM has recently been found in the Southern Ocean (Skinner et al. 2010). The degassing of this carbon with low ∆14C would then reduce ∆14Catm throughout the deglaciation. We have further developed the CLIMBER-2 model to include a cosmogenic production of 14C as well as an interactive atmospheric 14C reservoir. We investigate the role of both the sinking of brine and cosmogenic production, alongside iron fertilization mechanisms, to explain changes in ∆14Catm during the last deglaciation. In our simulations, not only is the sinking of brine mechanism consistent with past ∆14C data, but it also explains most of the differences in pCO2atm and ∆14Catm between the LGM and preindustrial times. Finally, this study represents the first time to our knowledge that a model experiment explains glacial-interglacial differences in pCO2atm, δ13C, and ∆14C together with a coherent LGM climate.

Do real balance effects invalidate the Taylor principle in closed and open economies?

This paper examines the determinacy implications of forecast-based monetary policy rules that set the interest rate in response to expected future inflation in a Neo-Wicksellian model that incorporates real balance effects. We show that the presence of such effects in closed economies restricts the ability of the Taylor principle to prevent indeterminacy of the rational expectations equilibrium. The problem is exacerbated in open economies, particularly if the policy rule reacts to consumer-price, rather than domestic-price, inflation. However, determinacy can be restored in both closed and open economies with the addition of monetary policy inertia.

The lifecycle of the North Atlantic storm track

The North Atlantic eddy-driven jet exhibits latitudinal variability, with evidence of three preferred latitudinal locations: south, middle and north. Here we examine the drivers of this variability and the variability of the associated storm track. We investigate the changes in the storm track characteristics for the three jet locations, and propose a mechanism by which enhanced storm track activity, as measured by upstream heat flux, is responsible for cyclical downstream latitudinal shifts in the jet. This mechanism is based on a nonlinear oscillator relationship between the enhanced meridional temperature gradient (and thus baroclinicity) and the meridional high-frequency (periods of shorter than 10 days) eddy heat flux. Such oscillations in baroclinicity and heat flux induce variability in eddy anisotropy which is associated with the changes in the dominant type of wave breaking and a different latitudinal deflection of the jet. Our results suggest that high heat flux is conducive to a northward deflection of the jet, whereas low heat flux is conducive to a more zonal jet. This jet deflecting effect was found to operate most prominently downstream of the storm track maximum, while the storm track and the jet remain anchored at a fixed latitudinal location at the beginning of the storm track. These cyclical changes in storm track characteristics can be viewed as different stages of the storm track’s spatio-temporal lifecycle.

Patterns of maximum body size evolution in Cenozoic land mammals: eco-evolutionary processes and abiotic forcing

There is accumulating evidence that macroevolutionary patterns of mammal evolution during the Cenozoic follow similar trajectories on different continents. This would suggest that such patterns are strongly determined by global abiotic factors, such as climate, or by basic eco-evolutionary processes such as filling of niches by specialization. The similarity of pattern would be expected to extend to the history of individual clades. Here, we investigate the temporal distribution of maximum size observed within individual orders globally and on separate continents. While the maximum size of individual orders of large land mammals show differences and comprise several families, the times at which orders reach their maximum size over time show strong congruence, peaking in the Middle Eocene, the Oligocene and the Plio-Pleistocene. The Eocene peak occurs when global temperature and land mammal diversity are high and is best explained as a result of niche expansion rather than abiotic forcing. Since the Eocene, there is a significant correlation between maximum size frequency and global temperature proxy. The Oligocene peak is not statistically significant and may in part be due to sampling issues. The peak in the Plio-Pleistocene occurs when global temperature and land mammal diversity are low, it is statistically the most robust one and it is best explained by global cooling. We conclude that the macroevolutionary patterns observed are a result of the interplay between eco-evolutionary processes and abiotic forcing

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).

Representation errors and retrievals in linear and nonlinear data assimilation

This article shows how one can formulate the representation problem starting from Bayes’ theorem. The purpose of this article is to raise awareness of the formal solutions,so that approximations can be placed in a proper context. The representation errors appear in the likelihood, and the different possibilities for the representation of reality in model and observations are discussed, including nonlinear representation probability density functions. Specifically, the assumptions needed in the usual procedure to add a representation error covariance to the error covariance of the observations are discussed,and it is shown that, when several sub-grid observations are present, their mean still has a representation error ; socalled ‘superobbing’ does not resolve the issue. Connection is made to the off-line or on-line retrieval problem, providing a new simple proof of the equivalence of assimilating linear retrievals and original observations. Furthermore, it is shown how nonlinear retrievals can be assimilated without loss of information. Finally we discuss how errors in the observation operator model can be treated consistently in the Bayesian framework, connecting to previous work in this area.

Thermospheric Control of the Auroral Source of O+Ions for the Magnetosphere

Linear theory, model ion-density profiles and MSIS neutral thermospheric predictions are used to investigate the stability of the auroral, topside ionosphere to oxygen cyclotron waves: variations of the critical height, above which the plasma is unstable, with field-aligned current, thermal ion density and exospheric temperature are considered. In addition, probabilities are assessed that interactions with neutral atomic gases prevent O+ ions from escaping into the magnetosphere after they have been transversely accelerated by these waves. The two studies are combined to give a rough estimate of the total O+ escape flux as a function of the field-aligned current density for an assumed rise in the perpendicular ion temperature. Charge exchange with neutral oxygen, not hydrogen, is shown to be the principle limitation to the escape of O+ ions, which occurs when the waves are driven unstable down to low altitudes. It is found that the largest observed field-aligned current densities can heat a maximum of about 5×1014 O+ ions m−2 to a threshold above which they are subsequently able to escape into the magnetosphere in the following 500s. Averaged over this period, this would constitute a flux of 1012 m−2 s−1 and in steady-state the peak outflow would then be limited to about 1013 m−2 s−1 by frictional drag on thermal O+ at lower altitudes. Maximum escape is at low plasma density unless the O+ scale height is very large. The outflow decreases with decreasing field-aligned current density and, to a lesser extent, with increasing exospheric temperature. Upward flowing ion events are evaluated as a source of O+ ions for the magnetosphere and as an explanation of the observed solar cycle variation of ring current O+ abundance.

Simple M-factor algorithm for improved estimation of the basic maximum usable frequency of radio waves reflected from the ionospheric F-region

The equations of Milsom are evaluated, giving the ground range and group delay of radio waves propagated via the horizontally stratified model ionosphere proposed by Bradley and Dudeney. Expressions for the ground range which allow for the effects of the underlying E- and F1-regions are used to evaluate the basic maximum usable frequency or M-factors for single F-layer hops. An algorithm for the rapid calculation of the M-factor at a given range is developed, and shown to be accurate to within 5%. The results reveal that the M(3000)F2-factor scaled from vertical-incidence ionograms using the standard URSI procedure can be up to 7.5% in error. A simple addition to the algorithm effects a correction to ionogram values to make these accurate to 0.5%.

Advances in the study of boundary value problems for nonlinear integrable PDEs

In this review I summarise some of the most significant advances of the last decade in the analysis and solution of boundary value problems for integrable partial differential equations in two independent variables. These equations arise widely in mathematical physics, and in order to model realistic applications, it is essential to consider bounded domain and inhomogeneous boundary conditions. I focus specifically on a general and widely applicable approach, usually referred to as the Unified Transform or Fokas Transform, that provides a substantial generalisation of the classical Inverse Scattering Transform. This approach preserves the conceptual efficiency and aesthetic appeal of the more classical transform approaches, but presents a distinctive and important difference. While the Inverse Scattering Transform follows the "separation of variables" philosophy, albeit in a nonlinear setting, the Unified Transform is a based on the idea of synthesis, rather than separation, of variables. I will outline the main ideas in the case of linear evolution equations, and then illustrate their generalisation to certain nonlinear cases of particular significance.

The elliptic sine-Gordon equation: a nonlinear elliptic integrable PDE

In this paper, we summarise this recent progress to underline the features specific to this nonlinear elliptic case, and we give a new classification of boundary conditions on the semistrip that satisfy a necessary condition for yielding a boundary value problem can be effectively linearised. This classification is based on formulation the equation in terms of an alternative Lax pair.