High-resolution sea wind hindcasts over the Mediterranean area

The goal of this study is to develop a high-resolution atmospheric hindcast over the Mediterranean area using the WRF-ARW model, focusing on offshore surface wind fields. In order to choose the most adequate model configuration, the study provides details on the calibration of the experimental saet-up through a sensitivity test considering the October–December 2001 period (the 2001 super-storm event in the West Mediterranean). A daily forecast outperforms the spectral technique of previous products and the boundary data from ERA-Interim reanalysis produces the most accurate estimates in terms of wind variability and hour-to-hour correspondence. According to the sensitivity test, two data sets of wind hindcast are produced: the SeaWind I (30-km horizontal resolution for a period of 60 years) and the SeaWind II (15-km horizontal resolution for 20 years). The validation of the resulting surface winds is undertaken considering two offshore observational datasets. On the one hand, hourly surface buoy stations are used to validate wind time series at specific locations; on the other hand, wind altimeter satellite observations are considered for spatial validation in the whole Mediterranean Sea. The results obtained from this validation process show a very good agreement with observations for the southern Europe region. Finally, SeaWind I and II are used to characterize offshore wind fields in the Mediterranean Sea. The statistical structure of sea surface wind is analyzed and the agreement with Weibull probability distribution is discussed. In addition, wind persistence and extreme wind speed (50 year return period) are characterized and relevant areas of wind power generation are described by estimating wind energy quantities.     

Simulation of Extreme Surges in the Taganrog Bay with Atmosphere and Ocean Circulation Models

The simulation the most extreme surges over the period of instrumental observations in the Taganrog Bay since 1881, the surges occurred on March 24, 2013 and September 24, 2014. The objective of the simulation is to study surge formation features and to reveal requirements for the accuracy of simulating atmospheric and oceanic circulation in the Sea of Azov. For this purpose, the Institute of Numerical Mathematics Ocean Model (INMOM) with the spatial resolution of ~4 km and ~250 m was used. The atmospheric forcing over the Black Sea region was specified using ERA-Interim reanalysis data and WRF model data with the spatial resolution of 80 and 10 km, respectively. It is shown that the quality of simulation of extreme surges in the Sea of Azov is more dependent on the quality of the input atmospheric forcing than on the spatial resolution of the ocean circulation model. The usage of WRF data as atmospheric forcing allows the more accurate simulation of extreme surges. However, the simulation of the extreme surge of 2014 overestimates, and simulations for the 2013 surge underestimate the surge level. Evidently, as the used version of INMOM does not take into account the coastal zone flooding, the maximum surge value is overestimated.     

Projected future changes in tropical cyclone-related wave climate in the North Atlantic

Tropical cyclones are a major hazard for numerous countries surrounding the tropical-to-subtropical North Atlantic sub-basin including the Caribbean Sea and Gulf of Mexico. Their intense winds, which can exceed 300 km h⁻¹, can cause serious damage, particularly along coastlines where the combined action of waves, currents and low atmospheric pressure leads to storm surge and coastal flooding. This work presents future projections of North Atlantic tropical cyclone-related wave climate. A new configuration of the ARPEGE-Climat global atmospheric model on a stretched grid reaching ~ 14 km resolution to the north-east of the eastern Caribbean is able to reproduce the distribution of tropical cyclone winds, including Category 5 hurricanes. Historical (1984–2013, 5 members) and future (2051–2080, 5 members) simulations with the IPCC RCP8.5 scenario are used to drive the MFWAM (Météo-France Wave Action Model) spectral wave model over the Atlantic basin during the hurricane season. An intermediate 50-km resolution grid is used to propagate mid-latitude swells into a higher 10-km resolution grid over the tropical cyclone main development region. Wave model performance is evaluated over the historical period with the ERA5 reanalysis and satellite altimetry data. Future projections exhibit a modest but widespread reduction in seasonal mean wave heights in response to weakening subtropical anticyclone, yet marked increases in tropical cyclone-related wind sea and extreme wave heights within a large region extending from the African coasts to the North American continent.

     

