Study of silicon-on-insulator multiple-gate MOS structures including band-gap engineering and self heating effects

The progresses of electron devices integration have proceeded for more than 40 years following the well–known Moore’s law, which states that the transistors density on chip doubles every 24 months. This trend has been possible due to the downsizing of the MOSFET dimensions (scaling); however, new issues and new challenges are arising, and the conventional ”bulk” architecture is becoming inadequate in order to face them. In order to overcome the limitations related to conventional structures, the researchers community is preparing different solutions, that need to be assessed. Possible solutions currently under scrutiny are represented by: • devices incorporating materials with properties different from those of silicon, for the channel and the source/drain regions; • new architectures as Silicon–On–Insulator (SOI) transistors: the body thickness of Ultra-Thin-Body SOI devices is a new design parameter, and it permits to keep under control Short–Channel–Effects without adopting high doping level in the channel. Among the solutions proposed in order to overcome the difficulties related to scaling, we can highlight heterojunctions at the channel edge, obtained by adopting for the source/drain regions materials with band–gap different from that of the channel material. This solution allows to increase the injection velocity of the particles travelling from the source into the channel, and therefore increase the performance of the transistor in terms of provided drain current. The first part of this thesis work addresses the use of heterojunctions in SOI transistors: chapter 3 outlines the basics of the heterojunctions theory and the adoption of such approach in older technologies as the heterojunction–bipolar–transistors; moreover the modifications introduced in the Monte Carlo code in order to simulate conduction band discontinuities are described, and the simulations performed on unidimensional simplified structures in order to validate them as well. Chapter 4 presents the results obtained from the Monte Carlo simulations performed on double–gate SOI transistors featuring conduction band offsets between the source and drain regions and the channel. In particular, attention has been focused on the drain current and to internal quantities as inversion charge, potential energy and carrier velocities. Both graded and abrupt discontinuities have been considered. The scaling of devices dimensions and the adoption of innovative architectures have consequences on the power dissipation as well. In SOI technologies the channel is thermally insulated from the underlying substrate by a SiO2 buried–oxide layer; this SiO2 layer features a thermal conductivity that is two orders of magnitude lower than the silicon one, and it impedes the dissipation of the heat generated in the active region. Moreover, the thermal conductivity of thin semiconductor films is much lower than that of silicon bulk, due to phonon confinement and boundary scattering. All these aspects cause severe self–heating effects, that detrimentally impact the carrier mobility and therefore the saturation drive current for high–performance transistors; as a consequence, thermal device design is becoming a fundamental part of integrated circuit engineering. The second part of this thesis discusses the problem of self–heating in SOI transistors. Chapter 5 describes the causes of heat generation and dissipation in SOI devices, and it provides a brief overview on the methods that have been proposed in order to model these phenomena. In order to understand how this problem impacts the performance of different SOI architectures, three–dimensional electro–thermal simulations have been applied to the analysis of SHE in planar single and double–gate SOI transistors as well as FinFET, featuring the same isothermal electrical characteristics. In chapter 6 the same simulation approach is extensively employed to study the impact of SHE on the performance of a FinFET representative of the high–performance transistor of the 45 nm technology node. Its effects on the ON–current, the maximum temperatures reached inside the device and the thermal resistance associated to the device itself, as well as the dependence of SHE on the main geometrical parameters have been analyzed. Furthermore, the consequences on self–heating of technological solutions such as raised S/D extensions regions or reduction of fin height are explored as well. Finally, conclusions are drawn in chapter 7.

Advanced modelling of time domain electromagnetic data with updated hydrogeological interpretations

Several countries have acquired, over the past decades, large amounts of area covering Airborne Electromagnetic data. Contribution of airborne geophysics has dramatically increased for both groundwater resource mapping and management proving how those systems are appropriate for large-scale and efficient groundwater surveying. We start with processing and inversion of two AEM dataset from two different systems collected over the Spiritwood Valley Aquifer area, Manitoba, Canada respectively, the AeroTEM III (commissioned by the Geological Survey of Canada in 2010) and the “Full waveform VTEM” dataset, collected and tested over the same survey area, during the fall 2011. We demonstrate that in the presence of multiple datasets, either AEM and ground data, due processing, inversion, post-processing, data integration and data calibration is the proper approach capable of providing reliable and consistent resistivity models. Our approach can be of interest to many end users, ranging from Geological Surveys, Universities to Private Companies, which are often proprietary of large geophysical databases to be interpreted for geological and\or hydrogeological purposes. In this study we deeply investigate the role of integration of several complimentary types of geophysical data collected over the same survey area. We show that data integration can improve inversions, reduce ambiguity and deliver high resolution results. We further attempt to use the final, most reliable output resistivity models as a solid basis for building a knowledge-driven 3D geological voxel-based model. A voxel approach allows a quantitative understanding of the hydrogeological setting of the area, and it can be further used to estimate the aquifers volumes (i.e. potential amount of groundwater resources) as well as hydrogeological flow model prediction. In addition, we investigated the impact of an AEM dataset towards hydrogeological mapping and 3D hydrogeological modeling, comparing it to having only a ground based TEM dataset and\or to having only boreholes data.

