ALBERT

Description: Predicting future thaw slump activity requires a sound understanding of the atmospheric drivers and geomorphic controls on mass wasting across a range of time scales. On sub-seasonal time scales, sparse measurements indicate that mass wasting at active slumps is often limited by the energy available for melting ground ice, but other factors such as rainfall or the formation of an insulating veneer may also be relevant. To study the sub-seasonal drivers, we derive topographic changes from single-pass radar interferometric data acquired by the TanDEM-X satellite (12 m resolution). The high vertical precision (around 30 cm), frequent observations (11 days) and large coverage (5000 km2) allow us to track volume losses as drivers such as the available energy change during summer in two study regions. We find that thaw slumps in the Tuktoyaktuk coastlands, Canada, are not energy limited in June, as they undergo limited mass wasting (height loss of around 0 cm/day) despite the ample available energy, indicating the widespread presence of an insulating snow or debris veneer. Later in summer, height losses generally increase (around 3 cm/day), but they do so in distinct ways. For many slumps, mass wasting tracks the available energy, a temporal pattern that is also observed at coastal yedoma cliffs on the Bykovsky Peninsula, Russia. However, the other two common temporal trajectories are asynchronous with the available energy, as they track strong precipitation events or show a sudden speed-up in late August, respectively. The observed temporal patterns are poorly related to slump characteristics like the slump area. The contrasting temporal behaviour of nearby thaw slumps highlights the importance of complex local and temporally varying controls on mass wasting.

Description: This paper investigates new observations from the poorly understood region between the Kara and Laptev Seas in the Eastern Arctic Ocean. We discuss relevant circulation features including riverine freshwater, Atlantic-derived water, and polynya-formed dense water, emphasize Vilkitsky Strait (VS) as an important Kara Sea gateway, and analyze the role of the adjacent 250 km-long submarine Vilkitsky Trough (VT) for the Arctic boundary current. Expeditions in 2013 and 2014 operated closely spaced hydrographic transects and 1 year-long oceanographic mooring near VT’s southern slope, and found persistent annually averaged flow of 0.2 m s21 toward the Nansen Basin. The flow is nearly barotropic from winter through early summer and becomes surface intensified with maximum velocities of 0.35 m s21 from August to October. Thermal wind shear is maximal above the southern flank at 30 m depth, in agreement with basinward flow above VT’s southern slope. The subsurface features a steep front separating warm (–0.58C) Atlantic-derived waters in central VT from cold (〈–1.58C) shelf waters, which episodically migrates across the trough indicated by current reversals and temperature fluctuations. Shelf-transformed waters dominate above VT’s slope, measuring near-freezing temperatures throughout the water column at salinities of 34–35. These dense waters are vigorously advected toward the Eurasian Basin and characterize VT as a conduit for near-freezing waters that could potentially supply the Arctic Ocean’s lower halocline, cool Atlantic water, and ventilate the deeper Arctic Ocean. Our observations from the northwest Laptev Sea highlight a topographically complex region with swift currents, several water masses, narrow fronts, polynyas, and topographically channeled storms.

Description: The mechanisms causing widespread flow acceleration of Jakobshavn Isbræ, West Greenland, remain unclear despite an abundance of observations and modeling studies. Here we simulate the glacier's evolution from 1985 to 2016 using a three-dimensional thermomechanical ice flow model. The model captures the timing and 90% of the observed changes by forcing the calving front. Basal drag in the trough is low, and lateral drag balances the ice stream's driving stress. The calving front position is the dominant control on changes of Jakobshavn Isbræ since the ice viscosity in the shear margins instantaneously drops in response to the stress perturbation caused by calving front retreat, which allows for widespread flow acceleration. Gradual shear margin warming contributes 5 to 10% to the total acceleration. Our simulations suggest that the glacier will contribute to eustatic sea level rise at a rate comparable to or higher than at present.

Description: Greenland's bed topography is a primary control on ice flow, grounding line migration, calving dynamics and subglacial drainage. Moreover, fjord bathymetry regulates the penetration of warm Atlantic Water (AW) that rapidly melts and undercuts Greenland's marine-terminating glaciers. Here, we present a new compilation of Greenland bed topography that assimilates seafloor bathymetry and ice thickness data through a mass conservation (MC) approach. A new 150-m horizontal resolution bed topography/bathymetric map of Greenland is constructed with seamless transitions at the ice/ocean interface, yielding major improvements over previous datasets, particularly in the marine-terminating sectors of northwest and southeast Greenland. Our map reveals the total sea level potential of the Greenland Ice Sheet is 7.42±0.05 m, which is 7 cm greater than previous estimates. Furthermore, it explains recent calving front response of numerous outlet glaciers and reveals new pathways by which AW can access glaciers with marine-based basins, thereby highlighting sectors of Greenland that are most vulnerable to future oceanic forcing.

Description: We present a simulation of Antarctic iceberg drift and melting that includes small, medium-sized, and giant tabular icebergs with a realistic size distribution. For the first time, an iceberg model is initialized with a set of nearly 7000 observed iceberg positions and sizes around Antarctica. The study highlights the necessity to account for larger and giant icebergs in order to obtain accurate melt climatologies. We simulate drift and lateral melt using iceberg-draft averaged ocean currents, temperature, and salinity. A new basal melting scheme, originally applied in ice shelf melting studies, uses in situ temperature, salinity, and relative velocities at an iceberg's bottom. Climatology estimates of Antarctic iceberg melting based on simulations of small (≤ 2.2 km), 'small-to-medium'-sized (≤ 10 km), and small-to-giant icebergs (including icebergs 〉 10 km) exhibit differential characteristics: successive inclusion of larger icebergs leads to a reduced seasonality of the iceberg meltwater flux and a shift of the mass input to the area north of 58 °S, while less meltwater is released into the coastal areas. This suggests that estimates of meltwater input solely based on the simulation of small icebergs introduce a systematic meridional bias; they underestimate the northward mass transport and are, thus, closer to the rather crude treatment of iceberg melting as coastal runoff in models without an interactive iceberg model. Future ocean simulations will benefit from the improved meridional distribution of iceberg melt, especially in climate change scenarios where the impact of iceberg melt is likely to increase due to increased calving from the Antarctic ice sheet.

Description: A suite of oxygenated volatile organic compounds (OVOCs – acetaldehyde, acetone, propanal, butanal and butanone) were measured concurrently in the surface water and atmosphere of the South China Sea and Sulu Sea in November 2011. A strong correlation was observed between all OVOC concentrations in the surface seawater along the entire cruise track, except for acetaldehyde, suggesting similar sources and sinks in the surface ocean. Additionally, several phytoplankton groups, such as haptophytes or pelagophytes, were also correlated to all OVOCs indicating that phytoplankton may be an important source for marine OVOCs in the South China and Sulu Seas. Humic and protein like fluorescent dissolved organic matter (FDOM) components seemed to be additional precursors for butanone and acetaldehyde. The atmospheric OVOC mixing ratios were relative high compared with literature values, suggesting the coastal region of North Borneo as a local hot spot for atmospheric OVOCs. The flux of atmospheric OVOCs was largely into the ocean for all 5 gases, with a few important exceptions near the coast of Borneo. The calculated amount of OVOCs entrained into the ocean seemed to be an important source of OVOCs to the surface ocean. When the fluxes were out of the ocean, marine OVOCs were found to be enough to control the local measured OVOC distribution in the atmosphere. Based on our model calculations, at least 0.4 ppb of marine derived acetone and butanone can reach the upper troposphere, where they may have an important influence on hydrogen oxide radical formation over the western Pacific Ocean.

Description: The Extrapolar SWIFT model is a fast ozone chemistry scheme for interactive calculation of the extrapolar stratospheric ozone layer in coupled general circulation models (GCMs). In contrast to the widely used prescribed ozone, the SWIFT ozone layer interacts with the model dynamics and can respond to atmospheric variability or climatological trends. The Extrapolar SWIFT model employs a repro-modelling approach, where algebraic functions are used to approximate the numerical output of a full stratospheric chemistry and transport model (ATLAS). The full model solves a coupled chemical differential equations system with 55 initial and boundary conditions (mixing ratio of various chemical species and atmospheric parameters). Hence the rate of change of ozone over 24  h is a function of 55 variables. Using covariances between these variables, we can find linear combinations in order to reduce the parameter space to the following nine basic variables: latitude, pressure altitude, temperature, local ozone column, mixing ratio of ozone and of the ozone depleting families (Cly, Bry, NOy and HOy). We will show that these 9 variables are sufficient to characterize the rate of change of ozone. An automated procedure fits a polynomial function of fourth degree to the rate of change of ozone obtained from several simulations with the ATLAS model. One polynomial function is determined per month which yields the rate of change of ozone over 24 h. A key aspect for the robustness of the Extrapolar SWIFT model is to include a wide range of stratospheric variability in the numerical output of the ATLAS model, also covering atmospheric states that will occur in a future climate (e.g. temperature and meridional circulation changes or reduction of stratospheric chlorine loading). For validation purposes, the Extrapolar SWIFT model has been integrated into the ATLAS model replacing the full stratospheric chemistry scheme. Simulations with SWIFT in ATLAS have proven that the systematic error is small and does not accumulate during the course of a simulation. In the context of a 10 year simulation, the ozone layer, simulated by SWIFT, shows a stable annual cycle, with inter-annual variations comparable to the ATLAS model. The application of Extrapolar SWIFT requires the evaluation of polynomial functions with 30–100 terms. Nowadays, computers can calculate such polynomial functions at thousands of model grid points in seconds. SWIFT provides the desired numerical efficiency and computes the ozone layer 104 times faster than the chemistry scheme in the ATLAS CTM.

