We document the development of the first version of the U.K. Earth System Model UKESM1. The model represents a major advance on its predecessor HadGEM2‐ES, with enhancements to all component models and new feedback mechanisms. These include a new core physical model with a well‐resolved stratosphere; terrestrial biogeochemistry with coupled carbon and nitrogen cycles and enhanced land management; tropospheric‐stratospheric chemistry allowing the holistic simulation of radiative forcing from ozone, methane, and nitrous oxide; two‐moment, five‐species, modal aerosol; and ocean biogeochemistry with two‐way coupling to the carbon cycle and atmospheric aerosols. The complexity of coupling between the ocean, land, and atmosphere physical climate and biogeochemical cycles in UKESM1 is unprecedented for an Earth system model. We describe in detail the process by which the coupled model was developed and tuned to achieve acceptable performance in key physical and Earth system quantities and discuss the challenges involved in mitigating biases in a model with complex connections between its components. Overall, the model performs well, with a stable pre‐industrial state and good agreement with observations in the latter period of its historical simulations. However, global mean surface temperature exhibits stronger‐than‐observed cooling from 1950 to 1970, followed by rapid warming from 1980 to 2014. Metrics from idealized simulations show a high climate sensitivity relative to previous generations of models: Equilibrium climate sensitivity is 5.4 K, transient climate response ranges from 2.68 to 2.85 K, and transient climate response to cumulative emissions is 2.49 to 2.66 K TtC⁻¹.

The responses of carbon dioxide (CO₂) and other climate variables to an emission pulse of CO₂ into the atmosphere are often used to compute the Global Warming Potential (GWP) and Global Temperature change Potential (GTP), to characterize the response timescales of Earth System models, and to build reduced-form models. In this carbon cycle-climate model intercomparison project, which spans the full model hierarchy, we quantify responses to emission pulses of different magnitudes injected under different conditions. The CO₂ response shows the known rapid decline in the first few decades followed by a millennium-scale tail. For a 100 Gt-C emission pulse added to a constant CO₂ concentration of 389 ppm, 25 ± 9% is still found in the atmosphere after 1000 yr; the ocean has absorbed 59 ± 12% and the land the remainder (16 ± 14%). The response in global mean surface air temperature is an increase by 0.20 ± 0.12 °C within the first twenty years; thereafter and until year 1000, temperature decreases only slightly, whereas ocean heat content and sea level continue to rise. Our best estimate for the Absolute Global Warming Potential, given by the time-integrated response in CO₂ at year 100 multiplied by its radiative efficiency, is 92.5 × 10⁻¹⁵ yr W m⁻² per kg-CO₂. This value very likely (5 to 95% confidence) lies within the range of (68 to 117) × 10⁻¹⁵ yr W m⁻² per kg-CO₂. Estimates for time-integrated response in CO₂ published in the IPCC First, Second, and Fourth Assessment and our multi-model best estimate all agree within 15% during the first 100 yr. The integrated CO₂ response, normalized by the pulse size, is lower for pre-industrial conditions, compared to present day, and lower for smaller pulses than larger pulses. In contrast, the response in temperature, sea level and ocean heat content is less sensitive to these choices. Although, choices in pulse size, background concentration, and model lead to uncertainties, the most important and subjective choice to determine AGWP of CO₂ and GWP is the time horizon.

In Paris, France, December 2015, the Conference of the Parties (COP) to the United Nations Framework Convention on Climate Change (UNFCCC) invited the Intergovernmental Panel on Climate Change (IPCC) to provide a special report in 2018 on the impacts of global warming of 1.5 °C above pre-industrial levels and related global greenhouse gas emission pathways. In Nairobi, Kenya, April 2016, the IPCC panel accepted the invitation. Here we describe the response devised within the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP) to provide tailored, cross-sectorally consistent impact projections to broaden the scientific basis for the report. The simulation protocol is designed to allow for (1) separation of the impacts of historical warming starting from pre-industrial conditions from impacts of other drivers such as historical land-use changes (based on pre-industrial and historical impact model simulations); (2) quantification of the impacts of additional warming up to 1.5 °C, including a potential overshoot and long-term impacts up to 2299, and comparison to higher levels of global mean temperature change (based on the low-emissions Representative Concentration Pathway RCP2.6 and a no-mitigation pathway RCP6.0) with socio-economic conditions fixed at 2005 levels; and (3) assessment of the climate effects based on the same climate scenarios while accounting for simultaneous changes in socio-economic conditions following the middle-of-the-road Shared Socioeconomic Pathway (SSP2, Fricko et al., 2016) and in particular differential bioenergy requirements associated with the transformation of the energy system to comply with RCP2.6 compared to RCP6.0. With the aim of providing the scientific basis for an aggregation of impacts across sectors and analysis of cross-sectoral interactions that may dampen or amplify sectoral impacts, the protocol is designed to facilitate consistent impact projections from a range of impact models across different sectors (global and regional hydrology, lakes, global crops, global vegetation, regional forests, global and regional marine ecosystems and fisheries, global and regional coastal infrastructure, energy supply and demand, temperature-related mortality, and global terrestrial biodiversity).