TALDICE-1 age scale of the Talos Dome deep ice core, East Antarctica

A new deep ice core drilling program, TALDICE, has been successfully handled by a European team at Talos Dome, in the Ross Sea sector of East Antarctica, down to 1620 m depth. Using stratigraphic markers and a new inverse method, we produce the first official chronology of the ice core, called TALDICE-1. We show that it notably improves an a priori chronology resulting from a one-dimensional ice flow model. It is in agreement with a posteriori controls of the resulting accumulation rate and thinning function along the core. An absolute uncertainty of only 300 yr is obtained over the course of the last deglaciation. This uncertainty remains lower than 600 yr over Marine Isotope Stage 3, back to 50 kyr BP. The phasing of the TALDICE ice core climate record with respect to the central East Antarctic plateau and Greenland records can thus be determined with a precision allowing for a discussion of the mechanisms at work at sub-millennial time scales.

Synchronisation of the EDML and EDC ice cores for the last 52 kyr by volcanic signature matching

A common time scale for the EPICA ice cores from Dome C (EDC) and Dronning Maud Land (EDML) has been established. Since the EDML core was not drilled on a dome, the development of the EDML1 time scale for the EPICA ice core drilled in Dronning Maud Land was based on the creation of a detailed stratigraphic link between EDML and EDC, which was dated by a simpler 1D ice-flow model. The synchronisation between the two EPICA ice cores was done through the identification of several common volcanic signatures. This paper describes the rigorous method, using the signature of volcanic sulfate, which was employed for the last 52 kyr of the record. We estimated the discrepancies between the modelled EDC and EDML glaciological age scales during the studied period, by evaluating the ratio R of the apparent duration of temporal intervals between pairs of isochrones. On average R ranges between 0.8 and 1.2 corresponding to an uncertainty of up to 20% in the estimate of the time duration in at least one of the two ice cores. Significant deviations of R up to 1.4–1.5 are observed between 18 and 28 kyr before present (BP), where present is defined as 1950. At this stage our approach does not allow us unequivocally to find out which of the models is affected by errors, but assuming that the thinning function at both sites and accumulation history at Dome C (which was drilled on a dome) are correct, this anomaly can be ascribed to a complex spatial accumulation variability (which may be different in the past compared to the present day) upstream of the EDML core.

Simulating dioxane transport in a heterogeneous glacial aquifer system (Washtenaw County, Michigan) using publicly available models and data

The primary challenge in groundwater and contaminant transport modeling is obtaining the data needed for constructing, calibrating and testing the models. Large amounts of data are necessary for describing the hydrostratigraphy in areas with complex geology. Increasingly states are making spatial data available that can be used for input to groundwater flow models. The appropriateness of this data for large-scale flow systems has not been tested. This study focuses on modeling a plume of 1,4-dioxane in a heterogeneous aquifer system in Scio Township, Washtenaw County, Michigan. The analysis consisted of: (1) characterization of hydrogeology of the area and construction of a conceptual model based on publicly available spatial data, (2) development and calibration of a regional flow model for the site, (3) conversion of the regional model to a more highly resolved local model, (4) simulation of the dioxane plume, and (5) evaluation of the model's ability to simulate field data and estimation of the possible dioxane sources and subsequent migration until maximum concentrations are at or below the Michigan Department of Environmental Quality's residential cleanup standard for groundwater (85 ppb). MODFLOW-2000 and MT3D programs were utilized to simulate the groundwater flow and the development and movement of the 1, 4-dioxane plume, respectively. MODFLOW simulates transient groundwater flow in a quasi-3-dimensional sense, subject to a variety of boundary conditions that can simulate recharge, pumping, and surface-/groundwater interactions. MT3D simulates solute advection with groundwater flow (using the flow solution from MODFLOW), dispersion, source/sink mixing, and chemical reaction of contaminants. This modeling approach was successful at simulating the groundwater flows by calibrating recharge and hydraulic conductivities. The plume transport was adequately simulated using literature dispersivity and sorption coefficients, although the plume geometries were not well constrained.

