Search NASA⌕ Search

SEARCH · Search NASA

Results for “Cluster Computing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22

Reinterpretation of the electronic absorption spectrum of the methylene amidogen radical (H 2 CN)

The odd electronic absorption spectrum of H 2 CN has been of great interest to experiment. Herein, this system is reported on extensively by applying theoretical methods to the ground and low-lying excited electronic states. Employing a large breadth of high-level ab initio computations, including coupled cluster [CCSD(T) and CCSDT(Q)] and multireference configuration interaction [MRCISD+Q] methods, we comprehensively demonstrate that the most recent experimental and theoretical interpretations of the electronic spectrum of H 2 CN are in error. The previous assignments of the two broad features in the spectrum as the origin 0$^{0}_{0}$ (~35 050 cm –1 ) and 4$^{2}_{0}$ (~35 600 cm –1 )$\widetilde{B}$ 2 A 1 ←$\widetilde{X}$ 2 B 2 transitions are both found to be incorrect. The presently reported transition energies suggest that the higher energy band near 35 600 cm –1 is the true origin band. Moreover, from the computed anharmonic vibrational frequencies of the $\widetilde{X}$ 2 B 2 and $\widetilde{B}$ 2 A 1 states, we demonstrate that this ~550 cm –1 band spacing cannot be attributed to a simple vibronic transition, as claimed by the 4$^{2}_{0}$ assignment. Possible alternative explanations for the appearance of the lower intensity band near 35 050 cm –1 are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

SPT-CL J2040–4451: An SZ-selected galaxy cluster at z = 1.478 with significant ongoing star formation

SPT-CL J2040–4451—spectroscopically confirmed at z = 1.478—is the highest-redshift galaxy cluster yet discovered via the Sunyaev-Zel'dovich effect. SPT-CL J2040–4451 was a candidate galaxy cluster identified in the first 720 deg{sup 2} of the South Pole Telescope Sunyaev-Zel'dovich (SPT-SZ) survey, and has been confirmed in follow-up imaging and spectroscopy. From multi-object spectroscopy with Magellan-I/Baade+IMACS we measure spectroscopic redshifts for 15 cluster member galaxies, all of which have strong [O II] λλ3727 emission. SPT-CL J2040–4451 has an SZ-measured mass of M {sub 500,} {sub SZ} = 3.2 ± 0.8 × 10{sup 14} M {sub ☉} h {sub 70}{sup −1}, corresponding to M {sub 200,} {sub SZ} = 5.8 ± 1.4 × 10{sup 14} M {sub ☉} h {sub 70}{sup −1}. The velocity dispersion measured entirely from blue star-forming members is σ {sub v} = 1500 ± 520 km s{sup –1}. The prevalence of star-forming cluster members (galaxies with >1.5 M {sub ☉} yr{sup –1}) implies that this massive, high-redshift cluster is experiencing a phase of active star formation, and supports recent results showing a marked increase in star formation occurring in galaxy clusters at z ≳ 1.4. We also compute the probability of finding a cluster as rare as this in the SPT-SZ survey to be >99%, indicating that its discovery is not in tension with the concordance ΛCDM cosmological model.

79 ASTRONOMY AND ASTROPHYSICS↗

Importance of elastic finite-size effects: Neutral defects in ionic compounds

Small system sizes are a well known source of error in DFT calculations, yet computational constraints frequently dictate the use of small supercells, often as small as 96 atoms in oxides and compound semiconductors. In ionic compounds, electrostatic finite size effects have been well characterised, but self-interaction of charge neutral defects is often discounted or assumed to follow an asymptotic behaviour and thus easily corrected with linear elastic theory. Here we show that elastic effect are also important in the description of defects in ionic compounds and can lead to qualitatively incorrect conclusions if inadequatly small supercells are used; moreover, the spurious self-interaction does not follow the behaviour predicted by linear elastic theory. Considering the exemplar cases of metal oxides with fluorite structure, we show that numerous previous studies, employing 96-atom supercells, misidentify the ground state structure of (charge neutral) Schottky defects. We show that the error is eliminated by employing larger cells (324, 768 and 1500 atoms), and careful analysis determines that elastic effects, not electrostatic, are responsible. The spurious self-interaction was also observed in non-oxide ionic compounds and irrespective of the computational method used, thereby resolving long standing discrepancies between DFT and force-field methods, previously attributed to the level of theory. The surprising magnitude of the elastic effects are a cautionary tale for defect calculations in ionic materials, particularly when employing computationally expensive methods (e.g. hybrid functionals) or when modelling large defect clusters. We propose two computationally practicable methods to test the magnitude of the elastic self-interaction in any ionic system. In commonly studies oxides, where electrostatic effects would be expected to be dominant, it is the elastic effects that dictate the need for larger supercells | greater than 96 atoms.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Studying CPU and memory utilization of applications on Fujitsu A64FX and Nvidia Grace Superchip

