Source data for Plos ONE publication Redefining thermal regimes to design reserves for coral reefs in the face of climate change

Reef managers cannot fight global warming through mitigation at local scale, but they can use information on thermal patterns to plan for reserve networks that maximize the probability of persistence of their reef system. Here we assess previous methods for the design of reserves for climate change and present a new approach to prioritize areas for conservation that leverages the most desirable properties of previous approaches. The new method moves the science of reserve design for climate change a step forwards by: (1) recognizing the role of seasonal acclimation in increasing the limits of environmental tolerance of corals and ameliorating the bleaching response; (2) including information from several bleaching events, which frequency is likely to increase in the future; (3) assessing relevant variability at country scales, where most management plans are carried out. We demonstrate the method in Honduras, where a reassessment of the marine spatial plan is in progress.

Rocks of the Western Kamchatka-Koryak continental margin volcanogenic belt: age, chemical and mineral compositions

An isotope-geochemical study of Eocene-Oligocene magmatic rocks from the Western Kamchatka-Koryak volcanogenic belt revealed lateral heterogeneity of mantle magma sources in its segments: Western Kamchatka, Central Koryak, and Northern Koryak ones. In the Western Kamchatka segment magmatic melts were generated from isotopically heterogeneous (depleted and/or insignificantly enriched) mantle sources significantly contaminated by quartz-feldspathic sialic sediments; higher 87Sr/86Sr (0.70429-0.70564) and lower 143Nd/144Nd [eNd(T) = 0.06-2.9] ratios in volcanic rocks from the Central Koryak segment presumably reflect contribution of an enriched mantle source; high positive eNd(T) and low 87Sr/86Sr ratios in magmatic rocks from the Northern Koryak segment area indicate their derivation from an isotopically depleted mantle source without significant contamination by sialic or mantle material enriched in radiogenic Sr and Nd. Significantly different contamination histories of Eocene-Oligocene mantle magmas in Kamchatka and Koryakia are related to their different thermal regimes: higher heat flow beneath Kamchatka led to crustal melting and contamination of mantle suprasubduction magmas by crustal melts. Cessation of suprasubduction volcanism in the Western Kamchatka segment of the continental margin belt was possibly related to accretion of the Achaivayam-Valagin terrane 40 Ma ago, whereas suprasubduction activity in the Koryak segment stopped due to closure of the Ukelayat basin in Oligocene.

Lithium sources in pore fluids of Peru slope

High Li concentrations, up to a maximum of 1155 µM are observed in the pore fluids of the Peru convergent margin slope sediments. At Ocean Drilling Program Sites 683 and 685 (ca. 9°S), the Li concentration depth gradients are twice as steep as at Site 682 and 688 (ca. 11°S). Within the sediments, the most important Li sources are from aluminosilicate minerals. Biogenic opal-A contains little Li and thus dilutes the Li concentration of the bulk sediments. The sediment compositions and the thermal regimes are similar at 9° and 11°S, suggesting there is an additional, non-sedimentary source for the observed high Li concentrations in the northern pore fluids. At 9°S, the 87Sr/86Sr ratios reach a maximum value of 0.709958. The observed radiogenic 87Sr/86Sr values in the pore fluids support the suggestion that the additional Li may derive from exchange reactions with underlying continental crust. The high concentrations of Li at 11°S may derive from basalt alteration at moderate to high temperatures, as suggested by the non-radiogenic 87Sr/86Sr ratios in these pore fluids, which reach a minimum value of 0.707218. Based on (1) Li concentrations in the pore fluids in slope sediments from Peru and several other margins, and (2) an approximate estimate of fluid flux from continental margins into the ocean, continental margins provide an estimated 1 to 3 * 10**10 moles Li/yr to the ocean. This source of oceanic Li, which has not been considered previously, is of the same order of magnitude as some estimates of hydrothermal and river Li fluxes and may have important consequences for the oceanic Li isotope budget. The sink is unknown for this newly discovered and possibly large Li source, but it may be more pervasive low-temperature alteration of oceanic basement than previously estimated, or burial of mineral phases, such as authigenic clay minerals, or metal oxyhydroxides which may be Li-rich.

Temperature, water level and bathymetry of thermokarst lakes in the continuous permafrost zone of northern Siberia - Lena River Delta, Siberia

Thermokarst lakes are typical features of the northern permafrost ecosystems, and play an important role in the thermal exchange between atmosphere and subsurface. The objective of this study is to describe the main thermal processes of the lakes and to quantify the heat exchange with the underlying sediments. The thermal regimes of five lakes located within the continuous permafrost zone of northern Siberia (Lena River Delta) were investigated using hourly water temperature and water level records covering a 3-year period (2009-2012), together with bathymetric survey data. The lakes included thermokarst lakes located on Holocene river terraces that may be connected to Lena River water during spring flooding, and a thermokarst lake located on deposits of the Pleistocene Ice Complex. Lakes were covered by ice up to 2 m thick that persisted for more than 7 months of the year, from October until about mid-June. Lake-bottom temperatures increased at the start of the ice-covered period due to upward-directed heat flux from the underlying thawed sediment. Prior to ice break-up, solar radiation effectively warmed the water beneath the ice cover and induced convective mixing. Ice break-up started at the beginning of June and lasted until the middle or end of June. Mixing occurred within the entire water column from the start of ice break-up and continued during the ice-free periods, as confirmed by the Wedderburn numbers, a quantitative measure of the balance between wind mixing and stratification that is important for describing the biogeochemical cycles of lakes. The lake thermal regime was modeled numerically using the FLake model. The model demonstrated good agreement with observations with regard to the mean lake temperature, with a good reproduction of the summer stratification during the ice-free period, but poor agreement during the ice-covered period. Modeled sensitivity to lake depth demonstrated that lakes in this climatic zone with mean depths > 5 m develop continuous stratification in summer for at least 1 month. The modeled vertical heat flux across the bottom sediment tends towards an annual mean of zero, with maximum downward fluxes of about 5 W/m**2 in summer and with heat released back into the water column at a rate of less than 1 W/m**2 during the ice-covered period. The lakes are shown to be efficient heat absorbers and effectively distribute the heat through mixing. Monthly bottom water temperatures during the ice-free period range up to 15 °C and are therefore higher than the associated monthly air or ground temperatures in the surrounding frozen permafrost landscape. The investigated lakes remain unfrozen at depth, with mean annual lake-bottom temperatures of between 2.7 and 4 °C.