Kinetic characterization for pretreatment of timber varieties and switchgrass using diluted acid hydrolysis

In recent years, growing attention has been devoted to the use of lignocellulosic biomass as a feedstock to produce renewable carbohydrates as a source of energy products, including liquid alternatives to fossil fuels. The benefits of developing woody biomass to ethanol technology are to increase the long-term national energy security, reduce fossil energy consumption, lower greenhouse gas emissions, use renewable rather than depletable resources, and create local jobs. Currently, research is driven by the need to reduce the cost of biomass-ethanol production. One of the preferred methods is to thermochemically pretreat the biomass material and subsequently, enzymatically hydrolyze the pretreated material to fermentable sugars that can then be converted to ethanol using specialized microorganisms. The goals of pretreatment are to remove the hemicellulose fraction from other biomass components, reduce bioconversion time, enhance enzymatic conversion of the cellulose fraction, and, hopefully, obtain a higher ethanol yield. The primary goal of this research is to obtain kinetic detailed data for dilute acid hydrolysis for several timber species from the Upper Peninsula of Michigan and switchgrass. These results will be used to identify optimum reaction conditions to maximize production of fermentable sugars and minimize production of non-fermentable byproducts. The structural carbohydrate analysis of the biomass species used in this project was performed using the procedure proposed by National Renewable Energy Laboratory (NREL). Subsequently, dilute acid-catalyzed hydrolysis of biomass, including aspen, basswood, balsam, red maple, and switchgrass, was studied at various temperatures, acid concentrations, and particle sizes in a 1-L well-mixed batch reactor (Parr Instruments, ii Model 4571). 25 g of biomass and 500 mL of diluted acid solution were added into a 1-L glass liner, and then put into the reactor. During the experiment, 5 mL samples were taken starting at 100°C at 3 min intervals until reaching the targeted temperature (160, 175, or 190°C), followed by 4 samples after achieving the desired temperature. The collected samples were then cooled in an ice bath immediately to stop the reaction. The cooled samples were filtered using 0.2 μm MILLIPORE membrane filter to remove suspended solids. The filtered samples were then analyzed using High Performance Liquid Chromatography (HPLC) with a Bio-Rad Aminex HPX-87P column, and refractive index detection to measure monomeric and polymeric sugars plus degradation byproducts. A first order reaction model was assumed and the kinetic parameters such as activation energy and pre-exponential factor from Arrhenius equation were obtained from a match between the model and experimental data. The reaction temperature increases linearly after 40 minutes during experiments. Xylose and other sugars were formed from hemicellulose hydrolysis over this heat up period until a maximum concentration was reached at the time near when the targeted temperature was reached. However, negligible amount of xylose byproducts and small concentrations of other soluble sugars, such as mannose, arabinose, and galactose were detected during this initial heat up period. Very little cellulose hydrolysis yielding glucose was observed during the initial heat up period. On the other hand, later in the reaction during the constant temperature period xylose was degraded to furfural. Glucose production from cellulose was increased during this constant temperature period at later time points in the reaction. The kinetic coefficient governing the generation of xylose from hemicellulose and the generation of furfural from xylose presented a coherent dependence on both temperature and acid concentration. However, no effect was observed in the particle size. There were three types of biomass used in this project; hardwood (aspen, basswood, and red maple), softwood (balsam), and a herbaceous crop (switchgrass). The activation energies and the pre-exponential factors of the timber species and switchgrass were in a range of 49 - 180 kJ/mol and from 7.5x104 - 2.6x1020 min-1, respectively, for the xylose formation model. In addition, for xylose degradation, the activation energies and the preexponential factors ranged from 130 - 170 kJ/mol and from 6.8x1013 - 3.7x1017 min-1, respectively. The results compare favorably with the literature values given by Ranganathan et al, 1985. Overall, up to 92 % of the xylose was able to generate from the dilute acid hydrolysis in this project.

Critical evaluation of diameter increment modelling in the Great Lakes region

