17 results
Search Results
Now showing 1 - 10 of 17
- ItemEvaluating stratospheric ozone and water vapour changes in CMIP6 models from 1850 to 2100(Katlenburg-Lindau : European Geosciences Union, 2021) Keeble, James; Hassler, Birgit; Banerjee, Antara; Checa-Garcia, Ramiro; Chiodo, Gabriel; Davis, Sean; Eyring, Veronika; Griffiths, Paul T.; Morgenstern, Olaf; Nowack, Peer; Zeng, Guang; Zhang, Jiankai; Bodeker, Greg; Burrows, Susannah; Cameron-Smith, Philip; Cugnet, David; Danek, Christopher; Deushi, Makoto; Horowitz, Larry W.; Kubin, Anne; Li, Lijuan; Lohmann, Gerrit; Michou, Martine; Mills, Michael J.; Nabat, Pierre; Olivié, Dirk; Park, Sungsu; Seland, Øyvind; Stoll, Jens; Wieners, Karl-Hermann; Wu, TongwenStratospheric ozone and water vapour are key components of the Earth system, and past and future changes to both have important impacts on global and regional climate. Here, we evaluate long-term changes in these species from the pre-industrial period (1850) to the end of the 21st century in Coupled Model Intercomparison Project phase 6 (CMIP6) models under a range of future emissions scenarios. There is good agreement between the CMIP multi-model mean and observations for total column ozone (TCO), although there is substantial variation between the individual CMIP6 models. For the CMIP6 multi-model mean, global mean TCO has increased from ∼300 DU in 1850 to ∼ 305 DU in 1960, before rapidly declining in the 1970s and 1980s following the use and emission of halogenated ozone-depleting substances (ODSs). TCO is projected to return to 1960s values by the middle of the 21st century under the SSP2-4.5, SSP3-7.0, SSP4-3.4, SSP4-6.0, and SSP5-8.5 scenarios, and under the SSP3-7.0 and SSP5-8.5 scenarios TCO values are projected to be ∼ 10 DU higher than the 1960s values by 2100. However, under the SSP1-1.9 and SSP1-1.6 scenarios, TCO is not projected to return to the 1960s values despite reductions in halogenated ODSs due to decreases in tropospheric ozone mixing ratios. This global pattern is similar to regional patterns, except in the tropics where TCO under most scenarios is not projected to return to 1960s values, either through reductions in tropospheric ozone under SSP1-1.9 and SSP1-2.6, or through reductions in lower stratospheric ozone resulting from an acceleration of the Brewer-Dobson circulation under other Shared Socioeconomic Pathways (SSPs). In contrast to TCO, there is poorer agreement between the CMIP6 multi-model mean and observed lower stratospheric water vapour mixing ratios, with the CMIP6 multi-model mean underestimating observed water vapour mixing ratios by ∼ 0.5 ppmv at 70 hPa. CMIP6 multi-model mean stratospheric water vapour mixing ratios in the tropical lower stratosphere have increased by ∼ 0.5 ppmv from the pre-industrial to the present-day period and are projected to increase further by the end of the 21st century. The largest increases (∼ 2 ppmv) are simulated under the future scenarios with the highest assumed forcing pathway (e.g. SSP5-8.5). Tropical lower stratospheric water vapour, and to a lesser extent TCO, shows large variations following explosive volcanic eruptions. © Author(s) 2021.
- ItemImportance of size representation and morphology in modelling optical properties of black carbon: comparison between laboratory measurements and model simulations(Katlenburg-Lindau : Copernicus, 2022) Romshoo, Baseerat; Pöhlker, Mira; Wiedensohler, Alfred; Pfeifer, Sascha; Saturno, Jorge; Nowak, Andreas; Ciupek, Krzysztof; Quincey, Paul; Vasilatou, Konstantina; Ess, Michaela N.; Gini, Maria; Eleftheriadis, Konstantinos; Robins, Chris; Gaie-Levrel, François; Müller, ThomasBlack carbon (BC) from incomplete combustion of biomass or fossil fuels is the strongest absorbing aerosol component in the atmosphere. Optical properties of BC are essential in climate models for quantification of their impact on radiative forcing. The global climate models, however, consider BC to be spherical particles, which causes uncertainties in their optical properties. Based on this, an increasing number of model-based studies provide databases and parameterization schemes for the optical properties of BC, using more realistic fractal aggregate morphologies. In this study, the reliability of the different modelling techniques of BC was investigated by comparing them to laboratory measurements. The modelling techniques were examined for bare BC particles in the first step and for BC particles with organic material in the second step. A total of six morphological representations of BC particles were compared, three each for spherical and fractal aggregate morphologies. In general, the aggregate representation performed well for modelling the particle light absorption coefficient σabs, single-scattering albedo SSA, and mass absorption cross-section MACBC for laboratory-generated BC particles with volume mean mobility diameters dp,V larger than 100nm. However, for modelling Ångström absorption exponent AAE, it was difficult to suggest a method due to size dependence, although the spherical assumption was in better agreement in some cases. The BC fractal aggregates are usually modelled using monodispersed particles, since their optical simulations are computationally expensive. In such studies, the modelled optical properties showed a 25% uncertainty in using the monodisperse size method. It is shown that using the polydisperse size distribution in combination with fractal aggregate morphology reduces the uncertainty in measured σabs to 10% for particles with dp,V between 60-160nm. Furthermore, the sensitivities of the BC optical properties to the various model input parameters such as the real and imaginary parts of the refractive index (mre and mim), the fractal dimension (Df), and the primary particle radius (app) of an aggregate were investigated. When the BC particle is small and rather fresh, the change in the Df had relatively little effect on the optical properties. There was, however, a significant relationship between app and the particle light scattering, which increased by a factor of up to 6 with increasing total particle size. The modelled optical properties of BC are well aligned with laboratory-measured values when the following assumptions are used in the fractal aggregate representation: mre between 1.6 and 2, mim between 0.50 and 1, Df from 1.7 to 1.9, and app between 10 and 14nm. Overall, this study provides experimental support for emphasizing the importance of an appropriate size representation (polydisperse size method) and an appropriate morphological representation for optical modelling and parameterization scheme development of BC.
- ItemComparison of particle number size distribution trends in ground measurements and climate models(Katlenburg-Lindau : EGU, 2022) Leinonen, Ville; Kokkola, Harri; Yli-Juuti, Taina; Mielonen, Tero; Kühn, Thomas; Nieminen, Tuomo; Heikkinen, Simo; Miinalainen, Tuuli; Bergman, Tommi; Carslaw, Ken; Decesari, Stefano; Fiebig, Markus; Hussein, Tareq; Kivekäs, Niku; Krejci, Radovan; Kulmala, Markku; Leskinen, Ari; Massling, Andreas; Mihalopoulos, Nikos; Mulcahy, Jane P.; Noe, Steffen M.; van Noije, Twan; O'Connor, Fiona M.; O'Dowd, Colin; Olivie, Dirk; Pernov, Jakob B.; Petäjä, Tuukka; Seland, Øyvind; Schulz, Michael; Scott, Catherine E.; Skov, Henrik; Swietlicki, Erik; Tuch, Thomas; Wiedensohler, Alfred; Virtanen, Annele; Mikkonen, SanttuDespite a large number of studies, out of all drivers of radiative forcing, the effect of aerosols has the largest uncertainty in global climate model radiative forcing estimates. There have been studies of aerosol optical properties in climate models, but the effects of particle number size distribution need a more thorough inspection. We investigated the trends and seasonality of particle number concentrations in nucleation, Aitken, and accumulation modes at 21 measurement sites in Europe and the Arctic. For 13 of those sites, with longer measurement time series, we compared the field observations with the results from five climate models, namely EC-Earth3, ECHAM-M7, ECHAM-SALSA, NorESM1.2, and UKESM1. This is the first extensive comparison of detailed aerosol size distribution trends between in situ observations from Europe and five earth system models (ESMs). We found that the trends of particle number concentrations were mostly consistent and decreasing in both measurements and models. However, for many sites, climate models showed weaker decreasing trends than the measurements. Seasonal variability in measured number concentrations, quantified by the ratio between maximum and minimum monthly number concentration, was typically stronger at northern measurement sites compared to other locations. Models had large differences in their seasonal representation, and they can be roughly divided into two categories: for EC-Earth and NorESM, the seasonal cycle was relatively similar for all sites, and for other models the pattern of seasonality varied between northern and southern sites. In addition, the variability in concentrations across sites varied between models, some having relatively similar concentrations for all sites, whereas others showed clear differences in concentrations between remote and urban sites. To conclude, although all of the model simulations had identical input data to describe anthropogenic mass emissions, trends in differently sized particles vary among the models due to assumptions in emission sizes and differences in how models treat size-dependent aerosol processes. The inter-model variability was largest in the accumulation mode, i.e. sizes which have implications for aerosol-cloud interactions. Our analysis also indicates that between models there is a large variation in efficiency of long-range transportation of aerosols to remote locations. The differences in model results are most likely due to the more complex effect of different processes instead of one specific feature (e.g. the representation of aerosol or emission size distributions). Hence, a more detailed characterization of microphysical processes and deposition processes affecting the long-range transport is needed to understand the model variability.
- ItemCRI-HOM: A novel chemical mechanism for simulating highly oxygenated organic molecules (HOMs) in global chemistry-aerosol-climate models(Katlenburg-Lindau : EGU, 2020) Weber, James; Archer-Nicholls, Scott; Griffiths, Paul; Berndt, Torsten; Jenkin, Michael; Gordon, Hamish; Knote, Christoph; Archibald, Alexander T.We present here results from a new mechanism, CRI-HOM, which we have developed to simulate the formation of highly oxygenated organic molecules (HOMs) from the gas-phase oxidation of-pinene, one of the most widely emitted biogenic volatile organic compounds (BVOCs) by mass. This concise scheme adds 12 species and 66 reactions to the Common Representative Intermediates (CRI) mechanism v2.2 Reduction 5 and enables the representation of semi-explicit HOM treatment suitable for long-term global chemistry-aerosol-climate modelling, within a comprehensive tropospheric chemical mechanism. The key features of the new mechanism are (i) representation of the autoxidation of peroxy radicals from the hydroxyl radical and ozone initiated reactions of-pinene, (ii) formation of multiple generations of peroxy radicals, (iii) formation of accretion products (dimers), and (iv) isoprene-driven suppression of accretion product formation, as observed in experiments. The mechanism has been constructed through optimisation against a series of flow tube laboratory experiments. The mechanism predicts a HOM yield of 2 %-4.5% under conditions of low to moderate NOx , in line with experimental observations, and reproduces qualitatively the decline in HOM yield and concentration at higher NOx levels. The mechanism gives a HOM yield that also increases with temperature, in line with observations, and our mechanism compares favourably to some of the limited observations of [HOM] observed in the boreal forest in Finland and in the southeast USA. The reproduction of isoprene-driven suppression of HOMs is a key step forward as it enables global climate models to capture the interaction between the major BVOC species, along with the potential climatic feedbacks. This suppression is demonstrated when the mechanism is used to simulate atmospheric profiles over the boreal forest and rainforest; different isoprene concentrations result in different [HOM] distributions, illustrating the importance of BVOC interactions in atmospheric composition and climate. Finally particle nucleation rates calculated from [HOM] in present-day and preindustrial atmospheres suggest that "sulfuric-acid-free"nucleation can compete effectively with other nucleation pathways in the boreal forest, particularly in the pre-industrial period, with important implications for the aerosol budget and radiative forcing. © Author(s) 2020.
- ItemHow the extreme 2019-2020 Australian wildfires affected global circulation and adjustments(Katlenburg-Lindau : EGU, 2023) Senf, Fabian; Heinold, Bernd; Kubin, Anne; Müller, Jason; Schrödner, Roland; Tegen, InaWildfires are a significant source of absorbing aerosols in the atmosphere. Extreme fires in particular, such as those during the 2019-2020 Australian wildfire season (Black Summer fires), can have considerable large-scale effects. In this context, the climate impact of extreme wildfires unfolds not only because of the emitted carbon dioxide but also due to smoke aerosol released up to an altitude of 17ĝ€¯km. The overall aerosol effects depend on a variety of factors, such as the amount emitted, the injection height, and the composition of the burned material, and is therefore subject to considerable uncertainty. In the present study, we address the global impact caused by the exceptionally strong and high-reaching smoke emissions from the Australian wildfires using simulations with a global aerosol-climate model. We show that the absorption of solar radiation by the black carbon contained in the emitted smoke led to a shortwave radiative forcing of more than +5ĝ€¯Wm-2 in the southern mid-latitudes of the lower stratosphere. Subsequent adjustment processes in the stratosphere slowed down the diabatically driven meridional circulation, thus redistributing the heating perturbation on a global scale. As a result of these stratospheric adjustments, a positive temperature perturbation developed in both hemispheres, leading to additional longwave radiation emitted back to space. According to the model results, this adjustment occurred in the stratosphere within the first 2 months after the event. At the top of the atmosphere (TOA), the net effective radiative forcing (ERF) averaged over the Southern Hemisphere was initially dominated by the instantaneous positive radiative forcing of about +0.5ĝ€¯Wm-2, for which the positive sign resulted mainly from the presence of clouds above the Southern Ocean. The longwave adjustments led to a compensation of the initially net positive TOA ERF, which is seen in the Southern Hemisphere, the tropics, and the northern mid-latitudes. The simulated changes in the lower stratosphere also affected the upper troposphere through a thermodynamic downward coupling. Subsequently, increased temperatures were also obtained in the upper troposphere, causing a global decrease in relative humidity, cirrus amount, and the ice water path of about 0.2ĝ€¯%. As a result, surface precipitation also decreased by a similar amount, which was accompanied by a weakening of the tropospheric circulation due to the given energetic constraints. In general, it appears that the radiative effects of smoke from single extreme wildfire events can lead to global impacts that affect the interplay of tropospheric and stratospheric budgets in complex ways. This emphasizes that future changes in extreme wildfires need to be included in projections of aerosol radiative forcing.
- ItemHarmonization of global land use change and management for the period 850–2100 (LUH2) for CMIP6(Katlenburg-Lindau : Copernicus, 2020) Hurtt, George C.; Chini, Louise; Sahajpal, Ritvik; Frolking, Steve; Bodirsky, Benjamin L.; Calvin, Katherine; Doelman, Jonathan C.; Fisk, Justin; Fujimori, Shinichiro; Klein Goldewijk, Kees; Hasegawa, Tomoko; Havlik, Peter; Heinimann, Andreas; Humpenöder, Florian; Jungclaus, Johan; Kaplan, Jed O.; Kennedy, Jennifer; Krisztin, Tamás; Lawrence, David; Lawrence, Peter; Ma, Lei; Mertz, Ole; Pongratz, Julia; Popp, Alexander; Poulter, Benjamin; Riahi, Keywan; Shevliakova, Elena; Stehfest, Elke; Thornton, Peter; Tubiello, Francesco N.; van Vuuren, Detlef P.; Zhang, XinHuman land use activities have resulted in large changes to the biogeochemical and biophysical properties of the Earth's surface, with consequences for climate and other ecosystem services. In the future, land use activities are likely to expand and/or intensify further to meet growing demands for food, fiber, and energy. As part of the World Climate Research Program Coupled Model Intercomparison Project (CMIP6), the international community has developed the next generation of advanced Earth system models (ESMs) to estimate the combined effects of human activities (e.g., land use and fossil fuel emissions) on the carbon–climate system. A new set of historical data based on the History of the Global Environment database (HYDE), and multiple alternative scenarios of the future (2015–2100) from Integrated Assessment Model (IAM) teams, is required as input for these models. With most ESM simulations for CMIP6 now completed, it is important to document the land use patterns used by those simulations. Here we present results from the Land-Use Harmonization 2 (LUH2) project, which smoothly connects updated historical reconstructions of land use with eight new future projections in the format required for ESMs. The harmonization strategy estimates the fractional land use patterns, underlying land use transitions, key agricultural management information, and resulting secondary lands annually, while minimizing the differences between the end of the historical reconstruction and IAM initial conditions and preserving changes depicted by the IAMs in the future. The new approach builds on a similar effort from CMIP5 and is now provided at higher resolution (0.25∘×0.25∘) over a longer time domain (850–2100, with extensions to 2300) with more detail (including multiple crop and pasture types and associated management practices) using more input datasets (including Landsat remote sensing data) and updated algorithms (wood harvest and shifting cultivation); it is assessed via a new diagnostic package. The new LUH2 products contain > 50 times the information content of the datasets used in CMIP5 and are designed to enable new and improved estimates of the combined effects of land use on the global carbon–climate system.
- ItemModelling mineral dust emissions and atmospheric dispersion with MADE3 in EMAC v2.54(Katlenburg-Lindau : Copernicus, 2020) Beer, Christof G.; Hendricks, Johannes; Righi, Mattia; Heinold, Bernd; Tegen, Ina; Groß, Silke; Sauer, Daniel; Walser, Adrian; Weinzierl, BernadettIt was hypothesized that using mineral dust emission climatologies in global chemistry climate models (GCCMs), i.e. prescribed monthly-mean dust emissions representative of a specific year, may lead to misrepresentations of strong dust burst events. This could result in a negative bias of model dust concentrations compared to observations for these episodes. Here, we apply the aerosol microphysics submodel MADE3 (Modal Aerosol Dynamics model for Europe, adapted for global applications, third generation) as part of the ECHAM/MESSy Atmospheric Chemistry (EMAC) general circulation model. We employ two different representations of mineral dust emissions for our model simulations: (i) a prescribed monthly-mean climatology of dust emissions representative of the year 2000 and (ii) an online dust parametrization which calculates wind-driven mineral dust emissions at every model time step. We evaluate model results for these two dust representations by comparison with observations of aerosol optical depth from ground-based station data. The model results show a better agreement with the observations for strong dust burst events when using the online dust representation compared to the prescribed dust emissions setup. Furthermore, we analyse the effect of increasing the vertical and horizontal model resolution on the mineral dust properties in our model. We compare results from simulations with T42L31 and T63L31 model resolution (2.8∘×2.8∘ and 1.9∘×1.9∘ in latitude and longitude, respectively; 31 vertical levels) with the reference setup (T42L19). The different model versions are evaluated against airborne in situ measurements performed during the SALTRACE mineral dust campaign (Saharan Aerosol Long-range Transport and Aerosol-Cloud Interaction Experiment, June–July 2013), i.e. observations of dust transported from the Sahara to the Caribbean. Results show that an increased horizontal and vertical model resolution is able to better represent the spatial distribution of airborne mineral dust, especially in the upper troposphere (above 400 hPa). Additionally, we analyse the effect of varying assumptions for the size distribution of emitted dust but find only a weak sensitivity concerning these changes. The results of this study will help to identify the model setup best suited for future studies and to further improve the representation of mineral dust particles in EMAC-MADE3.
- ItemReduced Complexity Model Intercomparison Project Phase 1: introduction and evaluation of global-mean temperature response(Katlenburg-Lindau : Copernicus, 2020) Nicholls, Zebedee R. J.; Meinshausen, Malte; Lewis, Jared; Gieseke, Robert; Dommenget, Dietmar; Dorheim, Kalyn; Fan, Chen-Shuo; Fuglestvedt, Jan S.; Gasser, Thomas; Golüke, Ulrich; Goodwin, Philip; Hartin, Corinne; Hope, Austin P.; Kriegler, Elmar; Leach, Nicholas J.; Marchegiani, Davide; McBride, Laura A.; Quilcaille, Yann; Rogelj, Joeri; Salawitch, Ross J.; Samset, Bjørn H.; Sandstad, Marit; Shiklomanov, Alexey N.; Skeie, Ragnhild B.; Smith, Christopher J.; Smith, Steve; Tanaka, Katsumasa; Tsutsui, Junichi; Xie, ZhiangReduced-complexity climate models (RCMs) are critical in the policy and decision making space, and are directly used within multiple Intergovernmental Panel on Climate Change (IPCC) reports to complement the results of more comprehensive Earth system models. To date, evaluation of RCMs has been limited to a few independent studies. Here we introduce a systematic evaluation of RCMs in the form of the Reduced Complexity Model Intercomparison Project (RCMIP). We expect RCMIP will extend over multiple phases, with Phase 1 being the first. In Phase 1, we focus on the RCMs' global-mean temperature responses, comparing them to observations, exploring the extent to which they emulate more complex models and considering how the relationship between temperature and cumulative emissions of CO2 varies across the RCMs. Our work uses experiments which mirror those found in the Coupled Model Intercomparison Project (CMIP), which focuses on complex Earth system and atmosphere–ocean general circulation models. Using both scenario-based and idealised experiments, we examine RCMs' global-mean temperature response under a range of forcings. We find that the RCMs can all reproduce the approximately 1 ∘C of warming since pre-industrial times, with varying representations of natural variability, volcanic eruptions and aerosols. We also find that RCMs can emulate the global-mean temperature response of CMIP models to within a root-mean-square error of 0.2 ∘C over a range of experiments. Furthermore, we find that, for the Representative Concentration Pathway (RCP) and Shared Socioeconomic Pathway (SSP)-based scenario pairs that share the same IPCC Fifth Assessment Report (AR5)-consistent stratospheric-adjusted radiative forcing, the RCMs indicate higher effective radiative forcings for the SSP-based scenarios and correspondingly higher temperatures when run with the same climate settings. In our idealised setup of RCMs with a climate sensitivity of 3 ∘C, the difference for the ssp585–rcp85 pair by 2100 is around 0.23∘C(±0.12 ∘C) due to a difference in effective radiative forcings between the two scenarios. Phase 1 demonstrates the utility of RCMIP's open-source infrastructure, paving the way for further phases of RCMIP to build on the research presented here and deepen our understanding of RCMs.
- ItemCoupling aerosols to (cirrus) clouds in the global EMAC-MADE3 aerosol–climate model(Katlenburg-Lindau : Copernicus, 2020) Righi, Mattia; Hendricks, Johannes; Lohmann, Ulrike; Beer, Christof Gerhard; Hahn, Valerian; Heinold, Bernd; Heller, Romy; Krämer, Martina; Ponater, Michael; Rolf, Christian; Tegen, Ina; Voigt, ChristianeA new cloud microphysical scheme including a detailed parameterization for aerosol-driven ice formation in cirrus clouds is implemented in the global ECHAM/MESSy Atmospheric Chemistry (EMAC) chemistry–climate model and coupled to the third generation of the Modal Aerosol Dynamics model for Europe adapted for global applications (MADE3) aerosol submodel. The new scheme is able to consistently simulate three regimes of stratiform clouds – liquid, mixed-, and ice-phase (cirrus) clouds – considering the activation of aerosol particles to form cloud droplets and the nucleation of ice crystals. In the cirrus regime, it allows for the competition between homogeneous and heterogeneous freezing for the available supersaturated water vapor, taking into account different types of ice-nucleating particles, whose specific ice-nucleating properties can be flexibly varied in the model setup. The new model configuration is tuned to find the optimal set of parameters that minimizes the model deviations with respect to observations. A detailed evaluation is also performed comparing the model results for standard cloud and radiation variables with a comprehensive set of observations from satellite retrievals and in situ measurements. The performance of EMAC-MADE3 in this new coupled configuration is in line with similar global coupled models and with other global aerosol models featuring ice cloud parameterizations. Some remaining discrepancies, namely a high positive bias in liquid water path in the Northern Hemisphere and overestimated (underestimated) cloud droplet number concentrations over the tropical oceans (in the extratropical regions), which are both a common problem in these kinds of models, need to be taken into account in future applications of the model. To further demonstrate the readiness of the new model system for application studies, an estimate of the anthropogenic aerosol effective radiative forcing (ERF) is provided, showing that EMAC-MADE3 simulates a relatively strong aerosol-induced cooling but within the range reported in the Intergovernmental Panel on Climate Change (IPCC) assessments.
- ItemInvestigating Mesozoic Climate Trends and Sensitivities With a Large Ensemble of Climate Model Simulations(Hoboken, NJ : Wiley, 2021) Landwehrs, Jan; Feulner, Georg; Petri, Stefan; Sames, Benjamin; Wagreich, MichaelThe Mesozoic era (∼252 to 66 million years ago) was a key interval in Earth's evolution toward its modern state, witnessing the breakup of the supercontinent Pangaea and significant biotic innovations like the early evolution of mammals. Plate tectonic dynamics drove a fundamental climatic transition from the early Mesozoic supercontinent toward the Late Cretaceous fragmented continental configuration. Here, key aspects of Mesozoic long-term environmental changes are assessed in a climate model ensemble framework. We analyze so far the most extended ensemble of equilibrium climate states simulated for evolving Mesozoic boundary conditions covering the period from 255 to 60 Ma in 5 Myr timesteps. Global mean temperatures are generally found to be elevated above the present and exhibit a baseline warming trend driven by rising sea levels and increasing solar luminosity. Warm (Triassic and mid-Cretaceous) and cool (Jurassic and end-Cretaceous) anomalies result from pCO2 changes indicated by different reconstructions. Seasonal and zonal temperature contrasts as well as continental aridity show an overall decrease from the Late Triassic-Early Jurassic to the Late Cretaceous. Meridional temperature gradients are reduced at higher global temperatures and less land area in the high latitudes. With systematic sensitivity experiments, the influence of paleogeography, sea level, vegetation patterns, pCO2, solar luminosity, and orbital configuration on these trends is investigated. For example, long-term seasonality trends are driven by paleogeography, but orbital cycles could have had similar-scale effects on shorter timescales. Global mean temperatures, continental humidity, and meridional temperature gradients are, however, also strongly affected by pCO2.