The development of parameterizations is a major task in the development of weather and climate models. Model improvement has been slow in the past decades, due to the difficulty of encompassing key physical processes into parameterizations, but also of calibrating or “tuning” the many free parameters involved in their formulation. Machine learning techniques have been recently used for speeding up the development process. While some studies propose to replace parameterizations by data‐driven neural networks, we rather advocate that keeping physical parameterizations is key for the reliability of climate projections. In this paper we propose to harness machine learning to improve physical parameterizations. In particular, we use Gaussian process‐based methods from uncertainty quantification to calibrate the model free parameters at a process level. To achieve this, we focus on the comparison of single‐column simulations and reference large‐eddy simulations over multiple boundary‐layer cases. Our method returns all values of the free parameters consistent with the references and any structural uncertainties, allowing a reduced domain of acceptable values to be considered when tuning the three‐dimensional (3D) global model. This tool allows to disentangle deficiencies due to poor parameter calibration from intrinsic limits rooted in the parameterization formulations. This paper describes the tool and the philosophy of tuning in single‐column mode. Part 2 shows how the results from our process‐based tuning can help in the 3D global model tuning.

This paper aims to understand the relationship between area level deprivation and monthly COVID-19 cases in England in response to government policy throughout 2020. The response variable is monthly reported COVID-19 cases at the Middle Super Output Area (MSOA) level by Public Health England, with Index of Multiple Deprivation (IMD), ethnicity (percentage of the population across 5 ethnicity categories) and the percentage of the population older than 70 years old and time as predictors. A GEE population-averaged panel-data model was employed to model trends in monthly COVID-19 cases with the population of each MSOA included as the exposure variable. Area level deprivation is significantly associated with COVID-19 cases from March 2020; however, this relationship is reversed in December 2020. Follow up analysis found that this reversal was maintained when controlling for the novel COVID-19 variant outbreak in the South East of England. This analysis indicates that changes in the role of deprivation and monthly reported COVID-19 over time cases may be linked to two government policies: (1) the premature easing of national restrictions in July 2020 when cases were still high in the most deprived areas in England and (2) the introduction of a regional tiered system in October predominantly in the North of England. The analysis adds to the evidence showing that deprivation is a key driver of COVID-19 outcomes and highlights the unintended negative impact of government policy.

Global assessments of air quality and health require comprehensive estimates of the exposures to air pollution that are experienced by populations in every country. However, there are many countries in which measurements from ground-based monitoring are sparse or non-existent, with quality-control and representativeness providing additional challenges. While ground-based monitoring provides a far from complete picture of global air quality, there are other sources of information that provide comprehensive coverage across the globe. The World Health Organization developed the Data Integration Model for Air Quality (DIMAQ) to combine information from ground measurements with that from other sources, such as atmospheric chemical transport models and estimates from remote sensing satellites in order to produce the information that is required for health burden assessment and the calculation of air pollution-related Sustainable Development Goals indicators. Here, we show an example of the use of DIMAQ with the Copernicus Atmosphere Monitoring Service Re-Analysis (CAMSRA) of atmospheric composition, which represents the best practices in meteorology and climate monitoring that were developed under the World Meteorological Organization’s Global Atmosphere Watch programme. Estimates of PM_2.5 from CAMSRA are integrated within the DIMAQ framework in order to produce high-resolution estimates of air pollution exposure that can be aggregated in a coherent fashion to produce country-level assessments of exposures.

Ice sheet models are used to study the deglaciation of North America at the end of the last ice age (past 21,000 years), so that we might understand whether and how existing ice sheets ma...

With the growing complexity of land surface models used to represent the terrestrial part of wider Earth system models, the need for sophisticated and robust parameter optimisation techniques is paramount. Quantifying parameter uncertainty is essential for both model development and more accurate projections. In this study, we assess the power of history matching by comparing results to the variational data assimilation approach commonly used in land surface models for parameter estimation. Although both approaches have different setups and goals, we can extract posterior parameter distributions from both methods and test the model–data fit of ensembles sampled from these distributions. Using a twin experiment, we test whether we can recover known parameter values. Through variational data assimilation, we closely match the observations. However, the known parameter values are not always contained in the posterior parameter distribution, highlighting the equifinality of the parameter space. In contrast, while more conservative, history matching still gives a reasonably good fit and provides more information about the model structure by allowing for non-Gaussian parameter distributions. Furthermore, the true parameters are contained in the posterior distributions. We then consider history matching's ability to ingest different metrics targeting different physical parts of the model, thus helping to reduce the parameter space further and improve the model–data fit. We find the best results when history matching is used with multiple metrics; not only is the model–data fit improved, but we also gain a deeper understanding of the model and how the different parameters constrain different parts of the seasonal cycle. We conclude by discussing the potential of history matching in future studies.

The development of parameterizations is a major task in the development of weather and climate models. Model improvement has been slow in the past decades, due to the difficulty of encompassing key physical processes into parameterizations, but also of calibrating or â\euro˜tuningâ\euro™ the many free parameters involved in their formulation. Machine learning techniques have been recently used for speeding up the development process. While some studies propose to replace parameterizations by data-driven neural networks, we rather advocate that keeping physical parameterizations is key for the reliability of climate projections. In this paper we propose to harness machine learning to improve physical parameterizations. In particular we use Gaussian process-based methods from uncertainty quantification to calibrate the model free parameters at a process level. To achieve this, we focus on the comparison of single-column simulations and reference large-eddy simulations over multiple boundary-layer cases. Our method returns all values of the free parameters consistent with the references and any structural uncertainties, allowing a reduced domain of acceptable values to be considered when tuning the 3D global model. This tool allows to disentangle deficiencies due to poor parameter calibration from intrinsic limits rooted in the parameterization formulations. This paper describes the tool and the philosophy of tuning in single-column mode. Part 2 shows how the results from our process-based tuning can help in the 3D global model tuning.