Description: The Southern Hemisphere Antarctic stratosphere experienced two noteworthy events in 2015: a significant injection of sulfur from the Calbuco volcanic eruption in Chile in April, and a record-large Antarctic ozone hole in October and November. Here, we quantify Calbuco's influence on stratospheric ozone depletion in austral spring 2015 using observations and an earth system model. We analyze ozonesondes, as well as data from the Microwave Limb Sounder. We employ the Community Earth System Model, version 1, with the Whole Atmosphere Community Climate Model (CESM1(WACCM)) in a specified dynamics setup, which includes calculations of volcanic effects. The Cloud Aerosol Lidar with Orthogonal Polarization data indicate enhanced volcanic liquid sulfate 532 nm backscatter values as far poleward as 68°S during October and November (in broad agreement with WACCM). Comparison of the location of the enhanced aerosols to ozone data supports the view that aerosols played a major role in increasing the ozone hole size, especially at pressure levels between 150 and 100 hPa. Ozonesonde vertical ozone profiles from the sites of Syowa, South Pole, and Neumayer, display the lowest individual October or November measurements at 150 hPa since the 1991 Mt. Pinatubo eruption period, with Davis showing similarly low values, but no available 1990s data. The analysis suggests that under the cold conditions ideal for ozone depletion, stratospheric volcanic aerosol particles from the moderate-magnitude eruption of Calbuco in 2015 greatly enhanced austral ozone depletion, particularly at 55–68°S, where liquid binary sulfate aerosols have a large influence on ozone concentrations.

Description: Ice retreat in the eastern Eurasian Arctic is a consequence of atmospheric and oceanic processes and regional feedback mechanisms acting on the ice cover, both in winter and summer. A correct representation of these processes in numerical models is important, since it will improve predictions of sea ice anomalies along the Northeast Passage and beyond. In this study, we highlight the importance of winter ice dynamics for local summer sea ice anomalies in thickness, volume and extent. By means of airborne sea ice thickness surveys made over pack ice areas in the south-eastern Laptev Sea, we show that years of offshore-directed sea ice transport have a thinning effect on the late-winter sea ice cover. To confirm the preconditioning effect of enhanced offshore advection in late winter on the summer sea ice cover, we perform a sensitivity study using a numerical model. Results verify that the preconditioning effect plays a bigger role for the regional ice extent. Furthermore, they indicate an increase in volume export from the Laptev Sea as a consequence of enhanced offshore advection, which has far-reaching consequences for the entire Arctic sea ice mass balance. Moreover we show that ice dynamics in winter not only preconditions local summer ice extent, but also accelerate fast-ice decay.

Description: Near-surface air temperatures close to 0°C were observed in situ over sea ice in the central Arctic during the last three winter seasons. Here we use in situ winter (December–March) temperature observations, such as those from Soviet North Pole drifting stations and ocean buoys, to determine how common Arctic winter warming events are. Observations of winter warming events exist over most of the Arctic Basin. Temperatures exceeding -5°C were observed during 〉30% of winters from 1954 to 2010 by North Pole drifting stations or ocean buoys. Using the ERA-Interim record (1979–2016), we show that the North Pole (NP) region typically experiences 10 warming events (T2m 〉 10°C) per winter, compared with only five in the Pacific Central Arctic (PCA). There is a positive trend in the overall duration of winter warming events for both the NP region (4.25 days/decade) and PCA (1.16 days/decade), due to an increased number of events of longer duration.

Description: Retrogressive thaw slumps (RTSs) are among the most active landforms in the Arctic; their number has increased significantly over the past decades. While processes initiating discrete RTSs are well identified, the major terrain controls on the development of coastal RTSs at a regional scale are not yet defined. Our research reveals the main geomorphic factors that determine the development of RTSs along a 238 km segment of the Yukon Coast, Canada. We (1) show the current extent of RTSs, (2) ascertain the factors controlling their activity and initiation, and (3) explain the spatial differences in the density and areal coverage of RTSs. We mapped and classified 287 RTSs using high-resolution satellite images acquired in 2011. We highlighted the main terrain controls over their development using univariate regression trees model. Coastal geomorphology influenced both the activity and initiation of RTSs: active RTSs and RTSs initiated after 1972 occurred primarily on terrains with slope angles greater than 3.9° and 5.9°, respectively. The density and areal coverage of RTSs were constrained by the volume and thickness of massive ice bodies. Differences in rates of coastal change along the coast did not affect the model. We infer that rates of coastal change averaged over a 39 year period are unable to reflect the complex relationship between RTSs and coastline dynamics. We emphasize the need for large-scale studies of RTSs to evaluate their impact on the ecosystem and to measure their contribution to the global carbon budget.

Description: The Mozambique Ridge, a prominent basement high in the southwestern Indian Ocean, consists of four major geomorphological segments associated with numerous phases of volcanic activity in the Lower Cretaceous. The nature and origin of the Mozambique Ridge have been intensely debated with one hypothesis suggesting a Large Igneous Province origin. High-resolution seismic reflection data reveal a large number of extrusion centres with a random distribution throughout the southern Mozambique Ridge and the nearby Transkei Rise. Intra-basement reflections emerge from the extrusion centres and are interpreted to represent massive lava flow sequences. Such lava flow sequences are characteristic of eruptions leading to the formation of continental and oceanic flood basalt provinces, hence supporting a Large Igneous Province origin of the Mozambique Ridge. We observe evidence for widespread post-sedimentary magmatic activity that we correlate with a southward propagation of the East African Rift System. Based on our volumetric analysis of the southern Mozambique Ridge we infer a rapid sequential emplacement between ~131 Ma and ~125 Ma, which is similar to the short formation periods of other Large Igneous Provinces like the Agulhas Plateau.

Description: Digital grids of basement age of the world’s oceans are essential for modern geodynamic and paleoceanographic studies. Any such grid is built using a plate kinematic model, whose accuracy and reli- ability directly influence the accuracy and reliability of the grid. We present a seafloor age grid for the South Atlantic based on a recent high-resolution plate kinematic model. The grid is built from a data set of points whose ages are defined in or for the plate kinematic model, incorporating breaks at tectonic boundaries like fracture zones where the age function is discontinuous. We compare predictions of the new grid and of a previously published one, which is based on an older plate kinematic model, to magnetic isochron pick data sets. The comparison shows the new grid to provide a more reliable depiction of seafloor age in the South Atlantic. Numerical estimates of the new grid’s uncertainty are determined by interpolation between (1) misfits at grid cells coinciding with magnetic isochron ages, (2) misfits implied by locational uncertainties in predicted isochrons propagated from uncertainties in the plate kinematic model, and (3) by the proxim- ities of cells to fracture zone traces or ridge-jump scars. Estimated total uncertainty is 〈10 My for 94% of the grid and 〈5 My for 72%, but much larger in areas where magnetic anomaly data are scarce (such as the Cretaceous Normal Superchron) and in the vicinity of long-offset fracture zones.

Description: The deep southern component water (SCW), comprising Lower Circumpolar Deep Water (LCDW) and Antarctic Bottom Water (AABW), is a major component of the global oceanic circulation. It has been suggested that the deep Atlantic water mass structure changed significantly during the last glacial/interglacial cycle. However, deep SCW source-proximal records remain sparse. Here we present three coherent deep SCW paleocurrent records from the deep Argentine continental margin shedding light on deep water circulation and deep SCW flow strength in the Southwest Atlantic since the Last Glacial Maximum (LGM). Based on increased sortable silt values, we propose enhanced deep SCW flow strength from 14 to 10 cal ka B.P. relative to the early deglacial/LGM and the Holocene. We propose a direct influence of deep northern component water (NCW) on deep SCW flow strength due to vertical narrowing of deep SCW spreading, concurrent with a migration of the high-energetic LCDW/AABW interface occupying our core sites. We suggest a shoaled NCW until 13 cal ka B.P., thereby providing space for deep SCW spreading that resulted in reduced carbonate preservation at our core sites. Increased carbonate content from 13 cal ka B.P. indicates that the NCW expanded changing deep water properties at our core sites in the deep Southwest Atlantic. However, southern sourced terrigenous sediments continued to be deposited at our core sites, suggesting that deep SCW flow was uninterrupted along the Argentine continental margin since the LGM.

Description: Measurements of late springtime nutrient concentrations in Arctic waters are relatively rare due to the extensive sea ice cover that makes sampling difficult. During the SUBICE cruise in May-June 2014, an extensive survey of hydrography and pre-bloom nutrient concentrations was conducted in the Chukchi Sea. Cold (〈 -1.5°C) winter water was prevalent throughout the Chukchi Sea shelf, and the water column was weakly stratified. Nitrate (NO3-) concentration averaged 12.6±1.92 µM in surface waters and 14.0±1.91 µM near the bottom and was significantly correlated with salinity. The highest NO3- concentrations were associated with winter water within the Central Channel flow path. NO3- concentrations were much reduced near the northern shelfbreak within the upper halocline waters of the Canada Basin and along the eastern side of the shelf near the Alaskan coast. Net community production (NCP), estimated as the difference in depth-integrated NO3- content between spring (this study) and summer (historical), varied from 28-38 g C m-2 a-1. This is much lower than previous NCP estimates using NO3- concentrations from the southeastern Bering Sea as a baseline. These results demonstrate the importance of using local profiles of NO3- measured as close to the beginning of the spring bloom as possible when estimating NCP.

