Genetic analysis of immunological traits in tilapia

The immunological response to handling stress of four tilapia species is evaluated.Polymorphism is examined in genes known to influence immune response in fish

Future Antarctic bed topography and its implications for ice sheet dynamics

The Antarctic bedrock is evolving as the solid Earth responds to the past and ongoing evolution of the ice sheet. A recently improved ice loading history suggests that the Antarctic Ice Sheet (AIS) has generally been losing its mass since the Last Glacial Maximum. In a sustained warming climate, the AIS is predicted to retreat at a greater pace, primarily via melting beneath the ice shelves. We employ the glacial isostatic adjustment (GIA) capability of the Ice Sheet System Model (ISSM) to combine these past and future ice loadings and provide the new solid Earth computations for the AIS. We find that past loading is relatively less important than future loading for the evolution of the future bed topography. Our computations predict that the West Antarctic Ice Sheet (WAIS) may uplift by a few meters and a few tens of meters at years AD 2100 and 2500, respectively, and that the East Antarctic Ice Sheet is likely to remain unchanged or subside minimally except around the Amery Ice Shelf. The Amundsen Sea Sector in particular is predicted to rise at the greatest rate; one hundred years of ice evolution in this region, for example, predicts that the coastline of Pine Island Bay will approach roughly 45 mm yr−1 in viscoelastic vertical motion. Of particular importance, we systematically demonstrate that the effect of a pervasive and large GIA uplift in the WAIS is generally associated with the flattening of reverse bed slope, reduction of local sea depth, and thus the extension of grounding line (GL) towards the continental shelf. Using the 3-D higher-order ice flow capability of ISSM, such a migration of GL is shown to inhibit the ice flow. This negative feedback between the ice sheet and the solid Earth may promote stability in marine portions of the ice sheet in the future

Representation of sharp rifts and faults mechanics in modeling ice shelf flow dynamics: Application to Brunt/Stancomb-Wills Ice Shelf, Antarctica

Ice shelves play a major role in buttressing ice sheet flow into the ocean, hence the importance of accurate numerical modeling of their stress regime. Commonly used ice flow models assume a continuous medium and are therefore complicated by the presence of rupture features (crevasses, rifts, and faults) that significantly affect the overall flow patterns. Here we apply contact mechanics and penalty methods to develop a new ice shelf flow model that captures the impact of rifts and faults on the rheology and stress distribution of ice shelves. The model achieves a best fit solution to satellite observations of ice shelf velocities to infer the following: (1) a spatial distribution of contact and friction points along detected faults and rifts, (2) a more realistic spatial pattern of ice shelf rheology, and (3) a better representation of the stress balance in the immediate vicinity of faults and rifts. Thus, applying the model to the Brunt/Stancomb-Wills Ice Shelf, Antarctica, we quantify the state of friction inside faults and the opening rates of rifts and obtain an ice shelf rheology that remains relatively constant everywhere else on the ice shelf. We further demonstrate that better stress representation has widespread application in examining aspects affecting ice shelf structure and dynamics including the extent of ice mélange in rifts and the change in fracture configurations. All are major applications for better insight into the important question of ice shelf stability

ISSM: Ice Sheet System Model

In order to have the capability to use satellite data from its own missions to inform future sea-level rise projections, JPL needed a full-fledged ice-sheet/iceshelf flow model, capable of modeling the mass balance of Antarctica and Greenland into the near future. ISSM was developed with such a goal in mind, as a massively parallelized, multi-purpose finite-element framework dedicated to ice-sheet modeling. ISSM features unstructured meshes (Tria in 2D, and Penta in 3D) along with corresponding finite elements for both types of meshes. Each finite element can carry out diagnostic, prognostic, transient, thermal 3D, surface, and bed slope simulations. Anisotropic meshing enables adaptation of meshes to a certain metric, and the 2D Shelfy-Stream, 3D Blatter/Pattyn, and 3D Full-Stokes formulations capture the bulk of the ice-flow physics. These elements can be coupled together, based on the Arlequin method, so that on a large scale model such as Antarctica, each type of finite element is used in the most efficient manner. For each finite element referenced above, ISSM implements an adjoint. This adjoint can be used to carry out model inversions of unknown model parameters, typically ice rheology and basal drag at the ice/bedrock interface, using a metric such as the observed InSAR surface velocity. This data assimilation capability is crucial to allow spinning up of ice flow models using available satellite data. ISSM relies on the PETSc library for its vectors, matrices, and solvers. This allows ISSM to run efficiently on any parallel platform, whether shared or distrib- ISSM: Ice Sheet System Model NASA's Jet Propulsion Laboratory, Pasadena, California uted. It can run on the largest clusters, and is fully scalable. This allows ISSM to tackle models the size of continents. ISSM is embedded into MATLAB and Python, both open scientific platforms. This improves its outreach within the science community. It is entirely written in C/C++, which gives it flexibility in its design, and the power/speed that C/C++ allows. ISSM is svn (subversion) hosted, on a JPL repository, to facilitate its development and maintenance. ISSM can also model propagation of rifts using contact mechanics and mesh splitting, and can interface to the Dakota software. To carry out sensitivity analysis, mesh partitioning algorithms are available, based on the Scotch, Chaco, and Metis partitioners that ensure equal area mesh partitions can be done, which are then usable for sampling and local reliability methods

