2009, Kawa, S.R., Stolarski, R.S., Newman, P.A., Douglass, A.R., Rex, M., Hofmann, D.J., Santee, M.L., Frieler, K.

The impact and significance of uncertainties in model calculations of stratospheric ozone loss resulting from known uncertainty in chemical kinetics parameters is evaluated in trajectory chemistry simulations for the Antarctic and Arctic polar vortices. The uncertainty in modeled ozone loss is derived from Monte Carlo scenario simulations varying the kinetic (reaction and photolysis rate) parameters within their estimated uncertainty bounds. Simulations of a typical winter/spring Antarctic vortex scenario and Match scenarios in the Arctic produce large uncertainty in ozone loss rates and integrated seasonal loss. The simulations clearly indicate that the dominant source of model uncertainty in polar ozone loss is uncertainty in the Cl2O 2 photolysis reaction, which arises from uncertainty in laboratory-measured molecular cross sections at atmospherically important wavelengths. This estimated uncertainty in JCl 2O2 from laboratory measurements seriously hinders our ability to model polar ozone loss within useful quantitative error limits. Atmospheric observations, however, suggest that the Cl2O2 photolysis uncertainty may be less than that derived from the lab data. Comparisons to Match, South Pole ozonesonde, and Aura Microwave Limb Sounder (MLS) data all show that the nominal recommended rate simulations agree with data within uncertainties when the Cl2O2 photolysis error is reduced by a factor of two, in line with previous in situ ClOx measurements. Comparisons to simulations using recent cross sections from Pope et al. (2007) are outside the constrained error bounds in each case. Other reactions producing significant sensitivity in polar ozone loss include BrO + ClO and its branching ratios. These uncertainties challenge our confidence in modeling polar ozone depletion and projecting future changes in response to changing halogen emissions and climate. Further laboratory, theoretical, and possibly atmospheric studies are needed.

2011, Boettle, M., Kropp, J.P., Reiber, L., Roithmeier, O., Rybski, D., Walther, C.

The assessment of coastal flood risks in a particular region requires the estimation of typical damages caused by storm surges of certain characteristics and annualities. Although the damage depends on a multitude of factors, including flow velocity, duration of flood, precaution, etc., the relationship between flood events and the corresponding average damages is usually described by a stage-damage function, which considers the maximum water level as the only damage influencing factor. Starting with different (microscale) building damage functions we elaborate a macroscopic damage function for the entire case study area Kalundborg (Denmark) on the basis of multiple coarse-graining methods and assumptions of the hydrological connectivity. We find that for small events, the macroscopic damage function mostly depends on the properties of the elevation model, while for large events it strongly depends on the assumed building damage function. In general, the damage in the case study increases exponentially up to a certain level and then less steep.

2013, Attwa, M., Günther, T.