ARM-based manycore CPU architectures are well-positioned to provide the rising memory throughput requirements of modern data intensive scientific applications in High Performance Computing (HPC). The Fujitsu A64FX CPU platform is based on the ARM v8.2A architecture, and is the processor of the flagship Japanese supercomputer - "Fugaku", which was previously ranked as the #1 supercomputer in the world according to the Top500 list. The Nvidia Grace superchip features 144 Neoverse V2 cores based on the ARMv9 architecture with 4x128b SVE2, providing exceptional computational power. The chip supports up to 480GB of memory, making it ideal for AI, machine learning, and scientific computing workloads. In this paper, we conduct a thorough performance exploration of a variety of parallel bandwidth-sensitive benchmarks and applications compiled with the native Fujitsu compiler on a Fugaku A64FX compute node and ARM (LLVM) Compiler on an NVIDIA Grace superchip compute node, engaging all the computational cores per cluster using OpenMP multithreading (assuming the cores can drive the available bandwidth). Our ultimate goals are to study the resource utilization of scientific applications and benchmarks on A64FX and Grace superchip, considering graph application scenarios ( GAP Benchmark suite) and eleven appli- cation proxies from the Rodinia heterogeneous benchmark suite (considering domains such as Data Mining, Bioinformatics, Fluid Dynamics, Pattern Recognition, etc.). Through exhaustive performance monitoring, we quantify the resource utilization of diverse OpenMP-based HPC applications on both the Fujitsu A64FX and the Nvidia Grace Superchip platforms.

benchmarking, Performance Analysis, High performan↗

Cosmological inference from an emulator based halo model. II. Joint analysis of galaxy-galaxy weak lensing and galaxy clustering from HSC-Y1 and SDSS

Here, we present high-fidelity cosmology results from a blinded joint analysis of galaxy-galaxy weak lensing (ΔΣ) and projected galaxy clustering (w p ) measured from the Hyper Suprime-Cam Year-1 (HSC-Y1) data and spectroscopic Sloan Digital Sky Survey (SDSS) galaxy catalogs in the redshift range 0.15 < z < 0.7. We define luminosity-limited samples of SDSS galaxies to serve as the tracers of w p in three spectroscopic redshift bins, and as the lens samples for ΔΣ. For the ΔΣ measurements, we select a single sample of 4×10 6 source galaxies over 140 deg 2 from HSC-Y1 with photometric redshifts (photo z) greater than 0.75, enabling a better handle of photo- z errors by comparing the ΔΣ amplitudes for the three lens redshift bins. The deep, high-quality HSC-Y1 data enable significant detections of the ΔΣ signals, with integrated signal-to-noise ratio S/N ~ 15 in the range 3 ≤ R/[h –1 Mpc] ≤ 30 for the three lens samples, despite the small area coverage. For cosmological parameter inference, we use an input galaxy-halo connection model built on the dark emulator package (which uses an ensemble set of high-resolution N-body simulations and enables fast, accurate computation of the clustering observables) with a halo occupation distribution that includes nuisance parameters to marginalize over modeling uncertainties. We model the ΔΣ and wp measurements on scales from R≃3 and 2h –1 Mpc , respectively, up to 30 h –1 Mpc (therefore excluding the baryon acoustic oscillations information) assuming a flat Λ CDM cosmology, marginalizing over about 20 nuisance parameters and demonstrating the robustness of our results to them. With various tests using mock catalogs described in Miyatake et al. [preceding paper, Phys. Rev. D 106, 083519 (2022)], we show that any bias in the clustering amplitude S 8 ≡ σ 8 (Ω m /0.3) 0.5 due to uncertainties in the galaxy-halo connection is less than ~ 50 % of the statistical uncertainty of S 8 , unless the assembly biaseffect is unexpectedly large. Our best-fit models have S 8 = 0.795$_{-0.042}^{+0.049}$ (mode and 68% credible interval) for the flat Λ CDM model; we find tighter constraints on the quantity S 8 (α = 0.17)≡σ 8 (Ω m /0.3) 0.17 = 0.745$_{-0.031}^{+0.039}$.