A high-resolution 44-year atmospheric hindcast for the Mediterranean Basin: contribution to the regional improvement of global reanalysis

A 44-year (1958–2001) high-resolution atmospheric hindcast for the whole Mediterranean Basin was performed within the EU-funded Hindcast of Dynamic Processes of the Ocean and Coastal Areas of Europe (HIPOCAS) Project. The long-term hindcasted data set, which comprises several atmospheric parameters at different levels, was produced by means of dynamical downscaling from the NCEP/NCAR global reanalysis using the atmospheric limited area model REMO. The REMO hindcast has been exhaustively validated. On that score, various hindcasted surface parameters, such as 10-m wind field, 2-m temperature and mean sea level pressure, have been compared to satellite data (ERS-1/2 scatterometer) and in-situ measurements from offshore stations. In addition, two ocean models (waves and sea level) have been forced with REMO hindcasted fields (mean sea level pressure and 10-m wind field). The validation of these ocean runs, performed through comparisons of simulated waves and sea level with oceanographic measurements, allows to evaluate "indirectly" the quality of the REMO hindcasted data used as atmospheric forcing. Once the quality of the hindcasted data was verified, the efficiency of the regional enhancement performed through dynamical downscaling on the NCEP global reanalysis was assessed. The regional improvement was evaluated through comparisons of REMO and NCEP performance in reproducing observations. The important improvement obtained in the characterization of extreme wind events is particularly remarkable.     

Modelling the main features of the Algarve coastal circulation during July 2004: A downscaling approach

This article presents a methodology for simulating the Algarve coastal circulation using realistic forcing (e.g. low-frequency circulation, tide, high-resolution atmospheric forcing). Low-frequency open boundary conditions are defined via a downscaling of the HYCOM-US operational solution for the Gulf of Cadiz. Atmospheric forcing is imposed using the MM5 high-resolution mesoscale model (9 km resolution near the coast). A 3-level nesting system based on the MOHID numerical system is implemented. The higher nesting level has a horizontal resolution of 0.02° along the Algarve coast. The methodology is first validated qualitatively. A comparison between the numerical results and the conceptual model of the circulation described in the literature is presented. A quantitative validation is also performed, based mainly on remote sensing data (sea surface temperature and altimetry) available for July 2004. The numerical system is able to reproduce many of the circulation features described in the literature (e.g. Azores current recirculation in the Gulf of Cadiz, the upwelling jet, Mediterranean Water undercurrent, Mediterranean Outflow splitting, generation of meddies) and observed with remote-sensing data (e.g. the signature in sea surface temperature (SST) during a regime of upwelling relaxation).     

Urban boundary layer analysis in the complex coastal terrain of Bilbao using Enviro-HIRLAM

This study analyses the atmospheric boundary layer over the Bilbao metropolitan area during summer (13–18 Jul 2009) and winter (20–29 Jan 2010) episodes using the Environment–High Resolution Limited Area Model (Enviro-HIRLAM) coupled with the building effect parameterisation (BEP). The main objectives of this study are: to evaluate the performance of the model to simulate the land–sea breezes over this complex terrain; to assess the simulations with the integration of an urban parameterisation in Enviro-HIRLAM and finally; and to analyse the urban–atmosphere interactions. Even if the hydrostraticity of the model is a limitation to simulate atmospheric flows over complex terrain, sensibility tests demonstrate that 2.4 km is the optimal horizontal resolution over Bilbao that allows at the same time: to obtain satisfactory reproducibility of the large-scale processes and to explore the urban effects at local scale. During the summer episode, a typical regime of diurnal sea breeze from the NW-N-NE direction and nocturnal valley breezes from the SE direction are observed over Bilbao. The urban heat island (UHI) phenomenon is developed in the city centre expanding to the suburbs from 22 to 10 local time (LT), covering an area of 130 km². The maximum UHI intensity, 1 °C, is reached at the end of the night (5 LT), and it is advected 12 km towards the sea by the land breezes. The urban boundary layer (UBL) height amplitude varies from 100 (night time) to 1,360 m (at 14 LT). During the winter episode, the land breeze dominates the atmospheric diffusion during the day and night time. The maximum UHI intensity, 1.7 °C, is observed at 01 LT. It is spread and remained over the city covering an area of 160 km², with a vertical extension of 33 m. The UBL reaches 780 m height at 16 LT the following day.     