Description: Climate trends in the Antarctic region remain poorly characterized, owing to the brevity and scarcity of direct climate observations and the large magnitude of interannual to decadal-scale climate variability. Here, within the framework of the PAGES Antarctica2k working group, we build an enlarged database of ice core water stable isotope records from Antarctica, consisting of 112 records. We produce both unweighted and weighted isotopic (δ18O) composites and temperature reconstructions since 0 CE, binned at 5- and 10-year resolution, for seven climatically distinct regions covering the Antarctic continent. Following earlier work of the Antarctica2k working group, we also produce composites and reconstructions for the broader regions of East Antarctica, West Antarctica and the whole continent. We use three methods for our temperature reconstructions: (i) a temperature scaling based on the δ18O–temperature relationship output from an ECHAM5-wiso model simulation nudged to ERA-Interim atmospheric reanalyses from 1979 to 2013, and adjusted for the West Antarctic Ice Sheet region to borehole temperature data, (ii) a temperature scaling of the isotopic normalized anomalies to the variance of the regional reanalysis temperature and (iii) a composite-plus-scaling approach used in a previous continent-scale reconstruction of Antarctic temperature since 1 CE but applied to the new Antarctic ice core database. Our new reconstructions confirm a significant cooling trend from 0 to 1900 CE across all Antarctic regions where records extend back into the 1st millennium, with the exception of the Wilkes Land coast and Weddell Sea coast regions. Within this long-term cooling trend from 0 to 1900 CE, we find that the warmest period occurs between 300 and 1000 CE, and the coldest interval occurs from 1200 to 1900 CE. Since 1900 CE, significant warming trends are identified for the West Antarctic Ice Sheet, the Dronning Maud Land coast and the Antarctic Peninsula regions, and these trends are robust across the distribution of records that contribute to the unweighted isotopic composites and also significant in the weighted temperature reconstructions. Only for the Antarctic Peninsula is this most recent century-scale trend unusual in the context of natural variability over the last 2000 years. However, projected warming of the Antarctic continent during the 21st century may soon see significant and unusual warming develop across other parts of the Antarctic continent. The extended Antarctica2k ice core isotope database developed by this working group opens up many avenues for developing a deeper understanding of the response of Antarctic climate to natural and anthropogenic climate forcings. The first long-term quantification of regional climate in Antarctica presented herein is a basis for data–model comparison and assessments of past, present and future driving factors of Antarctic climate.

Description: The Eger Rift is an active element of the European Cenozoic Rift System associated with intense Cenozoic intraplate alkaline volcanism and system of sedimentary basins. The intracontinental Cheb Basin at its western part displays geodynamic activity with fluid emanations, persistent seismicity, Cenozoic volcanism, and neotectonic crustal movements at the intersections of major intraplate faults. In this paper, we study detailed geometry of the crust/mantle boundary and its possible origin in the western Eger Rift. We review existing seismic and seismological studies, provide new interpretation of the reflection profile 9HR, and supplement it by new results from local seismicity. We identify significant lateral variations of the high-velocity lower crust and relate them to the distribution and chemical status of mantle-derived fluids and to xenolith studies from corresponding depths. New interpretation based on combined seismic and isotope study points to a local-scale magmatic emplacement at the base of the continental crust within a new rift environment. This concept of magmatic underplating is supported by detecting two types of the lower crust: a high-velocity lower crust with pronounced reflectivity and a high-velocity reflection-free lower crust. The character of the underplated material enables to differentiate timing and tectonic setting of two episodes with different times of origin of underplating events. The lower crust with high reflectivity evidences magmatic underplating west of the Eger Rift of the Late Variscan age. The reflection-free lower crust together with a strong reflector at its top at depths of ~28–30 km forms a magma body indicating magmatic underplating of the late Cenozoic (middle and upper Miocene) to recent. Spatial and temporal relations to recent geodynamic processes suggest active magmatic underplating in the intracontinental setting.

Description: Antarctica is the largest reservoir of ice on Earth. Understanding its ice sheet dynamics is crucial to unraveling past global climate change and making robust climatic and sea level predictions. Of the basic parameters that shape and control ice flow, the most poorly known is geothermal heat flux. Direct observations of heat flux are difficult to obtain in Antarctica, and until now continent-wide heat flux maps have only been derived from low-resolution satellite magnetic and seismological data. We present a high-resolution heat flux map and associated uncertainty derived from spectral analysis of the most advanced continental compilation of airborne magnetic data. Small-scale spatial variability and features consistent with known geology are better reproduced than in previous models, between 36% and 50%. Our high-resolution heat flux map and its uncertainty distribution provide an important new boundary condition to be used in studies on future subglacial hydrology, ice sheet dynamics, and sea level change.

Description: The evidence from both data and models indicates that specific equilibrium climate sensitivity S[X]—the global annual mean surface temperature change (ΔTg) as a response to a change in radiative forcing X (ΔR[X])—is state dependent. Such a state dependency implies that the best fit in the scatterplot of ΔTg versus ΔR[X] is not a linear regression but can be some nonlinear or even nonsmooth function. While for the conventional linear case the slope (gradient) of the regression is correctly interpreted as the specific equilibrium climate sensitivity S[X], the interpretation is not straightforward in the nonlinear case. We here explain how such a state-dependent scatterplot needs to be interpreted and provide a theoretical understanding—or generalization—how to quantify S[X] in the nonlinear case. Finally, from data covering the last 2.1 Myr we show that—due to state dependency—the specific equilibrium climate sensitivity which considers radiative forcing of CO2 and land ice sheet (LI) albedo, math formula, is larger during interglacial states than during glacial conditions by more than a factor 2.

Description: We present vertically resolved observations of aerosol composition during pristine summertime Arctic background conditions. The methansulfonic acid (MSA)-to-sulfate ratio peaked near the surface (mean 0.10), indicating a contribution from ocean-derived biogenic sulfur. Similarly, the organic aerosol (OA)-to-sulfate ratio increased toward the surface (mean 2.0). Both MSA-to-sulfate and OA-to-sulfate ratios were significantly correlated with FLEXPART-WRF-predicted air mass residence time over open water, indicating marine-influenced OA. External mixing of sea salt aerosol from a larger number fraction of organic, sulfate, and amine-containing particles, together with low wind speeds (median 4.7 m s−1), suggests a role for secondary organic aerosol formation. Cloud condensation nuclei concentrations were nearly constant (∼120 cm−3) when the OA fraction was 〈60% and increased to 350 cm−3 when the organic fraction was larger and residence times over open water were longer. Our observations illustrate the importance of marine-influenced OA under Arctic background conditions, which are likely to change as the Arctic transitions to larger areas of open water.

Description: Arctic sea ice has displayed significant thinning as well as an increase in drift speed in recent years. Taken together this suggests an associated rise in sea ice deformation rate. A winter and spring expedition to the sea ice covered region north of Svalbard–the Norwegian young sea ICE2015 expedition (N-ICE2015)—gave an opportunity to deploy extensive buoy arrays and to monitor the deformation of the first-year and secondyear ice now common in the majority of the Arctic Basin. During the 5 month long expedition, the ice cover underwent several strong deformation events, including a powerful storm in early February that damaged the ice cover irreversibly. The values of total deformation measured during N-ICE2015 exceed previously measured values in the Arctic Basin at similar scales: At 100 km scale, N-ICE2015 values averaged above 0.1 d-1, compared to rates of 0.08 d-1 or less for previous buoy arrays. The exponent of the power law between the deformation length scale and total deformation developed over the season from 0.37 to 0.54 with an abrupt increase immediately after the early February storm, indicating a weakened ice cover with more free drift of the sea ice floes. Our results point to a general increase in deformation associated with the younger and thinner Arctic sea ice and to a potentially destructive role of winter storms.

Description: We investigate the impact of different CO2 levels and different subarctic gateway configurations on the surface temperatures during the latest Cretaceous using the Earth System Model COSMOS. The simulated temperatures are compared with the surface temperature reconstructions based on a recent compilation of the latest Cretaceous proxies. In our numerical experiments, the CO2 level ranges from 1 to 6 times the preindustrial (PI) CO2 level of 280 ppm. On a global scale, the most reasonable match between modeling and proxy data is obtained for the experiments with 3 to 5 × PI CO2 concentrations. However, the simulated low- (high-) latitude temperatures are too high (low) as compared to the proxy data. The moderate CO2 levels scenarios might be more realistic, if we take into account proxy data and the dead zone effect criterion. Furthermore, we test if the model-data discrepancies can be caused by too simplistic proxy-data interpretations. This is distinctly seen at high latitudes, where most proxies are biased toward summer temperatures. Additional sensitivity experiments with different ocean gateway configurations and constant CO2 level indicate only minor surface temperatures changes (〈~1°C) on a global scale, with higher values (up to ~8°C) on a regional scale. These findings imply that modeled and reconstructed temperature gradients are to a large degree only qualitatively comparable, providing challenges for the interpretation of proxy data and/or model sensitivity. With respect to the latter, our results suggest that an assessment of greenhouse worlds is best constrained by temperatures in the midlatitudes.