79 ASTRONOMY AND ASTROPHYSICS↗

Multiscale modeling of fission gas behavior in U 3 Si 2 under LWR conditions

In this work, we present a model of fission gas behavior in U 3 Si 2 under light water reactor (LWR) conditions for application in engineering fuel performance codes. The model includes components for intra-granular and inter-granular behavior of fission gases. The intra-granular component is based on cluster dynamics and computes the evolution of intra-granular fission gas bubbles and swelling coupled to gas diffusion to grain boundaries. The inter-granular component describes the evolution of grain-boundary fission gas bubbles coupled to fission gas release. Given the lack of experimental data for U 3 Si 2 under LWR conditions, the model is informed with parameters calculated via atomistic simulations. In particular, we derive fission gas diffusivities through density functional theory calculations, and the re-solution rate of fission gas atoms from intra-granular bubbles through binary collision approximation calculations. Further, the developed model is applied to the simulation of an experiment for U 3 Si 2 irradiated under LWR conditions available from the literature. Results point out a credible representation of fission gas swelling and release in U 3 Si 2 . Finally, we perform a sensitivity analysis for the various model parameters. Based on the sensitivity analysis, indications are derived that can help in addressing future research on the characterization of the physical parameters relative to fission gas behavior in U 3 Si 2 . The developed model is intended to provide a suitable infrastructure for the engineering scale calculation of fission gas behavior in U 3 Si 2 that exploits a multiscale approach to fill the experimental data gap and can be progressively improved as new lower-length scale calculations and validation data become available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Conservation of Dynamics Associated with Biological Function in an Enzyme Superfamily

Enzyme superfamily members that share common chemical and/or biological functions also share common features. While the role of structure is well characterized, the link between enzyme function and dynamics is not well understood. We present a systematic characterization of intrinsic dynamics of over 20 members of the pancreatic-type RNase superfamily, which share a common structural fold. This study is motivated by the fact that the range of chemical activity as well as molecular motions of RNase homologs spans over 105 folds. Dynamics was characterized using a combination of nuclear magnetic resonance experiments and computer simulations. Phylogenetic clustering led to the grouping of sequences into functionally distinct subfamilies. Detailed characterization of the diverse RNases showed conserved dynamical traits for enzymes within subfamilies. These results suggest that selective pressure for the conservation of dynamical behavior, among other factors, may be linked to the distinct chemical and biological functions in an enzyme superfamily.

Biochemistry & Molecular Biology↗

Vehicle energy consumption estimation using large scale simulations and machine learning methods

Vehicle energy consumption and powertrain operations for future vehicle powertrain technologies are predicted through full-vehicle simulations using Autonomie, a simulation tool developed at Argonne National Laboratory. Over 1 million vehicles, with different powertrain types (conventional and numerous electrified powertrains) and component technologies (advanced engines, batteries, electrical machines, light weighting) are simulated for outputs including vehicle fuel economy and cost. Despite the ability of Autonomic to predict fuel economy given various combination of vehicle inputs, the full-vehicle simulation is very time-consuming. It takes approximately 84 h to simulate for 33,060 vehicles even using distributed computing on a cluster of 128 worker nodes. Moreover, Autonomic is unable to populate a continuous output space through simulation since it only yields outputs with discrete values instead of an input-output relation. This paper proposes a novel large-scale learning and prediction process (LSLPP) via machine learning approaches to answer the need for a continuous space of vehicle fuel economy and a more efficient simulation process. By learning in a supervised fashion to discover structures in data, the machine learning approaches seek to train a model that can explain the data with which we can perform prediction over unseen data. As such, we are able to accomplish analytic continuation on the current discrete space of simulation results. In addition, we present a random sampling -based numerosity reduction algorithm incorporated in LSLPP that reduces simulation runs and improves efficiency of data procurement. We analytically derive a theoretical bound in order to guarantee the effectiveness of our algorithm. We also compare random sampling to stratified sampling and suggest sampling strategies. A user-friendly tool is developed to support the process. Our experimental results confirm that the proposed LSLPP is able to greatly accelerate prediction and analysis of vehicle fuel economy. The output of this study is used by the US government to evaluate the impact of its R&D funding on energy security and CO 2 emission as well as for the setting of Corporate Average Fuel Economy (CAFE) standards.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Concordant Mode Approach (CMA): Vibrational Analysis of New and Upgraded Intermolecular Benchmarks for Noncovalent Bonding

