Aortic morphometry at endograft position as assessed by 3D image analysis affects risk of type I endoleak formation after TEVAR

PURPOSE The purpose of this study was to identify morphologic factors affecting type I endoleak formation and bird-beak configuration after thoracic endovascular aortic repair (TEVAR). METHODS Computed tomography (CT) data of 57 patients (40 males; median age, 66 years) undergoing TEVAR for thoracic aortic aneurysm (34 TAA, 19 TAAA) or penetrating aortic ulcer (n = 4) between 2001 and 2010 were retrospectively reviewed. In 28 patients, the Gore TAG® stent-graft was used, followed by the Medtronic Valiant® in 16 cases, the Medtronic Talent® in 8, and the Cook Zenith® in 5 cases. Proximal landing zone (PLZ) was in zone 1 in 13, zone 2 in 13, zone 3 in 23, and zone 4 in 8 patients. In 14 patients (25%), the procedure was urgent or emergent. In each case, pre- and postoperative CT angiography was analyzed using a dedicated image processing workstation and complimentary in-house developed software based on a 3D cylindrical intensity model to calculate aortic arch angulation and conicity of the landing zones (LZ). RESULTS Primary type Ia endoleak rate was 12% (7/57) and subsequent re-intervention rate was 86% (6/7). Left subclavian artery (LSA) coverage (p = 0.036) and conicity of the PLZ (5.9 vs. 2.6 mm; p = 0.016) were significantly associated with an increased type Ia endoleak rate. Bird-beak configuration was observed in 16 patients (28%) and was associated with a smaller radius of the aortic arch curvature (42 vs. 65 mm; p = 0.049). Type Ia endoleak was not associated with a bird-beak configuration (p = 0.388). Primary type Ib endoleak rate was 7% (4/57) and subsequent re-intervention rate was 100%. Conicity of the distal LZ was associated with an increased type Ib endoleak rate (8.3 vs. 2.6 mm; p = 0.038). CONCLUSIONS CT-based 3D aortic morphometry helps to identify risk factors of type I endoleak formation and bird-beak configuration during TEVAR. These factors were LSA coverage and conicity within the landing zones for type I endoleak formation and steep aortic angulation for bird-beak configuration.

Engineering an in vitro air-blood barrier by 3D bioprinting

Intensive efforts in recent years to develop and commercialize in vitro alternatives in the field of risk assessment have yielded new promising two- and three dimensional (3D) cell culture models. Nevertheless, a realistic 3D in vitro alveolar model is not available yet. Here we report on the biofabrication of the human air-blood tissue barrier analogue composed of an endothelial cell, basement membrane and epithelial cell layer by using a bioprinting technology. In contrary to the manual method, we demonstrate that this technique enables automatized and reproducible creation of thinner and more homogeneous cell layers, which is required for an optimal air-blood tissue barrier. This bioprinting platform will offer an excellent tool to engineer an advanced 3D lung model for high-throughput screening for safety assessment and drug efficacy testing.

Gene expression profile in newborn rat lungs after two days of recovery of mechanical ventilation.

