Modeling oxygen isotopes in the Pliocene: Large-scale features over the land and ocean

The first isotope-enabled general circulation model (GCM) simulations of the Pliocene are used to discuss the interpretation of δ18O measurements for a warm climate. The model suggests that spatial patterns of Pliocene ocean surface δ18O ( inline image) were similar to those of the preindustrial period; however, Arctic and coastal regions were relatively depleted, while South Atlantic and Mediterranean regions were relatively enriched. Modeled inline image anomalies are closely related to modeled salinity anomalies, which supports using inline image as a paleosalinity proxy. Modeled Pliocene precipitation δ18O ( inline image) was enriched relative to the preindustrial values (but with depletion of <2‰ over some tropical regions). While usually modest (<4‰), the enrichment can reach 25‰ over ice sheet regions. In the tropics inline image anomalies are related to precipitation amount anomalies, although there is usually a spatial offset between the two. This offset suggests that the location of precipitation change is more uncertain than the amplitude when interpreting inline image. At high latitudes inline image anomalies relate to temperature anomalies; however, the relationship is neither linear nor spatially coincident: a large inline image signal does not always translate to a large temperature signal. These results suggest that isotope modeling can lead to enhanced synergy between climate models and climate proxy data. The model can relate proxy data to climate in a physically based way even when the relationship is complex and nonlocal. The δ18O-climate relationships, identified here from a GCM, could not be determined from transfer functions or simple models

Can We Detect Changes in Amazon Forest Structure Using Measurements of the Isotopic Composition of Precipitation?

Large‐scale (>500 km) spatial gradients of precipitation oxygen isotope‐ratios (δ18Op) hold information about the hydrological cycle. They result from the interplay between rainout and evapotranspiration along air‐parcel paths, but these counteracting effects are difficult to disentangle complicating quantification of the effect of land cover change on δ18Op. We show that disentangling can qualitatively be achieved using climate model simulations with a land‐derived precipitation tracer for tropical South America. We then either vary land cover as observed since 1870 or by replacing Amazon forests with bare land to determine the resulting signals. Our results indicate that effects of historically changing land cover on annual mean δ18O isotope‐ratio gradients are small and unlikely detectable, although there is a noticeable signal during the dry season. Furthermore, the effect of changes in water recycling on Amazon δ18Op in paleo‐records may have been overestimated and need reinterpretation

The HadCM3 contribution to PlioMIP phase 2

We present the UK's input into the Pliocene Model Intercomparison Project phase 2 (PlioMIP2) using the Hadley Centre Climate Model version 3 (HadCM3). The 400 ppm CO2 Pliocene experiment has a mean annual surface air temperature that is 2.9 ∘C warmer than the pre-industrial and a polar amplification of between 1.7 and 2.2 times the global mean warming. The Pliocene Research Interpretation and Synoptic Mapping (PRISM4) enhanced Pliocene palaeogeography accounts for a warming of 1.4 ∘C, whilst the CO2 increase from 280 to 400 ppm leads to a further 1.5 ∘C of warming. Climate sensitivity is 3.5 ∘C for the pre-industrial and 2.9 ∘C for the Pliocene. Precipitation change between the pre-industrial and Pliocene is complex, with geographic and land surface changes primarily modifying the geographical extent of mean annual precipitation. Sea ice fraction and areal extent are reduced during the Pliocene, particularly in the Southern Hemisphere, although they persist through summer in both hemispheres. The Pliocene palaeogeography drives a more intense Pacific and Atlantic meridional overturning circulation (AMOC). This intensification of AMOC is coincident with more widespread deep convection in the North Atlantic. We conclude by examining additional sensitivity experiments and confirm that the choice of total solar insolation (1361 vs. 1365 Wm−2) and orbital configuration (modern vs. 3.205 Ma) does not significantly influence the anomaly-type analysis in use by the Pliocene community

Simulating stable carbon isotopes in the ocean component of the FAMOUS general circulation model with MOSES1 (XOAVI)

Ocean circulation and the marine carbon cycle can be indirectly inferred from stable and radiogenic carbon isotope ratios (δ13C and Δ14C, respectively), measured directly in the water column, or recorded in geological archives such as sedimentary microfossils and corals. However, interpreting these records is non-trivial because they reflect a complex interplay between physical and biogeochemical processes. By directly simulating multiple isotopic tracer fields within numerical models, we can improve our understanding of the processes that control large-scale isotope distributions and interpolate the spatiotemporal gaps in both modern and palaeo datasets. We have added the stable isotope 13C to the ocean component of the FAMOUS coupled atmosphere–ocean general circulation model, which is a valuable tool for simulating complex feedbacks between different Earth system processes on decadal to multi-millennial timescales. We tested three different biological fractionation parameterisations to account for the uncertainty associated with equilibrium fractionation during photosynthesis and used sensitivity experiments to quantify the effects of fractionation during air–sea gas exchange and primary productivity on the simulated δ13CDIC distributions. Following a 10 000-year pre-industrial spin-up, we simulated the Suess effect (the isotopic imprint of anthropogenic fossil fuel burning) to assess the performance of the model in replicating modern observations. Our implementation captures the large-scale structure and range of δ13CDIC observations in the surface ocean, but the simulated values are too high at all depths, which we infer is due to biases in the biological pump. In the first instance, the new 13C tracer will therefore be useful for recalibrating both the physical and biogeochemical components of FAMOUS