DNA BARCODING OF FISH SPECIES FROM THE MEDITERRANEAN COAST OF ISRAEL

Accurately-classified genomic data in the Barcode of Life Data System (BOLD) database is vital to the protection and conservation of marine biodiversity in the Mediterranean Sea. The taxonomic classifications of 468 fish of 50 Mediterranean species were analyzed using the BOLD Identifier tool for variation in the cytochrome oxidase subunit I (COI) mitochondrial gene. Within species, nucleotide maximum composite likelihood was low with a mean of 0.0044±0.0008. Three presumptive species had significantly higher values e.g., Arnoglossus spp. (0.07), Torquigener flavimaculosus (0.013) and Boops boops (0.028). However, samples of Arnoglossus species were sub-classified into two groups that were finally identified as two different species e.g., Arnoglossus laterna and Arnoglossus thori. For the different species, BLAST searches against the BOLD database using our DNA barcoding data as the query sequences designated the most similar targets into groups. For each analyzed species, the similarity of the first and second threshold groups ranged from 95 to 99% and from 83 to 98%, respectively. Sequence based classification for the first threshold group was concordant with morphology-based identification. However, for 34 analyzed species (68%) overlaps of species between the two threshold groups hampered classification. Tree-based phylogeny analysis detected more than one cluster in the first threshold group for 22 out of 50 species, representing genetic subgroups and geographic origins. There was a tendency for higher conservation and lower number of clusters in the Lessepsian (Red Sea) migrant versus indigenous species

High-resolution ice-thickness mapping in South Greenland

Airborne radar sounding is difficult in South Greenland because of the presence of englacial water, which prevents the signal from reaching the bed. Data coverage remains suboptimal for traditional methods of ice-thickness and bed mapping that rely on geostatistical techniques, such as kriging, because important features are missing. Here we apply two alternative approaches of highresolution (̃300 m) ice-thickness mapping, that are based on the conservation of mass, to two regions of South Greenland: (1) Qooqqup Sermia and Kiattuut Sermiat, and (2) Ikertivaq. These two algorithms solve optimization problems, for which the conservation of mass is either enforced as a hard constraint, or as a soft constraint. For the first region, very few measurements are available but there is no gap in ice motion data, whereas for Ikertivaq, more ice-thickness measurements are available, but there are gaps in ice motion data. We show that mass-conservation algorithms can be used as validation tools for radar sounding. We also show that it is preferable to apply mass conservation as a hard constraint, rather than a soft constraint, as it better preserves elongated features, such as glacial valleys and ridges

Complex Basal Thermal Transition Near the Onset of Petermann Glacier, Greenland

The basal thermal regime of ice sheets exerts a strong control on ice‐sheet stability and the onset of rapidly streaming flow. However, the nature of this thermal transition where sliding initiates is largely unconstrained by geophysical observations. In the Greenland Ice Sheet, topographic troughs or elevated geothermal heat fluxes typically define the onset of outlet glaciers. In contrast, Petermann Glacier in Northern Greenland does not have any distinct bed troughs or localized geothermal heating associated with its onset, making it an ideal site to investigate the basal thermal state and examine its role in the onset of Petermann Glacier. Here we use radar bed reflectivity and an ice‐sheet thermomechanical model to examine the basal thermal regime beneath Petermann Glacier. Our results reveal a complex thermal transition near the onset of Petermann Glacier. As the bed shifts from largely frozen to largely thawed with increasing distances from the ice divide, our results show that this thermal transition happens through alternating bands of frozen and thawed bands. The complex thermal state across the onset region suggests that lateral meltwater injection and local meltwater production determine the location of Petermann Glacier. Given the lack of topographic pinning at the onset location, the upstream margin of Petermann is vulnerable to migrate depending on a combination of advective cooling and meltwater supply from the interior of the ice sheet