Properties of adjoint sea ice sensitivities to atmospheric forcing and implications for the causes of the long term trend of Arctic sea ice

Based on adjoint sensitivities of the coupled Massachusetts Institute of Technology ocean–sea ice circulation model, the potential influence of thermodynamic atmospheric forcing on the interannual variability of the September sea ice area (AREA) and volume (VOLUME) in the Arctic is investigated for the three periods 1980–1989, 1990–1999 and 2000–2009. Sensitivities suggest that only large forcing anomalies prior to the spring melting onset in May can influence the September sea ice characteristics while even small changes in the atmospheric variables during subsequent months can significantly influence September sea ice state. Specifically, AREA close to the ice edge in the Arctic seas is highly sensitive to thermodynamic atmospheric forcing changes from June to July. In contrast, VOLUME is highly sensitive to atmospheric temperature changes occurring during the same period over the central parts of the Arctic Ocean. A comparison of the sea ice conditions and sensitivities during three different periods reveals that, due to the strong decline of sea ice concentration and sea ice thickness, sea ice area became substantially more sensitive to the same amplitude thermodynamic atmospheric forcing anomalies during 2000–2009 relative to the earlier periods. To obtain a quantitative estimate of changes that can be expected from existing atmospheric trends, adjoint sensitivities are multiplied by monthly temperature differences between 1980s and two following decades. Strongest contributions of surface atmospheric temperature differences to AREA and VOLUME changes are observed during May and September. The strongest contribution from the downward long-wave heat flux to AREA changes occurs in September and to VOLUME changes in July–August. About 62 % of the AREA decrease simulated by the model can be explained by summing all contributions to the thermodynamic atmospheric forcing. The changing sea ice state (sensitivity) is found to enhance the decline and accounts for about one third of the explained reduction. For the VOLUME decrease, the explained fraction of the decrease is only about 37 %.     

Thirty-two-year ocean–atmosphere coupled downscaling of global reanalysis over the Intra-American Seas

This study examines the oceanic and atmospheric variability over the Intra-American Seas (IAS) from a 32-year integration of a 15-km coupled regional climate model consisting of the Regional Spectral Model (RSM) for the atmosphere and the Regional Ocean Modeling System (ROMS) for the ocean. It is forced at the lateral boundaries by National Centers for Environmental Prediction-Department of Energy (NCEP-DOE R-2) atmospheric global reanalysis and Simplified Ocean Data Assimilation global oceanic reanalysis. This coupled downscaling integration is a free run without any heat flux correction and is referred as the Regional Ocean–Atmosphere coupled downscaling of global Reanalysis over the Intra-American Seas (ROARS). The paper examines the fidelity of ROARS with respect to independent observations that are both satellite based and in situ. In order to provide a perspective on the fidelity of the ROARS simulation, we also compare it with the Climate Forecast System Reanalysis (CFSR), a modern global ocean–atmosphere reanalysis product. Our analysis reveals that ROARS exhibits reasonable climatology and interannual variability over the IAS region, with climatological SST errors less than 1 °C except along the coastlines. The anomaly correlation of the monthly SST and precipitation anomalies in ROARS are well over 0.5 over the Gulf of Mexico, Caribbean Sea, Western Atlantic and Eastern Pacific Oceans. A highlight of the ROARS simulation is its resolution of the loop current and the episodic eddy events off of it. This is rather poorly simulated in the CFSR. This is also reflected in the simulated, albeit, higher variance of the sea surface height in ROARS and the lack of any variability in the sea surface height of the CFSR over the IAS. However the anomaly correlations of the monthly heat content anomalies of ROARS are comparatively lower, especially over the Gulf of Mexico and the Caribbean Sea. This is a result of ROARS exhibiting a bias of underestimation (overestimation) of high (low) clouds. ROARS like CFSR is also able to capture the Caribbean Low Level Jet and its seasonal variability reasonably well.     

