A confined–unconfined aquifer model for subglacial hydrology and its application to the Northeast Greenland Ice Stream

Subglacial hydrology plays an important role in ice sheet dynamics as it determines the sliding velocity. It also drives freshwater into the ocean, leading to undercutting of calving fronts by plumes. Modeling subglacial water has been a challenge for decades. Only recently have new approaches been developed such as representing subglacial channels and thin water sheets by separate layers of variable hydraulic conductivity. We extend this concept by modeling a confined–unconfined aquifer system (CUAS) in a single layer of an equivalent porous medium (EPM). The advantage of this formulation is that it prevents unphysical values of pressure at reasonable computational cost. We performed sensitivity tests to investigate the effect of different model parameters. The strongest influence of model parameters was detected in terms of governing the opening and closure of the system. Furthermore, we applied the model to the Northeast Greenland Ice Stream, where an efficient system independent of seasonal input was identified about 500&thinsp;km downstream from the ice divide. Using the effective pressure from the hydrology model, the Ice Sheet System Model (ISSM) showed considerable improvements in modeled velocities in the coastal region.</p

Groundwater data improve modelling of headwater stream CO₂ outgassing with a stable DIC isotope approach

A large portion of terrestrially derived carbon outgasses as carbon dioxide (CO2) from streams and rivers to the atmosphere. Particularly, the amount of CO2 outgassing from small headwater streams is highly uncertain. Conservative estimates suggest that they contribute 36 % (i.e. 0.93 petagrams (Pg) C yr−1) of total CO2 outgassing from all fluvial ecosystems on the globe. In this study, stream pCO2, dissolved inorganic carbon (DIC), and δ13CDIC data were used to determine CO2 outgassing from an acidic headwater stream in the Uhlířská catchment (Czech Republic). This stream drains a catchment with silicate bedrock. The applied stable isotope model is based on the principle that the 13C ∕ 12C ratio of its sources and the intensity of CO2 outgassing control the isotope ratio of DIC in stream water. It avoids the use of the gas transfer velocity parameter (k), which is highly variable and mostly difficult to constrain. Model results indicate that CO2 outgassing contributed more than 80 % to the annual stream inorganic carbon loss in the Uhlířská catchment. This translated to a CO2 outgassing rate from the stream of 34.9 kg C m−2 yr−1 when normalised to the stream surface area. Large temporal variations with maximum values shortly before spring snowmelt and in summer emphasise the need for investigations at higher temporal resolution. We improved the model uncertainty by incorporating groundwater data to better constrain the isotope compositions of initial DIC. Due to the large global abundance of acidic, humic-rich headwaters, we underline the importance of this integral approach for global applications

Scalability and some optimization of the Finite-volumE Sea ice–Ocean Model, Version 2.0 (FESOM2)

A study of the scalability of the Finite-volumE Sea ice–Ocean circulation Model, Version 2.0 (FESOM2), the first mature global model of its kind formulated on unstructured meshes, is presented. This study includes an analysis of the main computational kernels with a special focus on bottlenecks in parallel scalability. Several model enhancements improving this scalability for large numbers of processes are described and tested. Model grids at different resolutions are used on four high-performance computing (HPC) systems with differing computational and communication hardware to demonstrate the model's scalability and throughput. Furthermore, strategies for improvements in parallel performance are presented and assessed. We show that, in terms of throughput, FESOM2 is on a par with state-of-the-art structured ocean models and, in a realistic eddy-resolving configuration (1/10∘ resolution), can achieve about 16 years per day on 14 000 cores. This suggests that unstructured-mesh models are becoming very competitive tools in high-resolution climate modeling. We show that the main bottlenecks of FESOM2 parallel scalability are the two-dimensional components of the model, namely the computations of the external (barotropic) mode and the sea-ice model. It is argued that these bottlenecks are shared with other general ocean circulation models

Multiscale modeling of coastal, shelf, and global ocean dynamics

In contemporary ocean science, modeling systems that integrate understanding of complex multiscale phenomena and utilize efficient numerics are paramount. Many of today's fundamental ocean science questions involve multiple scales and multiple dynamics. A new generation of modeling systems would allow to study such questions quantitatively by being less restrictive dynamically and more efficient numerically than more traditional systems