Polar Region Bathymetry: Critical Knowledge for the Prediction of Global Sea Level Rise

The ocean and the marine parts of the cryosphere interact directly with, and are affected by, the seafloor and its primary properties of depth (bathymetry) and shape (morphology) in many ways. Bottom currents are largely constrained by undersea terrain with consequences for both regional and global heat transport. Deep ocean mixing is controlled by seafloor roughness, and the bathymetry directly influences where marine outlet glaciers are susceptible to the inflow relatively warm subsurface waters - an issue of great importance for ice-sheet discharge, i.e., the loss of mass from calving and undersea melting. Mass loss from glaciers and the Greenland and Antarctic ice sheets, is among the primary drivers of global sea-level rise, together now contributing more to sea-level rise than the thermal expansion of the ocean. Recent research suggests that the upper bounds of predicted sea-level rise by the year 2100 under the scenarios presented in IPCC’s Special Report on the Ocean and Cryosphere in a Changing Climate (SROCCC) likely are conservative because of the many unknowns regarding ice dynamics. In this paper we highlight the poorly mapped seafloor in the Polar regions as a critical knowledge gap that needs to be filled to move marine cryosphere science forward and produce improved understanding of the factors impacting ice-discharge and, with that, improved predictions of, among other things, global sea-level. We analyze the bathymetric data coverage in the Arctic Ocean specifically and use the results to discuss challenges that must be overcome to map the most remotely located areas in the Polar regions in general.

Temporal Distribution of Arsenic and Metals in Soil From King George Island, Antarctica

The polar regions are vulnerable to impacts caused by local and global pollution. The Antarctic continent has been considered an environment that has remained little affected by human activities. Direct exposure to contaminants may occur in areas continuously occupied by research stations for several decades. Admiralty Bay on the southeast coast of King George Island, has potential for being affected by human activities due research stations operating in the area, including the Brazilian Commandant Ferraz Antarctic Station (CFAS). The levels of metals and arsenic were determined in soils collected near CFAS (points 5, 6, 7, and 9), Base G and at two points distant from the CFAS: Refuge II and Hennequin. Samples were collected after the fire in CFAS occurred in February 2012, up to December 2018 to assess the environmental impacts in the area. Al and As were related with Base G. Refuge II and Hennequin can be considered as control points for this region. As a consequence of the accident, the increased levels for Cd, Cu, Pb, and Zn, especially at point 9 (inside the CFAS) and in the soil surrounding the CFAS in 2013. The results from 2016 to 2018 demonstrated a reduction in levels of all studied metals near CFAS, which may be related to the leaching of metals into Admiralty Bay; it is thus, being important the continue monitoring soil, sediments, and Antarctic biota.

Temperature, carbon dioxide and methane

1) Globally-representative monthly rates of change of atmospheric carbon dioxide and methane are compared with global rates of change of sea ice and with Arctic and Antarctic air temperatures. 2) Carbon dioxide is very strongly correlated with sea ice dynamics, with the carbon dioxide rate at Mauna Loa lagging sea ice extent rate by 7 months. 3) Methane is very strongly correlated with sea ice dynamics, with the global (and Mauna Loa) methane rate lagging sea ice extent rate by 5 months. 4) Sea ice melt rate peaks in very tight synchrony with temperature in each Hemisphere. 5) The very high synchrony of the two gases is most parsimoniously explained by a common causality acting in both Hemispheres. 6) Time lags between variables indicate primary drivers of the gas dynamics are due to solar action on the polar regions, not mid-latitudes as is conventionally believed. 7) Results are consistent with a proposed role of a high-latitude temperature-dependent abiotic variable such as sea ice in the annual cycles of carbon dioxide and methane. 8) If sea ice does not drive the net flux of these gases, it is a highly precise proxy for whatever does. 9) Potential mechanisms should be investigated urgently.

Cartesian Meshing Spherical Earth (CMSE): A Code Package to Incorporate the Spherical Earth in SPECFEM3D Cartesian Simulations