Meteorological drought in the Hadley Centre global climate model is assessed using the Palmer Drought Severity Index (PDSI), a commonly used drought index. At interannual time scales, for the majority of the land surface, the model captures the observed relationship between the El Niño–Southern Oscillation and regions of relative wetness and dryness represented by high and low values of the PDSI respectively. At decadal time scales, on a global basis, the model reproduces the observed drying trend (decreasing PDSI) since 1952. An optimal detection analysis shows that there is a significant influence of anthropogenic emissions of greenhouse gasses and sulphate aerosols in the production of this drying trend. On a regional basis, the specific regions of wetting and drying are not always accurately simulated. In this paper, present-day drought events are defined as continuous time periods where the PDSI is less than the 20th percentile of the PDSI distribution between 1952 and 1998 (i.e., on average 20% of the land surface is in drought at any one time). Overall, the model predicts slightly less frequent but longer events than are observed. Future projections of drought in the twenty-first century made using the Special Report on Emissions Scenarios (SRES) A2 emission scenario show regions of strong wetting and drying with a net overall global drying trend. For example, the proportion of the land surface in extreme drought is predicted to increase from 1% for the present day to 30% by the end of the twenty-first century.

Climate change contributes to the increased frequency and intensity of wildfires globally, with significant impacts on society and the environment. However, our understanding of the global distribution of extreme fires remains skewed, primarily influenced by media coverage and regionalised research efforts. This inaugural State of Wildfires report systematically analyses fire activity worldwide, identifying extreme events from the March 2023–February 2024 fire season. We assess the causes, predictability, and attribution of these events to climate change and land use and forecast future risks under different climate scenarios. During the 2023–2024 fire season, 3.9×10⁶ km² burned globally, slightly below the average of previous seasons, but fire carbon (C) emissions were 16 % above average, totalling 2.4 Pg C. Global fire C emissions were increased by record emissions in Canadian boreal forests (over 9 times the average) and reduced by low emissions from African savannahs. Notable events included record-breaking fire extent and emissions in Canada, the largest recorded wildfire in the European Union (Greece), drought-driven fires in western Amazonia and northern parts of South America, and deadly fires in Hawaii (100 deaths) and Chile (131 deaths). Over 232 000 people were evacuated in Canada alone, highlighting the severity of human impact. Our analyses revealed that multiple drivers were needed to cause areas of extreme fire activity. In Canada and Greece, a combination of high fire weather and an abundance of dry fuels increased the probability of fires, whereas burned area anomalies were weaker in regions with lower fuel loads and higher direct suppression, particularly in Canada. Fire weather prediction in Canada showed a mild anomalous signal 1 to 2 months in advance, whereas events in Greece and Amazonia had shorter predictability horizons. Attribution analyses indicated that modelled anomalies in burned area were up to 40 %, 18 %, and 50 % higher due to climate change in Canada, Greece, and western Amazonia during the 2023–2024 fire season, respectively. Meanwhile, the probability of extreme fire seasons of these magnitudes has increased significantly due to anthropogenic climate change, with a 2.9–3.6-fold increase in likelihood of high fire weather in Canada and a 20.0–28.5-fold increase in Amazonia. By the end of the century, events of similar magnitude to 2023 in Canada are projected to occur 6.3–10.8 times more frequently under a medium–high emission scenario (SSP370). This report represents our first annual effort to catalogue extreme wildfire events, explain their occurrence, and predict future risks. By consolidating state-of-the-art wildfire science and delivering key insights relevant to policymakers, disaster management services, firefighting agencies, and land managers, we aim to enhance society's resilience to wildfires and promote advances in preparedness, mitigation, and adaptation. New datasets presented in this work are available from https://doi.org/10.5281/zenodo.11400539 (Jones et al., 2024) and https://doi.org/10.5281/zenodo.11420742 (Kelley et al., 2024a).