BACKGROUND Preterm infants having immature lungs often require respiratory support, potentially leading to bronchopulmonary dysplasia (BPD). Conventional BPD rodent models based on mechanical ventilation (MV) present outcome measured at the end of the ventilation period. A reversible intubation and ventilation model in newborn rats recently allowed discovering that different sets of genes modified their expression related to time after MV. In a newborn rat model, the expression profile 48 h after MV was analyzed with gene arrays to detect potentially interesting candidates with an impact on BPD development. METHODS Rat pups were injected P4-5 with 2 mg/kg lipopolysaccharide (LPS). One day later, MV with 21 or 60% oxygen was applied during 6 h. Animals were sacrified 48 h after end of ventilation. Affymetrix gene arrays assessed the total gene expression profile in lung tissue. RESULTS In fully treated animals (LPS + MV + 60% O(2)) vs. controls, 271 genes changed expression significantly. All modified genes could be classified in six pathways: tissue remodeling/wound repair, immune system and inflammatory response, hematopoiesis, vasodilatation, and oxidative stress. Major alterations were found in the MMP and complement system. CONCLUSION MMPs and complement factors play a central role in several of the pathways identified and may represent interesting targets for BPD treatment/prevention.Bronchopulmonary dysplasia (BPD) is a chronic lung disease occurring in ~30% of preterm infants born less than 30 wk of gestation (1). Its main risk factors include lung immaturity due to preterm delivery, mechanical ventilation (MV), oxygen toxicity, chorioamnionitis, and sepsis. The main feature is an arrest of alveolar and capillary formation (2). Models trying to decipher genes involved in the pathophysiology of BPD are mainly based on MV and oxygen application to young mammals with immature lungs of different species (3). In newborn rodent models, analyses of lung structure and gene and protein expression are performed for practical reasons directly at the end of MV (4,5,6). However, later appearing changes of gene expression might also have an impact on lung development and the evolution towards BPD and cannot be discovered by such models. Recently, we developed a newborn rat model of MV using an atraumatic (orotracheal) intubation technique that allows the weaning of the newborn animal off anesthesia and MV, the extubation to spontaneous breathing, and therefore allows the evaluation of effects of MV after a ventilation-free period of recovery (7). Indeed, applying this concept of atraumatic intubation by direct laryngoscopy, we recently were able to show significant differences between gene expression changes appearing directly after MV compared to those measured after a ventilation-free interval of 48 h. Immediately after MV, inflammation-related genes showed a transitory modified expression, while another set of more structurally related genes changed their expression only after a delay of 2 d (7). Lung structure, analyzed by conventional 2D histology and also by 3D reconstruction using synchrotron x-ray tomographic microscopy revealed, 48 h after end of MV, a reduced complexity of lung architecture compared to the nonventilated rat lungs, similar to the typical findings in BPD. To extend these observations about late gene expression modifications, we performed with a similar model a full gene expression profile of lung tissue 48 h after the end of MV with either room air or 60% oxygen. Essentially, we measured changes in the expression of genes related to the MMPs and complement system which played a role in many of the six identified mostly affected pathways.

3D-modeling of Mercury's solar wind sputtered surface-exosphere environment

The efficiency of sputtered refractory elements by H+ and He++ solar wind ions from Mercury's surface and their contribution to the exosphere are studied for various solar wind conditions. A 3D solar wind-planetary interaction hybrid model is used for the evaluation of precipitation maps of the sputter agents on Mercury's surface. By assuming a global mineralogical surface composition, the related sputter yields are calculated by means of the 2013 SRIM code and are coupled with a 3D exosphere model. Because of Mercury's magnetic field, for quiet and nominal solar wind conditions the plasma can only precipitate around the polar areas, while for extreme solar events (fast solar wind, coronal mass ejections, interplanetary magnetic clouds) the solar wind plasma has access to the entire dayside. In that case the release of particles form the planet's surface can result in an exosphere density increase of more than one order of magnitude. The corresponding escape rates are also about an order of magnitude higher. Moreover, the amount of He++ ions in the precipitating solar plasma flow enhances also the release of sputtered elements from the surface in the exosphere. A comparison of our model results with MESSENGER observations of sputtered Mg and Ca elements in the exosphere shows a reasonable quantitative agreement. (C) 2015 Elsevier Ltd. All rights reserved.

Modern diatom, cladocera, chironomid, and chrysophyte cyst assemblages as quantitative indicators for the reconstruction of past environmental conditions in the Alps. II. Nutrients