Description: State-of-the-art Arctic Ocean mean sea surface (MSS) models and global geoid models (GGMs) are used to support sea ice freeboard estimation from satellite altimeters, as well as in oceanographic studies such as mapping sea level anomalies and mean dynamic ocean topography. However, errors in a given model in the high frequency domain, primarily due to unresolved gravity features, can result in errors in the estimated along-track freeboard. These errors are exacerbated in areas with a sparse lead distribution in consolidated ice pack conditions. Additionally model errors can impact ocean geostrophic currents, derived from satellite altimeter data, while remaining biases in these models may impact longer-term, multi-sensor oceanographic time-series of sea level change in the Arctic. This study focuses on an assessment of five state-of-the-art Arctic MSS models (UCL13/04, DTU15/13/10) and a commonly used GGM (EGM2008). We describe errors due to unresolved gravity features, inter-satellite biases, and remaining satellite orbit errors, and their impact on the derivation of sea ice freeboard. The latest MSS models, incorporating CryoSat-2 sea surface height measurements, show improved definition of gravity features, such as the Gakkel Ridge. The standard deviation between models ranges 0.03-0.25 m. The impact of remaining MSS/GGM errors on freeboard retrieval can reach several decimeters in parts of the Arctic. While the maximum observed freeboard difference found in the central Arctic was 0.59 m (UCL13 MSS minus EGM2008 GGM), the standard deviation in freeboard differences is 0.03-0.06 m.

Description: Species flocks (SFs) fascinate evolutionary biologists who wonder whether such striking diversification can be driven by normal evolutionary processes. Multiple definitions of SFs have hindered the study of their origins. Previous studies identified a monophyletic taxon as a SF if it displays high speciosity in an area in which it is endemic (criterion 1), high ecological diversity among species (criterion 2), and if it dominates the habitat in terms of biomass (criterion 3); we used these criteria in our analyses. Our starting hypothesis is that normal evolutionary processes may provide a sufficient explanation for most SFs. We thus clearly separate each criterion and identify which biological (intrinsic) and environmental (extrinsic) traits are most favourable to their realization. The first part focuses on evolutionary processes. We highlight that some popular putative causes of SFs, such as key innovations or ecological speciation, are neither necessary nor sufficient to fulfill some or all of the three criteria. Initial differentiation mechanisms are diverse and difficult to identify a posteriori because a primary differentiation of one type (genetic, ecological or geographical) often promotes other types of differentiation. Furthermore, the criteria are not independent: positive feedbacks between speciosity and ecological diversity among species are expected whatever the initial cause of differentiation, and ecological diversity should enhance habitat dominance at the clade level. We then identify intrinsic and extrinsic factors that favour each criterion. Low dispersal emerges as a convincing driver of speciosity. Except for a genomic architecture favouring ecological speciation, for which assessment is difficult, high effective population sizes are the single intrinsic factor that directly enhances speciosity, ecological diversity and habitat dominance. No extrinsic factor appeared to enhance all criteria simultaneously but a combination of factors (insularity, fragmentation and environmental stability) may favour the three criteria, although the effect is indirect for habitat dominance. We then apply this analytical framework to Antarctic marine environments by analysing data from 18 speciose clades belonging to echinoderms (five unrelated clades), notothenioid fishes (five clades) and peracarid crustaceans (eight clades). Antarctic shelf environments and history appear favourable to endemicity and speciosity, but not to ecological specialization. Two main patterns are distinguished among taxa. (i) In echinoderms, many brooding, species‐rich and endemic clades are reported, but without remarkable ecological diversity or habitat dominance. In these taxa, loss of the larval stage is probably a consequence of past Antarctic environmental factors, and brooding is suggested to be responsible for enhanced allopatric speciation (via dispersal limitation). (ii) In notothenioids and peracarids, many clades fulfill all three SF criteria. This could result from unusual features in fish and crustaceans: chromosome instability and key innovations (antifreeze proteins) in notothenioids, ecological opportunity in peracarids, and a genomic architecture favouring ecological speciation in both groups. Therefore, the data do not support our starting point that normal evolutionary factors or processes drive SFs because in these two groups uncommon intrinsic features or ecological opportunity provide the best explanation. The utility of the three‐criterion SF concept is therefore questioned and guidelines are given for future studies.

Description: Our study followed the seasonal cycling of soluble (SFe), colloidal (CFe), dissolved (DFe), total dissolvable (TDFe), labile particulate (LPFe) and total particulate (TPFe) iron in the Celtic Sea (NE Atlantic Ocean). Preferential uptake of SFe occurred during the spring bloom, preceding the removal of CFe. Uptake and export of Fe during the spring bloom, coupled with a reduction in vertical exchange, led to Fe deplete surface waters (〈0.2 nM DFe; 0.11 nM LPFe, 0.45 nM TDFe, 1.84 nM TPFe) during summer stratification. Below the seasonal thermocline, DFe concentrations increased from spring to autumn, mirroring NO3- and consistent with supply from remineralised sinking organic material, and cycled independently of particulate Fe over seasonal timescales. These results demonstrate that summer Fe availability is comparable to the seasonally Fe limited Ross Sea shelf, and therefore is likely low enough to affect phytoplankton growth and species composition.

Description: Northwestern Namibia, at the landfall of the Walvis Ridge, was affected by the Tristan da Cunha mantle plume during continental rupture between Africa and South America, as evidenced by the presence of the Etendeka continental flood basalts. Here we use data from a passive-source seismological network to investigate the upper mantle structure and to elucidate the Cretaceous mantle plume-lithosphere interaction. Receiver functions reveal an interface associated with a negative velocity contrast within the lithosphere at an average depth of 80 km. We interpret this interface as the relic of the lithosphereasthenosphere boundary (LAB) formed during the Mesozoic by interaction of the Tristan da Cunha plume head with the pre-existing lithosphere. The velocity contrast might be explained by stagnated and ‘‘frozen’’ melts beneath an intensively depleted and dehydrated peridotitic mantle. The present-day LAB is poorly visible with converted waves, indicating a gradual impedance contrast. Beneath much of the study area, converted phases of the 410 and 660 km mantle transition zone discontinuities arrive 1.5 s earlier than in the landward plume-unaffected continental interior, suggesting high velocities in the upper mantle caused by a thick lithosphere. This indicates that after lithospheric thinning during continental breakup, the lithosphere has increased in thickness during the last 132 Myr. Thermal cooling of the continental lithosphere alone cannot produce the lithospheric thickness required here. We propose that the remnant plume material, which has a higher seismic velocity than the ambient mantle due to melt depletion and dehydration, significantly contributed to the thickening of the mantle lithosphere.

Description: Snow on sea ice alters the properties of the underlying ice cover as well as associated physical and biological processes at the interfaces between atmosphere, sea ice, and ocean. The Antarctic snow cover persists during most of the year and contributes significantly to the sea-ice mass due to the widespread surface flooding and related snow-ice formation. Snow also enhances the sea-ice surface reflectivity of incoming shortwave radiation and determines therefore the amount of light being reflected, absorbed, and transmitted to the upper ocean. Here, we present results of a case study of spectral solar radiation measurements under Antarctic pack ice with an instrumented Remotely Operated Vehicle in the Weddell Sea in 2013. In order to identify the key variables controlling the spatial distribution of the under-ice light regime, we exploit under-ice optical measurements in combination with simultaneous characterization of surface properties, such as sea-ice thickness and snow depth. Our results reveal that the distribution of flooded and nonflooded sea-ice areas dominates the spatial scales of under-ice light variability for areas smaller than 100 m-by-100 m. However, the heterogeneous and highly metamorphous snow on Antarctic pack ice obscures a direct correlation between the under-ice light field and snow depth. Compared to the Arctic, light levels under Antarctic pack ice are extremely low during spring (〈0.1%). This is mostly a result of the distinctly different dominant sea ice and snow properties with seasonal snow cover (including strong surface melt and summer melt ponds) in the Arctic and a year-round snow cover and widespread surface flooding in the Southern Ocean.

Description: Snow thickness on sea ice is a largely under-sampled parameter, yet of importance for the sea-ice mass balance and for satellite based sea-ice thickness estimates and thus our general understanding of global ice-volume change. Traditional direct thickness measurements with meter sticks can provide accurate but only spot information, referred to as "needles" due to their pinpoint focus and information while airborne and satellite remote sensing snow products, referred to as "the haystack" have large uncertainties due to their scale. We demonstrate the remarkable accuracy and applicability of ground penetrating radar (GPR) snow-thickness measurements by comparing them with in-situ, meter-stick data from two field campaigns to Antarctica in late winter/early spring. The efficiency and millimeter- to centimeter accuracy of GPR enables practitioners to acquire extensive, semi-regional data with the potential to upscale "needles" to "the haystack" and to potentially calibrate satellite remote sensing products that we confirm to derive roughly 30 % of the in-situ thickness. We find the radar wave propagation velocity in snow to be rather constant (+/- 6%), encouraging regional snow-thickness surveys. Snow thinner than 10 cm is under the detection limit with the off-the-shelf GPR setup utilized in our study.