Downscaling over Vietnam using the stretched-grid CCAM: verification of the mean and interannual variability of rainfall

Rainfall over Vietnam is highly variable from north to south, due to the interaction of the monsoonal winds with the terrain. There is high rainfall from April to September, and little rainfall from October to March (except along the central Vietnam coast). In order to study the ability of the Commonwealth Scientific and Industrial Research Organisation stretched-grid Conformal Cubic Atmospheric Model (CCAM) to capture the climatic and interannual variability of rainfall, downscaled simulations at approximately 20 km horizontal resolution over the region were produced for the period 1979–2001. A scale-selective digital filter was used to force the winds, temperature and sea-level pressure from the ERA-Interim reanalysis for length scales greater than about 700 km. For wind and temperature, the forcing is applied for pressure-sigma levels above about 0.9. ERA-Interim sea surface temperatures were used over the oceans. The simulations were primarily validated against the gridded Asian Precipitation Highly Resolved Observational Data Integration Toward Evaluation of the Water Resources rainfall dataset and station observations using standard statistical methods. It was found that CCAM reproduces well the amount and spatial variability of rainfall, with an area-averaged bias for the entire study domain of less than 1 mm day^?1; CCAM is also able to capture the rainfall pattern under different El Niño Southern Oscillation phases reasonably well for the dry season. For interannual variability, the simulation generally performed better for North and Central Vietnam than for South Vietnam, where rainfall variability was overestimated.     

Heat balance and eddies in the Peru-Chile current system

The Peru-Chile current System (PCS) is a region of persistent biases in global climate models. It has strong coastal upwelling, alongshore boundary currents, and mesoscale eddies. These oceanic phenomena provide essential heat transport to maintain a cool oceanic surface underneath the prevalent atmospheric stratus cloud deck, through a combination of mean circulation and eddy flux. We demonstrate these behaviors in a regional, quasi-equilibrium oceanic model that adequately resolves the mesoscale eddies with climatological forcing. The key result is that the atmospheric heating is large (>50 W m^?2) over a substantial strip >500 km wide off the coast of Peru, and the balancing lateral oceanic flux is much larger than provided by the offshore Ekman flux alone. The atmospheric heating is weaker and the coastally influenced strip is narrower off Chile, but again the Ekman flux is not sufficient for heat balance. The eddy contribution to the oceanic flux is substantial. Analysis of eddy properties shows strong surface temperature fronts and associated large vorticity, especially off Peru. Cyclonic eddies moderately dominate the surface layer, and anticyclonic eddies, originating from the nearshore poleward Peru-Chile Undercurrent (PCUC), dominate the subsurface, especially off Chile. The sensitivity of the PCS heat balance to equatorial intra-seasonal oscillations is found to be small. We demonstrate that forcing the regional model with a representative, coarse-resolution global reanalysis wind product has dramatic and deleterious consequences for the oceanic circulation and climate heat balance, the eddy heat flux in particular.     

Evaluating a modified point-based method to downscale cell-based climate variable data to high-resolution grids

To address the demand for high spatial resolution gridded climate data, we have advanced the Daymet point-based interpolation algorithm for downscaling global, coarsely gridded data with additional output variables. The updated algorithm, High-Resolution Climate Downscaler (HRCD), performs very good downscaling of daily, global, historical reanalysis data from 1° input resolution to 2.5 arcmin output resolution for day length, downward longwave radiation, pressure, maximum and minimum temperature, and vapor pressure deficit. It gives good results for monthly and yearly cumulative precipitation and fair results for wind speed distributions and modeled downward shortwave radiation. Over complex terrain, 2.5 arcmin resolution is likely too low and aggregating it up to 15 arcmin preserves accuracy. HRCD performs comparably to existing daily and monthly US datasets but with a global extent for nine daily climate variables spanning 1948–2006. Furthermore, HRCD can readily be applied to other gridded climate datasets.     

