Exploiting spectrally resolved satellite observations to assess the performance of climate models

The accurate representation of the Earth Radiation Budget by General Circulation Models (GCMs) is a fundamental requirement to provide reliable historical and future climate simulations. In this study, we found reasonable agreement between the integrated energy fluxes at the top of the atmosphere simulated by 34 state-of-the-art climate models and the observations provided by the Cloud and Earth Radiant Energy System (CERES) mission on a global scale, but large regional biases have been detected throughout the globe. Furthermore, we highlighted that a good agreement between simulated and observed integrated Outgoing Longwave Radiation (OLR) fluxes may be obtained from the cancellation of opposite-in-sign systematic errors, localized in different spectral ranges. To avoid this and to understand the causes of these biases, we compared the observed Earth emission spectra, measured by the Infrared Atmospheric Sounding Interferometer (IASI) in the period 2008-2016, with the synthetic radiances computed on the basis of the atmospheric fields provided by the EC-Earth GCM. To this purpose, the fast σ-IASI radiative transfer model was used, after its validation and implementation in EC-Earth. From the comparison between observed and simulated spectral radiances, a positive temperature bias in the stratosphere and a negative temperature bias in the middle troposphere, as well as a dry bias of the water vapor concentration in the upper troposphere, have been identified in the EC-Earth climate model. The analysis has been performed in clear-sky conditions, but the feasibility of its extension in the presence of clouds, whose impact on the radiation represents the greatest source of uncertainty in climate models, has also been proven. Finally, the analysis of simulated and observed OLR trends indicated good agreement and provided detailed information on the spectral fingerprints of the evolution of the main climate variables.

Evapotranspiration: Present and Future Challenges

The thesis first explored and evaluated some of the most used models that were developed to account for the effect of CO2 on evapotranspiration. This review depicts the complexity of the modeling procedure and underlines the advantages and shortcomings of each model. Then, the projected climate change in the near future (2021-2050) in different locations in Emilia-Romagna (Italy) was studied, with an emphasis on the opposite effect of an increase in both air temperature and CO2 levels on ETo. The case study used reanalysis data as a surrogate to historical weather stations measurements and an ensemble of regional climate models (RCMs) for the future projections. Results show that higher CO2 levels moderated the increase in ETo that accompanies an increase in air temperature, taking in consideration the change in other weather variables i.e. solar radiation, wind speed and dew point temperature. The outcomes of this study show that considering the CO2 fertilization effect when calculating reference evapotranspiration might give a more realistic estimation of water use efficiency and irrigation requirements in Emilia-Romagna and a better analysis of the future availability and distribution of water resources in the region. Finally, data from a model forecasting reference evapotranspiration (FRET) and the different variables involved in its calculation for the state of California (USA) were compared with similar data from the regional weather station network (CIMIS) to evaluate their accuracy and reliability. The evaluation was done in locations with different microclimates and included also sample irrigation schedules developed using FRET ETo. The obtained results demonstrate that FRET ETo forecasts are a viable alternative to traditional ETo measurements with some differences depending on the climatic condition of the location considered in this study. This implies that FRET could be replicated in other areas with similar climate settings.

The role of the Rocky Mountains in shaping the atmospheric mean state and the response to tropical forcing in idealised simulations.