Description: New two-year long records from three moorings, located at 76°S along the eastern flank and shelf of the Filchner Trough, give insight in the seasonal cycle of hydrography to a region where Modified Warm Deep Water (MWDW) enters the southern Weddell Sea continental shelf, possibly reaching the Filchner Ronne Ice Shelf, the biggest ice shelf (by volume) in Antarctica. A persistent northward flow of Ice Shelf Water (ISW) is found along the slope of the trough at 600 m depth, while the data on the shelf indicate a seasonal cycle, characterized by four phases. A distinct warm inflow period (separated in two phases), with maximum temperatures of −1°C, appears to be related to the seasonal heaving of the Antarctic Slope Front thermocline along the continental shelf break further north and a seasonal extension of the ISW layer onto the Eastern Shelf. The density gradients between the ISW in the trough and the MWDW on the adjacent shelf drive the southward flow during these phases. A flow reversal is found in winter, ceasing the southward flow along the eastern flank of the trough. Weaker density gradients between the trough and the shelf during winter allow a westward flow, partly driven by a N-S density gradient, existing across the Eastern Shelf during this time. From spring through to summer the ISW layer in the trough extends onto the eastern shelf where it occupies the bottom layer at our moorings and it is associated with northward flow.

Description: The Norwegian young sea ICE (N-ICE2015) expedition was designed to investigate the atmosphere-snow-ice-ocean interactions in the young and thin sea ice regime north of Svalbard. Radiosondes were launched twice daily during the expedition from January to June 2015. Here we use these upper air measurements to study the multiple cyclonic events observed during N-ICE2015 with respect to changes in the vertical thermodynamic structure, moisture content, and boundary layer characteristics. We provide statistics of temperature inversion characteristics, static stability, and boundary layer extent. During winter, when radiative cooling is most effective, we find the strongest impact of synoptic cyclones. Changes to thermodynamic characteristics of the boundary layer are associated with transitions between the radiatively “clear” and “opaque” atmospheric states. In spring, radiative fluxes warm the surface leading to lifted temperature inversions and a statically unstable boundary layer. Further, we compare the N-ICE2015 static stability distributions to corresponding profiles from ERA-Interim reanalysis, from the closest land station in the Arctic North Atlantic sector, Ny-Ålesund, and to soundings from the SHEBA expedition (1997/1998). We find similar stability characteristics for N-ICE2015 and SHEBA throughout the troposphere, despite differences in location, sea ice thickness, and snow cover. For Ny-Ålesund, we observe similar characteristics above 1000 m, while the topography and ice-free fjord surrounding Ny-Ålesund generate great differences below. The long-term radiosonde record (1993–2014) from Ny-Ålesund indicates that during the N-ICE2015 spring period, temperatures were close to the climatological mean, while the lowest 3000 m were 1–3∘C warmer than the climatology during winter.

Description: Focused research on the Pine Island and Thwaites glaciers, which drain the West Antarctic Ice Shelf (WAIS) into the Amundsen Sea Embayment (ASE), revealed strong signs of instability in recent decades that result from variety of reasons, such as inflow of warmer ocean currents and reverse bedrock topogra- phy, and has been established as the Marine Ice Sheet Instability hypothesis. Geothermal heat flux (GHF) is a poorly constrained parameter in Antarctica and suspected to affect basal conditions of ice sheets, i.e., basal melting and subglacial hydrology. Thermomechanical models demonstrate the influential boundary condition of geothermal heat flux for (paleo) ice sheet stability. Due to a complex tectonic and magmatic history of West Antarctica, the region is suspected to exhibit strong heterogeneous geothermal heat flux variations. We present an approach to investigate ranges of realistic heat fluxes in the ASE by different methods, discuss direct observations, and 3-D numerical models that incorporate boundary conditions derived from various geophysical studies, including our new Depth to the Bottom of the Magnetic Source (DBMS) estimates. Our in situ temperature measurements at 26 sites in the ASE more than triples the number of direct GHF observations in West Antarctica. We demonstrate by our numerical 3-D models that GHF spatially varies from 68 up to 110 mW m-2.

Description: The Barents Sea is one of the fastest changing regions of the Arctic, and has experienced the strongest decline in winter-time sea-ice area in the Arctic, at -23 +/- 4% per decade. Sea-ice thickness in the Barents Sea is not well studied. We present two previously unpublished helicopter-borne electromagnetic (HEM) ice thickness measurements from the northwestern Barents Sea acquired in March 2003 and 2014. The HEM data are compared to ice thickness calculated from ice draft measured by ULS deployed between 1994 and 1996. These data show that ice thickness varies greatly from year to year; influenced by the thermodynamic and dynamic processes that govern local formation vs long-range advection. In a year with a large inflow of sea-ice from the Arctic Basin, the Barents Sea ice cover is dominated by thick multiyear ice; as was the case in 2003 and 1995. In a year with an ice cover that was mainly grown in situ, the ice will be thin and mechanically unstable; as was the case in 2014. The HEM data allow us to explore the spatial and temporal variability in ice thickness. In 2003 the dominant ice class was more than 2 years old; and modal sea-ice thickness varied regionally from 0.6 to 1.4 m, with the thinner ice being either first-year ice, or multiyear ice which had come into contact with warm Atlantic water. In 2014 the ice cover was predominantly locally grown ice less than 1 month old (regional modes of 0.5–0.8 m). These two situations represent two extremes of a range of possible ice thickness distributions that can present very different conditions for shipping traffic; or have a different impact on heat transport from ocean to atmosphere.

Description: Antarctic pack ice serves as habitat for microalgae which contribute to Southern Ocean primary production and serve as important food source for pelagic herbivores. Ice algal biomass is highly patchy and remains severely undersampled by classical methods such as spatially restricted ice coring surveys. Here we provide an unprecedented view of ice algal biomass distribution, mapped (as chlorophyll a) in a 100 m by 100 m area of a Weddell Sea pack ice floe, using under-ice irradiance measurements taken with an instrumented remotely operated vehicle. We identified significant correlations (p 〈 0.001) between algal biomass and concomitant in situ surface measurements of snow depth, ice thickness, and estimated sea ice freeboard levels using a statistical model. The model’s explanatory power (r2 = 0.30) indicates that these parameters alone may provide a first basis for spatial prediction of ice algal biomass, but parameterization of additional determinants is needed to inform more robust upscaling efforts.

Description: The global ocean state for the modern age and for the Last Glacial Maximum (LGM) was dynamically reconstructed with a sophisticated data assimilation technique. A substantial amount of data including global seawater temperature, salinity (only for the modern estimate), and the isotopic composition of oxygen and carbon (only in the Atlantic for the LGM) were integrated into an ocean general circulation model with the help of the adjoint method, thereby the model was optimized to reconstruct plausible continuous fields of tracers, overturning circulation and water mass distribution. The adjoint-based LGM state estimation of this study represents the state of the art in terms of the length of forward model runs, the number of observations assimilated, and the model domain. Compared to the modern state, the reconstructed continuous sea-surface temperature field for the LGM shows a global-mean cooling of 2.2 K, and the reconstructed LGM ocean has a more vigorous Atlantic meridional overturning circulation, shallower North Atlantic Deep Water (NADW) equivalent, stronger stratification, and more saline deep water.

Description: Nioghalvfjerdsfjorden Glacier (79NG) is the largest of three marine-terminating outlet glaciers draining the Northeast Greenland Ice Stream. To understand how Atlantic waters supply waters in the cavity beneath the floating 79NG, we analyze historic and recent bathymetric, hydrographic, and velocity observations obtained on the Northeast Greenland continental shelf. The bathymetry is characterized by a trough system, consisting of the Westwind Trough and the Norske Trough in the northern and southern part of the continental shelf, respectively. Atlantic waters recirculating in Fram Strait cross the shelf break and enter the trough system at its southeastern inlet toward the inner shelf. Warm Atlantic Intermediate Water (AIW) present below 200 m in the Norske Trough shows large contributions of the recirculating Atlantic water. We found that the bathymetry is sufficiently deep to provide a direct subsurface pathway for warm AIW between the shelf break and the 79NG cavity via the Norske Trough. Likewise, based on the hydrographic data, we show that the Norske Trough supplies AIW warmer than 1°C to the 79NG, which is not present in the Westwind Trough. Our moored and lowered velocity measurements indicate that a boundary current carries warm AIW along the northeastern slope of Norske Trough toward the 79NG. We suggest that anomalies in Atlantic water temperatures in Fram Strait could reach 79NG within less than 1.5 years, thereby modifying the glacier's basal melt rates.

Description: The Central Yakutian permafrost landscape is rapidly being modified by land use and global warming, but small-scale thermokarst process variability and hydrological conditions are poorly understood. We analyze lake-area changes and thaw subsidence of young thermokarst lakes on ice-complex deposits (yedoma lakes) in comparison to residual lakes in alas basins during the last 70 years for a local study site and we record regional lake size and distribution on different ice-rich permafrost terraces using satellite and historical airborne imagery. Statistical analysis of climatic and ground-temperature data identified driving factors of yedoma- and alas-lake changes. Overall, lake area is larger today than in 1944 but alas-lake levels have oscillated greatly over 70 years, with a mean alas-lake-radius change rate of 1.663.0 m/yr. Anthropogenic disturbance and forest degradation initiated, and climate forced rapid, continuous yedoma-lake growth. The mean yedoma lake-radius change rate equals 1.261.0 m/yr over the whole observation period. Mean thaw subsidence below yedoma lakes is 6.261.4 cm/yr. Multiple regression analysis suggests that winter precipitation, winter temperature, and active-layer properties are primary controllers of area changes in both lake types; summer weather and permafrost conditions additionally influence yedoma-lake growth rates. The main controlling factors of alas-lake changes are unclear due to larger catchment areas and subsurface hydrological conditions. Increasing thermokarst activity is currently linked to older terraces with higher ground-ice contents, but thermokarst activity will likely stay high and wet conditions will persist within the near future in Central Yakutian alas basins.