The SPECFEM3D_Cartesian code package is widely used in simulating seismic wave propagation on local and regional scales due to its computational efficiency compared with the one-chunk version of the SPECFEM3D_Globe code. In SPECFEM3D_Cartesian, the built-in meshing tool maps a spherically curved cube to a rectangular cube using the Universal Transverse Mercator projection (UTM). Meanwhile, the geodetic east, north, and up directions are assigned as the local x–y–z directions. This causes coordinate orientation issues in simulating waveform propagation in regions larger than 6° × 6° or near the Earth’s polar regions. In this study, we introduce a new code package, named Cartesian Meshing Spherical Earth (CMSE), that can accurately mesh the 3D geometry of the Earth’s surface under the Cartesian coordinate frame, while retaining the geodetic directions. To benchmark our new package, we calculate the residual amplitude of the CMSE synthetics with respect to the reference synthetics calculated by SPECFEM3D_Globe. In the regional scale simulations with an area of 1300 km × 1300 km, we find a maximum of 5% amplitude residual for the SPECFEM3D_Cartesian synthetics using the mesh generated by the CMSE, much smaller than the maximum amplitude residual of 100% for the synthetics based on its built-in meshing tool. Therefore, our new meshing tool CMSE overcomes the limitations of the internal mesher used by SPECFEM3D_Cartesian and can be used for more accurate waveform simulations in larger regions beyond one UTM zone. Furthermore, CMSE can deal with regions at the south and north poles that cannot be handled by the UTM projection. Although other external code packages can be used to mesh the curvature of the Earth, the advantage of the CMSE code is that it is open-source, easy to use, and fully integrated with SPECFEM3D_Cartesian.

Science goals and new mission concepts for future exploration of Titan’s atmosphere, geology and habitability: titan POlar scout/orbitEr and in situ lake lander and DrONe explorer (POSEIDON)

In response to ESA’s “Voyage 2050” announcement of opportunity, we propose an ambitious L-class mission to explore one of the most exciting bodies in the Solar System, Saturn’s largest moon Titan. Titan, a “world with two oceans”, is an organic-rich body with interior-surface-atmosphere interactions that are comparable in complexity to the Earth. Titan is also one of the few places in the Solar System with habitability potential. Titan’s remarkable nature was only partly revealed by the Cassini-Huygens mission and still holds mysteries requiring a complete exploration using a variety of vehicles and instruments. The proposed mission concept POSEIDON (Titan POlar Scout/orbitEr and In situ lake lander DrONe explorer) would perform joint orbital and in situ investigations of Titan. It is designed to build on and exceed the scope and scientific/technological accomplishments of Cassini-Huygens, exploring Titan in ways that were not previously possible, in particular through full close-up and in situ coverage over long periods of time. In the proposed mission architecture, POSEIDON consists of two major elements: a spacecraft with a large set of instruments that would orbit Titan, preferably in a low-eccentricity polar orbit, and a suite of in situ investigation components, i.e. a lake lander, a “heavy” drone (possibly amphibious) and/or a fleet of mini-drones, dedicated to the exploration of the polar regions. The ideal arrival time at Titan would be slightly before the next northern Spring equinox (2039), as equinoxes are the most active periods to monitor still largely unknown atmospheric and surface seasonal changes. The exploration of Titan’s northern latitudes with an orbiter and in situ element(s) would be highly complementary in terms of timing (with possible mission timing overlap), locations, and science goals with the upcoming NASA New Frontiers Dragonfly mission that will provide in situ exploration of Titan’s equatorial regions, in the mid-2030s.

Physical processes behind interactions of microplastic particles with natural ice

Microplastic particles (MPs, <5 mm) are found in marine ice in larger quantities than in seawater, however, the distribution pattern within the ice cores is not consistent. To get insights into the most general physical processes behind interactions of ice and plastic particles in cool natural environments, information from academic and applied research is integrated and verified against available field observations. Non-polar molecules of common-market plastics are hydrophobic, so MPs are weak ice nucleators, are repelled from water and ice, and concentrate within air bubbles and brine channels. A large difference in thermal properties of ice and plastics favours concentration of MPs at the ice surface during freeze/thaw cycles. Under low environmental temperatures, falling in polar regions below the glass / brittle-ductile transition temperatures of the common-use plastics, they become brittle. This might partially explain the absence of floating macroplastics in polar waters. Freshwater freezes at the temperature well below that of its maximum density, so the water column is stably stratified, and MPs eventually concentrate at the ice surface and in air bubbles. In contrast, below growing sea ice, mechanisms of suspension freezing under conditions of (thermal plus haline) convection should permanently entangle MPs into ice. During further sea ice growth and aging, MPs are repelled from water and ice into air bubbles, brine channels, and to the upper/lower boundaries of the ice column. Sea ice permeability, especially while melting periods, can re-distribute sub-millimeter MPs through the brine channels, thus potentially introducing the variability of contamination with time. In accord with field observations, analysis reveals several competing factors that influence the distribution of MPs in sea ice. A thorough sampling of the upper ice surface, prevention of brine leakage while sampling and handling, considering the ice structure while segmenting the ice core – these steps may be advantageous for further understanding the pattern of plastic contamination in natural ice.

Effect of Space Fractional Parameter on Nonlinear Ion Acoustic Shock Wave Excitation in an Unmagnetized Relativistic Plasma