Asian monsoons and aridification response to Paleogene sea retreat and Neogene westerly shielding indicated by seasonality in Paratethys oysters

Asian climate patterns, characterised by highly seasonal monsoons and continentality, are thought to originate in the Eocene epoch (56 to 34 million years ago – Ma) in response to global climate, Tibetan Plateau uplift and the disappearance of the giant Proto-Paratethys sea formerly extending over Eurasia. The influence of this sea on Asian climate has hitherto not been constrained by proxy records despite being recognised as a major driver by climate models. We report here strongly seasonal records preserved in annual lamina of Eocene oysters from the Proto-Paratethys with sedimentological and numerical data showing that monsoons were not dampened by the sea and that aridification was modulated by westerly moisture sourced from the sea. Hot and arid summers despite the presence of the sea suggest a strong anticyclonic zone at Central Asian latitudes and an orographic effect from the emerging Tibetan Plateau. Westerly moisture precipitating during cold and wetter winters appear to have decreased in two steps. First in response to the late Eocene (34–37 Ma) sea retreat; second by the orogeny of the Tian Shan and Pamir ranges shielding the westerlies after 25 Ma. Paleogene sea retreat and Neogene westerly shielding thus provide two successive mechanisms forcing coeval Asian desertification and biotic crises

Cell killing and resistance in pre-operative breast cancer chemotherapy

<p>Abstract</p> <p>Background</p> <p>Despite the recent development of technologies giving detailed images of tumours <it>in vivo</it>, direct or indirect ways to measure how many cells are actually killed by a treatment or are resistant to it are still beyond our reach.</p> <p>Methods</p> <p>We designed a simple model of tumour progression during treatment, based on descriptions of the key phenomena of proliferation, quiescence, cell killing and resistance, and giving as output the macroscopically measurable tumour volume and growth fraction. The model was applied to a database of the time course of volumes of breast cancer in patients undergoing pre-operative chemotherapy, for which the initial estimate of proliferating cells by the measure of the percentage of Ki67-positive cells was available.</p> <p>Results</p> <p>The analysis recognises different patterns of response to treatment. In one subgroup of patients the fitting implied drug resistance. In another subgroup there was a shift to higher sensitivity during the therapy. In the subgroup of patients where killing of cycling cells had the highest score, the drugs showed variable efficacy against quiescent cells.</p> <p>Conclusion</p> <p>The approach was feasible, providing items of information not otherwise available. Additional data, particularly sequential Ki67 measures, could be added to the system, potentially reducing uncertainty in estimates of parameter values.</p

Antarctic last interglacial isotope peak in response to sea ice retreat not ice-sheet collapse

Several studies have suggested that sea-level rise during the last interglacial implies retreat of the West Antarctic Ice Sheet (WAIS). The prevalent hypothesis is that the retreat coincided with the peak Antarctic temperature and stable water isotope values from 128,000 years ago (128 ka); very early in the last interglacial. Here, by analysing climate model simulations of last interglacial WAIS loss featuring water isotopes, we show instead that the isotopic response to WAIS loss is in opposition to the isotopic evidence at 128 ka. Instead, a reduction in winter sea ice area of 65±7% fully explains the 128 ka ice core evidence. Our finding of a marked retreat of the sea ice at 128 ka demonstrates the sensitivity of Antarctic sea ice extent to climate warming

A combined model reduction algorithm for controlled biochemical systems