Surface sediments from 68 small lakes in the Alps and 9 well-dated sediment core samples that cover a gradient of total phosphorus (TP) concentrations of 6 to 520 μg TP l-1 were studied for diatom, chrysophyte cyst, cladocera, and chironomid assemblages. Inference models for mean circulation log10 TP were developed for diatoms, chironomids, and benthic cladocera using weighted-averaging partial least squares. After screening for outliers, the final transfer functions have coefficients of determination (r2, as assessed by cross-validation, of 0.79 (diatoms), 0.68 (chironomids), and 0.49 (benthic cladocera). Planktonic cladocera and chrysophytes show very weak relationships to TP and no TP inference models were developed for these biota. Diatoms showed the best relationship with TP, whereas the other biota all have large secondary gradients, suggesting that variables other than TP have a strong influence on their composition and abundance. Comparison with other diatom – TP inference models shows that our model has high predictive power and a low root mean squared error of prediction, as assessed by cross-validation.

Modeling Techniques for the High-Resolution Interpretation of Cryo-Electron Microscopy Reconstructions

Essential biological processes are governed by organized, dynamic interactions between multiple biomolecular systems. Complexes are thus formed to enable the biological function and get dissembled as the process is completed. Examples of such processes include the translation of the messenger RNA into protein by the ribosome, the folding of proteins by chaperonins or the entry of viruses in host cells. Understanding these fundamental processes by characterizing the molecular mechanisms that enable then, would allow the (better) design of therapies and drugs. Such molecular mechanisms may be revealed trough the structural elucidation of the biomolecular assemblies at the core of these processes. Various experimental techniques may be applied to investigate the molecular architecture of biomolecular assemblies. High-resolution techniques, such as X-ray crystallography, may solve the atomic structure of the system, but are typically constrained to biomolecules of reduced flexibility and dimensions. In particular, X-ray crystallography requires the sample to form a three dimensional (3D) crystal lattice which is technically di‑cult, if not impossible, to obtain, especially for large, dynamic systems. Often these techniques solve the structure of the different constituent components within the assembly, but encounter difficulties when investigating the entire system. On the other hand, imaging techniques, such as cryo-electron microscopy (cryo-EM), are able to depict large systems in near-native environment, without requiring the formation of crystals. The structures solved by cryo-EM cover a wide range of resolutions, from very low level of detail where only the overall shape of the system is visible, to high-resolution that approach, but not yet reach, atomic level of detail. In this dissertation, several modeling methods are introduced to either integrate cryo-EM datasets with structural data from X-ray crystallography, or to directly interpret the cryo-EM reconstruction. Such computational techniques were developed with the goal of creating an atomic model for the cryo-EM data. The low-resolution reconstructions lack the level of detail to permit a direct atomic interpretation, i.e. one cannot reliably locate the atoms or amino-acid residues within the structure obtained by cryo-EM. Thereby one needs to consider additional information, for example, structural data from other sources such as X-ray crystallography, in order to enable such a high-resolution interpretation. Modeling techniques are thus developed to integrate the structural data from the different biophysical sources, examples including the work described in the manuscript I and II of this dissertation. At intermediate and high-resolution, cryo-EM reconstructions depict consistent 3D folds such as tubular features which in general correspond to alpha-helices. Such features can be annotated and later on used to build the atomic model of the system, see manuscript III as alternative. Three manuscripts are presented as part of the PhD dissertation, each introducing a computational technique that facilitates the interpretation of cryo-EM reconstructions. The first manuscript is an application paper that describes a heuristics to generate the atomic model for the protein envelope of the Rift Valley fever virus. The second manuscript introduces the evolutionary tabu search strategies to enable the integration of multiple component atomic structures with the cryo-EM map of their assembly. Finally, the third manuscript develops further the latter technique and apply it to annotate consistent 3D patterns in intermediate-resolution cryo-EM reconstructions. The first manuscript, titled An assembly model for Rift Valley fever virus, was submitted for publication in the Journal of Molecular Biology. The cryo-EM structure of the Rift Valley fever virus was previously solved at 27Å-resolution by Dr. Freiberg and collaborators. Such reconstruction shows the overall shape of the virus envelope, yet the reduced level of detail prevents the direct atomic interpretation. High-resolution structures are not yet available for the entire virus nor for the two different component glycoproteins that form its envelope. However, homology models may be generated for these glycoproteins based on similar structures that are available at atomic resolutions. The manuscript presents the steps required to identify an atomic model of the entire virus envelope, based on the low-resolution cryo-EM map of the envelope and the homology models of the two glycoproteins. Starting with the results of the exhaustive search to place the two glycoproteins, the model is built iterative by running multiple multi-body refinements to hierarchically generate models for the different regions of the envelope. The generated atomic model is supported by prior knowledge regarding virus biology and contains valuable information about the molecular architecture of the system. It provides the basis for further investigations seeking to reveal different processes in which the virus is involved such as assembly or fusion. The second manuscript was recently published in the of Journal of Structural Biology (doi:10.1016/j.jsb.2009.12.028) under the title Evolutionary tabu search strategies for the simultaneous registration of multiple atomic structures in cryo-EM reconstructions. This manuscript introduces the evolutionary tabu search strategies applied to enable a multi-body registration. This technique is a hybrid approach that combines a genetic algorithm with a tabu search strategy to promote the proper exploration of the high-dimensional search space. Similar to the Rift Valley fever virus, it is common that the structure of a large multi-component assembly is available at low-resolution from cryo-EM, while high-resolution structures are solved for the different components but lack for the entire system. Evolutionary tabu search strategies enable the building of an atomic model for the entire system by considering simultaneously the different components. Such registration indirectly introduces spatial constrains as all components need to be placed within the assembly, enabling the proper docked in the low-resolution map of the entire assembly. Along with the method description, the manuscript covers the validation, presenting the benefit of the technique in both synthetic and experimental test cases. Such approach successfully docked multiple components up to resolutions of 40Å. The third manuscript is entitled Evolutionary Bidirectional Expansion for the Annotation of Alpha Helices in Electron Cryo-Microscopy Reconstructions and was submitted for publication in the Journal of Structural Biology. The modeling approach described in this manuscript applies the evolutionary tabu search strategies in combination with the bidirectional expansion to annotate secondary structure elements in intermediate resolution cryo-EM reconstructions. In particular, secondary structure elements such as alpha helices show consistent patterns in cryo-EM data, and are visible as rod-like patterns of high density. The evolutionary tabu search strategy is applied to identify the placement of the different alpha helices, while the bidirectional expansion characterizes their length and curvature. The manuscript presents the validation of the approach at resolutions ranging between 6 and 14Å, a level of detail where alpha helices are visible. Up to resolution of 12 Å, the method measures sensitivities between 70-100% as estimated in experimental test cases, i.e. 70-100% of the alpha-helices were correctly predicted in an automatic manner in the experimental data. The three manuscripts presented in this PhD dissertation cover different computation methods for the integration and interpretation of cryo-EM reconstructions. The methods were developed in the molecular modeling software Sculptor (http://sculptor.biomachina.org) and are available for the scientific community interested in the multi-resolution modeling of cryo-EM data. The work spans a wide range of resolution covering multi-body refinement and registration at low-resolution along with annotation of consistent patterns at high-resolution. Such methods are essential for the modeling of cryo-EM data, and may be applied in other fields where similar spatial problems are encountered, such as medical imaging.

Age determination and sea-surface reconstruction of sediment cores from the Norwegian Sea

Polycystine radiolarians are used to reconstruct summer sea surface temperatures (SSSTs) for the Late Pleistocene-Holocene (600-13,400 14C years BP) in the Norwegian Sea. At 13,200 14C years BP, the SSST was close to the average Holocene SSST (~12°C). It then gradually dropped to 7.1°C in the Younger Dryas. Near the Younger Dryas-Holocene transition (~10,000 14C years BP), the SSST increased 5°C in about 530 years. Four abrupt cooling events, with temperature drops of up to 2.1°C, are recognized during the Holocene: at 9340, 7100 ("8200 calendar years event"), 6400 and 1650 14C years BP. Radiolarian SSSTs and the isotopic signal from the GISP2 ice core are strongly coupled, stressing the importance of the Norwegian Sea as a mediator of heat/precipitation exchange between the North Atlantic, the atmosphere, and the Greenland ice sheet. Radiolarian and diatom-derived SSSTs display similarities, with the former not showing the recently reported Holocene cooling trend.

Planktonic foraminifera and sea surface temperture reconstruction for Red Sea sediment cores

In order to assess how insolation-driven climate change superimposed on sea level rise and millennial events influenced the Red Sea during the Holocene, we present new paleoceanographic records from two sediment cores to develop a comprehensive reconstruction of Holocene circulation dynamics in the basin. We show that the recovery of the planktonic foraminiferal fauna after the Younger Dryas was completed earlier in the northern than in the central Red Sea, implying significant changes in the hydrological balance of the northern Red Sea region during the deglaciation. In the early part of the Holocene, the environment of the Red Sea closely followed the development of the Indian summer monsoon and was dominated by a circulation mode similar to the current summer circulation, with low productivity throughout the central and northern Red Sea. The climatic signal during the late Holocene is dominated by a faunal transient event centered around 2.4 ka BP. Its timing corresponds to that of North Atlantic Bond event 2 and to a widespread regionally recorded dry period. This faunal transient is characterized by a more productive foraminiferal fauna and can be explained by an intensification of the winter circulation mode and high evaporation. The modern distribution pattern of planktonic foraminifera, reflecting the prevailing circulation system, was established after 1.7 ka BP.

Oxygen isotopes, IRD and SST reconstruction for the interval MIS31-19 of IODP Site 306-U1314

Early and Mid-Pleistocene climate, ocean hydrography and ice sheet dynamics have been reconstructed using a high-resolution data set (planktonic and benthic d18O time series, faunal-based sea surface temperature (SST) reconstructions and ice-rafted debris (IRD)) record from a high-deposition-rate sedimentary succession recovered at the Gardar Drift formation in the subpolar North Atlantic (Integrated Ocean Drilling Program Leg 306, Site U1314). Our sedimentary record spans from late in Marine Isotope Stage (MIS) 31 to MIS 19 (1069-779 ka). Different trends of the benthic and planktonic oxygen isotopes, SST and IRD records before and after MIS 25 (~940 ka) evidence the large increase in Northern Hemisphere ice-volume, linked to the cyclicity change from the 41-kyr to the 100-kyr that occurred during the Mid-Pleistocene Transition (MPT). Beside longer glacial-interglacial (G-IG) variability, millennial-scale fluctuations were a pervasive feature across our study. Negative excursions in the benthic d18O time series observed at the times of IRD events may be related to glacio-eustatic changes due to ice sheets retreats and/or to changes in deep hydrography. Time series analysis on surface water proxies (IRD, SST and planktonic d18O) of the interval between MIS 31 to MIS 26 shows that the timing of these millennial-scale climate changes are related to half-precessional (10 kyr) components of the insolation forcing, which are interpreted as cross-equatorial heat transport toward high latitudes during both equinox insolation maxima at the equator.

Mean dynamic topography and oceanographic parameters estimated from an inverse model and satellite geodesy, with link to model result in one single NetCDF file (392 MB), including the inverse data error covariance

Geostrophic surface velocities can be derived from the gradients of the mean dynamic topography-the difference between the mean sea surface and the geoid. Therefore, independently observed mean dynamic topography data are valuable input parameters and constraints for ocean circulation models. For a successful fit to observational dynamic topography data, not only the mean dynamic topography on the particular ocean model grid is required, but also information about its inverse covariance matrix. The calculation of the mean dynamic topography from satellite-based gravity field models and altimetric sea surface height measurements, however, is not straightforward. For this purpose, we previously developed an integrated approach to combining these two different observation groups in a consistent way without using the common filter approaches (Becker et al. in J Geodyn 59(60):99-110, 2012, doi:10.1016/j.jog.2011.07.0069; Becker in Konsistente Kombination von Schwerefeld, Altimetrie und hydrographischen Daten zur Modellierung der dynamischen Ozeantopographie, 2012, http://nbn-resolving.de/nbn:de:hbz:5n-29199). Within this combination method, the full spectral range of the observations is considered. Further, it allows the direct determination of the normal equations (i.e., the inverse of the error covariance matrix) of the mean dynamic topography on arbitrary grids, which is one of the requirements for ocean data assimilation. In this paper, we report progress through selection and improved processing of altimetric data sets. We focus on the preprocessing steps of along-track altimetry data from Jason-1 and Envisat to obtain a mean sea surface profile. During this procedure, a rigorous variance propagation is accomplished, so that, for the first time, the full covariance matrix of the mean sea surface is available. The combination of the mean profile and a combined GRACE/GOCE gravity field model yields a mean dynamic topography model for the North Atlantic Ocean that is characterized by a defined set of assumptions. We show that including the geodetically derived mean dynamic topography with the full error structure in a 3D stationary inverse ocean model improves modeled oceanographic features over previous estimates.

Model results of ocean physical field of the reference state and the targets in NetCDF format

In a feasibility study, the potential of proxy data for the temperature and salinity during the Last Glacial Maximum (LGM, about 19 000 to 23 000 years before present) in constraining the strength of the Atlantic meridional overturning circulation (AMOC) with a general ocean circulation model was explored. The proxy data were simulated by drawing data from four different model simulations at the ocean sediment core locations of the Multiproxy Approach for the Reconstruction of the Glacial Ocean surface (MARGO) project, and perturbing these data with realistic noise estimates. The results suggest that our method has the potential to provide estimates of the past strength of the AMOC even from sparse data, but in general, paleo-sea-surface temperature data without additional prior knowledge about the ocean state during the LGM is not adequate to constrain the model. On the one hand, additional data in the deep-ocean and salinity data are shown to be highly important in estimating the LGM circulation. On the other hand, increasing the amount of surface data alone does not appear to be enough for better estimates. Finally, better initial guesses to start the state estimation procedure would greatly improve the performance of the method. Indeed, with a sufficiently good first guess, just the sea-surface temperature data from the MARGO project promise to be sufficient for reliable estimates of the strength of the AMOC.

Planktonic foraminifera record and SST reconstruction for the last 130 ka in the Western South Atlantic

This work reconstructs Late Quaternary paleoceanographic changes in the western South Atlantic Ocean based on sedimentary core GL-77, recovered from the lower continental slope in the Campos basin, offshore SE Brazil. The studied interval comprises the last 130 ka. Changes in sea surface temperature (SST) and paleoproductivity were estimated using the total planktonic foraminiferal fauna and oxygen isotope analyses. The age model was based on the oxygen isotope record, biostratigraphic datums and AMS 14C dating. It was observed that the Pleistocene/Holocene transition occurs within Globorotalia menardii Biozone Y, and is not coeval with the base of Biozone Z. The range between summer and winter SST estimates is larger during the glacial period compared to interglacials. Three peaks of low SST around 70, 50 - 45 and 20 ka coincided with periods of enhanced SE trade winds. Despite faunal differences between the last interglacial (MIS 5e) and the Holocene, our SST estimates suggest that SSTs did not differ significantly between these intervals.

Northwards and eastwards velocities extracted from a historical run (for the year 2000) of the UVic Earth System Climate Model of Weaver et al. (2001)

This study provides a theoretical assessment of the potential bias due to differential lateral transport on multi-proxy studies based on a range of marine microfossils. Microfossils preserved in marine sediments are at the centre of numerous proxies for paleoenvironmental reconstructions. The precision of proxies is based on the assumption that they accurately represent the overlying watercolumn properties and faunas. Here we assess the possibility of a syn-depositional bias in sediment assemblages caused by horizontal drift in the water column, due to differential settling velocities of sedimenting particles based on their shape, size and density, and due to differences in current velocities. Specifically we calculate the post-mortem lateral transport undergone by planktic foraminifera and a range of other biological proxy carriers (diatoms, radiolaria and fecal pellets transporting coccolithophores) in several regions with high current velocities. We find that lateral transport of different planktic foraminiferal species is minimal due to high settling velocities. No significant shape- or size-dependent sorting occurs before reaching the sediment, making planktic foraminiferal ideal proxy carriers. In contrast, diatoms, radiolaria and fecal pellets can be transported up to 500km in some areas. For example in the Agulhas current, transport can lead to differences of up to 2°C in temperature reconstructions between different proxies in response to settling velocities. Therefore, sediment samples are likely to contain different proportions of local and imported particles, decreasing the precision of proxies based on these groups and the accuracy of the temperature reconstruction.

Paleoproductivity reconstruction for the eastern equatorial Pacific

Productivity at six core locations in the eastern equatorial Pacific (EEP) was reconstructed with a benthic foraminiferal transfer function. The core records show strong regionality, especially where affected by Peru margin upwelling of deeper Equatorial Undercurrent Water (EUC) (originally coming from the subantarctic). This "Peru margin" record differs from that seen along the equator where divergence leads to shallow upwelling, and it is generally inverse to that seen in cores outside the areas of equatorial upwelling. Principal components analysis shows that the main productivity pattern correlates well to the global oxygen isotope record and has lowest values during isotope stages 2 and 4. In addition to this, equatorial cores show a higher frequency pattern of variation which becomes much more pronounced during MIS 3 and 2. The reconstructions based on benthic foraminifera were tested against those from nonaccumulation rate based inorganic chemical proxies of export production. These were found to correlate well in the region influenced by Peru upwelling, and also to share common features for sites along the equator. All the nonaccumulation rate based paleotracers are consistent with one another and differ from accumulation rate derived proxies. The differences between the two classes of paleotracers may result from uncertainties in calculating actual biogenic fluxes since 230Th-normalized results conform more to those we obtained. Analysis of planktonic carbon isotope values for the EEP, and their comparison to the record of the Pacific subantarctic, indicates that the subantarctic contribution to the EUC was reduced during MIS 3 and 2.

Model-based changes in global annual mean surface temperature change (Delta T_g) and radiative forcing due to land ice albedo changes (Delta R_[LI]) over the last 5 Myr, supplementary material

It is still an open question how equilibrium warming in response to increasing radiative forcing - the specific equilibrium climate sensitivity S - depends on background climate. We here present palaeodata-based evidence on the state dependency of S, by using CO2 proxy data together with a 3-D ice-sheet-model-based reconstruction of land ice albedo over the last 5 million years (Myr). We find that the land ice albedo forcing depends non-linearly on the background climate, while any non-linearity of CO2 radiative forcing depends on the CO2 data set used. This non-linearity has not, so far, been accounted for in similar approaches due to previously more simplistic approximations, in which land ice albedo radiative forcing was a linear function of sea level change. The latitudinal dependency of ice-sheet area changes is important for the non-linearity between land ice albedo and sea level. In our set-up, in which the radiative forcing of CO2 and of the land ice albedo (LI) is combined, we find a state dependence in the calculated specific equilibrium climate sensitivity, S[CO2,LI], for most of the Pleistocene (last 2.1 Myr). During Pleistocene intermediate glaciated climates and interglacial periods, S[CO2,LI] is on average ~ 45 % larger than during Pleistocene full glacial conditions. In the Pliocene part of our analysis (2.6-5 Myr BP) the CO2 data uncertainties prevent a well-supported calculation for S[CO2,LI], but our analysis suggests that during times without a large land ice area in the Northern Hemisphere (e.g. before 2.82 Myr BP), the specific equilibrium climate sensitivity, S[CO2,LI], was smaller than during interglacials of the Pleistocene. We thus find support for a previously proposed state change in the climate system with the widespread appearance of northern hemispheric ice sheets. This study points for the first time to a so far overlooked non-linearity in the land ice albedo radiative forcing, which is important for similar palaeodata-based approaches to calculate climate sensitivity. However, the implications of this study for a suggested warming under CO2 doubling are not yet entirely clear since the details of necessary corrections for other slow feedbacks are not fully known and the uncertainties that exist in the ice-sheet simulations and global temperature reconstructions are large.