Description: © The Author(s), 2017. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Molecular Biology and Evolution 34 (2017): 1890-1901, doi:10.1093/molbev/msx125.

Description: The highly conserved ADAR enzymes, found in all multicellular metazoans, catalyze the editing of mRNA transcripts by the deamination of adenosines to inosines. This type of editing has two general outcomes: site specific editing, which frequently leads to recoding, and clustered editing, which is usually found in transcribed genomic repeats. Here, for the first time, we looked for both editing of isolated sites and clustered, non-specific sites in a basal metazoan, the coral Acropora millepora during spawning event, in order to reveal its editing pattern. We found that the coral editome resembles the mammalian one: it contains more than 500,000 sites, virtually all of which are clustered in non-coding regions that are enriched for predicted dsRNA structures. RNA editing levels were increased during spawning and increased further still in newly released gametes. This may suggest that editing plays a role in introducing variability in coral gametes.

Description: This work was supported by the Australian Research Council (to PK), the European Research Council (grant 311257), the I-CORE Program of the Planning and Budgeting Committee in Israel (grants 41/11 and 1796/12), and the Israel Science Foundation (1380/14).

Description: Author Posting. © The Authors, 2016. This article is posted here by permission of Oxford University Press for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 208 (2017): 1026-1042, doi:10.1093/gji/ggw435.

Description: In recent years, marine controlled source electromagnetics (CSEM) has found increasing use in hydrocarbon exploration due to its ability to detect thin resistive zones beneath the seafloor. It is the purpose of this paper to evaluate the physics of CSEM for an ocean whose electrical thickness is comparable to or much thinner than that of the overburden using the in-line configuration through examination of the elliptically polarized seafloor electric field, the time-averaged energy flow depicted by the real part of the complex Poynting vector, energy dissipation through Joule heating and the Fréchet derivatives of the seafloor field with respect to the subseafloor conductivity that is assumed to be isotropic. The deep water (ocean layer electrically much thicker than the overburden) seafloor EM response for a model containing a resistive reservoir layer has a greater amplitude and reduced phase as a function of offset compared to that for a half-space, or a stronger and faster response. For an ocean whose electrical thickness is comparable to or much smaller than that of the overburden, the electric field displays a greater amplitude and reduced phase at small offsets, shifting to a stronger amplitude and increased phase at intermediate offsets and a weaker amplitude and enhanced phase at long offsets, or a stronger and faster response that first changes to stronger and slower, and then transitions to weaker and slower. These transitions can be understood by visualizing the energy flow throughout the structure caused by the competing influences of the dipole source and guided energy flow in the reservoir layer, and the air interaction caused by coupling of the entire subseafloor resistivity structure with the sea surface. A stronger and faster response occurs when guided energy flow is dominant, while a weaker and slower response occurs when the air interaction is dominant. However, at intermediate offsets for some models, the air interaction can partially or fully reverse the direction of energy flux in the reservoir layer toward rather than away from the source, resulting in a stronger and slower response. The Fréchet derivatives are dominated by preferential sensitivity to the reservoir layer conductivity for all water depths except at high frequencies, but also display a shift with offset from the galvanic to the inductive mode in the underburden and overburden due to the interplay of guided energy flow and the air interaction. This means that the sensitivity to the horizontal conductivity is almost as strong as to the vertical component in the shallow parts of the subsurface, and in fact is stronger than the vertical sensitivity deeper down. However, the sensitivity to horizontal conductivity is still weak compared to the vertical component within thin resistive regions. The horizontal sensitivity is gradually decreased when the water becomes deep. These observations in part explain the success of shallow towed CSEM using only measurements of the in-line component of the electric field.

Description: © The Author(s), 2017. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in PLoS One 12 (2017): e0173350, doi:10.1371/journal.pone.0173350.

Description: Scleractinian coral are experiencing unprecedented rates of mortality due to increases in sea surface temperatures in response to global climate change. Some coral species however, survive high temperature events due to a reduced susceptibility to bleaching. We investigated the relationship between bleaching susceptibility and expression of five metabolically related genes of Symbiodinium spp. from the coral Porites astreoides originating from an inshore and offshore reef in the Florida Keys. The acclimatization potential of Symbiodinium spp. to changing temperature regimes was also measured via a two-year reciprocal transplant between the sites. Offshore coral fragments displayed significantly higher expression in Symbiodinium spp. genes PCNA, SCP2, G3PDH, PCP and psaE than their inshore counterparts (p〈0.05), a pattern consistent with increased bleaching susceptibility in offshore corals. Additionally, gene expression patterns in Symbiodinium spp. from site of origin were conserved throughout the two-year reciprocal transplant, indicating acclimatization did not occur within this multi-season time frame. Further, laboratory experiments were used to investigate the influence of acute high temperature (32°C for eight hours) and disease (lipopolysaccharide of Serratia marcescens) on the five metabolically related symbiont genes from the same offshore and inshore P. astreoides fragments. Gene expression did not differ between reef fragments, or as a consequence of acute exposure to heat or heat and disease, contrasting to results found in the field. Gene expression reported here indicates functional variation in populations of Symbiodinium spp. associated with P. astreoides in the Florida Keys, and is likely a result of localized adaptation. However, gene expression patterns observed in the lab imply that functional variation in zooxanthellae observed under conditions of chronic moderate stress is lost under the acute extreme conditions studied here.

Description: Funding for this research was provided by Coastal Preservation network (KBS, BHS).

Description: © The Author(s), 2017. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in PLoS One 12 (2017): e0179318, doi: 10.1371/journal.pone.0179318.

Description: Trace metals are essential for health but toxic when present in excess. The maintenance of trace metals at physiologic levels reflects both import and export by cells and absorption and excretion by organs. The mechanism by which this maintenance is achieved in vertebrate organisms is incompletely understood. To explore this, we chose zebrafish as our model organism, as they are amenable to both pharmacologic and genetic manipulation and comprise an ideal system for genetic screens and toxicological studies. To characterize trace metal content in developing zebrafish, we measured levels of three trace elements, copper, zinc, and manganese, from the oocyte stage to 30 days post-fertilization using inductively coupled plasma mass spectrometry. Our results indicate that metal levels are stable until zebrafish can acquire metals from the environment and imply that the early embryo relies on maternal contribution of metals to the oocyte. We also measured metal levels in bodies and yolks of embryos reared in presence and absence of the copper chelator neocuproine. All three metals exhibited different relative abundances between yolks and bodies of embryos. While neocuproine treatment led to an expected phenotype of copper deficiency, total copper levels were unaffected, indicating that measurement of total metal levels does not equate with measurement of biologically active metal levels. Overall, our data not only can be used in the design and execution of genetic, physiologic, and toxicologic studies but also has implications for the understanding of vertebrate metal homeostasis.

Description: This work was supported by the National Institutes of Health, R00 DK84122.

Description: © The Author(s), 2017. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Bioscience 67 (2017): 760–768, doi:10.1093/biosci/bix059.

Description: As the sampling frequency and resolution of Earth observation imagery increase, there are growing opportunities for novel applications in population monitoring. New methods are required to apply established analytical approaches to data collected from new observation platforms (e.g., satellites and unmanned aerial vehicles). Here, we present a method that estimates regional seasonal abundances for an understudied and growing population of gray seals (Halichoerus grypus) in southeastern Massachusetts, using opportunistic observations in Google Earth imagery. Abundance estimates are derived from digital aerial survey counts by adapting established correction-based analyses with telemetry behavioral observation to quantify survey biases. The result is a first regional understanding of gray seal abundance in the northeast US through opportunistic Earth observation imagery and repurposed animal telemetry data. As species observation data from Earth observation imagery become more ubiquitous, such methods provide a robust, adaptable, and cost-effective solution to monitoring animal colonies and understanding species abundances.

Description: We would like to thank generous support from International Fund for Animal Welfare, the Bureau of Ocean Energy, and the Oak Foundation for funding support for the telemetry devices.

Description: Author Posting. © The Authors, 2017. This article is posted here by permission of Oxford University Press for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 211 (2017): 1046–1061, doi:10.1093/gji/ggx360.

Description: In recent years, marine controlled source electromagnetics (CSEM) has found increasing use in hydrocarbon exploration due to its ability to detect thin resistive zones beneath the seafloor. It is the purpose of this paper to evaluate the physics of CSEM for an ocean whose electrical thickness is comparable to or much thinner than that of the overburden using the in-line configuration through examination of the elliptically-polarized seafloor electric field, the time-averaged energy flow depicted by the real part of the complex Poynting vector, energy dissipation through Joule heating and the Fréchet derivatives of the seafloor field with respect to the sub-seafloor conductivity that is assumed to be transversely anisotropic, with a vertical-to-horizontal resistivity ratio of 3:1. For an ocean whose electrical thickness is comparable to that of the overburden, the seafloor electromagnetic response for a model containing a resistive reservoir layer has a greater amplitude and reduced phase as a function of offset compared to that for a halfspace, or a stronger and faster response, and displays little to no evidence for the air interaction. For an ocean whose electrical thickness is much smaller than that of the overburden, the electric field displays a greater amplitude and reduced phase at small offsets, shifting to a stronger amplitude and increased phase at intermediate offsets, and a weaker amplitude and enhanced phase at long offsets, or a stronger and faster response that first changes to stronger and slower, and then transitions to weaker and slower. By comparison to the isotropic case with the same horizontal conductivity, transverse anisotropy stretches the Poynting vector and the electric field response from a thin resistive layer to much longer offsets. These phenomena can be understood by visualizing the energy flow throughout the structure caused by the competing influences of the dipole source and guided energy flow in the reservoir layer, and the air interaction caused by coupling of the entire sub-seafloor resistivity structure with the sea surface. The Fréchet derivatives are dominated by preferential sensitivity to the vertical conductivity in the reservoir layer and overburden at short offsets. The horizontal conductivity Fréchet derivatives are weaker than to comparable to the vertical derivatives at long offsets in the substrate. This means that the sensitivity to the horizontal conductivity is present in the shallow parts of the subsurface. In the presence of transverse anisotropy, it is necessary to go to higher frequencies to sense the horizontal conductivity in the overburden as compared to an isotropic model with the same horizontal conductivity. These observations in part explain the success of shallow towed CSEM using only measurements of the in-line component of the electric field.