Simulations of forest stand dynamics in a modelling framework including Forest Vegetation Simulator (FVS) are diameter driven, thus the diameter or basal area increment model needs a special attention. This dissertation critically evaluates diameter or basal area increment models and modelling approaches in the context of the Great Lakes region of the United States and Canada. A set of related studies are presented that critically evaluate the sub-model for change in individual tree basal diameter used in the Forest Vegetation Simulator (FVS), a dominant forestry model in the Great Lakes region. Various historical implementations of the STEMS (Stand and Tree Evaluation and Modeling System) family of diameter increment models, including the current public release of the Lake States variant of FVS (LS-FVS), were tested for the 30 most common tree species using data from the Michigan Forest Inventory and Analysis (FIA) program. The results showed that current public release of the LS-FVS diameter increment model over-predicts 10-year diameter increment by 17% on average. Also the study affirms that a simple adjustment factor as a function of a single predictor, dbh (diameter at breast height) used in the past versions, provides an inadequate correction of model prediction bias. In order to re-engineer the basal diameter increment model, the historical, conceptual and philosophical differences among the individual tree increment model families and their modelling approaches were analyzed and discussed. Two underlying conceptual approaches toward diameter or basal area increment modelling have been often used: the potential-modifier (POTMOD) and composite (COMP) approaches, which are exemplified by the STEMS/TWIGS and Prognosis models, respectively. It is argued that both approaches essentially use a similar base function and neither is conceptually different from a biological perspective, even though they look different in their model forms. No matter what modelling approach is used, the base function is the foundation of an increment model. Two base functions – gamma and Box-Lucas – were identified as candidate base functions for forestry applications. The results of a comparative analysis of empirical fits showed that quality of fit is essentially similar, and both are sufficiently detailed and flexible for forestry applications. The choice of either base function in order to model diameter or basal area increment is dependent upon personal preference; however, the gamma base function may be preferred over the Box-Lucas, as it fits the periodic increment data in both a linear and nonlinear composite model form. Finally, the utility of site index as a predictor variable has been criticized, as it has been widely used in models for complex, mixed species forest stands though not well suited for this purpose. An alternative to site index in an increment model was explored, using site index and a combination of climate variables and Forest Ecosystem Classification (FEC) ecosites and data from the Province of Ontario, Canada. The results showed that a combination of climate and FEC ecosites variables can replace site index in the diameter increment model.

COMPUTATIONAL ANALYSIS OF THE AEROELASTIC DYNAMICS OF A WIND TURBINE BLADE WITH VARIABLE-SPEED STALL CONTROL

The accuracy of simulating the aerodynamics and structural properties of the blades is crucial in the wind-turbine technology. Hence the models used to implement these features need to be very precise and their level of detailing needs to be high. With the variety of blade designs being developed the models should be versatile enough to adapt to the changes required by every design. We are going to implement a combination of numerical models which are associated with the structural and the aerodynamic part of the simulation using the computational power of a parallel HPC cluster. The structural part models the heterogeneous internal structure of the beam based on a novel implementation of the Generalized Timoshenko Beam Model Technique.. Using this technique the 3-D structure of the blade is reduced into a 1-D beam which is asymptotically equivalent. This reduces the computational cost of the model without compromising its accuracy. This structural model interacts with the Flow model which is a modified version of the Blade Element Momentum Theory. The modified version of the BEM accounts for the large deflections of the blade and also considers the pre-defined structure of the blade. The coning, sweeping of the blade, tilt of the nacelle and the twist of the sections along the blade length are all computed by the model which aren’t considered in the classical BEM theory. Each of these two models provides feedback to the other and the interactive computations lead to more accurate outputs. We successfully implemented the computational models to analyze and simulate the structural and aerodynamic aspects of the blades. The interactive nature of these models and their ability to recompute data using the feedback from each other makes this code more efficient than the commercial codes available. In this thesis we start off with the verification of these models by testing it on the well-known benchmark blade for the NREL-5MW Reference Wind Turbine, an alternative fixed-speed stall-controlled blade design proposed by Delft University, and a novel alternative design that we proposed for a variable-speed stall-controlled turbine, which offers the potential for more uniform power control and improved annual energy production.. To optimize the power output of the stall-controlled blade we modify the existing designs and study their behavior using the aforementioned aero elastic model.

Coordinated Lot-sizing and Dynamic Prizing under a Supplier All-units Quantity Discount

