J. F. González-Rouco

Simulation and inversion of borehole temperature profiles in surrogate climates: Spatial distribution and surface coupling

A heat-conduction forward model driven by ground surface temperature from three 1000-year climate simulations with the state-of-the-art ECHO-g model has been used to simulate underground temperature perturbation profiles. An inversion approach has been applied to reconstruct ground surface temperature histories from the simulated profiles and to compare them with the climate model temperatures. Results support the skill of borehole inversion methods to retrieve long-term temperature trends, and the robustness of using the present-day borehole network for reconstructing SAT variations.

Natural and anthropogenic modes of surface temperature variations in the last thousand years

The spatial patterns of surface air-temperature variations in the period 1000 to 2100, simulated with the ECHO-G atmosphere-ocean coupled model, are analyzed. The model was driven by solar, volcanic and greenhouse gas forcing. The leading mode of temperature variability in the preindustrial period represents an almost global coherent variation of temperatures, with larger amplitudes over the continents and Northern Hemisphere. This mode also describes a large part of the spatial structure of the warming simulated in the 21st century. However, in the 21st century, regional departures from this spatial structure are also present and can be ascribed to atmospheric circulation responses to anthropogenic forcing in the last decades of the 21st century.

Agreement between Observed Rainfall Trends and Climate Change Simulations in the Southwest of Europe

Abstract The lowest spatial scale at which current climate models are considered to be skillful is on the order of 1000 km because of resolution and computer capabilities. The estimation of the regional changes caused by anthropogenic emissions of greenhouse gases and aerosols therefore is problematic. Here a statistical downscaling scheme is used to study the relationship between large-scale sea level pressure and regional precipitation in southwestern Europe, both in observed data and in outputs from a general circulation model (GCM) forced with increasing levels of greenhouse gases and sulfate aerosols. The results indicate that the GCM does reproduce the main aspects of the large- to local-scale coupled variability. Furthermore, these large- to local-scale relationships remain stable in the scenario simulations. The GCM runs predict increases of advection of oceanic air masses to the Iberian Peninsula that will produce a slight decrease of precipitation amounts in the north coast and the opposite effe...

Sensitivity of the MM5 mesoscale model to physical parameterizations for regional climate studies: Annual cycle

We present an analysis of the sensitivity to different physical parameterizations of a high-resolution simulation of the MM5 mesoscale model over the Iberian Peninsula. Several (16) 5-year runs of the MM5 model with varying parameterizations of microphysics, cumulus, planetary boundary layer and longwave radiation have been carried out. The results have been extensively compared with observational precipitation and surface temperature data. The parameterization uncertainty has also been compared with that related to the boundary conditions and the varying observational data sets. The annual cycles of precipitation and surface temperature are well reproduced. The summer season presents the largest deviations, with a 5 K cold bias in the southeast and noticeable precipitation errors over mountain areas. The cold bias seems to be related to the surface, probably because of the excessive moisture availability of the five-layer soil scheme used. No parameterization combination was found to perform best in simulating both precipitation and surface temperature in every season and subregion. The Kain-Fritsch cumulus scheme was found to produce unrealistically high summer precipitation. The longwave radiation parameterizations tested were found to have little impact on our target variables. Other factors, such as the choice of boundary conditions, have an impact on the results as large as the selection of parameterizations. The range of variability in the MM5 physics ensemble is of the same order of magnitude as the observational uncertainty, except in summer, when it is larger and probably related to the inaccuracy of the model to reproduce the summer precipitation over the area.

Spatial performance of four climate field reconstruction methods targeting the Common Era

The spatial skill of four climate field reconstruction (CFR) methods is investigated using pseudoproxy experiments (PPEs) based on two millennial-length general circulation model simulations. Results indicate that presently available global and hemispheric CFRs for the Common Era likely suffer from spatial uncertainties not previously characterized. No individual method produced CFRs with universally superior spatial error statistics, making it difficult to advocate for one method over another. Northern Hemisphere means are shown to be insufficient for evaluating spatial skill, indicating that the spatial performance of future CFRs should be rigorously tested for dependence on proxy type and location, target data and employed methodologies. Observed model-dependent methodological performance also indicates that CFR methods must be tested across multiple models and conclusions from PPEs should be carefully evaluated against the spatial statistics of real-world climatic fields. Citation: Smerdon, J. E., A. Kaplan, E. Zorita, J. F. Gonzalez-Rouco, and M. N. Evans (2011), Spatial performance of four climate field reconstruction methods targeting the Common Era, Geophys. Res. Lett., 38, L11705, doi: 10.1029/2011GL047372.

Megadroughts in Southwestern North America in ECHO-G Millennial Simulations and Their Comparison to Proxy Drought Reconstructions*

AbstractSimulated hydroclimate variability in millennium-length forced transient and control simulations from the ECHAM and the global Hamburg Ocean Primitive Equation (ECHO-G) coupled atmosphere–ocean general circulation model (AOGCM) is analyzed and compared to 1000 years of reconstructed Palmer drought severity index (PDSI) variability from the North American Drought Atlas (NADA). The ability of the model to simulate megadroughts in the North American southwest is evaluated. (NASW: 25°–42.5°N, 125°–105°W). Megadroughts in the ECHO-G AOGCM are found to be similar in duration and magnitude to those estimated from the NADA. The droughts in the forced simulation are not, however, temporally synchronous with those in the paleoclimate record, nor are there significant differences between the drought features simulated in the forced and control runs. These results indicate that model-simulated megadroughts can result from internal variability of the modeled climate system rather than as a response to changes ...