High resolution re-analysis for the Baltic Sea region during 1965–2005 period

Regional reanalysis database BaltAn65+ comprising meteorological data for Baltic Sea region for the time period 1965–2005 is described. For data assimilation and hindcasts, the numerical weather prediction model HIRLAM 7.1.4 is applied, with 11 km horizontal and 60-layer vertical resolution. Reanalysis includes three-dimensional weather analysis data. Standard surface observations and meteorological soundings together with ship and buoy measurements from WMO observational network are used in analysis. Boundary fields are obtained from ECMWF ERA-40 global re-analysis. The BaltAn65+ can be considered as a regional refinement of ERA-40 for Baltic Sea region, providing the historical weather and climate data with enhanced spatial resolution, which is main motivation for creation of this novel reanalysis database.     

Scale dependency of regional climate modeling of current and future climate extremes in Germany

A 26-year simulation (1980–2005) was performed with the Weather Research and Forecast (WRF) model over the Volta Basin in West Africa. This was to investigate the ability of a climate version of WRF to reproduce present day temperature and precipitation over the Volta Basin. The ERA-Interim reanalysis and one realization of the ECHAM6 global circulation model (GCM) data were dynamically downscaled using two nested domains within the WRF model. The outer domain had a horizontal resolution of 50 km and covered the whole of West Africa while the inner domain had a horizontal resolution of 10 km. It was observed that biases in the respective forcing data were carried over to the RCM, but also the RCM itself contributed to the mean bias of the model. Also, the biases in the 50-km domain were transferred unchanged, especially in the case of temperature, to the 10-km domain, but, for precipitation, the higher-resolution simulations increased the mean bias in some cases. While in general, WRF underestimated temperature in both the outer (mean biases of ?1.6 and ?2.3 K for ERA-Interim and ECHAM6, respectively) and the inner (mean biases of ?0.9 K for the reanalysis and ?1.8 K for the GCM) domains, WRF slightly underestimated precipitation in the coarser domain but overestimated precipitation in the finer domain over the Volta Basin. The performance of the GCM, in general, is good, particularly for temperature with mean bias of ?0.7 K over the outer domain. However, for precipitation, the added value of the RCM cannot be overlooked, especially over the whole West African region on the annual time scale (mean biases of ?3% for WRF and ?8% for ECHAM6). Over the whole Volta Basin and the Soudano-Sahel for the month of April and spring (MAM) rainfall, respectively, mean bias close to 0% was simulated. Biases in the interannual variability in both temperature and precipitation over the basin were smaller in the WRF than the ECHAM6. High spatial pattern correlations between 0.7 and 0.8 were achieved for the autumn precipitation and low spatial correlation in the range of 0.0 and 0.2 for the winter season precipitation over the whole basin and all the three belts over the basin.     

Intra-annual variation of atmospheric static stability in the Mediterranean region: a 60-year climatology

The seasonal characteristics of atmospheric static stability in the Mediterranean region are examined, for the 60-year period 1948–2007 and for the four 15-year sub-periods 1948–1962, 1963–1977, 1978–1992 and 1993–2007. S-Mode and T-Mode Factor Analysis are applied to the mean 5-day values of K static stability index over the Mediterranean region. Three dominant modes are revealed for both, the intra-annual variation and the spatial distribution of K-index. It is found that these modes are connected to the seasonal characteristics of the main atmospheric circulation systems affecting the region and the thermal properties of the Earth’s surface (land or sea). The differences among the results of the four sub-periods partially reflect the inter-decadal variations of the strength of the above factors.     

Gravity wave breaking in easterly flow over Greenland and associated low level barrier- and reverse tip-jets

A first evidence of severe turbulence in the lower stratosphere during easterly tropospheric flow over Greenland is presented. A numerical simulation shows the turbulence to be associated with gravity wave breaking and that simulating with a horizontal resolution of 3 km gives substantially greater and more realistic turbulence than at a 9 km horizontal resolution. It is concluded that real-time simulations at high resolutions would improve aviation forecasts. As the atmospheric flow impinges on South-Greenland a barrier jet, a reverse tip jet and amplified mountain waves with secondary wave breaking are generated at the same time.     

