We provide research findings on the physics of aerosol and droplet dispersion relevant to the hypothesized aerosol transmission of SARS-CoV-2 during the current pandemic. We utilize physics-based modeling at different levels of complexity, along with previous literature on coronaviruses, to investigate the possibility of airborne transmission. The previous literature, our 0D-3D simulations by various physics-based models, and theoretical calculations, indicate that the typical size range of speech and cough originated droplets ( ) allows lingering in the air for ) so that they could be inhaled. Consistent with the previous literature, numerical evidence on the rapid drying process of even large droplets, up to sizes , into droplet nuclei/aerosols is provided. Based on the literature and the public media sources, we provide evidence that the individuals, who have been tested positive on COVID-19, could have been exposed to aerosols/droplet nuclei by inhaling them in significant numbers e.g. . By 3D scale-resolving computational fluid dynamics (CFD) simulations, we give various examples on the transport and dilution of aerosols ( ) over distances in generic environments. We study susceptible and infected individuals in generic public places by Monte-Carlo modelling. The developed model takes into account the locally varying aerosol concentration levels which the susceptible accumulate via inhalation. The introduced concept, ’exposure time’ to virus containing aerosols is proposed to complement the traditional ’safety distance’ thinking. We show that the exposure time to inhale aerosols could range from to or even to depending on the situation. The Monte-Carlo simulations, along with the theory, provide clear quantitative insight to the exposure time in different public indoor environments.

Abstract. The Fourth IPCC Assessment Report concluded that ice sheet flow models, in their current state, were unable to provide accurate forecast for the increase of polar ice sheet discharge and the associated contribution to sea level rise. Since then, the glaciological community has undertaken a huge effort to develop and improve a new generation of ice flow models, and as a result a significant number of new ice sheet models have emerged. Among them is the parallel finite-element model Elmer/Ice, based on the opensource multi-physics code Elmer. It was one of the first fullStokes models used to make projections for the evolution of the whole Greenland ice sheet for the coming two centuries. Originally developed to solve local ice flow problems of high mechanical and physical complexity, Elmer/Ice has today reached the maturity to solve larger-scale problems, earning the status of an ice sheet model. Here, we summarise almost 10 yr of development performed by different groups. Elmer/Ice solves the full-Stokes equations, for isotropic but also anisotropic ice rheology, resolves the grounding line dynamics as a contact problem, and contains various basal friction laws. Derived fields, like the age of the ice, the strain rate or stress, can also be computed. Elmer/Ice includes two recently proposed inverse methods to infer badly known parameters. Elmer is a highly parallelised code thanks to recent developments and the implementation of a block preconditioned solver for the Stokes system. In this paper, all these components are presented in detail, as well as the numerical performance of the Stokes solver and developments planned for the future.

[1] Sliding velocity and basal drag are strongly influenced by changes in subglacial water pressure or subglacial water storage associated with opening and closing of water cavities in the lee of bedrock obstacles. To better understand this influence, finite-element simulations of ice flowing past bedrock obstacles with cavity formation are carried out for different synthetic periodic bedrock shapes. In the numerical model, the cavity roof is treated as an unknown free surface and is part of the solution. As an improvement over earlier studies, the cases of nonlinear ice rheology and infinite bedrock slopes are treated. Our results show that the relationship between basal drag and sliding velocity, the friction law, can be easily extended from linear to nonlinear ice rheology and is bounded even for bedrocks with locally infinite slopes. Combining our results with earlier works by others, a phenomenological friction law is proposed that includes three independent parameters that depend only on the bedrock geometry. This formulation yields an upper bound of the basal drag for finite sliding velocity and a decrease in the basal drag at low effective pressure or high velocity. This law should dramatically alter results of models of temperate glaciers and should also have important repercussions on models of glacier surges and fast glacier flows.

Iceberg calving accounts for around half of all mass loss from both the Greenland and Antarctic ice sheets. The diverse nature of calving and its complex links to both internal dynamics and climate make it challenging to incorporate into models of glaciers and ice sheets. Here we present results from a new open‐source 3‐D full‐Stokes calving model developed in Elmer/Ice. The calving model implements the crevasse depth criterion, which states that calving occurs when surface and basal crevasses penetrate the full thickness of the glacier. The model also implements a new 3‐D rediscretization approach and a time‐evolution scheme which allow the calving front to evolve realistically through time. We test the model in an application to Store Glacier, one of the largest outlet glaciers in West Greenland, and find that it realistically simulates the seasonal advance and retreat when two principal environmental forcings are applied. These forcings are (1) submarine melting in distributed and concentrated forms and (2) ice mélange buttressing. We find that ice mélange buttressing is primarily responsible for Store Glacier's seasonal advance and retreat. Distributed submarine melting prevents the glacier from forming a permanent floating tongue, while concentrated plume melting has a disproportionately large and potentially destabilizing effect on the calving front position. Our results also highlight the importance of basal topography, which exerts a strong control on calving, explaining why Store Glacier has remained stable during a period when neighboring glaciers have undergone prolonged interannual retreat.

The Fourth IPCC Assessment Report concluded that ice-sheet flow models are unable to forecast the current increase of polar ice sheet discharge and the associated contribution to sea-level rise. Since then, the glaciological community has undertaken a huge effort to develop and improve a new generation of ice-flow models, and as a result, a significant number of new ice-sheet models have emerged. Among them is the parallel finite-element model Elmer/Ice, based on the open-source multi-physics code Elmer. It was one of the first full-Stokes models used to make projections for the evolution of the whole Greenland ice sheet for the coming two centuries. Originally developed to solve local ice flow problems of high mechanical and physical complexity, Elmer/Ice has today reached the maturity to solve larger scale problems, earning the status of an ice-sheet model. Here, we summarise almost 10 yr of development performed by different groups. We present the components already included in Elmer/Ice, its numerical performance, selected applications, as well as developments planned for the future