The role of the land-surface model for climate change projections over the Iberian Peninsula

The importance of land-surface processes within Regional Climate Models for accurately reproducing the present-day climate is well known. However, their role when projecting future climate is still poorly reported. Hence, this work assesses the influence of the land-surface processes, particularly the contribution of soil moisture, when projecting future changes for temperature, precipitation and wind over a complex area as the Iberian Peninsula, which, in addition, shows great sensitivity to climate change. The main signals are found for the summer season, when the results indicate a strengthening in the increases projected for both mean temperature and temperature variability as a consequence of the future intensification of the positive soil moisture-temperature feedback. The more severe warming over the inner dry Iberian Peninsula further implies an intensification of the Iberian thermal low and, thus, of the cyclonic circulation. Furthermore, the land-atmosphere coupling leads to the projection of a wider future daily temperature range, since maximum temperatures are more affected than minima, a feature absent in non-coupled simulations. Regarding variability, the areas where the land-atmosphere coupling introduces larger changes are those where the reduction in the soil moisture content is more dramatic in future simulations, i.e., the so-called transitional zones. As regards precipitation, weaker positive signals for convective precipitation and more intense negative signals for non-convective precipitation are obtained as a result of the soil moisture-atmosphere interactions. These results highlight the crucial contribution of soil moisture to climate change projections and suggest its plausible key role for future projections of extreme events.

A Pacific Centennial Oscillation Predicted by Coupled GCMs

Internal climate variability at the centennial time scale is investigated using long control integrations from three state-of-the-art global coupled general circulation models. In the absence of external forcing, all three models produce centennial variability in the mean zonal sea surface temperature (SST) and sea level pressure (SLP) gradients in the equatorial Pacific with counterparts in the extratropics. The centennial pattern in the tropical Pacific is dissimilar to that of the interannual El Nino-Southern Oscillation (ENSO), in that the most prominent expression in temperature is found beneath the surface of the western Pacific warm pool. Some global repercussions nevertheless are analogous, such as a hemispherically symmetric atmospheric wave pattern of alternating highs and lows. Centennial variability in western equatorial Pacific SST is a result of the strong asymmetry of interannual ocean heat content anomalies, while the eastern equatorial Pacific exhibits a lagged, Bjerknes-like response to temperature and convection in the west. The extratropical counterpart is shown to be a flux-driven response to the hemispherically symmetric circulation anomalies emanating from the tropical Pacific. Significant centennial-length trends in the zonal SST and SLP gradients rivaling those estimated from observations and model simulations forced with increasing CO_2 appear to be inherent features of the internal climate dynamics simulated by all three models. Unforced variability and trends on the centennial time scale therefore need to be addressed in estimated uncertainties, beyond more traditional signal-to-noise estimates that do not account for natural variability on the centennial time scale.

Temperature sensitivity to the land-surface model in MM5 climate simulations over the Iberian Peninsula

Three different Land Surface Models have been used in three high resolution climate simulations performed with the mesoscale model MM5 over the Iberian Peninsula. The main difference among them lies in the soil moisture treatment, which is dynamically modelled by only two of them (Noah and Pleim & Xiu models), while in the simplest model (Simple Five-Layers) it is fixed to climatological values. The simulated period covers 1958-2002, using the ERA40 reanalysis data as driving conditions. Focusing on near-surface air temperature, this work evaluates the skill of each simulation in reproducing mean values and temporal variability, by comparing the simulations with observed temperature series. When the simplest simulation was analyzed, the greatest discrepances were observed for the summer season, when both, the mean values and the temporal variability of the temperature series, were badly underestimated. These weaknesses are largely overcome in the other two simulations (performed by coupling a more advanced soil model to MM5), and there was greater concordance between the simulated and observed spatial patterns. The influence of a dynamic soil moisture parameterization and, therefore, a more realistic simulation of the latent and sensible heat fluxes between the land and the atmosphere, helps to explain these results.

North Atlantic atmospheric circulation and surface wind in the Northeast of the Iberian Peninsula: uncertainty and long term downscaled variability

The variability and predictability of the surface wind field at the regional scale is explored over a complex terrain region in the northeastern Iberian Peninsula by means of a downscaling technique based on Canonical Correlation Analysis. More than a decade of observations (1992–2005) allows for calibrating and validating a statistical method that elicits the main associations between the large scale atmospheric circulation over the North Atlantic and Mediterranean areas and the regional wind field. In an initial step the downscaling model is designed by selecting parameter values from practise. To a large extent, the variability of the wind at monthly timescales is found to be governed by the large scale circulation modulated by the particular orographic features of the area. The sensitivity of the downscaling methodology to the selection of the model parameter values is explored, in a second step, by performing a systematic sampling of the parameters space, avoiding a heuristic selection. This provides a metric for the uncertainty associated with the various possible model configurations. The uncertainties associated with the model configuration are considerably dependent on the spatial variability of the wind. While the sampling of the parameters space in the model set up moderately impact estimations during the calibration period, the regional wind variability is very sensitive to the parameters selection at longer timescales. This fact illustrates that downscaling exercises based on a single configuration of parameters should be interpreted with extreme caution. The downscaling model is used to extend the estimations several centuries to the past using long datasets of sea level pressure, thereby illustrating the large temporal variability of the regional wind field from interannual to multicentennial timescales. The analysis does not evidence long term trends throughout the twentieth century, however anomalous episodes of high/low wind speeds are identified.