Description: This work was supported at WHOI by an Independent Research and Development award, and by the Walter A. and Hope Noyes Smith Chair for Excellence in Oceanography.

Description: A major uncertainty in determining the mass balance of the Antarctic ice sheet from measurements of satellite gravimetry, and to a lesser extent satellite altimetry, is the poorly known correction for the ongoing deformation of the solid Earth caused by glacial isostatic adjustment (GIA). Although much progress has been made in consistently modeling the ice-sheet evolution throughout the last glacial cycle, as well as the induced bedrock deformation caused by these load changes, forward models of GIA remain ambiguous due to the lack of observational constraints on the ice sheet's past extent and thickness and mantle rheology beneath the continent. As an alternative to forward-modeling GIA, we estimate GIA from multiple space-geodetic observations: Gravity Recovery and Climate Experiment (GRACE), Envisat/ICESat and Global Positioning System (GPS). Making use of the different sensitivities of the respective satellite observations to current and past surface-mass (ice mass) change and solid Earth processes, we estimate GIA based on viscoelastic response functions to disc load forcing. We calculate and distribute the viscoelastic response functions according to estimates of the variability of lithosphere thickness and mantle viscosity in Antarctica. We compare our GIA estimate with published GIA corrections and evaluate its impact in determining the ice-mass balance in Antarctica from GRACE and satellite altimetry. Particular focus is applied to the Amundsen Sea Sector in West Antarctica, where uplift rates of several centimetres per year have been measured by GPS. We show that most of this uplift is caused by the rapid viscoelastic response to recent ice-load changes, enabled by the presence of a low-viscosity upper mantle in West Antarctica. This paper presents the second and final contributions summarizing the work carried out within a European Space Agency funded study, REGINA (www.regina-science.eu).

Description: © The Author(s), 2017. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Genome Biology and Evolution 9 (2017): 659-676, doi:10.1093/gbe/evx023.

Description: Understanding and predicting the fate of populations in changing environments require knowledge about the mechanisms that support phenotypic plasticity and the adaptive value and evolutionary fate of genetic variation within populations. Atlantic killifish (Fundulus heteroclitus) exhibit extensive phenotypic plasticity that supports large population sizes in highly fluctuating estuarine environments. Populations have also evolved diverse local adaptations. To yield insights into the genomic variation that supports their adaptability, we sequenced a reference genome and 48 additional whole genomes from a wild population. Evolution of genes associated with cell cycle regulation and apoptosis is accelerated along the killifish lineage, which is likely tied to adaptations for life in highly variable estuarine environments. Genome-wide standing genetic variation, including nucleotide diversity and copy number variation, is extremely high. The highest diversity genes are those associated with immune function and olfaction, whereas genes under greatest evolutionary constraint are those associated with neurological, developmental, and cytoskeletal functions. Reduced genetic variation is detected for tight junction proteins, which in killifish regulate paracellular permeability that supports their extreme physiological flexibility. Low-diversity genes engage in more regulatory interactions than high-diversity genes, consistent with the influence of pleiotropic constraint on molecular evolution. High genetic variation is crucial for continued persistence of species given the pace of contemporary environmental change. Killifish populations harbor among the highest levels of nucleotide diversity yet reported for a vertebrate species, and thus may serve as a useful model system for studying evolutionary potential in variable and changing environments.

Description: This work was primarily supported by a grant from the National Science Foundation (collaborative research grants DEB-1265282, DEB-1120512, DEB-1120013, DEB-1120263, DEB-1120333, DEB-1120398 to J.K.C., D.L.C., M.E.H., S.I.K., M.F.O., J.R.S., W.W., and A.W.). Further support was provided by the National Institute of Environmental Health Sciences (1R01ES021934-01 to A.W., P42ES7373 to T.H.H., P42ES007381 to M.E.H., and R01ES019324 to J.R.S.), the National Institute of General Medical Sciences (P20GM103423 and P20GM104318 to B.L.K.), and the National Science Foundation (DBI-0640462 and XSEDE-MCB100147 to D.G.).

Description: © The Author(s), 2017. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in PLoS ONE 12 (2017): e0184849, doi:10.1371/journal.pone.0184849.

Description: Diatoms are important components of marine ecosystems and contribute greatly to the world's primary production. Despite their important roles in ecosystems, the molecular basis of how diatoms cope with oxidative stress caused by nutrient fluctuations remains largely unknown. Here, an isobaric tags for relative and absolute quantitation (iTRAQ) proteomic method was coupled with a series of physiological and biochemical techniques to explore oxidative stress- and cell fate decision-related cellular and metabolic responses of the diatom Thalassiosira pseudonana to nitrate (N) and inorganic phosphate (P) stresses. A total of 1151 proteins were detected; 122 and 56 were significantly differentially expressed from control under N- and P-limited conditions, respectively. In N-limited cells, responsive proteins were related to reactive oxygen species (ROS) accumulation, oxidative stress responses and cell death, corresponding to a significant decrease in photosynthetic efficiency, marked intracellular ROS accumulation, and caspase-mediated programmed cell death activation. None of these responses were identified in P-limited cells; however, a significant up-regulation of alkaline phosphatase proteins was observed, which could be the major contributor for P-limited cells to cope with ambient P deficiency. These findings demonstrate that fundamentally different metabolic responses and cellular regulations are employed by the diatom in response to different nutrient stresses and to keep the cells viable.

Description: This study was funded by the National Natural Science Foundation of China (41576138, 41076080, 41576138) to Dr. Jun-Rong Liang; the Woods Hole Center for Oceans and Human Health, National Science Foundation (OCE-1314642) to Dr. DonaldM Anderson; the National Institute of Environmental Health Sciences (1-P01-ES021923- 01) to Dr. DonaldM Anderson; and the ERC Advanced Award Diatomite and ANR project DiaDomOil to Dr. Chris Bowler.

Description: © The Author(s), 2016. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in PLoS ONE 11 (2016): e0155049, doi:10.1371/journal.pone.0155049.

Description: Many reefs have shifted from coral and fish dominated habitats to less productive macroalgal dominated habitats, and current research is investigating means of reversing this phase shift. In the tropical Pacific, overfished reefs with inadequate herbivory can become dominated by the brown alga Sargassum polycystum. This alga suppresses recruitment and survival of corals and fishes, thus limiting the potential for reef recovery. Here we investigate the mechanisms that reinforce S. polycystum dominance and show that in addition to negatively affecting other species, this species acts in a self-reinforcing manner, positively promoting survival and growth of conspecifics. We found that survival and growth of both recruit-sized and mature S. polycystum fronds were higher within Sargassum beds than outside the beds and these results were found in both protected and fished reefs. Much of this benefit resulted from reduced herbivory within the Sargassum beds, but adult fronds also grew ~50% more within the beds even when herbivory did not appear to be occurring, suggesting some physiological advantage despite the intraspecific crowding. Thus via positive feedbacks, S. polycystum enhances its own growth and resistance to herbivores, facilitating its dominance (perhaps also expansion) and thus its resilience on degraded reefs. This may be a key feedback mechanism suppressing the recovery of coral communities in reefs dominated by macroalgal beds.

Description: Financial support came to MEH from the National Science Foundation (OCE 0929119), the National Institutes of Health (2 U19 TW007401-10), and the Teasley Endowment to the Georgia Institute of Technology.

Description: © The Author(s), 2017. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Conservation Physiology 5 (2017): cox061, doi:10.1093/conphys/cox061.

Description: Recent studies have demonstrated that some hormones are present in baleen powder from bowhead (Balaena mysticetus) and North Atlantic right (Eubalaena glacialis) whales. To test the potential generalizability of this technique for studies of stress and reproduction in large whales, we sought to determine whether all major classes of steroid and thyroid hormones are detectable in baleen, and whether these hormones are detectable in other mysticetes. Powdered baleen samples were recovered from single specimens of North Atlantic right, bowhead, blue (Balaenoptera [B.]musculus), sei (B. borealis), minke (B. acutorostrata), fin (B. physalus), humpback (Megaptera novaeangliae) and gray (Eschrichtius robustus) whales. Hormones were extracted with a methanol vortex method, after which we tested all species with commercial enzyme immunoassays (EIAs, Arbor Assays) for progesterone, testosterone, 17β-estradiol, cortisol, corticosterone, aldosterone, thyroxine and tri-iodothyronine, representing a wide array of steroid and thyroid hormones of interest for whale physiology research. In total, 64 parallelism tests (8 species × 8 hormones) were evaluated to verify good binding affinity of the assay antibodies to hormones in baleen. We also tested assay accuracy, although available sample volume limited this test to progesterone, testosterone and cortisol. All tested hormones were detectable in baleen powder of all species, and all assays passed parallelism and accuracy tests. Although only single individuals were tested, the consistent detectability of all hormones in all species indicates that baleen hormone analysis is likely applicable to a broad range of mysticetes, and that the EIA kits tested here perform well with baleen extract. Quantification of hormones in baleen may be a suitable technique with which to explore questions that have historically been difficult to address in large whales, including pregnancy and inter-calving interval, age of sexual maturation, timing and duration of seasonal reproductive cycles, adrenal physiology and metabolic rate.