Understanding the natural and forced variability of the atmospheric general circulation and its drivers is one of the grand challenges in climate science. It is of paramount importance to understand to what extent the systematic error of climate models affects the processes driving such variability. This is done by performing a set of simulations (ROCK experiments) with an intermediate complexity atmospheric model (SPEEDY), in which the Rocky Mountains orography is increased or decreased to influence the structure of the North Pacific jet stream. For each of these modified-orography experiments, the climatic response to idealized sea surface temperature anomalies of varying intensity in the El Niño Southern Oscillation (ENSO) region is studied. ROCK experiments are characterized by variations in the Pacific jet stream intensity whose extension encompasses the spread of the systematic error found in Coupled Model Intercomparison Project (CMIP6) models. When forced with ENSO-like idealised anomalies, they exhibit a non-negligible sensitivity in the response pattern over the Pacific North American region, indicating that the model mean state can affect the model response to ENSO. It is found that the classical Rossby wave train response to ENSO is more meridionally oriented when the Pacific jet stream is weaker and more zonally oriented with a stronger jet. Rossby wave linear theory suggests that a stronger jet implies a stronger waveguide, which traps Rossby waves at a lower latitude, favouring a zonal propagation of Rossby waves. The shape of the dynamical response to ENSO affects the ENSO impacts on surface temperature and precipitation over Central and North America. A comparison of the SPEEDY results with CMIP6 models suggests a wider applicability of the results to more resources-demanding climate general circulation models (GCMs), opening up to future works focusing on the relationship between Pacific jet misrepresentation and response to external forcing in fully-fledged GCMs.

Advanced representation of the ocean/sea ice dynamics at high latitudes

Sea ice is a fundamental element of global climate system, with numerous impacts on the polar environment. The ongoing drastic changes in the Earth’s sea ice cover highlight the necessity of monitoring the polar regions and systematically evaluating the quality of different numerical products. The main objective of this thesis is to improve our knowledge of the representation of Arctic and Antarctic sea ice using comprehensive global ocean reanalyses and coupled climate models. The dissertation will explore (i) the Antarctic marginal ice zone (MIZ) and pack ice area in the ensemble mean of four global ocean reanalyses called GREP; (ii) historical representation of the Arctic and Antarctic sea ice state in HighResMIP models; (iii) the future evolution of Arctic sea ice in HighResMIP models. Global ocean reanalyses and GREP are found to adequately capture interannual and seasonal variability in both pack ice and MIZ areas at hemispheric and regional scales. The advantage of the ensemble-mean approach is proved as GREP smooths the strengths and weaknesses of single systems and provides the most consistent and reliable estimates. This work is intended to encourage the use of GREP in a wide range of applications. The analysis of sea ice representation in the coupled climate models shows no systematic impact of the increased horizontal resolution. We argue that a few minor improvements in sea ice representation with the enhanced horizontal resolution are presumably not worth the major effort of costly computations. The thesis highlights the critical importance to distinguish the MIZ from consolidated pack ice both for investigating changes in sea ice distribution and evaluating the product’s performance. Considering that the MIZ is predicted to dominate the Arctic sea ice cover, the model physics parameterizations and sea ice rheology might require modifications. The results of the work can be useful for modelling community.

Creep and damage models for the service behaviour of structural members

Deformability is often a crucial to the conception of many civil-engineering structural elements. Also, design is all the more burdensome if both long- and short-term deformability has to be considered. In this thesis, long- and short-term deformability has been studied from the material and the structural modelling point of view. Moreover, two materials have been handled: pultruded composites and concrete. A new finite element model for thin-walled beams has been introduced. As a main assumption, cross-sections rigid are considered rigid in their plane; this hypothesis replaces that of the classical beam theory of plane cross-sections in the deformed state. That also allows reducing the total number of degrees of freedom, and therefore making analysis faster compared with twodimensional finite elements. Longitudinal direction warping is left free, allowing describing phenomena such as the shear lag. The new finite-element model has been first applied to concrete thin-walled beams (such as roof high span girders or bridge girders) subject to instantaneous service loadings. Concrete in his cracked state has been considered through a smeared crack model for beams under bending. At a second stage, the FE-model has been extended to the viscoelastic field and applied to pultruded composite beams under sustained loadings. The generalized Maxwell model has been adopted. As far as materials are concerned, long-term creep tests have been carried out on pultruded specimens. Both tension and shear tests have been executed. Some specimen has been strengthened with carbon fibre plies to reduce short- and long- term deformability. Tests have been done in a climate room and specimens kept 2 years under constant load in time. As for concrete, a model for tertiary creep has been proposed. The basic idea is to couple the UMLV linear creep model with a damage model in order to describe nonlinearity. An effective strain tensor, weighting the total and the elasto-damaged strain tensors, controls damage evolution through the damage loading function. Creep strains are related to the effective stresses (defined by damage models) and so associated to the intact material.