High-resolution modelling of the Antarctic surface mass balance, application for the twentieth, twenty first and twenty second centuries

About 75 % of the Antarctic surface mass gain occurs over areas below 2,000 m asl, which cover 40 % of the grounded ice-sheet. As the topography is complex in many of these regions, surface mass balance modelling is highly dependent on horizontal resolution, and studying the impact of Antarctica on the future rise in sea level requires physical approaches. We have developed a computationally efficient, physical downscaling model for high-resolution (15 km) long-term surface mass balance (SMB) projections. Here, we present results of this model, called SMHiL (surface mass balance high-resolution downscaling), which was forced with the LMDZ4 atmospheric general circulation model to assess Antarctic SMB variability in the twenty first and the twenty second centuries under two different scenarios. The higher resolution of SMHiL better reproduces the geographical patterns of SMB and increase significantly the averaged SMB over the grounded ice-sheet for the end of the twentieth century. A comparison with more than 3200 quality-controlled field data shows that LMDZ4 and SMHiL reproduce the observed values equally well. Nevertheless, field data below 2,000 m asl are too scarce to efficiently show the added value of SMHiL and measuring the SMB in these undocumented areas should be a future scientific priority. Our results suggest that running LMDZ4 at a finer resolution (15 km) may give a future increase in SMB in Antarctica that is about 30 % higher than by using its standard resolution (60 km) due to the higher increase in precipitation in coastal areas at 15 km. However, a part (～15 %) of these discrepancies could be an artefact from SMHiL since it neglects the foehn effect and likely overestimates the precipitation increase. Future changes in the Antarctic SMB at low elevations will result from the competition between higher snow accumulation and runoff. For this reason, developing downscaling models is crucial to represent processes in sufficient detail and correctly model the SMB in coastal areas.     

The role of horizontal resolution in simulating drivers of the global hydrological cycle

The role of atmospheric general circulation model (AGCM) horizontal resolution in representing the global energy budget and hydrological cycle is assessed, with the aim of improving the understanding of model uncertainties in simulating the hydrological cycle. We use two AGCMs from the UK Met Office Hadley Centre: HadGEM1-A at resolutions ranging from 270 to 60 km, and HadGEM3-A ranging from 135 to 25 km. The models exhibit a stable hydrological cycle, although too intense compared to reanalyses and observations. This over-intensity is explained by excess surface shortwave radiation, a common error in general circulation models (GCMs). This result is insensitive to resolution. However, as resolution is increased, precipitation decreases over the ocean and increases over the land. This is associated with an increase in atmospheric moisture transport from ocean to land, which changes the partitioning of moisture fluxes that contribute to precipitation over land from less local to more non-local moisture sources. The results start to converge at 60-km resolution, which underlines the excessive reliance of the mean hydrological cycle on physical parametrization (local unresolved processes) versus model dynamics (large-scale resolved processes) in coarser HadGEM1 and HadGEM3 GCMs. This finding may be valid for other GCMs, showing the necessity to analyze other chains of GCMs that may become available in the future with such a range of horizontal resolutions. Our finding supports the hypothesis that heterogeneity in model parametrization is one of the underlying causes of model disagreement in the Coupled Model Intercomparison Project (CMIP) exercises.     

Simulating the diurnal cycle of rainfall in global climate models: resolution versus parameterization

