Search

Deoxygenation affects many continental shelf seas across the world today and results in increasing areas of hypoxia (dissolved oxygen concentration ([O2])

Supplemental material, Supplementary_Information_1_and_2 for Mid to late-Holocene sea-surface temperature variability off north-eastern Newfoundland and its linkage to the North Atlantic Oscillation by Lisa C Orme, Arto Miettinen, Marit-Solveig Seidenkrantz, Kirsi Tuominen, Christof Pearce, Dmitry V Divine, Mimmi Oksman and Antoon Kuijpers in The Holocene

The chemical composition of foraminiferal calcite is widely used for studying past environmental conditions and biogeochemistry. However, high rates of microbial organic matter degradation and abundant dissolved metal sources in sediments and pore waters may impede the application of foraminifera-based proxies due to formation of secondary carbonates or other authigenic minerals on and/or inside of foraminiferal tests. Secondary carbonate precipitation severely alters the bulk foraminiferal geochemistry and can be difficult to eliminate through standard foraminiferal trace element cleaning procedures. We present results showing the mineral composition and formation sequence of diagenetic coatings on foraminiferal tests formed under extreme anoxic conditions in the Baltic Sea's deepest basin (the Landsort Deep, IODP Exp. 347, Site M0063). Our study focuses primarily on the diagenetic carbonates present on and in the tests of the low-oxygen tolerant benthic foraminiferal species Elphidium selseyensis and Elphidium clavatum. We applied various geochemical and imaging methods to ascertain the diagenesis processes and the authigenic mineral formation sequence on foraminifera. The authigenic carbonates were enriched in Mn, Mg, Fe and Ba, depending on the environmental redox conditions when the authigenic carbonates were precipitated. Concentrations of redox-sensitive elements such as Mn and Fe were particularly elevated in bottom waters and sedimentary pore waters under oxygen-depleted conditions in the Landsort Deep, resulting in formation of carbonates with highly elevated Mn and Fe contents. In addition to Mn- and Fe carbonates, several other types of authigenic minerals also formed on and in the foraminiferal chambers, including authigenic calcite, and non-carbonate accessory minerals. The formation sequence reveals the redox sensitivities of different elements and the preferential sequence of oxidants used by the microbes during organic matter oxidation and secondary redox reactions. This study provides a case study of extreme early diagenesis of foraminiferal calcite and may serve as a valuable guide when interpreting foraminiferal trace element records from low oxygen environments.

Marine fjords with active glacier outlets are hot spots for organic matter burial in the sediments and subsequent microbial mineralization, and will be increasingly important as climate warming causes more rapid glacial melt. Here, we investigated controls on microbial community assembly in sub-arctic glacier-influenced (GI) and non-glacier-influenced (NGI) marine sediments in the Godthåbsfjord region, south-western Greenland. We used a correlative approach integrating 16S rRNA gene and dissimilatory sulfite reductase (dsrB) amplicon sequence data over six meters of depth with biogeochemistry, sulfur-cycling activities, and sediment ages. GI sediments were characterized by comparably high sedimentation rates and had ‘young’ sediment ages of Desulfobacteraceaemembers, correlated positively with sulfate reduction rates, indicating their active contributions to sulfur-cycling processes. In contrast, other surface community members, such asDesulfatiglans, AtribacteriaandChloroflexi, survived the slow sediment burial at NGI stations and dominated in the deepest sediment layers. These taxa are typical for the energy-limited marine deep biosphere and their relative abundances correlated positively with sediment age. In conclusion, our data suggests that high rates of sediment accumulation caused by glacier runoff and associated changes in biogeochemistry, promote persistence of sulfur-cycling activity and burial of a larger fraction of the surface microbial community into the deep subsurface.Contribution to the Field StatementIn most coastal marine sediments organic matter turnover and total energy flux are highest at the surface and decrease significantly with increasing sediment depth, causing depth-dependent changes in the microbial community composition. Glacial runoff in arctic and subarctic fjords alters the composition of the microbial community at the surface, mainly due to different availabilities of organic matter and metals. Here we show that glacial runoff also modifies microbial community assembly with sediment depth. Sediment age was a key driver of microbial community composition in six-meter-long marine sediment cores from the Godthåbsfjord region, south-western Greenland. High sedimentation rates at glacier-influenced sediment stations enabled a complex community of sulfur-cycling-associated microorganisms to continuously thrive at high relative abundances from the surface into the sediment subsurface. These communities consisted of putative fermenters, sulfate reducers and sulfur oxidizers, which likely depended on high metal concentrations in the relatively young, glacier-influenced sediments. In non-glacier-influenced sediments with lower sedimentation rates, these sulfur-cycling-associated microorganisms were only present near the surface. With increasing sediment depth these surface microorganisms were largely replaced by other surface microorganisms that positively correlated with sediment age and belong to known taxa of the energy-limited, marine deep biosphere.

The East Asian winter monsoon (EAWM) significantly impacts living conditions in a large part of Asia, and therefore, it is important to understand its major driving mechanisms. Winter sea surface temperature (SSTW) and circulation in the southern Okinawa Trough are today both primarily controlled by the EAWM. Here we present a new SSTW reconstruction for the last millennium based on a diatom record from sediment core MD05-2908, from the continental slope of the southern Okinawa Trough off northeastern Taiwan. Our reconstruction indicates that SSTW varied between 14.1 and 19.6°C over the past 1,000 years. Changes in SSTW in the southern Okinawa Trough correspond closely to the index of warm winters based on historical documents from the East Asian monsoon domain. This implies that our SSTW record can be used to reconstruct EAWM variability during the last millennium. Comparisons with the reconstructed winter Arctic Oscillation (AO, developed from historical snow anomaly events in Eastern Asia) and Arctic sea ice cover reveal a significant positive correlation between the EAWM and AO during the time interval from 1000–1400 Common Era (C.E.), coinciding with reduced sea ice cover. However, there is no significant correlation with increased sea ice cover during the interval from 1400 to 1700 C.E. This suggests that the reduction in Arctic sea ice may periodically have played a role in strengthening the relationship between the EAWM and the AO during the last millennium and that the current and future reduction in Arctic sea ice may have significant consequences for the EAWM.

Climate during the Last Interglacial period (LIG, Marine Isotope Stage 5e) was on average warmer than the present, with a higher global sea level but also more unstable conditions. Today, the Baltic Sea interacts strongly with conditions in the North Atlantic region, and this interaction was likely even stronger during the LIG. We here present a reconstruction of seawater conditions during the LIG based on benthic foraminiferal geochemistry (stable isotopes and trace elements) and compare these records with modern marine monitoring data to evaluate seasonal hydrographic conditions in the western and southern Baltic Sea during the first half of the LIG (130–123 ka BP). Our reconstructions reflect the evolution of seasonal temperature and salinity, rather than annual mean conditions. The spring LIG bottom water temperatures in the Skagerrak and Kattegat were ∼2–3 °C higher compared to the modern spring bottom water temperatures. During the LIG, there was an increase in seasonality in bottom water temperature (progressively warmer summers and cooler springs) in the southern Baltic Sea, which can be linked to seasonal insolation changes. Moreover, our data suggest a decreased gradient of bottom water salinity along a transect through the Skagerrak-Kattegat-Danish Straits-southern Baltic Sea, supporting previous investigations inferring a stronger ocean-water influx into the Baltic Sea during the LIG than at present.

We present a reconstruction of Holocene sea-ice variability from sediment core GA306-GC3, from the Holsteinsborg Dyb off West Greenland, which provides an index of palaeoceanographic and palaeoenvironmental conditions within this climatically sensitive region during the last 6700 yr. The reconstructed sea-ice record, combined with previously published proxy data, suggests that relatively warm conditions with reduced sea-ice extent prior to 5000 cal. yr BP were associated with the Holocene Thermal Maximum. Subsequent cooling and extensive sea-ice cover between ca. 5000 and 4000 cal. yr BP was followed by even colder conditions and persistent sea-ice cover during the Neoglacial cooling phase, particularly after ca. 1500 cal. yr BP. There is a positive correlation between West Greenland sea-ice cover and solar activity over the past 5000 yr, but the correlation is much weaker prior to 5000 cal. yr BP. In addition, there is a strong link between West Greenland sea ice and changes in the abundance of arctic benthic foraminifera species, related to different water masses of the West Greenland Current during the entire interval, even prior to 5000 cal. yr BP. Our findings indicate that sea-ice variability off West Greenland was driven not only by solar activity, but also by ocean circulation (the strength of cold Polar water from the East Greenland Current and warm Atlantic water from the Irminger Current).

Using qualitative and quantitative analyses of planktonic foraminifera assemblages, δ18O measurements, Mg/Ca–temperature relationship ( Globigerinoides ruber pink sensu stricto) data and sea-surface temperatures (SSTs) derived from artificial neural network (ANN) transfer functions, a reconstruction was made of upper water mass variability in the Anegada–Jungfern Passage (AJP), northeastern Caribbean Sea, over the last c. 11,100 years. Our record is based on the study of two marine sediment cores and reveals three main circulation stages during the Holocene. In the early-Holocene (11,100–6300 cal. yr BP), ANN-based SST estimations indicate moderately cooler than present winter condition in the NE Caribbean with SSTs of ~25.5°C. These conditions presumably reflect advection of well-mixed upper water masses from the Guyana upwelling area associated with a strong Atlantic Meridional Overturning Circulation (AMOC) and enhanced trade wind activity, linked to a more northerly location of the inter-tropical convergence zone (ITCZ). Between 6300 and 3700 cal. yr BP, a relative warming of winter SSTs (~26.5°C) was probably related to a weaker circulation and upper water mass mixing because of less intense trade wind activity, as the ITCZ moved southwards. From 3700 cal. yr BP to the present, the region was characterised by small seasonal SST variations and generally stable winter and summer SSTs. The data suggest a minor shift in the (sub)surface inflow pattern during the last 2000 years, possibly related to changes in Northern Hemisphere large-scale atmospheric circulation also observed at higher latitude. The ANN-based temperature pattern is supported by fluctuations in the Mg/Ca-derived temperature record, although temperature maxima derived from the Mg/Ca ratio appear anomalously high. On a Holocene timescale, we conclude that the northeastern Caribbean SST and circulation regime have been mainly dependent on the position of the ITCZ, which, in turn, is controlled by changes in hemispheric solar insolation.

We here reconstruct a glacial and climate history of arctic NW Siberia for the last ~600,000 years, based on the stratigraphy and chronology of 35 studied river sections on the southern Taimyr Peninsula. From this stratigraphic mosaic we have identified four glacial events, marked by tills/glaciotectonics, which are intercalated with mainly marine sediments deposited in proglacial settings during transitions from glacial conditions into subsequent interglacials/interstadials. The traces of early shelf-based Kara Sea Ice Sheet (KSIS) glaciations in marine isotope stages (MIS) 12–14 and 8 are sparsely preserved, but these ice advances are suggested to have terminated far south into the central Siberian uplands, as also was the case with the younger Taz glaciation (MIS 6). The inception phase of the latter glaciation was complex, with ice advancing into a proglacial marine basin both from the south (Putorana – Anabar uplands) and the north. The deglaciation leading into the Karginsky interglacial (MIS 5e) was marked by the development of the southernmost ice-marginal zones (IMZs) on the Taimyr lowlands – the Urdakh and Sampesa IMZs. The most recent (late Pleistocene) glacial cycle is recorded by three successively smaller KSIS advances from the Kara Sea shelf onto Taimyr, of which only the first, during Early Zyryanka (MIS 5d), reached south of the Byrranga Mountains, with its maximum extent marked by the Jangoda – Syntabul – Severokokorsky IMZ. Retreat of the ice margin during MIS 5c-b was accompanied by deposition of glaciomarine sediment in the proglacial basin and deposition of large successions of delta sediments in the foothills of the Byrranga Mountains, reaching ≥100 m above present sea level. The region north of the Byrranga Mountains was subjected to two subsequent KSIS glaciations, during MIS 4 and MIS 2, while the area south of the Byrranga Mountains transitioned to a terrestrial environment from the Middle into the Lower Zyryanka, as evidenced by deposition of fluvial, aeolian and ice-complex (Yedoma) sediments.

International audience

Supplemental material, Supplementary_Information for Holocene water mass changes in the Labrador Current by Annalena Antonia Lochte, Janne Repschläger, Marit-Solveig Seidenkrantz, Markus Kienast, Thomas Blanz and Ralph R Schneider in The Holocene

This paper presents new marine evidence of Lateglacial and Holocene environmental changes in the western part of Limfjorden, and provides a review of the geological history/development of this part of northern Jylland, Denmark. Lateglacial clay without fossils is widespread in the region and is probably a glaciolacustrine deposit. Limfjorden began to form as a strait in the Early Holocene due to rising relative sea level and the oldest marine shells are dated to c. 9300 cal. years BP. We propose a new relative sealevel curve for the region based on new and published data, which appear to confirm that the relative sea-level change was not extremely rapid, which was suggested earlier. During the Mid-Holocene a wide connection existed from the western part of Limfjorden to the North Sea in the west and more narrow connections existed between Limfjorden and Skagerrak in the north. The marine fauna included several species that indicate warmer and more salty waters than at present. Gradually, the connections to the North Sea and Skagerrak closed due to long-shore sediment transport and deposition of aeolian sand combined with a fall in the relative sea level during the Middle- to Late Holocene. During the Viking Age, 800–1050 CE (Common Era), the western connection to the North Sea was still open, but around 1200 CE it was closed by a coastal sandy barrier and the western part of Limfjorden became brackish. The coastal barrier was flooded on several occasions but soon formed again. After 1825 CE the western connection from Limfjorden to the North Sea has been maintained artificially.

Although geological and modelling evidence indicate that the last glacial inception in North America was in NE Canada, little is known about the glacial response of the nearby western Greenland Ice Sheet (GIS) during the glacial advance of marine oxygen isotope stage 4 (MIS4). Our multi-proxy study of a marine sediment core collected about 60 km southwest of the Outer Hellefisk Moraines demonstrates that in the southern Davis Strait region the most extreme Greenland shelf glaciation of the last glacial cycle occurred during MIS 4, with another prominent glacial advance at 37–33 kyr BP. During those periods the GIS likely reached the Outer Hellefisk Moraines in this area. Except for these two periods, our data suggest significant advection of relatively warm Irminger Sea Water by the West Greenland Current since MIS 4. This advection likely limited the extent of the MIS2 glaciation on the SW Greenland shelf. Decreased precipitation over southwestern Greenland predicted by atmospheric models as a downstream effect of a much larger MIS2 Laurentide Ice Sheet may have played an additional role.

The mass loss from the Greenland Ice Sheet has increased over the past 2 decades. Marine-terminating glaciers contribute significantly to this mass loss due to increased melting and ice discharge. Periods of rapid retreat of these tidewater glaciers have been linked to the concurrent inflow of warm Atlantic-sourced waters. However, little is known about the variability of these Atlantic-derived waters within the fjords, due to a lack of multi-annual in situ measurements. Thus, to better understand the potential role of ocean warming on glacier retreat, reconstructions that characterize the variability of Atlantic water inflow to the fjords are required. Here, we investigate foraminiferal assemblages in a sediment core from Upernavik Fjord, West Greenland, in which the major ice stream Upernavik Isstrøm terminates. We conclude that the foraminiferal assemblage is predominantly controlled by changes in bottom water composition and provide a reconstruction of Atlantic water inflow to Upernavik Fjord, spanning the period 1925–2012. This reconstruction reveals peak Atlantic water influx during the 1930s and again after 2000, a pattern that is comparable to the Atlantic Multidecadal Oscillation (AMO). The comparison of these results to historical observations of front positions of Upernavik Isstrøm reveals that inflow of warm Atlantic-derived waters likely contributed to high retreat rates in the 1930s and after 2000. However, moderate retreat rates of Upernavik Isstrøm also prevailed in the 1960s and 1970s, showing that glacier retreat continued despite a reduced Atlantic water inflow, albeit at a lower rate. Considering the link between bottom water variability and the AMO in Upernavik Fjord, and the fact that a persistent negative phase of the AMO is expected for the next decade, Atlantic water inflow into the fjord may decrease in the coming decade, potentially minimizing or stabilizing the retreat of Upernavik Isstrøm during this time interval.

Benthic foraminiferal populations were studied in a shallow bay of San Salvador Island, the Bahamas. Surface sediments and marine macrophytes were collected from 14 sample sites along a 500 m transect at Grahams Harbour to investigate the foraminiferal assemblage in each microhabitat and to test the link between dead foraminiferal test accumulation patterns and living epiphytic and sedimentary foraminiferal assemblages, macrophyte distribution, and environmental gradients. The analyses include grain size measurements, macrophyte biomass quantification, and qualitative and quantitative studies of benthic foraminifera. The foraminifera found attached to macrophytes differed between macrophyte habitats. However, a correlation between these living communities and the dead assemblages in the sediments at the same sites could not be observed. Principal component analysis (PCA) and redundancy analysis (RDA) suggest that the presence of the macroalgae Halimeda explains 16 % of the residual faunal variation in the dead foraminiferal assemblage after the effects of sorting according to fall speed are partialled out. The RDA also reflects a positive correlation between foraminifera larger than 1.0 mm in diameter and the 0.25–0.5 mm sediment grain size, indicating sedimentological processes as the main factor controlling the sedimentary epiphytic foraminiferal assemblages. These sedimentary processes overprint most effects of ecological features or macrophyte-specific association.

Pockmarks or crater-like depressions formed due to fluid venting, is a well-known geological phenomenon. On a global scale, pockmarks expel a significant amount of fluids and greenhouse gasses to the marine environment and the atmosphere. Despite pockmarks being well-studied, the longevity and history of activity are rarely addressed. In this study, we analyse more than twenty nearshore seafloor depressions of varying morphology and size in the shallow waters of the western Limfjord in northern Denmark. By using high-resolution seismo-acoustic and sediment core data, we characterize the morphology and sediments of the pockmarks. We further use 14C-dating and pore water analyses to address the timing and activity of the depressions. Our results show that the depressions coincide with submarine groundwater discharge areas and hence represent freshwater pockmarks. The pockmarks most likely formed initially due to two eruptive events of thermogenic light oil/gas ca. 2000 and 900 yrs. BP, both triggered by relative sea-level fall. Freshwater seepage continues to sustain and develop the pockmark morphology and we show that pockmarks grow by collapse of pockmark walls and amalgamation of discrete erosive seepage sites at the bottom. The episodic expulsions of thermogenic fluids and the continuous seepage of freshwater probably prolonged the longevity of the studied pockmarks and further demonstrates the complex nature of pockmark formation and activity. While thermogenic and biogenic gas escape to the surface, on the global scale, is relevant for the greenhouse gas budget, the outflow of groundwater which may or may not be polluted is relevant with respect to the geochemistry, biota and productivity of the receiving marine environment especially in shallow water settings.

Evaluating and understanding bio-environmental relationships at different spatial and temporal scales is important for accurately assessing marine productivity in relation to oceanographic conditions in changing Arctic ecosystems. In this paper, we compare four different bio-environmental datasets, in order to shed new light on the complex marine ecosystem of the Baffin Bay-Labrador Sea region, as well as the historical links between environmental changes and the ecosystem. Satellite dataset was used to describe the regional distribution of chlorophyll a, sea ice concentration (SIC), sea surface temperature (SST) and sea surface salinity (SSS). Chlorophyll a data were also analysed together with shrimp assessment dataset to examine the relationship between marine pelagic productivity and bottom feeders off the Greenland coast. A microfossil dataset (diatom valves) from 57 surface sediment stations was used to describe diatom concentration and diversity across the study region. It was also used in combination with satellite data from the same sites to analyse species-environment relationships and develop a diatom-based transfer function for reconstructing past changes in SIC, SST and SSS. Subsequently, we applied the diatom transfer function on three short sediment cores from off SW and W Greenland spanning the last few centuries. A monitoring dataset of sea ice observations and instrumental measurements of water temperature and salinity, together with atmospheric data (air temperature and wind conditions) from Greenland, were used to develop a more complete picture of 20th and 21st century changes in ocean-climate conditions. Our comparative analyses show clear spatial distribution patterns of marine phytoplankton and identify key ecological groups linked to oceanic conditions. In addition, we find a significant correlation between the different trophic levels, culminating in a ‘productivity boom’ in 2001 CE. On a longer temporal scale, the transfer function reconstructions reveal a good correlation with the satellite, monitoring and observational trends, as well as a consistent, long-term decrease in sea ice linked to ocean-climate variability.

Bhutan, located in the Himalayas in the South Asian monsoon region, has extremely high variation in elevation, climatic conditions, and land cover despite its small geographical area, as well as great biodiversity. This paper provides the first comprehensive description of climatic conditions in Bhutan. It assesses the spatial variation of temperature and precipitation across the country and evaluates the causes for this variation based on daily data from 70 meteorological stations that have been recording data for time spans ranging from 3 to 21 years. Temperature and precipitation show contrasting spatial variation, with temperature primarily affected by elevation and precipitation by latitude. Models were developed using mixed linear regression models to predict seasonal and annual mean temperature and precipitation based on geographical location. Using linear regression we found that temperatures changed by about 0.5°C for every 100 m of change in elevation, with lapse rates being highest in February, March, and November and lowest from June to August. The lapse rate was highest for minimum temperatures and lowest for maximum temperatures, with the greatest difference during winter. The spatial distribution of precipitation was mainly controlled by latitude, having a quadratic relationship, with the highest rates in the southern foothills of the Himalayan range and the lowest at midlatitudes. The land cover is affected by topography and local climate, with variations in temperature being a main deciding factor for vegetation types; most human settlements and associated land uses are concentrated at lower elevations.

Arctic sea ice represents an important component of the climate system, and the present reduction of sea ice in the Arctic is of major concern. Despite its importance, little is known about past changes in sea-ice cover and the underlying forcing mechanisms. Here, we use diatom assemblages from a marine sediment core collected from the West Greenland shelf to reconstruct changes in sea-ice cover over the last millennium. The proxy-based reconstruction demonstrates a generally strong link between changes in sea-ice cover and solar variability during the last millennium. Weaker (or stronger) solar forcing may result in the increase (or decrease) in sea-ice cover west of Greenland. In addition, model simulations show that variations in solar activity not only affect local sea-ice formation, but also control the sea-ice transport from the Arctic Ocean through a sea-ice–ocean–atmosphere feedback mechanism. The role of solar forcing, however, appears to have been more ambiguous during an interval around AD 1500, after the transition from the Medieval Climate Anomaly to the Little Ice Age, likely to be driven by a range of factors.

Late- and postglacial geological evolution of the southern Kattegat connection to the Great Belt was investigated from high-resolution seismic data and radiocarbon-dated sediment cores in order to elucidate the Ancylus Lake drainage/Littorina Sea transgression. It was found that glacial deposits form the acoustic basement and are covered by Lateglacial (LG) marine sediments and postglacial (PG; Holocene) material. The LG deposits form a highstand systems tract, whereas the PG deposits cover a full depositional sequence, consisting of a lowstand systems tract (PG I), transgressive systems tract (PG II; subdivided into three parasequences) and finally a highstand systems tract (PG III). PG I sand deposits (11.7–10.8 cal. ka BP) are found in a major western channel and in a secondary eastern channel. PG II (10.8–9.8 cal. ka BP) consists of estuarine and coastal deposits linked to an estuary located at the mouth of the channels. Both channels drained fresh water from south to north. The PG III, that is younger than 9.8 cal. ka BP, represents the threshold marine flooding at the southeastern branch of the palaeo-Great Belt channel. At 9.3 cal. ka BP, fully marine conditions were established, shortly before the flooding of the threshold to the northern part of the Great Belt. These early Holocene spits and sand bars are preserved as features on the present seabed, probably as a result of the rapid sea-level rise that led to back-stepping of the early Holocene palaeo-coast system. This study shows no evidence of major erosion or delta deposition linked to the emptying of the Ancylus Lake, which suggests that continuous water flow from the south characterized the area, without any major drainage event of the Ancylus Lake impacting the southwestern Kattegat.

The last deglaciation was marked by intervals of rapid climatic fluctuations accompanied by glacial advances and retreats along the eastern margin of the Laurentide Ice Sheet. One of these climatic events, the Younger Dryas cold reversal, was accompanied by a detrital carbonate event referred to as “Heinrich event 0” (H0) that deposited ice-rafted debris and especially detrital carbonate-rich (DC) sediment layers in the western and southern Labrador Sea. The precise age, duration, source, and geographical distribution of the H0 DC event, however, are not entirely clear. A high-sedimentation rate sequence cored off southern Newfoundland yielded an age of ~11.5–11.3 calibrated kyr B.P. for this layer, thus pointing to its deposition directly following the Younger Dryas termination, likely associated with rapid ice retreat and warming at the onset of the Holocene. At the study site, the H0 layer contains increased concentrations of detrital carbonates (calcite/dolomite ratio ~ 2.5:1) and a mature biomarker composition, similar to that found in Heinrich Layers 1, 2, 4, and 5, which together suggest an origin from Hudson Strait. Grain size analysis indicates that the H0-sediments were transported to the study site mostly through shallow plumes along the Labrador Shelf, with some hyperpycnal dispersal, deeper, along the slope and rise. Our data thus point to a large meltwater release caused by retreat of the ice margin in the Hudson Strait area. The relvatively short duration of the H0 layer, and its timing right after a major climatic transition, combined with its ubiquitous regional occurrence makes it valuable for correlating different sediment records from the western Labrador Sea.

International audience; The Last Interglacial (LIG) represents an invaluable case study to investigate the response of components of the Earth system to global warming. However, the scarcity of absolute age constraints in most archives leads to extensive use of various stratigraphic alignments to different reference chronologies. This feature sets limitations to the accuracy of the stratigraphic assignment of the climatic sequence of events across the globe during the LIG. Here, we review the strengths and limitations of the methods that are commonly used to date or develop chronologies in various climatic archives for the time span (∼140–100 ka) encompassing the penultimate deglaciation, the LIG and the glacial inception. Climatic hypotheses underlying record alignment strategies and the interpretation of tracers are explicitly described. Quantitative estimates of the associated absolute and relative age uncertainties are provided.Recommendations are subsequently formulated on how best to define absolute and relative chronologies. Future climato-stratigraphic alignments should provide (1) a clear statement of climate hypotheses involved, (2) a detailed understanding of environmental parameters controlling selected tracers and (3) a careful evaluation of the synchronicity of aligned paleoclimatic records. We underscore the need to (1) systematically report quantitative estimates of relative and absolute age uncertainties, (2) assess the coherence of chronologies when comparing different records, and (3) integrate these uncertainties in paleoclimatic interpretations and comparisons with climate simulations.Finally, we provide a sequence of major climatic events with associated age uncertainties for the period 140–105 ka, which should serve as a new benchmark to disentangle mechanisms of the Earth system's response to orbital forcing and evaluate transient climate simulations.

Palynomorph and foraminiferal assemblages have been studied from the Upper Pliensbachian to Bathonian of a borehole section on the island of Anholt in the Kattegat, situated near the eastern margin of the Norwegian-Danish Basin. Palynomorphs were recorded throughout the succession and have been used for both biostratigraphical and palaeoenvironmental assessments. Foraminifera were recorded from only the lower part of the succession, where they proved useful for interpreting palaeoenvironments. Four palynomorph and three foraminiferal zones have been established.The interval encompasses a palaeoenvironmental transition from a marine, inner shelf setting to mainly terrestrial conditions. The Upper Pliensbachian-Toarcian boundary marks the beginning of a major regression, which continued through the Toarcian and Aalenian. Hence, it took place significantly earlier at Anholt than in the centre of the Norwegian-Danish Basin, where a lowering of sea level did not occur until the late Toarcian.The Jurassic succession on Anholt spans the Fjerritslev and Haldager Sand formations; biostratigraphical data indicate that the Lower-Middle Jurassic boundary is here located within the uppermost part of the Fjerritslev Formation. The Fjerritslev and Haldager Sand formational transition was previously considered to coincide with the Lower-Middle Jurassic boundary.

Rising global sea level caused by melting ice sheets poses a major challenge ina persistently warming climate. The Greenland Ice Sheet (GrIS) is among the maincontributors, and in order to make accurate predictions of future ice retreat and sealevel rise, it is imperative to understand how the ice sheet responded to global warming in the past. Reconstructions of relative sea level (RSL) are a key constraint in models of past ice sheet fluctuations, however, high-precision data has until now been sparse in North Greenland. In this study, we present a RSL reconstruction for Finderup Land, North Greenland based on five isolation lakes located between 19.6 and 81.2 m a.s.l. The transition between marine and lacustrine sediments has been identified using XRF, lithological interpretation, and foraminiferal analysis. Age constraints are based on 14C dating of foraminifera and paleomagnetic age correlation. Our results show that Finderup Land was ice free by 10.8 0.2 cal ka BP with a subsequent rapid RSL fall occurring from 9.5 0.2 to 8.0 cal ka BP, at which point the RSL started to approach present level. Furthermore, we establish the marine limit to be minimum at 81.2 m a.s.l. We compare our data to modeled RSL predictions for the area and our results indicate a faster RSL fall, which in turn reflects that the ice retreat was more rapid than estimated and possibly, that the ice sheet in North and Northeast Greenland was larger than previous estimates suggest.

Dinoflagellate cyst assemblages in a marine sediment core from the Iceland Basin were used to carry out qualitative and quantitative assessments of upper ocean conditions in the northern branch of the Subpolar Gyre (SPG) during the Holocene. The data revealed that the early Holocene, i.e. prior to similar to 6.8 kyr B.P., experienced elevated summer sea surface temperatures (SST). In contrast, winter SST was lower compared to modern conditions. This stronger-than-present seasonality revealed by our data thus shows the influence of the early Holocene insolation forcing. Higher seasonality also dominated the reconstructed sea surface salinity (SSS) in the same period, with lower salinity during both summers and winters. The lower SSS suggests freshwater advection into the SPG, and the dinoflagellate cyst content from this interval indicates a dominant sourcing from the Canadian Labrador coast, and thus the eastern margin of the Laurentide Ice Sheet. Within this interval of low SSS, the onset of the Holocene Climate Optimum stands out by its distinct "marine" assemblage components and associated high summer SST and relative SSS increase between 9.3 and 7.8 kyr B.P. A prominent shift in the overall assemblage composition occurs at 6.8 kyr B.P., with lower seasonality in sea surface conditions from that time onwards. Increased summer and winter SSS suggest that this change is linked to the halted influence of meltwater from the Laurentide Ice Sheet, a westward shift of the subpolar front and a contracted SPG allowing an easier northward movement of southern-sourced North Atlantic waters, contributing to and consequently sustained by the spin-up of the overturning circulation. The subpolar front appears to gradually have moved slightly eastward again between 2.8 and 1.9 kyr B.P. towards its modern-day position, with a subtle increase of elements characterized by a typical East Greenland Current signature in the Iceland Basin. Finally, spectral analyses indicate that a range of periodicities may be present in the new SSS and SST records, most notably a 2000-yr periodicity in SSS changes, although none of the periodicities are dominant in both records or throughout the time interval covered by the records.

An edited version of this paper was published by AGU. Copyright 2018 American Geophysical Union. Limoges, A., Ribeiro, S., Weckström, K., Heikkilä, M., Zamelczyk, K., Andersen, T.J, ... Seidenkrantz, M. (2018). Linking the modern distribution of biogenic proxies in High Arctic Greenland shelf sediments to sea ice, primary production and Arctic-Atlantic inflow. Journal of Geophysical Research: Biogeosciences, 123(3), 760-786. https://doi.org/10.1002/2017JG003840. To view the published open abstract, go to https://doi.org/10.1002/2017JG003840. The eastern north coast of Greenland is considered to be highly sensitive to the ongoing Arctic warming, but there is a general lack of data on modern conditions and in particular on the modern distribution of climate and environmental proxies to provide a baseline and context for studies on past variability. Here we present a detailed investigation of 11 biogenic proxies preserved in surface sediments from the remote High Arctic Wandel Sea shelf, the entrance to the Independence, Hagen, and Danmark fjords. The composition of organic matter (organic carbon, C:N ratios, δ13C, δ15N, biogenic silica, and IP25) and microfossil assemblages revealed an overall low primary production dominated by benthic diatoms, especially at the shallow sites. While the benthic and planktic foraminiferal assemblages underline the intrusion of chilled Atlantic waters into the deeper parts of the study area, the distribution of organic‐walled dinoflagellate cysts is controlled by the local bathymetry and sea ice conditions. The distribution of the dinoflagellate cyst Polarella glacialis matches that of seasonal sea ice and the specific biomarker IP25, highlighting the potential of this species for paleo sea ice studies. The information inferred from our multiproxy study has important implications for the interpretation of the biogenic‐proxy signal preserved in sediments from circum‐Arctic fjords and shelf regions and can serve as a baseline for future studies. This is the first study of its kind in this area.

Atmospheric circulation may change with future climate change in response to modification of meridional temperature gradients, but the potential influence on ocean circulation is as yet unclear. Over the mid-late Holocene, atmospheric circulation in the North Atlantic region has fluctuated on millennial timescales; therefore, the ocean response to these changes can be investigated using the paleoceanographic records that have been developed in the north-eastern subpolar North Atlantic. Here, we present a diatom-based sea surface temperature reconstruction from the Iceland Basin, south of Iceland; the reconstruction shows the warmest temperatures of the record at 6.1–4 ka BP, cooler temperatures at 4-2 ka BP and warmer temperatures thereafter. Inter-record comparisons indicate that the cold period at c. 4-2 ka BP may have resulted from a strengthened East Greenland Current and/or melting of the Greenland ice sheet, in response to a negative North Atlantic Oscillation. The findings highlight that atmospheric circulation changes are likely to cause pronounced variations in the latitudinal exchange of heat, which may have consequences for deep-water formation and global ocean circulation.

Atlantic Water (AW) advection plays an important role in climatic, oceanographic and environmental conditions in the eastern Arctic. Situated along the only deep connection between the Atlantic and the Arctic oceans, the Svalbard Archipelago is an ideal location to reconstruct the past AW advection history and document its linkage with local glacier dynamics, as illustrated in the present study of a 275 cm long sedimentary record from Woodfjorden (northern Spitsbergen; water depth: 171 m) spanning the last ∼ 15 500 years. Sedimentological, micropalaeontological and geochemical analyses were used to reconstruct changes in marine environmental conditions, sea ice cover and glacier activity. Data illustrate a partial break-up of the Svalbard–Barents Sea Ice Sheet from Heinrich Stadial 1 onwards (until ∼ 14.6 ka). During the Bølling–Allerød ( ∼ 14.6–12.7 ka), AW penetrated as a bottom water mass into the fjord system and contributed significantly to the destabilization of local glaciers. During the Younger Dryas ( ∼ 12.7–11.7 ka), it intruded into intermediate waters while evidence for a glacier advance is lacking. A short-term deepening of the halocline occurred at the very end of this interval. During the early Holocene ( ∼ 11.7–7.8 ka), mild conditions led to glacier retreat, a reduced sea ice cover and increasing sea surface temperatures, with a brief interruption during the Preboreal Oscillation ( ∼ 11.1–10.8 ka). Due to a ∼ 6000-year gap, the mid-Holocene is not recorded in this sediment core. During the late Holocene ( ∼ 1.8–0.4 ka), a slightly reduced AW inflow and lower sea surface temperatures compared to the early Holocene are reconstructed. Glaciers, which previously retreated to the shallower inner parts of the Woodfjorden system, likely advanced during the late Holocene. In particular, topographic control in concert with the reduced summer insolation partly decoupled glacier dynamics from AW advection during this recent interval.

Benthic foraminiferal assemblages supported by selected geochemical data from three marine sediment cores collected in Placentia Bay, SE Newfoundland, are used to construct an ~13,000-year-long record of regional oceanographic changes in the SW Labrador Sea. The area is located in the boundary zone between the cold, ice-loaded Labrador Current (LC) in the north and the warm Gulf Stream (GS) waters to the south. After the Younger Dryas termination, the influence of GS-derived water increased and was further strengthened at 10.7 cal. kyr BP through enhanced northward flow of Atlantic water via the Slopewater Current. A short-term event of increased terrestrial input and water column stratification at 8.4 cal. kyr BP was likely linked to the distal drainage of glacial Lake Agassiz. After 7.3 cal. kyr BP, a stronger LC weakened the inflow of warmer subsurface waters from the GS. This may be explained by extensive meltwater release from ice sheets in Arctic Canada and is concurrent with a general shift in oceanographic conditions in the Labrador Sea region. Around 4.0 cal. kyr BP, conditions became more stable with a slight increase in salinity, indicating a decrease in meltwater transported via the LC. The Northern Hemisphere neoglacial cooling around 2.8 cal. kyr BP was characterized off SE Newfoundland by a further stabilization of the current system, dominated by the LC with some continued influx of GS water.

Mounting evidence from proxy records suggests that variations in solar activity have played a significant role in triggering past climate changes. However, the mechanisms for sun-climate links remain a topic of debate. Here we present a high-resolution summer sea-surface temperature (SST) record covering the past 9300 yr from a site located at the present-day boundary between polar and Atlantic surface-water masses. The record is age constrained via the identification of 15 independently dated tephra markers from terrestrial archives, circumventing marine reservoir age variability problems. Our results indicate a close link between solar activity and SSTs in the northern North Atlantic during the past 4000 yr; they suggest that the climate system in this area is more susceptible to the influence of solar variations during cool periods with less vigorous ocean circulation. Furthermore, the high-resolution SST record indicates that climate in the North Atlantic regions follows solar activity variations on multidecadal to centennial time scales.

Køge Bugt, in southeast Greenland, hosts three of the largest glaciers of the Greenland Ice Sheet; these have been major contributors to ice loss in the last two decades. Despite its importance, the Holocene history of this area has not been investigated. We present a 9100 year sediment core record of glaciological and oceanographic changes from analysis of foraminiferal assemblages, the abundance of ice-rafted debris, and sortable silt grain size data. Results show that ice-rafted debris accumulated constantly throughout the core; this demonstrates that glaciers in Køge Bugt remained in tidewater settings throughout the last 9100 years. This observation constrains maximum Holocene glacier retreat here to less than 6 km from present-day positions. Retreat was minimal despite oceanic and climatic conditions during the early-Holocene that were at least as warm as the present-day. The limited Holocene retreat of glaciers in Køge Bugt was controlled by the subglacial topography of the area; the steeply sloping bed allowed glaciers here to stabilise during retreat. These findings underscore the need to account for individual glacier geometry when predicting future behaviour. We anticipate that glaciers in Køge Bugt will remain in stable configurations in the near-future, despite the predicted continuation of atmospheric and oceanic warming.

Gennem sin historie har Jorden gennemlevet mange forandringer i klimaet. Hvordan kan vi forklare disse klimavariationer, og hvordan kan man i dag skelne mellem naturlige og menneskeskabte klimaforandringer, når nu Jordens klima ikke er stabilt? For at kunne forstå dette, skal vi forstå det komplicerede samspil mellem processer, der styrer klimaet

The study of active microbial populations in deep, energy-limited marine sediments has extended our knowledge of the limits of life on Earth. Typically, microbial activity in the deep biosphere is calculated by transport-reaction modelling of pore water solutes or from experimental measurements involving radiotracers. Here we modelled microbial activity from the degree of D:L-aspartic acid racemization in microbial necromass (remains of dead microbial biomass) in sediments up to ten million years old. This recently developed approach (D:L-amino acid modelling) does not require incubation experiments and is highly sensitive in stable, low-activity environments. We applied for the first time newly established constraints on several important input parameters of the D:L-amino acid model, such as a higher aspartic acid racemization rate constant and a lower cell-specific carbon content of sub-seafloor microorganisms. Our model results show that the pool of necromass amino acids is turned over by microbial activity every few thousand years, while the turnover times of vegetative cells are in the order of years to decades. Notably, microbial turnover times in million-year-old sediment from the Peru Margin are up to 100-fold shorter than previous estimates, highlighting the influence of microbial activities on element cycling over geologic time scales.

Source at https://doi.org/10.1007/s13280-016-0894-2. In order to establish a baseline for proxy-based reconstructions for the Young Sound–Tyrolerfjord system (Northeast Greenland), we analysed the spatial distribution of primary production and sea ice proxies in surface sediments from the fjord, against monitoring data from the Greenland Ecosystem Monitoring Programme. Clear spatial gradients in organic carbon and biogenic silica contents reflected marine influence, nutrient availability and river-induced turbidity, in good agreement with in situ measurements. The sea ice proxy IPz25 was detected at all sites but at low concentrations, indicating that IP25 records from fjords need to be carefully considered and not directly compared to marine settings. The sea ice-associated biomarker HBI III revealed an open-water signature, with highest concentrations near the mid-July ice edge. This proxy evaluation is an important step towards reliable palaeoenvironmental reconstructions that will, ultimately, contribute to better predictions for this High Arctic ecosystem in a warming climate.

Reproducible climate reconstructions of the Common Era (1 CE to present) are key to placing industrial-era warming into the context of natural climatic variability. Here we present a community-sourced database of temperature-sensitive proxy records from the PAGES2k initiative. The database gathers 692 records from 648 locations, including all continental regions and major ocean basins. The records are from trees, ice, sediment, corals, speleothems, documentary evidence, and other archives. They range in length from 50 to 2000 years, with a median of 547 years, while temporal resolution ranges from biweekly to centennial. Nearly half of the proxy time series are significantly correlated with HadCRUT4.2 surface temperature over the period 1850-2014. Global temperature composites show a remarkable degree of coherence between high- and low-resolution archives, with broadly similar patterns across archive types, terrestrial versus marine locations, and screening criteria. The database is suited to investigations of global and regional temperature variability over the Common Era, and is shared in the Linked Paleo Data (LiPD) format, including serializations in Matlab, R and Python., PAGES, a core project of Future Earth, is supported by the U.S. and Swiss National Science Foundations. Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government. Some of this work was conducted as part of the North America 2k Working Group supported by the John Wesley Powell Center for Analysis and Synthesis, funded by the U.S. Geological Survey. B. Bauer, W. Gross, and E. Gille (NOAA National Centers for Environmental Information) are gratefully acknowledged for helping assemble the data citations and creating the NCEI versions of the PAGES 2k data records. We thank all the investigators whose commitment to data sharing enables the open science ethos embodied by this project.

In total, 39 diatom species belonging to 22 genera were identified and photographed from Holocene marine sediment cores of Eastern Newfoundland bays. The cores were retrieved from Placentia Bay in the southeast and Bonavista Bay in the northeast of Newfoundland. The study area lies at the meeting point of the cold Labrador Current and the warmer Gulf Stream and is sensitive to changes in ocean circulation. It is thus an ideal location for paleoceanographic reconstructions based on the analysis of (sub)fossil diatoms, for which a good taxonomic framework is essential. The studied sediment cores encompass different age ranges and together represent the entire Holocene epoch. The most abundant species belonged to the genera Thalassiosira, represented by 10 different species, Fragilariopsis (four species), Detonula, Thalassionema and Odontella (one species each). This study provides an annotated list of the most commonly identified diatom species and references to further literature. All taxa are documented b...

A diatom-based sea-ice concentration (SIC) transfer function was developed by using 72 surface samples from west of Greenland and around Iceland, and validated against associated modem SIC. Canonical correspondence analysis on surface sediment diatoms and monthly average of SIC indicated that April SIC is the most important environmental factor controlling the distribution of diatoms in the area, justifying the development of a diatombased SIC transfer function. The agreement between reconstructed SIC based on diatoms from West Greenland and the satellite and modelled sea-ice data during the last -75 yr suggests that the diatom-based SIC reconstruction is reliable for studying the palaeoceanography off West Greenland.Relatively warm conditions with a strong influence of the Irminger Current (IC) were indicated for the early part of the record (-5000-3860 cal. yr BP), corresponding in time to the latest part of the Holocene Thermal Maximum. Between 3860 and 1510 cal. yr BP, April SIC oscillated around the mean value (55%) and during the time interval 1510-1120 cal. yr BP and after 650 cal. yr BP was above the mean, indicating more extensive sea-ice cover in Disko Bugt.Agreement between reconstructed April SIC and changes in the diatom species suggests that the sea-ice condition in Disko Bugt was strongly influenced by variations in the relative strength of two components of the West Greenland Current, i.e. the cold East Greenland Current and the relatively warm IC. Further analysis of the reconstructed SIC record suggests that solar radiation may be an important forcing mechanism behind the historic sea-ice changes. (C) 2014 Elsevier B.V. All rights reserved.

We present a new reconstruction of summer sea-surface salinity (SSS) over the past 15 000 years based on a diatom record from piston core 17940, located on the northern slope of the South China Sea (SCS). The reconstructed diatom-based summer SSS values for the modern period are in accord with instrumental observations of summer SSS in the area. Here, the modern summer SSS is primarily controlled by river runoff, in particular from the Pearl River. The reconstruction presented in this study shows that the summer SSS varied between 33.3 and 34.2 psu over the past 15 000 years. The long-term summer SSS trend closely followed the trend of the orbitally controlled solar insolation at 20°N, suggesting that orbital forcing was the dominant driver of changes in summer SSS in this area. Comparisons to speleothem δ18O data and studies of surface hydrography in the region suggest that changes in solar insolation affected the summer SSS through changes in the East Asian Monsoon and sea-level changes associated with the last deglaciation. Univariate spectral analyses indicate that centennial-scale oscillatory variations in summer SSS were superimposed on the long-term trend. During the deglacial period (c. 12 000–9000 cal. a BP), the dominant periodicity was centred around 230–250 years, whereas a ∼350-year oscillation dominated in the period 2200–4500 cal. a BP. The balance of evidence suggests that these centennial-scale changes in summer SSS may have been driven by solar-induced changes in the East Asian Monsoon, but further evidence is needed to firmly establish this relationship.

Sea ice is an important component of the Earth system with complex dynamics imperfectly documented from direct observations, which are primarily limited to the last 40 years. Whereas large amplitude variations of sea ice have been recorded, especially in the Arctic, with a strikingly fast decrease in recent years partly attributed to the impact of anthropogenic climate changes, little is known about the natural variability of the sea ice cover at multi-decadal to multi-millennial time scales. Hence, there is a need to establish longer sea ice time series to document the full range of sea ice variations under natural forcings. To do this, several approaches based on biogenic or geochemical proxies have been developed from marine, ice core and coastal records. The status of the sea ice proxies has been discussed by the Sea Ice Proxy (SIP) working group endorsed by PAGES during a first workshop held at GEOTOP in Montréal. The present volume contains a set of papers addressing various sea ice proxies and their application to large scale sea ice reconstruction. Here we summarize the contents of the volume, including a table of various proxies available in marine sediments and ice cores, with their possibilities and limitations.

A high-resolution diatom record from site MD05-2908 in the Southern Okinawa Trough, East China Sea, reveals pronounced multidecadal- to centennial-scale palaeoceanographic changes throughout the last millennium. Summer sea-surface salinity (SSS) was reconstructed using a weighted averaging partial least squares diatom-based training set. The reconstructed SSS shows slightly decreasing values during the period AD 9051930 with considerable fluctuations superimposed on this general trend. Relatively high-salinity conditions during the interval AD 9051450 probably suggest a low flood frequency in north-eastern Taiwan. Furthermore, the high SSS values are associated with a strong and stable influence of the Kuroshio Current on the Southern Okinawa Trough during the Medieval Climate Anomaly. The period AD 14501930 is characterized by three low-salinity intervals (AD 14501500, AD 16251725 and AD 17701880) separated by periods of relatively high salinity. The low SSS intervals indicate increased freshwater discharge into the Southern Okinawa Trough during the Little Ice Age, probably as a result of higher flood frequencies in north-eastern Taiwan. Spectral and wavelet analyses suggest that this pattern was linked to multidecadal variations in summer SSS, presumably associated with the Pacific Decadal Oscillation. Copyright (c) 2012 John Wiley & Sons, Ltd.

The deglacial history and oceanography of Uummannaq Trough, central West Greenland continental shelf, was investigated using foraminiferal, sedimentological, and bathymetric records together with a radiocarbon chronology, providing a timeline for the retreat of glacial ice after the Last Glacial Maximum (LGM). To map ice stream retreat, data were collected from cores from the outer (JR175-VC45 and JR175-VC43) and inner (JR175-VC42) Uummannaq Trough. A large ice stream, fed by confluent glaciers draining the interior of the Greenland Ice Sheet, extended across the outer shelf during the LGM and was in retreat by 15.0 cal kyr BP. Foraminiferal data indicate that the ‘warm’ West Greenland Current (WGC) was established prior to 14.0 cal kyr BP, which is the hitherto earliest record of Atlantic Water found on the West Greenland shelf. For each of the cores, foraminifera indicate that ice sheet retreat was followed quickly by incursion of the WGC, suggesting that the warm water may have enhanced ice retreat. Prior to the Younger Dryas cold event, the radiocarbon chronology indicates that the ice sheet retreated to the mid-shelf, where it subsequently stabilised and formed a large grounding-zone wedge (GZW). After the Younger Dryas, around 11.5 cal kyr BP, the ice retreated rapidly from the GZW and into the fjords.

Bhutan is located in the Himalayas with extreme variation in elevation, climatic conditions and land use. The high dependency of the economy on agriculture and natural resources emphasizes the importance of understanding inter- and intra-seasonal variation in water balance linked to monsoonal precipitation and evapotranspiration. We used data from 71 meteorological stations across Bhutan, each encompassing data series ranging from a few years till two decades. The temperature-based Hargreaves (H) equation for reference evapotranspiration (ETo) calculation was calibrated in reference to the FAO Penman-Monteith (PM) equation giving good estimates of average monthly ETo. Various published versions of the H equations consistently computed average ETo around 32% higher than the average ETo computed by PM. The difference could be largely attributed to higher mean relative humidity and a lower average wind speed compared to global average conditions. The spatial variation in PM ETo is linked to variation in sunshine hours during the summer season and to temperature in other seasons. Seasonal and annual ETo was mainly affected by elevation and latitude, which is linked to temperature and sunshine duration. Comparison with the Southern Oscillation Index (SOI) provided a positive correlation of both precipitation and water balance to SOI, while there was a negative correlation between ETo and SOI. When testing our models for predicting ETo and water balances versus the global CRU gridded data, our models performed clearly better than the global CRU gridded data for Bhutan, indicating that local higher-resolution data and models are preferred in mountainous regions. A positive water balance is found in across east-west belts in southern and northern Bhutan, where mainly broadleaf forest types dominate, whereas an either small or negative water balance in the east-west belt of central Bhutan cause dry conditions in this region, which is dominated by coniferous forests.