The Concordant Mode Approach (CMA) is a novel method that offers tremendous potential for increasing the system size and the level of theory attainable in quantum chemical computations of molecular vibrational frequencies. To investigate the extension of CMA to intermolecular vibrations, computations with coupled cluster singles and doubles with perturbative triples theory [CCSD(T)] using two augmented correlation-consistent polarized-valence triple-ζ basis sets (aug-cc-pVTZ or h-aug-ccpVTZ) were performed on 17 prototypical loosely bound complexes of hydrogen-bonded, dispersion, and mixed character. These Level A results provide new and upgraded benchmarks for noncovalent bonding and a severe test for CMA vibrational analyses. The Level A target frequencies were recovered remarkably well using second-order Møller−Plesset perturbation theory (MP2) with h-aug-cc-pVTZ for generating the underlying (Level B) normal modes of the CMA scheme. Employing this Level B within the lowest-rung CMA-0A method reproduces the 435 benchmark frequencies with a mean absolute error (MAE) of 0.23 cm −1 and a corresponding standard deviation (σ) of 0.84 cm −1 ; strikingly, the corresponding subset of 106 interfragment frequencies exhibits MAE = 0.34 cm −1 and σ = 0.90 cm −1 . Subsequent application of the higher-rung CMA-2A scheme eliminates all outliers and reduces the overall MAE to a minuscule 0.08 cm −1 with the inclusion of only 3.0% of the off-diagonal couplings not accounted for by CMA-0A. Accordingly, the highly efficient CMA methodology proves to be robust even for vibrations on flat potential energy surfaces.

Aromatic compounds↗

Crystallization of glass-forming liquids: Specific surface energy

A generalization of the Stefan-Skapski-Turnbull relation for the melt-crystal specific interfacial energy is developed in terms of the generalized Gibbs approach extending its standard formulation to thermodynamic non-equilibrium states. With respect to crystal nucleation, this relation is required in order to determine the parameters of the critical crystal clusters being a prerequisite for the computation of the work of critical cluster formation. As one of its consequences, a relation for the dependence of the specific surface energy of critical clusters on temperature and pressure is derived applicable for small and moderate deviations from liquid-crystal macroscopic equilibrium states. Employing the Stefan-Skapski-Turnbull relation, general expressions for the size and the work of formation of critical crystal clusters are formulated. The resulting expressions are much more complex as compared to the respective relations obtained via the classical Gibbs theory. Latter relations are retained as limiting cases of these more general expressions for moderate undercoolings. By this reason, the formulated, here, general relations for the specification of the critical cluster size and the work of critical cluster formation give a key for an appropriate interpretation of a variety of crystallization phenomena occurring at large undercoolings which cannot be understood in terms of the Gibbs’ classical treatment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

VIBES: a workflow for annotating and visualizing viral sequences integrated into bacterial genomes

Abstract Bacteriophages are viruses that infect bacteria. Many bacteriophages integrate their genomes into the bacterial chromosome and become prophages. Prophages may substantially burden or benefit host bacteria fitness, acting in some cases as parasites and in others as mutualists. Some prophages have been demonstrated to increase host virulence. The increasing ease of bacterial genome sequencing provides an opportunity to deeply explore prophage prevalence and insertion sites. Here we present VIBES (Viral Integrations in Bacterial genomES), a workflow intended to automate prophage annotation in complete bacterial genome sequences. VIBES provides additional context to prophage annotations by annotating bacterial genes and viral proteins in user-provided bacterial and viral genomes. The VIBES pipeline is implemented as a Nextflow-driven workflow, providing a simple, unified interface for execution on local, cluster and cloud computing environments. For each step of the pipeline, a container including all necessary software dependencies is provided. VIBES produces results in simple tab-separated format and generates intuitive and interactive visualizations for data exploration. Despite VIBES’s primary emphasis on prophage annotation, its generic alignment-based design allows it to be deployed as a general-purpose sequence similarity search manager. We demonstrate the utility of the VIBES prophage annotation workflow by searching for 178 Pf phage genomes across 1072 Pseudomonas spp. genomes.