We consider an economic order quantity model where the supplier offers an all-units quantity discount and a price sensitive customer demand. We compare a decentralized decision framework where selling price and replenishment policy are determined independently to simultaneous decision making. Constant and dynamic pricing are distinguished. We derive structural properties and develop algorithms that determine the optimal pricing and replenishment policy and show how quantity discounts not only influence the purchasing strategy but also the pricing policy. A sensitivity analysis indicates the impact of the fixed-holding cost ratio, the discount policy, and the customers' price sensitivity on the optimal decisions.

Tree species and diversity effects on soil water seepage in a tropical plantation

Plant diversity has been shown to influence the water cycle of forest ecosystems by differences in water consumption and the associated effects on groundwater recharge. However, the effects of biodiversity on soil water fluxes remain poorly understood for native tree species plantations in the tropics. Therefore, we estimated soil water fluxes and assessed the effects of tree species and diversity on these fluxes in an experimental native tree species plantation in Sardinilla (Panama). The study was conducted during the wet season 2008 on plots of monocultures and mixtures of three or six tree species. Rainfall and soil water content were measured and evapotranspiration was estimated with the Penman-Monteith equation. Soil water fluxes were estimated using a simple soil water budget model considering water input, output, and soil water and groundwater storage changes and in addition, were simulated using the physically based one-dimensional water flow model Hydrus-1D. In general, the Hydrus simulation did not reflect the observed pressure heads, in that modeled pressure heads were higher compared to measured ones. On the other hand, the results of the water balance equation (WBE) reproduced observed water use patterns well. In monocultures, the downward fluxes through the 200 cm-depth plane were highest below Hura crepitans (6.13 mm day−1) and lowest below Luehea seemannii (5.18 mm day−1). The average seepage rate in monocultures (±SE) was 5.66 ± 0.18 mm day−1, and therefore, significantly higher than below six-species mixtures (5.49 ± 0.04 mm day−1) according to overyielding analyses. The three-species mixtures had an average seepage rate of 5.63 ± 0.12 mm day−1 and their values did not differ significantly from the average values of the corresponding species in monocultures. Seepage rates were driven by the transpiration of the varying biomass among the plots (r = 0.61, p = 0.017). Thus, a mixture of trees with different growth rates resulted in moderate seepage rates compared to monocultures of either fast growing or slow growing tree species. Our results demonstrate that tree-species specific biomass production and tree diversity are important controls of seepage rates in the Sardinilla plantation during the wet season.

On new maximal supergravity and its BPS domain-walls

We revise the SU(3)-invariant sector of N  = 8 supergravity with dyonic SO(8) gaugings. By using the embedding tensor formalism, analytic expressions for the scalar potential, superpotential(s) and fermion mass terms are obtained as a function of the electromagnetic phase ω and the scalars in the theory. Equipped with these results, we explore non-supersymmetric AdS critical points at ω ≠ 0 for which perturbative stability could not be analysed before. The ω-dependent superpotential is then used to derive first-order flow equations and obtain new BPS domain-wall solutions at ω ≠ 0. We numerically look at steepest-descent paths motivated by the (conjectured) RG flows.

Where to find 1.5 million yr old ice for the IPICS "Oldest-Ice" ice core

The recovery of a 1.5 million yr long ice core from Antarctica represents a keystone of our understanding of Quaternary climate, the progression of glaciation over this time period and the role of greenhouse gas cycles in this progression. Here we tackle the question of where such ice may still be found in the Antarctic ice sheet. We can show that such old ice is most likely to exist in the plateau area of the East Antarctic ice sheet (EAIS) without stratigraphic disturbance and should be able to be recovered after careful pre-site selection studies. Based on a simple ice and heat flow model and glaciological observations, we conclude that positions in the vicinity of major domes and saddle position on the East Antarctic Plateau will most likely have such old ice in store and represent the best study areas for dedicated reconnaissance studies in the near future. In contrast to previous ice core drill site selections, however, we strongly suggest significantly reduced ice thickness to avoid bottom melting. For example for the geothermal heat flux and accumulation conditions at Dome C, an ice thickness lower than but close to about 2500 m would be required to find 1.5 Myr old ice (i.e., more than 700 m less than at the current EPICA Dome C drill site). Within this constraint, the resolution of an Oldest-Ice record and the distance of such old ice to the bedrock should be maximized to avoid ice flow disturbances, for example, by finding locations with minimum geothermal heat flux. As the geothermal heat flux is largely unknown for the EAIS, this parameter has to be carefully determined beforehand. In addition, detailed bedrock topography and ice flow history has to be reconstructed for candidates of an Oldest-Ice ice coring site. Finally, we argue strongly for rapid access drilling before any full, deep ice coring activity commences to bring datable samples to the surface and to allow an age check of the oldest ice.