High-resolution dynamical downscaling products for the North Slope of Alaska and surrounding areas, links to model result files

Climatic changes are most pronounced in northern high latitude regions. Yet, there is a paucity of observational data, both spatially and temporally, such that regional-scale dynamics are not fully captured, limiting our ability to make reliable projections. In this study, a group of dynamical downscaling products were created for the period 1950 to 2100 to better understand climate change and its impacts on hydrology, permafrost, and ecosystems at a resolution suitable for northern Alaska. An ERA-interim reanalysis dataset and the Community Earth System Model (CESM) served as the forcing mechanisms in this dynamical downscaling framework, and the Weather Research & Forecast (WRF) model, embedded with an optimization for the Arctic (Polar WRF), served as the Regional Climate Model (RCM). This downscaled output consists of multiple climatic variables (precipitation, temperature, wind speed, dew point temperature, and surface air pressure) for a 10 km grid spacing at three-hour intervals. The modeling products were evaluated and calibrated using a bias-correction approach. The ERA-interim forced WRF (ERA-WRF) produced reasonable climatic variables as a result, yielding a more closely correlated temperature field than precipitation field when long-term monthly climatology was compared with its forcing and observational data. A linear scaling method then further corrected the bias, based on ERA-interim monthly climatology, and bias-corrected ERA-WRF fields were applied as a reference for calibration of both the historical and the projected CESM forced WRF (CESM-WRF) products. Biases, such as, a cold temperature bias during summer and a warm temperature bias during winter as well as a wet bias for annual precipitation that CESM holds over northern Alaska persisted in CESM-WRF runs. The linear scaling of CESM-WRF eventually produced high-resolution downscaling products for the Alaskan North Slope for hydrological and ecological research, together with the calibrated ERA-WRF run, and its capability extends far beyond that. Other climatic research has been proposed, including exploration of historical and projected climatic extreme events and their possible connections to low-frequency sea-atmospheric oscillations, as well as near-surface permafrost degradation and ice regime shifts of lakes. These dynamically downscaled, bias corrected climatic datasets provide improved spatial and temporal resolution data necessary for ongoing modeling efforts in northern Alaska focused on reconstructing and projecting hydrologic changes, ecosystem processes and responses, and permafrost thermal regimes. The dynamical downscaling methods presented in this study can also be used to create more suitable model input datasets for other sub-regions of the Arctic.

Helium and oxygen isotope ratios of ODP Sites 193-1188 and 193-1189

Fluid mixing processes and thermal regimes within the Snowcap and Roman Ruins vent sites of the PACMANUS hydrothermal system, Papua New Guinea, were investigated using 3He/4He ratios from fluid inclusions within pyrite and anhydrite and the d18O signature of anhydrite. Depressed 3He/4He ratios of 0.2-6.91RA appear to be caused by significant atmospheric diffusive exchange, whilst He-Ne diffusive fractionation precludes correction using 20Ne. 40Ar/36Ar ratios of 295-310 are elevated above seawater, indicating the majority of argon is seawater derived but with a magmatic component. d18O anhydrite ratios are 6.5 per mil to 11 per mil for Snowcap and 6.4 per mil to 11.9 per mil for Roman Ruins. Using oxygen isotope fractionation factors for the anhydrite-water system, the temperatures calculated assuming isotopic equilibrium at depth are up to 100 °C cooler than fluid inclusion trapping temperatures. It is likely that anhydrite is precipitated rapidly, preventing d18O equilibration. By comparing new d18O values for anhydrite with corresponding published 87Sr/86Sr ratios, seawater is inferred to penetrate deep into the Snowcap system with little conductive heating. A simple fluid mixing model has been constructed whereby the differing venting styles can be explained by a plumbing system at depth which favors delivery of end-member hydrothermal fluid to the high temperature sites.

Theoretical model simulations for variations in thermal reaction norms

Thermal reaction norms for growth rates of six Emiliania huxleyi isolates originating from the central Atlantic (Azores, Portugal) and five isolates from the coastal North Atlantic (Bergen, Norway) were assessed. We used the template mode of variation model to decompose variations in growth rates into modes of biological interest: vertical shift, horizontal shift, and generalist-specialist variation. In line with the actual habitat conditions, isolates from Bergen (Bergen population) grew well at lower temperatures, and isolates from the Azores (Azores population) performed better at higher temperatures. The optimum growth temperature of the Azores population was significantly higher than that of the Bergen population. Neutral genetic differentiation was found between populations by microsatellite analysis. These findings indicate that E. huxleyi populations are adapted to local temperature regimes. Next to between-population variation, we also found variation within populations. Genotype-by-environment interactions resulted in the most pronounced phenotypic differences when isolates were exposed to temperatures outside the range they naturally encounter. Variation in thermal reaction norms between and within populations emphasizes the importance of using more than one isolate when studying the consequences of global change on marine phytoplankton. Phenotypic plasticity and standing genetic variation will be important in determining the potential of natural E. huxleyi populations to cope with global climate change.