In this work, the model equation with space fractional-order (FO) is used to investigate the nonlinear ion acoustic shock wave excitations (NIASWEs) in an unmagnetized collisionless weakly relativistic plasma having inertial relativistic ions fluid with viscous effects, inertial-less non-thermal electrons and inertial-less Boltzmann positrons. To do it, the Korteweg-de Vries Burgers equation (KdVBE) is derived from the considered fluid model equations by implementing the standard reductive perturbation method. Accordingly, such equation is converted to space fractional KdVBE via Agrawal’s variational principle with the help of the beta fractional derivative and its properties. The exact analytical solutions of KdVBE with space FO are determined via the modified Kudryashov method. The influence of space fractional and other related plasma parameters on NIASWEs are investigated. The outcomes would be useful to understand the nature of shocks with the presence of non-local or local space in many astrophysical and space environments (especially in the relativistic wind of pulsar magnetosphere, polar regions of neutron stars, etc.) and further laboratory verification.

OClO as observed by TROPOMI: a comparison with meteorological parameters and polar stratospheric cloud observations

Chlorine dioxide (OClO) is a by-product of the ozone-depleting halogen chemistry in the stratosphere. Although it is rapidly photolysed at low solar zenith angles (SZAs), it plays an important role as an indicator of the chlorine activation in polar regions during polar winter and spring at twilight conditions because of the nearly linear dependence of its formation on chlorine oxide (ClO). Here, we compare slant column densities (SCDs) of chlorine dioxide (OClO) retrieved by means of differential optical absorption spectroscopy (DOAS) from spectra measured by the TROPOspheric Monitoring Instrument (TROPOMI) with meteorological data for both Antarctic and Arctic regions for the first three winters in each of the hemispheres (November 2017–October 2020). TROPOMI, a UV–Vis–NIR–SWIR instrument on board of the Sentinel-5P satellite, monitors the Earth's atmosphere in a near-polar orbit at an unprecedented spatial resolution and signal-to-noise ratio and provides daily global coverage at the Equator and thus even more frequent observations at polar regions. The observed OClO SCDs are generally well correlated with the meteorological conditions in the polar winter stratosphere; for example, the chlorine activation signal appears as a sharp gradient in the time series of the OClO SCDs once the temperature drops to values well below the nitric acid trihydrate (NAT) existence temperature (TNAT). Also a relation of enhanced OClO values at lee sides of mountains can be observed at the beginning of the winters, indicating a possible effect of lee waves on chlorine activation. The dataset is also compared with CALIPSO Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) polar stratospheric cloud (PSC) observations. In general, OClO SCDs coincide well with CALIOP measurements for which PSCs are detected. Very high OClO levels are observed for the northern hemispheric winter 2019/20, with an extraordinarily long period with a stable polar vortex being even close to the values found for southern hemispheric winters. An extraordinary winter in the Southern Hemisphere was also observed in 2019, with a minor sudden stratospheric warming at the beginning of September. In this winter, similar OClO values were measured in comparison to the previous (usual) winter till that event but with a OClO deactivation that was 1–2 weeks earlier.

Stream Boundary Detection of a Hyper-Arid, Polar Region Using a U-Net Architecture: Taylor Valley, Antarctica

Convolutional neural networks (CNNs) are becoming an increasingly popular approach for classification mapping of large complex regions where manual data collection is too time consuming. Stream boundaries in hyper-arid polar regions such as the McMurdo Dry Valleys (MDVs) in Antarctica are difficult to locate because they have little hydraulic flow throughout the short summer months. This paper utilizes a U-Net CNN to map stream boundaries from lidar derived rasters in Taylor Valley located within the MDVs, covering ∼770 km2. The training dataset consists of 217 (300 × 300 m2) well-distributed tiles of manually classified stream boundaries with diverse geometries (straight, sinuous, meandering, and braided) throughout the valley. The U-Net CNN is trained on elevation, slope, lidar intensity returns, and flow accumulation rasters. These features were used for detection of stream boundaries by providing potential topographic cues such as inflection points at stream boundaries and reflective properties of streams such as linear patterns of wetted soil, water, or ice. Various combinations of these features were analyzed based on performance. The test set performance revealed that elevation and slope had the highest performance of the feature combinations. The test set performance analysis revealed that the CNN model trained with elevation independently received a precision, recall, and F1 score of 0.94±0.05, 0.95±0.04, and 0.94±0.04 respectively, while slope received 0.96±0.03, 0.93±0.04, and 0.94±0.04, respectively. The performance of the test set revealed higher stream boundary prediction accuracies along the coast, while inland performance varied. Meandering streams had the highest stream boundary prediction performance on the test set compared to the other stream geometries tested here because meandering streams are further evolved and have more distinguishable breaks in slope, indicating stream boundaries. These methods provide a novel approach for mapping stream boundaries semi-automatically in complex regions such as hyper-arid environments over larger scales than is possible for current methods.