The Coupled Model Intercomparison Project (CMIP) is one of the biggest international efforts aimed at better understanding the past, present, and future of climate changes in a multi-model context. A total of 21 model intercomparison projects (MIPs) were endorsed in its sixth phase (CMIP6), which included 190 different experiments that were used to simulate 40 000 years and produced around 40 PB of data in total. This paper presents the main findings obtained from the CPMIP (the Computational Performance Model Intercomparison Project), a collection of a common set of metrics, specifically designed for assessing climate model performance. These metrics were exclusively collected from the production runs of experiments used in CMIP6 and primarily from institutions within the IS-ENES3 consortium. The document presents the full set of CPMIP metrics per institution and experiment, including a detailed analysis and discussion of each of the measurements. During the analysis, we found a positive correlation between the core hours needed, the complexity of the models, and the resolution used. Likewise, we show that between 5 %–15 % of the execution cost is spent in the coupling between independent components, and it only gets worse by increasing the number of resources. From the data, it is clear that queue times have a great impact on the actual speed achieved and have a huge variability across different institutions, ranging from none to up to 78 % execution overhead. Furthermore, our evaluation shows that the estimated carbon footprint of running such big simulations within the IS-ENES3 consortium is 1692 t of CO₂ equivalent. As a result of the collection, we contribute to the creation of a comprehensive database for future community reference, establishing a benchmark for evaluation and facilitating the multi-model, multi-platform comparisons crucial for understanding climate modelling performance. Given the diverse range of applications, configurations, and hardware utilised, further work is required for the standardisation and formulation of general rules. The paper concludes with recommendations for future exercises aimed at addressing the encountered challenges which will facilitate more collections of a similar nature.

We review how the international modelling community, encompassing integrated assessment models, global and regional Earth system and climate models, and impact models, has worked together over the past few decades to advance understanding of Earth system change and its impacts on society and the environment and thereby support international climate policy. We go on to recommend a number of priority research areas for the coming decade, a timescale that encompasses a number of newly starting international modelling activities, as well as the IPCC Seventh Assessment Report (AR7) and the second UNFCCC Global Stocktake. Progress in these priority areas will significantly advance our understanding of Earth system change and its impacts, increasing the quality and utility of science support to climate policy. We emphasize the need for continued improvement in our understanding of, and ability to simulate, the coupled Earth system and the impacts of Earth system change. There is an urgent need to investigate plausible pathways and emission scenarios that realize the Paris climate targets – for example, pathways that overshoot 1.5 or 2 °C global warming, before returning to these levels at some later date. Earth system models need to be capable of thoroughly assessing such warming overshoots – in particular, the efficacy of mitigation measures, such as negative CO₂ emissions, in reducing atmospheric CO₂ and driving global cooling. An improved assessment of the long-term consequences of stabilizing climate at 1.5 or 2 °C above pre-industrial temperatures is also required. We recommend Earth system models run overshoot scenarios in CO₂-emission mode to more fully represent coupled climate–carbon-cycle feedbacks and, wherever possible, interactively simulate other key Earth system phenomena at risk of rapid change during overshoot. Regional downscaling and impact models should use forcing data from these simulations, so impact and regional climate projections cover a more complete range of potential responses to a warming overshoot. An accurate simulation of the observed, historical record remains a fundamental requirement of models, as does accurate simulation of key metrics, such as the effective climate sensitivity and the transient climate response to cumulative carbon emissions. For adaptation, a key demand is improved guidance on potential changes in climate extremes and the modes of variability these extremes develop within. Such improvements will most likely be realized through a combination of increased model resolution, improvement of key model parameterizations, and enhanced representation of important Earth system processes, combined with targeted use of new artificial intelligence (AI) and machine learning (ML) techniques. We propose a deeper collaboration across such efforts over the coming decade. With respect to sampling future uncertainty, increased collaboration between approaches that emphasize large model ensembles and those focussed on statistical emulation is required. We recommend an increased focus on high-impact–low-likelihood (HILL) outcomes – in particular, the risk and consequences of exceeding critical tipping points during a warming overshoot and the potential impacts arising from this. For a comprehensive assessment of the impacts of Earth system change, including impacts arising directly as a result of climate mitigation actions, it is important that spatially detailed, disaggregated information used to generate future scenarios in integrated assessment models be available for use in impact models. Conversely, there is a need to develop methods that enable potential societal responses to projected Earth system change to be incorporated into scenario development. The new models, simulations, data, and scientific advances proposed in this article will not be possible without long-term development and maintenance of a robust, globally connected infrastructure ecosystem. This system must be easily accessible and useable by modelling communities across the world, allowing the global research community to be fully engaged in developing and delivering new scientific knowledge to support international climate policy.

A catalog containing details of the highest-energy cosmic rays recorded through the detection of extensive air showers at the Pierre Auger Observatory is presented with the aim of opening the data to detailed examination. Descriptions of the 100 showers created by the highest-energy particles recorded between 2004 January 1 and 2020 December 31 are given for cosmic rays that have energies in the range 78–166 EeV. Details are also given on a further nine very energetic events that have been used in the calibration procedure adopted to determine the energy of each primary. A sky plot of the arrival directions of the most energetic particles is shown. No interpretations of the data are offered.

The flux of ultra-high energy cosmic rays reaching Earth above the ankle energy (5 EeV) can be described as a mixture of nuclei injected by extragalactic sources with very hard spectra and a low rigidity cutoff. Extragalactic magnetic fields existing between the Earth and the closest sources can affect the observed CR spectrum by reducing the flux of low-rigidity particles reaching Earth. We perform a combined fit of the spectrum and distributions of depth of shower maximum measured with the Pierre Auger Observatory including the effect of this magnetic horizon in the propagation of UHECRs in the intergalactic space. We find that, within a specific range of the various experimental and phenomenological systematics, the magnetic horizon effect can be relevant for turbulent magnetic field strengths in the local neighbourhood in which the closest sources lie of order B_rms ≃ (50–100) nG (20 Mpc/d_s)( 100 kpc/L_coh)^1/2, with d_s the typical intersource separation and L_coh the magnetic field coherence length. When this is the case, the inferred slope of the source spectrum becomes softer and can be closer to the expectations of diffusive shock acceleration, i.e., ∝ E^-2. An additional cosmic-ray population with higher source density and softer spectra, presumably also extragalactic and dominating the cosmic-ray flux at EeV energies, is also required to reproduce the overall spectrum and composition results for all energies down to 0.6 EeV.

This paper presents the message passing interface (MPI)-based parallelization of the three-dimensional hydrodynamic model SHYFEM (System of HydrodYnamic Finite Element Modules). The original sequential version of the code was parallelized in order to reduce the execution time of high-resolution configurations using state-of-the-art high-performance computing (HPC) systems. A distributed memory approach was used, based on the MPI. Optimized numerical libraries were used to partition the unstructured grid (with a focus on load balancing) and to solve the sparse linear system of equations in parallel in the case of semi-to-fully implicit time stepping. The parallel implementation of the model was validated by comparing the outputs with those obtained from the sequential version. The performance assessment demonstrates a good level of scalability with a realistic configuration used as benchmark.