Description: This work was supported by (1) the Center for Bioengineering Innovation at Northern Arizona University and (2) the New England Aquarium.

Description: © The Author(s), 2017. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in PLoS One 12 (2017): e0188257, doi:10.1371/journal.pone.0188257.

Description: Preservation of three-dimensional structure in the gut is necessary in order to analyze the spatial organization of the gut microbiota and gut luminal contents. In this study, we evaluated preparation methods for mouse gut with the goal of preserving micron-scale spatial structure while performing fluorescence imaging assays. Our evaluation of embedding methods showed that commonly used media such as Tissue-Tek Optimal Cutting Temperature (OCT) compound, paraffin, and polyester waxes resulted in redistribution of luminal contents. By contrast, a hydrophilic methacrylate resin, Technovit H8100, preserved three-dimensional organization. Our mouse intestinal preparation protocol optimized using the Technovit H8100 embedding method was compatible with microbial fluorescence in situ hybridization (FISH) and other labeling techniques, including immunostaining and staining with both wheat germ agglutinin (WGA) and 4', 6-diamidino-2-phenylindole (DAPI). Mucus could be visualized whether the sample was fixed with paraformaldehyde (PFA) or with Carnoy’s fixative. The protocol optimized in this study enabled simultaneous visualization of micron-scale spatial patterns formed by microbial cells in the mouse intestines along with biogeographical landmarks such as host-derived mucus and food particles.

Description: Funding provided by National Science Foundation (www.nsf.gov) grant 1650141 to J.L.M.W. and National Institutes of Health National Institute of Dental and Craniofacial Research (www.nidcr.nih.gov) grant DE022586 to G.G.B.

Description: This is an open access article, free of all copyright. The definitive version was published in PLoS ONE 12 (2017): e0172839, doi:10.1371/journal.pone.0172839.

Description: Marine animals, such as turtles, seabirds and pelagic fishes, are observed to travel and congregate around eddies in the open ocean. Mesoscale eddies, large swirling ocean vortices with radius scales of approximately 50–100 km, provide environmental variability that can structure these populations. In this study, we investigate the use of mesoscale eddies by 24 individual juvenile loggerhead sea turtles (Caretta caretta) in the Brazil-Malvinas Confluence region. The influence of eddies on turtles is assessed by collocating the turtle trajectories to the tracks of mesoscale eddies identified in maps of sea level anomaly. Juvenile loggerhead sea turtles are significantly more likely to be located in the interiors of anticyclones in this region. The distribution of surface drifters in eddy interiors reveals no significant association with the interiors of cyclones or anticyclones, suggesting higher prevalence of turtles in anticyclones is a result of their behavior. In the southern portion of the Brazil-Malvinas Confluence region, turtle swimming speed is significantly slower in the interiors of anticyclones, when compared to the periphery, suggesting that these turtles are possibly feeding on prey items associated with anomalously low near-surface chlorophyll concentrations observed in those features.

Description: A Woods Hole Postdoctoral Fellowship awarded to PG and NASA grant NNX13AE47G funded this work. CB is grateful for support from NOAA and the NSF-GRFP (Grant No. 1314109). DJM gratefully acknowledges support from NASA and NSF.

Description: This is an open access article, free of all copyright. The definitive version was published in PLoS ONE 12 (2017): e0170962, doi:10.1371/journal.pone.0170962.

Description: Bio-logging tags are an important tool for the study of cetaceans, but superficial tags inevitably increase hydrodynamic loading. Substantial forces can be generated by tags on fast-swimming animals, potentially affecting behavior and energetics or promoting early tag removal. Streamlined forms have been used to reduce loading, but these designs can accelerate flow over the top of the tag. This non-axisymmetric flow results in large lift forces (normal to the animal) that become the dominant force component at high speeds. In order to reduce lift and minimize total hydrodynamic loading this work presents a new tag design (Model A) that incorporates a hydrodynamic body, a channel to reduce fluid speed differences above and below the housing and wing to redirect flow to counter lift. Additionally, three derivatives of the Model A design were used to examine the contribution of individual flow control features to overall performance. Hydrodynamic loadings of four models were compared using computational fluid dynamics (CFD). The Model A design eliminated all lift force and generated up to ~30 N of downward force in simulated 6 m/s aligned flow. The simulations were validated using particle image velocimetry (PIV) to experimentally characterize the flow around the tag design. The results of these experiments confirm the trends predicted by the simulations and demonstrate the potential benefit of flow control elements for the reduction of tag induced forces on the animal.

Description: This project was funded by the National Oceanographic Partnership Program [National Science Foundation via the Office of Naval Research N00014-11-1-0113]. C. Spencer Garborg was supported by a Grove City College Swezey Student Fellowship to Erik Anderson. Mark Johnson was funded by a Marie Curie-Sklodowska grant from the European Union.

Description: © The Author(s), 2017. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Annual Review of Marine Science 9 (2017): 173-203, doi:10.1146/annurev-marine-010816-060733.

Description: The events that followed the Tohoku earthquake and tsunami on March 11, 2011, included the loss of power and overheating at the Fukushima Daiichi nuclear power plants, which led to extensive releases of radioactive gases, volatiles, and liquids, particularly to the coastal ocean. The fate of these radionuclides depends in large part on their oceanic geochemistry, physical processes, and biological uptake. Whereas radioactivity on land can be resampled and its distribution mapped, releases to the marine environment are harder to characterize owing to variability in ocean currents and the general challenges of sampling at sea. Five years later, it is appropriate to review what happened in terms of the sources, transport, and fate of these radionuclides in the ocean. In addition to the oceanic behavior of these contaminants, this review considers the potential health effects and societal impacts.

Description: K.B. was supported in part by the Gordon and Betty Moore Foundation and the Deerbrook Charitable Trust. P.M. was supported in part by the Generalitat de Catalunya through MERS (grant 2014 SGR 1356), the European Commission 7th Framework COMET-FRAME project (grant agreement 604974), and the Ministerio de Economía y Competitividad of Spain (project CTM2011-15152-E). S.C. was supported in part by the French program Investissement d'Avenir run by the National Research Agency (AMORAD project, grant ANR-11-RSNR-0002). D.O. was supported in part by the Center for Environmental Radioactivity (NFR Centers of Excellence grant 223268/F50). J.N.S. was supported in part by the Marine Environmental Observation, Prediction, and Response Network.

Description: © The Author(s), 2017. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in PLoS ONE 12 (2017): e0186156, doi:10.1371/journal.pone.0186156.

Description: Bowhead whales (Balaena mysticetus) have a nearly circumpolar distribution, and occasionally occupy warmer shallow coastal areas during summertime that may facilitate molting. However, relatively little is known about the occurrence of molting and associated behaviors in bowhead whales. We opportunistically observed whales in Cumberland Sound, Nunavut, Canada with skin irregularities consistent with molting during August 2014, and collected a skin sample from a biopsied whale that revealed loose epidermis and sloughing. During August 2016, we flew a small unmanned aerial system (sUAS) over whales to take video and still images to: 1) determine unique individuals; 2) estimate the proportion of the body of unique individuals that exhibited sloughing skin; 3) determine the presence or absence of superficial lines representative of rock-rubbing behavior; and 4) measure body lengths to infer age-class. The still images revealed that all individuals (n = 81 whales) were sloughing skin, and that nearly 40% of them had mottled skin over more than two-thirds of their bodies. The video images captured bowhead whales rubbing on large rocks in shallow, coastal areas—likely to facilitate molting. Molting and rock rubbing appears to be pervasive during late summer for whales in the eastern Canadian Arctic.

Description: Fieldwork was funded by: Fisheries and Oceans Canada (Emerging Fisheries 41436-810-120-4D875), Nunavut Wildlife Research Trust Fund (project 3-13-29 Bowhead Whale Movement's and Ecology), Molson Foundation (318366), Ocean Tracking Network (NSERC NETGP 375118-08), and ArcticNet Centre of Excellence (317588) awarded to S.H. Ferguson and World Wildlife Fund Canada (Arctic Species Conservation Fund) awarded to W.R. Koski and S.M.E. Fortune. Additional field work support was provided by: U.S. Department of the Interior, Minerals Management Service (MMS; now Bureau of Ocean Energy Management), through Inter-agency Agreement No. M08PG20021 with the U.S. Department of Commerce, National Oceanic and Atmospheric Administration, as part of the MMS Alaska Environmental Studies Program awarded to M.F. Baumgartner. Personnel support was provided by: Natural Sciences and Engineering Research Council Canadian Graduate Scholarship, the W. Garfield Weston Award for Northern Research, University of British Columbia Affiliated Fellowship and Northern Scientific Training Program (Canadian Polar Commission) awarded to S.M.E. Fortune