Background: Systems Biology continues to produce increasingly large models of complex biochemical reaction networks. In applications requiring, for example, parameter estimation, the use of agent-based modelling approaches, or real-time simulation, this growing model complexity can present a significant hurdle. Often, however, not all portions of a model are of equal interest in a given setting. In such situations methods of model reduction offer one possible approach for addressing the issue of complexity by seeking to eliminate those portions of a pathway that can be shown to have the least effect upon the properties of interest. Methods: In this paper a model reduction algorithm bringing together the complementary aspects of proper lumping and empirical balanced truncation is presented. Additional contributions include the development of a criterion for the selection of state-variable elimination via conservation analysis and use of an ‘averaged’ lumping inverse. This combined algorithm is highly automatable and of particular applicability in the context of ‘controlled’ biochemical networks. Results: The algorithm is demonstrated here via application to two examples; an 11 dimensional model of bacterial chemotaxis in Escherichia coli and a 99 dimensional model of extracellular regulatory kinase activation (ERK) mediated via the epidermal growth factor (EGF) and nerve growth factor (NGF) receptor pathways. In the case of the chemotaxis model the algorithm was able to reduce the model to 2 state-variables producing a maximal relative error between the dynamics of the original and reduced models of only 2.8% whilst yielding a 26 fold speed up in simulation time. For the ERK activation model the algorithm was able to reduce the system to 7 state-variables, incurring a maximal relative error of 4.8%, and producing an approximately 10 fold speed up in the rate of simulation. Indices of controllability and observability are additionally developed and demonstrated throughout the paper. These provide insight into the relative importance of individual reactants in mediating a biochemical system’s input-output response even for highly complex networks. Conclusions: Through application, this paper demonstrates that combined model reduction methods can produce a significant simplification of complex Systems Biology models whilst retaining a high degree of predictive accuracy. In particular, it is shown that by combining the methods of proper lumping and empirical balanced truncation it is often possible to produce more accurate reductions than can be obtained by the use of either method in isolation

Evaluation of Arctic warming in mid-Pliocene climate simulations

Palaeoclimate simulations improve our understanding of the climate, inform us about the performance of climate models in a different climate scenario, and help to identify robust features of the climate system. Here, we analyse Arctic warming in an ensemble of 16 simulations of the mid-Pliocene Warm Period (mPWP), derived from the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2). The PlioMIP2 ensemble simulates Arctic (60–90∘ N) annual mean surface air temperature (SAT) increases of 3.7 to 11.6 ∘C compared to the pre-industrial period, with a multi-model mean (MMM) increase of 7.2 ∘C. The Arctic warming amplification ratio relative to global SAT anomalies in the ensemble ranges from 1.8 to 3.1 (MMM is 2.3). Sea ice extent anomalies range from −3.0 to −10.4×10^{6} km^{2}, with a MMM anomaly of −5.6×10^{6} km^{2}, which constitutes a decrease of 53 % compared to the pre-industrial period. The majority (11 out of 16) of models simulate summer sea-ice-free conditions (≤1×10^{6} km^{2}) in their mPWP simulation. The ensemble tends to underestimate SAT in the Arctic when compared to available reconstructions, although the degree of underestimation varies strongly between the simulations. The simulations with the highest Arctic SAT anomalies tend to match the proxy dataset in its current form better. The ensemble shows some agreement with reconstructions of sea ice, particularly with regard to seasonal sea ice. Large uncertainties limit the confidence that can be placed in the findings and the compatibility of the different proxy datasets. We show that while reducing uncertainties in the reconstructions could decrease the SAT data–model discord substantially, further improvements are likely to be found in enhanced boundary conditions or model physics. Lastly, we compare the Arctic warming in the mPWP to projections of future Arctic warming and find that the PlioMIP2 ensemble simulates greater Arctic amplification than CMIP5 future climate simulations and an increase instead of a decrease in Atlantic Meridional Overturning Circulation (AMOC) strength compared to pre-industrial period. The results highlight the importance of slow feedbacks in equilibrium climate simulations, and that caution must be taken when using simulations of the mPWP as an analogue for future climate change

A return to large-scale features of Pliocene climate: the Pliocene Model Intercomparison Project Phase 2

The Pliocene epoch has great potential to improve our understanding of the long-term climatic and environmental consequences of an atmospheric CO2 concentration near ~ 400 parts per million by volume. Here we present the large-scale features of Pliocene climate as simulated by a new ensemble of climate models of varying complexity and spatial resolution and based on new reconstructions of boundary conditions (the Pliocene Model Intercomparison Project Phase 2; PlioMIP2). As a global annual average, modelled surface air temperatures increase by between 1.4 and 4.7 °C relative to pre-industrial with a multi-model mean value of 2.8 °C. Annual mean total precipitation rates increase by 6 % (range: 2 %–13 %). On average, surface air temperature (SAT) increases are 1.3 °C greater over the land than over the oceans, and there is a clear pattern of polar amplification with warming polewards of 60° N and 60° S exceeding the global mean warming by a factor of 2.4. In the Atlantic and Pacific Oceans, meridional temperature gradients are reduced, while tropical zonal gradients remain largely unchanged. Although there are some modelling constraints, there is a statistically significant relationship between a model's climate response associated with a doubling in CO2 (Equilibrium Climate Sensitivity; ECS) and its simulated Pliocene surface temperature response. The mean ensemble earth system response to doubling of CO2 (including ice sheet feedbacks) is approximately 50 % greater than ECS, consistent with results from the PlioMIP1 ensemble. Proxy-derived estimates of Pliocene sea-surface temperatures are used to assess model estimates of ECS and indicate a range in ECS from 2.5 to 4.3 °C. This result is in general accord with the range in ECS presented by previous IPCC Assessment Reports