59 BASIC BIOLOGICAL SCIENCES↗

Δ isobars and nuclear saturation

In this paper, we construct a nuclear interaction in chiral effective field theory with explicit inclusion of the $\mathrm{{\Delta}}$-isobar $\mathrm{{\Delta}}(1232)$ degree of freedom at all orders up to next-to-next-to-leading order (NNLO). We use pion-nucleon (${\pi}N$) low-energy constants (LECs) from a Roy-Steiner analysis of ${\pi}N$ scattering data, optimize the LECs in the contact potentials up to NNLO to reproduce low-energy nucleon-nucleon scattering phase shifts, and constrain the three-nucleon interaction at NNLO to reproduce the binding energy and point-proton radius of $^{4}\mathrm{He}$. For heavier nuclei we use the coupled-cluster method to compute binding energies, radii, and neutron skins. We find that radii and binding energies are much improved for interactions with explicit inclusion of $\mathrm{{\Delta}}(1232)$, while $\mathrm{{\Delta}}$-less interactions produce nuclei that are not bound with respect to breakup into ${\alpha}$ particles. Finally, the saturation of nuclear matter is significantly improved, and its symmetry energy is consistent with empirical estimates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Cast of Thousands: How the IDEAS Productivity Project Has Advanced Software Productivity and Sustainability

Computational and data-enabled science and engineering are revolutionizing advances throughout science and society, at all scales of computing. For example, teams in the U.S. Department of Energy’s Exascale Computing Project have been tackling new frontiers in modeling, simulation, and analysis by exploiting unprecedented exascale computing capabilities—building an advanced software ecosystem that supports next-generation applications and addresses disruptive changes in computer architectures. However, concerns are growing about the productivity of the developers of scientific software. Members of the Interoperable Design of Extreme-scale Application Software project serve as catalysts to address these challenges through fostering software communities, incubating and curating methodologies and resources, and disseminating knowledge to advance developer productivity and software sustainability. This article discusses how these synergistic activities are advancing scientific discovery—mitigating technical risks by building a firmer foundation for reproducible, sustainable science at all scales of computing, from laptops to clusters to exascale and beyond.

97 MATHEMATICS AND COMPUTING↗

panhandle

A project to provide user activity monitoring for High Performance Computing systems and clusters. The goal is to provide effective user activity monitoring with minimal performance impact on the host running this service.

McGee, David [@LANL @USMC @DoD]↗

The interaction of RNA helicase DDX3 with HIV-1 Rev-CRM1-RanGTP complex during the HIV replication cycle

Molecular traffic between the nucleus and the cytoplasm is regulated by the nuclear pore complex (NPC), which acts as a highly selective channel perforating the nuclear envelope in eukaryotic cells. The human immunodeficiency virus (HIV) exploits the nucleocytoplasmic pathway to export its RNA transcripts across the NPC to the cytoplasm. Despite extensive study on the HIV life cycle and the many drugs developed to target this cycle, no current drugs have been successful in targeting the critical process of viral nuclear export, even though HIV’s reliance on a single host protein, CRM1, to export its unspliced and partially spliced RNA transcripts makes it a tempting target. Due to recent findings implicating a DEAD-box helicase, DDX3, in HIV replication and a member of the export complex, it has become an appealing target for anti-HIV drug inhibition. In the present research, we have applied a hybrid computational protocol to analyze protein-protein interactions in the HIV mRNA export cycle. This method is based on molecular docking followed by molecular dynamics simulation and accompanied by approximate free energy calculation (MM/GBSA), computational alanine scanning, clustering, and evolutionary analysis. We highlight here some of the most likely binding modes and interfacial residues between DDX3 and CRM1 both in the absence and presence of RanGTP. This work shows that although DDX3 can bind to free CRM1, addition of RanGTP leads to more concentrated distribution of binding modes and stronger binding between CRM1 and RanGTP.

59 BASIC BIOLOGICAL SCIENCES↗

Using containers to speed up development, to run integration tests and to teach about distributed systems