Ice Sheet Model Intercomparison Project (ISMIP6) contribution to CMIP6

Reducing the uncertainty in the past, present, and future contribution of ice sheets to sea-level change requires a coordinated effort between the climate and glaciology communities. The Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6) is the primary activity within the Coupled Model Intercomparison Project – phase 6 (CMIP6) focusing on the Greenland and Antarctic ice sheets. In this paper, we describe the framework for ISMIP6 and its relationship with other activities within CMIP6. The ISMIP6 experimental design relies on CMIP6 climate models and includes, for the first time within CMIP, coupled ice-sheet–climate models as well as standalone ice-sheet models. To facilitate analysis of the multi-model ensemble and to generate a set of standard climate inputs for standalone ice-sheet models, ISMIP6 defines a protocol for all variables related to ice sheets. ISMIP6 will provide a basis for investigating the feedbacks, impacts, and sea-level changes associated with dynamic ice sheets and for quantifying the uncertainty in ice-sheet-sourced global sea-level change

Retreat of Thwaites Glacier, West Antarctica, over the next 100 years using various ice flow models, ice shelf melt scenarios and basal friction laws

Thwaites Glacier (TG), West Antarctica, has experienced rapid, potentially irreversible grounding line retreat and mass loss in response to enhanced ice shelf melting. Results from recent numerical models suggest a large spread in the evolution of the glacier in the coming decades to a century. It is therefore important to investigate how different approximations of the ice stress balance, parameterizations of basal friction and ice shelf melt parameterizations may affect projections. Here, we simulate the evolution of TG using ice sheet models of varying levels of complexity, different basal friction laws and ice shelf melt to quantify their effect on the projections. We find that the grounding line retreat and its sensitivity to ice shelf melt are enhanced when a full-Stokes model is used, a Budd friction is used and ice shelf melt is applied on partially floating elements. Initial conditions also impact the model results. Yet, all simulations suggest a rapid, sustained retreat of the glacier along the same preferred pathway. The fastest retreat rate occurs on the eastern side of the glacier, and the slowest retreat occurs across a subglacial ridge on the western side. All the simulations indicate that TG will undergo an accelerated retreat once the glacier retreats past the western subglacial ridge. Combining all the simulations, we find that the uncertainty of the projections is small in the first 30 years, with a cumulative contribution to sea level rise of 5&thinsp;mm, similar to the current rate. After 30 years, the contribution to sea level depends on the model configurations, with differences up to 300&thinsp;% over the next 100 years, ranging from 14 to 42&thinsp;mm.</p

Recruiting Diverse Smokers: Enrollment Yields and Cost

To help tobacco control research better include vulnerable populations, we sought to identify effective ways to recruit diverse smokers. In 2014–2015, we recruited 2149 adult cigarette smokers in California and North Carolina, United States, to participate in a randomized trial of pictorial cigarette pack warnings. The most effective means of recruiting smokers were the classified advertising website Craigslist (28% of participants), word of mouth (23%), Facebook (16%), and flyers or postcards (14%). Low-income and African American smokers were more likely to respond to interpersonal contact (including staff in-person recruitment and word of mouth) than were high-income and non-African American smokers (all p < 0.05). Hispanic and gay, lesbian, and bisexual smokers were more likely to be recruited by Craigslist than non-Hispanic and straight smokers (both p < 0.05). Of the recruitment methods requiring cost, the cheapest was Craigslist ($3–7 per smoker). The most expensive methods were newspaper ads in California ($375 per smoker) and staff in-person recruiting in North Carolina ($180 per smoker). Successfully recruiting diverse smokers requires using multiple methods including interpersonal, online, and other media. Craigslist and word of mouth are especially useful and low-cost ways to recruit diverse smokers