Atmospheric radiative transfer in multiple scattering conditions. Application to NWP models.

High spectral resolution radiative transfer (RT) codes are essential tools in the study of the radiative energy transfer in the Earth atmosphere and a support for the development of parameterizations for fast RT codes used in climate and weather prediction models. Cirrus clouds cover permanently 30% of the Earth's surface, representing an important contribution to the Earth-atmosphere radiation balance. The work has been focussed on the development of the RT model LBLMS. The model, widely tested in the infra-red spectral range, has been extended to the short wave spectrum and it has been used in comparison with airborne and satellite measurements to study the optical properties of cirrus clouds. A new database of single scattering properties has been developed for mid latitude cirrus clouds. Ice clouds are treated as a mixture of ice crystals with various habits. The optical properties of the mixture are tested in comparison to radiometric measurements in selected case studies. Finally, a parameterization of the mixture for application to weather prediction and global circulation models has been developed. The bulk optical properties of ice crystals are parameterized as functions of the effective dimension of measured particle size distributions that are representative of mid latitude cirrus clouds. Tests with the Limited Area Weather Prediction model COSMO have shown the impact of the new parameterization with respect to cirrus cloud optical properties based on ice spheres.

Climate impact studies of sediment transport on the Adriatic coastal zone

The study of the impact of climate change on the environment has been based, until very recently, on an global approach, whose interest from a local point of view is very limited. This thesis, on the contrary, has treated the study of the impact of climate change in the Adriatic Sea basin following a twofold strategy of regionalization and integration of numerical models in order to reproduce the present and future scenarios of the system through a more and more realistic and solid approach. In particular the focus of the study was on the impact on the physical environment and on the sediment transport in the basin. This latter is a very new and original issue, to our knowledge still uninvestigated. The study case of the coastal area of Montenegro was particularly studied, since it is characterized by an important supply of sediment through the Buna/Bojana river, second most important in the Adriatic basin in terms of flow. To do this, a methodology to introduce the tidal processes in a baroclinic primitive equations Ocean General Circulation Model was applied and tidal processes were successfully reproduced in the Adriatic Sea, analyzing also the impacts they have on the mean general circulation, on salt and heat transport and on mixing and stratification of the water column in the different seasons of the year. The new hydrodynamical model has been further coupled with a wave model and with a river and sea sediment transport model, showing good results in the reproduction of sediment transport processes. Finally this complex coupled platform was integrated in the period 2001-2030 under the A1B scenario of IPCC, and the impact of climate change on the physical system and on sediment transport was preliminarily evaluated.

Sea-Level climate variability in the Mediterranean Sea

Sea-level variability is characterized by multiple interacting factors described in the Fourth Assessment Report (Bindoff et al., 2007) of the Intergovernmental Panel on Climate Change (IPCC) that act over wide spectra of temporal and spatial scales. In Church et al. (2010) sea-level variability and changes are defined as manifestations of climate variability and change. The European Environmental Agency (EEA) defines sea level as one of most important indicators for monitoring climate change, as it integrates the response of different components of the Earths system and is also affected by anthropogenic contributions (EEA, 2011). The balance between the different sea-level contributions represents an important source of uncertainty, involving stochastic processes that are very difficult to describe and understand in detail, to the point that they are defined as an enigma in Munk (2002). Sea-level rate estimates are affected by all these uncertainties, in particular if we look at possible responses to sea-level contributions to future climate. At the regional scale, lateral fluxes also contribute to sea-level variability, adding complexity to sea-level dynamics. The research strategy adopted in this work to approach such an interesting and challenging topic has been to develop an objective methodology to study sea-level variability at different temporal and spatial scales, applicable in each part of the Mediterranean basin in particular, and in the global ocean in general, using all the best calibrated sources of data (for the Mediterranean): in-situ, remote-sensig and numerical models data. The global objective of this work was to achieve a deep understanding of all of the components of the sea-level signal contributing to sea-level variability, tendency and trend and to quantify them.