The effects of horizontal resolution and the treatment of convection on simulation of the diurnal cycle of precipitation during boreal summer are analyzed in several innovative weather and climate model integrations. The simulations include: season-long integrations of the Non-hydrostatic Icosahedral Atmospheric Model (NICAM) with explicit clouds and convection; year-long integrations of the operational Integrated Forecast System (IFS) from the European Centre for Medium-range Weather Forecasts at three resolutions (125, 39 and 16 km); seasonal simulations of the same model at 10 km resolution; and seasonal simulations of the National Center for Atmospheric Research (NCAR) low-resolution climate model with and without an embedded two-dimensional cloud-resolving model in each grid box. NICAM with explicit convection simulates best the phase of the diurnal cycle, as well as many regional features such as rainfall triggered by advancing sea breezes or high topography. However, NICAM greatly overestimates mean rainfall and the magnitude of the diurnal cycle. Introduction of an embedded cloud model within the NCAR model significantly improves global statistics of the seasonal mean and diurnal cycle of rainfall, as well as many regional features. However, errors often remain larger than for the other higher-resolution models. Increasing resolution alone has little impact on the timing of daily rainfall in IFS with parameterized convection, yet the amplitude of the diurnal cycle does improve along with the representation of mean rainfall. Variations during the day in atmospheric prognostic fields appear quite similar among models, suggesting that the distinctive treatments of model physics account for the differences in representing the diurnal cycle of precipitation.     

Storm surges in the Western Baltic Sea: the present and a possible future

Globally-coupled climate models are generally capable of reproducing the observed trends in the globally averaged atmospheric temperature or mean sea level. However, the global models do not perform as well on regional/local scales. Here, we present results from four 100-year ocean model experiments for the Western Baltic Sea. In order to simulate storm surges in this region, we have used the General Estuarine Transport Model (GETM) as a high-resolution local model (spatial resolution ≈ 1?km), nested into a regional atmospheric and regional oceanic model in a fully baroclinic downscaling approach. The downscaling is based on the global model ECHAM5/MPI-OM. The projections are imbedded into two greenhouse-gas emission scenarios, A1B and B1, for the period 2000–2100, each with two realisations. Two control runs from 1960 to 2000 are used for validation. We use this modelling system to statistically reproduce the present distribution of surge extremes. The usage of the high-resolution local model leads to an improvement in surge heights of at least 10% compared to the driving model. To quantify uncertainties associated with climate projections, we investigate the impact of enhanced wind velocities and changes in mean sea levels. The analysis revealed a linear dependence of surge height and mean sea level, although the slope parameter is spatially varying. Furthermore, the modelling system is used to project possible changes within the next century. The results show that the sea level rise has greater potential to increase surge levels than does increased wind speed. The simulations further indicate that the changes in storm surge height in the scenarios can be consistently explained by the increase in mean sea level and variation in wind speed.     

Moisture transport from the Atlantic to the Pacific basin and its response to North Atlantic cooling and global warming

Atmospheric moisture transport from the Atlantic to the Pacific basin plays an important role in regulating North Atlantic salinity and thus the strength of the thermohaline circulation. Potential changes in the strength of this moisture transport are investigated for two different climate-change scenarios: North Atlantic cooling representative of Heinrich events, and increased greenhouse gas (GHG) forcing. The effect of North Atlantic cooling is studied using a coupled regional model with comparatively high resolution that successfully simulates Central American gap winds and other important aspects of the region. Cooler North Atlantic sea surface temperature (SST) in this model leads to a regional decrease of atmospheric moisture but also to an increase in wind speed across Central America via an anomalous pressure gradient. The latter effect dominates, resulting in a 0.13 Sv (1 Sv = 10⁶ m³ s^?1) increase in overall moisture transport to the Pacific basin. In fresh water forcing simulations with four different general circulation models, the wind speed effect is also present but not strong enough to completely offset the effect of moisture decrease except in one model. The influence of GHG forcing is studied using simulations from the Intergovernmental Panel on Climate Change archive. In these simulations atmospheric moisture increases globally, resulting in an increase of moisture transport by 0.25 Sv from the Atlantic to Pacific. Thus, in both scenarios, moisture transport changes act to stabilize the thermohaline circulation. The notion that the Andes effectively block moisture transport from the Atlantic to the Pacific basin is not supported by the simulations and atmospheric reanalyses examined here. This indicates that such a blocking effect does not exist or else that higher resolution is needed to adequately represent the steep orography of the Andes.