Performance analysis of a miniature turbine generator for intracorporeal energy harvesting

Replacement intervals of implantable medical devices are commonly dictated by battery life. Therefore, intracorporeal energy harvesting has the potential to reduce the number of surgical interventions by extending the life cycle of active devices. Given the accumulated experience with intravascular devices such as stents, heart valves, and cardiac assist devices, the idea to harvest a small fraction of the hydraulic energy available in the cardiovascular circulation is revisited. The aim of this article is to explore the technical feasibility of harvesting 1 mW electric power using a miniature hydrodynamic turbine powered by about 1% of the cardiac output flow in a peripheral artery. To this end, numerical modelling of the fluid mechanics and experimental verification of the overall performance of a 1:1 scale friction turbine are performed in vitro. The numerical flow model is validated for a range of turbine configurations and flow conditions (up to 250 mL/min) in terms of hydromechanic efficiency; up to 15% could be achieved with the nonoptimized configurations of the study. Although this article does not entail the clinical feasibility of intravascular turbines in terms of hemocompatibility and impact on the circulatory system, the numerical model does provide first estimates of the mechanical shear forces relevant to blood trauma and platelet activation. It is concluded that the time-integrated shear stress exposure is significantly lower than in cardiac assist devices due to lower flow velocities and predominantly laminar flow.

Northern Hemisphere temperature reconstruction during the last millennium using multiple annual proxies

Previous studies have either exclusively used annual tree-ring data or have combined tree-ring series with other, lower temporal resolution proxy series. Both approaches can lead to significant uncertainties, as tree-rings may underestimate the amplitude of past temperature variations, and the validity of non-annual records cannot be clearly assessed. In this study, we assembled 45 published Northern Hemisphere (NH) temperature proxy records covering the past millennium, each of which satisfied 3 essential criteria: the series must be of annual resolution, span at least a thousand years, and represent an explicit temperature signal. Suitable climate archives included ice cores, varved lake sediments, tree-rings and speleothems. We reconstructed the average annual land temperature series for the NH over the last millennium by applying 3 different reconstruction techniques: (1) principal components (PC) plus second-order autoregressive model (AR2), (2) composite plus scale (CPS) and (3) regularized errors-in-variables approach (EIV). Our reconstruction is in excellent agreement with 6 climate model simulations (including the first 5 models derived from the fifth phase of the Coupled Model Intercomparison Project (CMIP5) and an earth system model of intermediate complexity (LOVECLIM), showing similar temperatures at multi-decadal timescales; however, all simulations appear to underestimate the temperature during the Medieval Warm Period (MWP). A comparison with other NH reconstructions shows that our results are consistent with earlier studies. These results indicate that well-validated annual proxy series should be used to minimize proxy-based artifacts, and that these proxy series contain sufficient information to reconstruct the low-frequency climate variability over the past millennium.

Spatial gradients of temperature, accumulation and delta18O-ice in Greenland over a series of Dansgaard-Oeschger events

Air and water stable isotope measurements from four Greenland deep ice cores (GRIP, GISP2, NGRIP and NEEM) are investigated over a series of Dansgaard–Oeschger events (DO 8, 9 and 10), which are representative of glacial millennial scale variability. Combined with firn modeling, air isotope data allow us to quantify abrupt temperature increases for each drill site (1σ = 0.6 °C for NEEM, GRIP and GISP2, 1.5 °C for NGRIP). Our data show that the magnitude of stadial–interstadial temperature increase is up to 2 °C larger in central and North Greenland than in northwest Greenland: i.e., for DO 8, a magnitude of +8.8 °C is inferred, which is significantly smaller than the +11.1 °C inferred at GISP2. The same spatial pattern is seen for accumulation increases. This pattern is coherent with climate simulations in response to reduced sea-ice extent in the Nordic seas. The temporal water isotope (δ18O)–temperature relationship varies between 0.3 and 0.6 (±0.08) ‰ °C−1 and is systematically larger at NEEM, possibly due to limited changes in precipitation seasonality compared to GISP2, GRIP or NGRIP. The gas age−ice age difference of warming events represented in water and air isotopes can only be modeled when assuming a 26% (NGRIP) to 40% (GRIP) lower accumulation than that derived from a Dansgaard–Johnsen ice flow model.