ALBERT

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: 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: 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: 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 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 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.