GlideinWMS is a workload manager provisioning resources for many experiments including CMS and DUNE. The software is distributed both as native packages and specialized production containers. Following an approach used in other communities like web development we built our workspaces, system-like containers to ease development and testing. Developers can change the source tree or check out a different branch and quickly reconfigure the services to see the effect of their changes. In this paper, we’ll talk about what differentiates workspaces from other containers. We’ll describe our base system composed of three containers. A one-node cluster including a compute element and a batch system. A GlideinWMS Factory controlling pilot jobs. And a scheduler and Frontend, to submit jobs and provision resources. Additional containers can be used for optional components. This system can easily run on a laptop and we’ll share our evaluation of different container runtimes, with an eye for ease of use and performance. Finally, we’ll talk about our experience as developers and with students. The GlideinWMS workspaces are easily integrated with IDEs like VS Code, simplifying debugging and allowing development and testing of the system also when offline. They simplified the training and onboarding of new team members and Summer interns. And they were useful in workshops where students could have first-hand experience with the mechanisms and components that, in production, run millions of jobs.

Mambelli, Marco↗

Using Containers to Speed Up Development, to Run Integration Tests and to Teach About Distributed Systems

GlideinWMS is a workload manager provisioning resources for many experiments, including CMS and DUNE. The software is distributed both as native packages and specialized production containers. Following an approach used in other communities like web development, we built our workspaces, system-like containers to ease development and testing. Developers can change the source tree or check out a different branch and quickly reconfigure the services to see the effect of their changes. In this paper, we will talk about what differentiates workspaces from other containers. We will describe our base system, composed of three containers: a one-node cluster including a compute element and a batch system, a GlideinWMS Factory controlling pilot jobs, and a scheduler and Frontend to submit jobs and provision resources. Additional containers can be used for optional components. This system can easily run on a laptop, and we will share our evaluation of different container runtimes, with an eye for ease of use and performance. Finally, we will talk about our experience as developers and with students. The GlideinWMS workspaces are easily integrated with IDEs like VS Code, simplifying debugging and allowing development and testing of the system even when offline. They simplified the training and onboarding of new team members and summer interns. And they were useful in workshops where students could have first-hand experience with the mechanisms and components that, in production, run millions of jobs.

Mambelli, Marco [Fermilab] (ORCID:0000000294892681↗

A new biogeochemical modelling framework (FLaMe-v1.0) for lake methane emissions on the regional scale: development and application to the European domain

This study presents a new physical-biogeochemical modelling framework for simulating lake methane (CH 4 ) emissions at regional scales. The new model, FLaMe-v1.0 (Fluxes of Lake Methane), rests on an innovative, computationally efficient lake clustering approach that enables the simulation of CH 4 emissions across a large number of lakes. Building on the Canadian Small Lake Model (CSLM) that simulates the lake physics, we develop a suite of biogeochemical modules to simulate transient dynamics of organic Carbon (C), Oxygen (O 2 ), and CH 4 . We first test the performance of FLaMe-v1.0 by analyzing physical and biogeochemical processes in two theoretical lakes with characteristics that can be considered representative for many lakes (an oligotrophic, deep lake driven by cold climate versus a eutrophic, shallow lake driven by warm climate). Next, we evaluate the model by comparing simulated and observed timeseries of CH 4 emissions in four well-surveyed lakes. We then apply FLaMe-v1.0 at the European scale to evaluate simulated diffusive and ebullitive lake CH 4 fluxes against in-situ measurements in both boreal and central European regions. Finally, we provide a first assessment of the spatio-temporal variability in CH 4 emissions from European lakes with a surface area comprised between 0.1–1000 km 2 (n= 108 407, total area = 1.33 × 105 km 2 ), indicating a total emission of 0.97 ± 0.23 Tg CH 4 yr −1 , with the uncertainty constrained by combining FLaMe-v1.0 and machine learning techniques. Moreover, 30 % and 70 % of these CH 4 emissions are through diffusive and ebullitive pathways, respectively. Annually averaged CH 4 emission rates per unit lake area during 2010–2016 have a South-to-North decreasing gradient, resulting in a mean over the European domain as 7.39 g CH 4 m −2 yr −1 . Our simulations reveal a strong seasonality (with ice-blocking effects accounted for) in European lake CH 4 emissions, with nearly ten times higher emissions during late summer than during winter. This pronounced seasonal variation highlights the importance of accounting for the sub-annual variability in CH 4 emissions to accurately constrain regional CH 4 budgets. In the future, FLaMe-v1.0 could be embedded into Earth System Models to investigate the feedback between climate warming and global lake CH 4 emissions.

Maisonnier, Manon [Free Univ. of Brussels (Belgium↗