Field and laboratory spectral induced polarization (SIP) measurements are integrated to characterize the hydrogeological conditions at the Schillerslage test site in Germany. The phase images are capable of monitoring thin peat layers within the sandy aquifers. However, the field results show limitations of decreasing resolution with depth. In comparison with the field inversion results, the SIP laboratory measurements show a certain shift in SIP response due to different compaction and sorting of the samples. The SIP data are analyzed to derive an empirical relationship for predicting the hydraulic conductivity (K). In particular, two significant but weak correlations between individual real resistivities (ρ') and relaxation times (τ), based on a Debye decomposition (DD) model, with measured K are found for the upper groundwater aquifer. The maximum relaxation time (τmax) and logarithmically weighted average relaxation time (τlw) show a better relation with K values than the median value τ50. A combined power law relation between individual ρ' and τ with K is developed with an expression of A · (ρ')B · (τlw)C, where A, B and C are determined using a least-squares fit between the measured and predicted K. The suggested approach with the calculated coefficients of the first aquifer is applied for the second. Results show good correlation with the measured K indicating that the derived relationship is superior to single phase angle models as Börner or Slater models.

2011, Jones, C.D., Hughes, J.K., Bellouin, N., Hardiman, S.C., Jones, G.S., Knight, J., Liddicoat, S., O'Connor, F.M., Andres, R.J., Bell, C., Boo, K.-O., Bozzo, A., Butchart, N., Cadule, P., Corbin, K.D., Doutriaux-Boucher, M., Friedlingstein, P., Gornall, J., Gray, L., Halloran, P.R., Hurtt, G., Ingram, W.J., Lamarque, J.-F., Law, R.M., Meinshausen, M., Osprey, S., Palin, E.J., Parsons, Chini, L., Raddatz, T., Sanderson, M.G., Sellar, A.A., Schurer, A., Valdes, P., Wood, N., Woodward, S., Yoshioka, M., Zerroukat, M.

The scientific understanding of the Earth's climate system, including thecentral question of how the climate system is likely to respond tohuman-induced perturbations, is comprehensively captured in GCMs and EarthSystem Models (ESM). Diagnosing the simulated climate response, andcomparing responses across different models, is crucially dependent ontransparent assumptions of how the GCM/ESM has been driven - especiallybecause the implementation can involve subjective decisions and may differbetween modelling groups performing the same experiment. This paper outlinesthe climate forcings and setup of the Met Office Hadley Centre ESM, HadGEM2-ES for the CMIP5 set of centennial experiments. We document theprescribed greenhouse gas concentrations, aerosol precursors, stratosphericand tropospheric ozone assumptions, as well as implementation of land-usechange and natural forcings for the HadGEM2-ES historical and futureexperiments following the Representative Concentration Pathways. Inaddition, we provide details of how HadGEM2-ES ensemble members wereinitialised from the control run and how the palaeoclimate and AMIPexperiments, as well as the "emission-driven" RCP experiments wereperformed.

2010, Baumgarten, G., Fiedler, J., Rapp, M.

Noctilucent clouds (NLC) in the polar summer mesopause region have been observed in Norway (69° N, 16° E) between 1998 and 2009 by 3-color lidar technique. Assuming a mono-modal Gaussian size distribution we deduce mean and width of the particle sizes throughout the clouds. We observe a quasi linear relationship between distribution width and mean of the particle size at the top of the clouds and a deviation from this behavior for particle sizes larger than 40 nm, most often in the lower part of the layer. The vertically integrated particle properties show that 65% of the data follows the linear relationship with a slope of 0.42±0.02 for mean particle sizes up to 40 nm. For the vertically resolved particle properties (Δz = Combining double low line 0.15 km) the slope is comparable and about 0.39±0.03. For particles larger than 40 nm the distribution width becomes nearly independent of particle size and even decreases in the lower part of the layer. We compare our observations to microphysical modeling of noctilucent clouds and find that the distribution width depends on turbulence, the time that turbulence can act (cloud age), and the sampling volume/time (atmospheric variability). The model results nicely reproduce the measurements and show that the observed slope can be explained by eddy diffusion profiles as observed from rocket measurements. © 2010 Author(s).

2013, Brocca, L., Liersch, S., Melone, F., Moramarco, T., Volk, M.

A framework for a comprehensive synthetic rainfall-runoff database was developed to study catchment response to a variety of rainfall events. The framework supports effective flood risk assessment and management and implements simple approaches. It consists of three flexible components, a rainfall generator, a continuous rainfallrunoff model, and a database management system. The system was developed and tested at two gauged river sections along the upper Tiber River (central Italy). One of the main questions was to investigate how simple such approaches can be applied without impairing the quality of the results. The rainfall-runoff model was used to simulate runoff on the basis of a large number of rainfall events. The resulting rainfallrunoff database stores pre-simulated events classified on the basis of the rainfall amount, initial wetness conditions and initial discharge. The real-time operational forecasts follow an analogue method that does not need new model simulations. However, the forecasts are based on the simulation results available in the rainfall-runoff database (for the specific class to which the forecast belongs). Therefore, the database can be used as an effective tool to assess possible streamflow scenarios assuming different rainfall volumes for the following days. The application to the study site shows that magnitudes of real flood events were appropriately captured by the database. Further work should be dedicated to introduce a component for taking account of the actual temporal distribution of rainfall events into the stochastic rainfall generator and to the use of different rainfall-runoff models to enhance the usability of the proposed procedure.

2013, Igel, J., Günther, T., Kuntzer, M.

Freshwater lenses, as important resource for drinking water, are sensitive to climate changes and sea level rise. To simulate this impact on the groundwater systems, hydraulic subsurface models have to be designed. Geophysical techniques can provide information for generating realistic models. The aim of our work is to show how ground-penetrating radar (GPR) investigations can contribute to such hydrological simulations. In the pilot area, Borkum island, GPR was used to map the shape of the groundwater table (GWT) and to characterise the aquifer. In total, 20 km of constant offset (CO) profiles were measured with centre frequencies of 80 and 200 MHz. Wave velocities were determined by common midpoint (CMP) measurements and vertical radar profiling (VRP) in a monitoring well. The 80 MHz CO data show a clear reflection at the groundwater table, whereas the reflection is weaker for the 200 MHz data. After correcting the GPR water tables for the capillary rise, they are in good accordance with the pressure heads of the observation wells in the area. In the centre of the island, the groundwater table is found up to 3.5 m above sea level, however it is lower towards the coastline and marshland. Some local depressions are observed in the region of dune valleys and around pumping stations of the local water supplier. GPR also reveals details within the sediments and highly-permeable aeolian sands can be distinguished from less-permeable marine sediments. Further, a silt loam layer below the water table could be mapped on a large area. The reflection characteristics indicates scattered erosion channels in this layer that cause it to be an aquitard with some leakage. GPR provides a high resolution map of the groundwater table and insight into the stratigraphy of the sediments and their hydraulic properties. This is valuable complementary information to the observation of sparsely distributed monitoring wells as input to hydraulic simulation.

2009, Schepanski, K., Tegen, I., MacKe, A.

We present a study of Saharan dust export towards the tropical North Atlantic using the regional dust emission, transport and deposition model LM-MUSCAT. Horizontal and vertical distribution of dust optical thickness, concentration, and dry and wet deposition rates are used to describe seasonality of dust export and deposition towards the eastern Atlantic for three typical months in different seasons. Deposition rates strongly depend on the vertical dust distribution, which differs with seasons. Furthermore the contribution of dust originating from the Bod́eĺe Depression to Saharan dust over the Atlantic is investigated. A maximum contribution of Bod́eĺe dust transported towards the Cape Verde Islands is evident in winter when the Bod́eĺe source area is most active and dominant with regard to activation frequency and dust emission. Limitations of using satellite retrievals to estimate dust deposition are highlighted.

2012, Sulzbacher, H., Wiederhold, H., Siemon, B., Grinat, M., Igel, J., Burschil, T., Günther, T., Hinsby, K.

A numerical, density dependent groundwater model is set up for the North Sea Island of Borkum to estimate climate change impacts on coastal aquifers and especially the situation of barrier islands in the Wadden Sea. The database includes information from boreholes, a seismic survey, a helicopter-borne electromagnetic (HEM) survey, monitoring of the freshwater-saltwater boundary by vertical electrode chains in two boreholes, measurements of groundwater table, pumping and slug tests, as well as water samples. Based on a statistical analysis of borehole columns, seismic sections and HEM, a hydrogeological model is set up. The groundwater model is developed using the finite-element programme FEFLOW. The density dependent groundwater model is calibrated on the basis of hydraulic, hydrological and geophysical data, in particular spatial HEM and local monitoring data. Verification runs with the calibrated model show good agreement between measured and computed hydraulic heads. A good agreement is also obtained between measured and computed density or total dissolved solids data for both the entire freshwater lens on a large scale and in the area of the well fields on a small scale. For simulating future changes in this coastal groundwater system until the end of the current century, we use the climate scenario A2, specified by the Intergovernmental Panel on Climate Change and, in particular, the data for the German North Sea coast. Simulation runs show proceeding salinisation with time beneath the well fields of the two waterworks Waterdelle and Ostland. The modelling study shows that the spreading of well fields is an appropriate protection measure against excessive salinisation of the water supply until the end of the current century.

2010, Lorenz, A., Held, H., Bauer, E., von Deimling, T.S.

Greenland ice-core data containing the 8.2 ka event are utilized by a model-data intercomparison within the Earth system model of intermediate complexity, CLIMBER-2.3 to investigate their potential for constraining the range of uncertain ocean diffusivity properties. Within a stochastic version of the model (Bauer et al. in Paleoceanography 19:PA3014, 2004) it has been possible to mimic the pronounced cooling of the 8.2 ka event with relatively good accuracy considering the timing of the event in comparison to other modelling exercises. When statistically inferring from the 8.2 ka event on diffusivity the technical difficulty arises to establish the related likelihood numerically per realisation of the uncertain model parameters: while mainstream uncertainty analyses can assume a quasi-Gaussian shape of likelihood, with weather fluctuating around a long term mean, the 8.2 ka event as a highly nonlinear effect precludes such an a priori assumption. As a result of this study the Bayesian Analysis leads to a sharp single-mode likelihood for ocean diffusivity parameters within CLIMBER-2.3. Depending on the prior distribution this likelihood leads to a reduction of uncertainty in ocean diffusivity parameters (e. g. for flat prior uncertainty in the vertical ocean diffusivity parameter is reduced by factor 2). These results highlight the potential of paleo data to constrain uncertain system properties and strongly suggest to make further steps with more complex models and richer data sets to harvest this potential. © The Author(s) 2009.