Using paleoclimate proxy-data to select optimal realisations in an ensemble of simulations of the climate of the past millennium

We present and describe in detail the advantages and limitations of a technique that combines in an optimal way model results and proxy-data time series in order to obtain states of the climate system consistent with model physics, reconstruction of past radiative forcing and proxy records. To achieve this goal, we select among an ensemble of simulations covering the last millennium performed with a low-resolution 3-D climate model the ones that minimise a cost function. This cost function measures the misfit between model results and proxy records. In the framework of the tests performed here, an ensemble of 30 to 40 simulations appears sufficient to reach reasonable correlations between model results and reconstructions, in configurations for which a small amount of data is available as well as in data-rich areas. Preliminary applications of the technique show that it can be used to provide reconstructions of past large-scale temperature changes, complementary to the ones obtained by statistical methods. Furthermore, as model results include a representation of atmospheric and oceanic circulations, it can be used to provide insights into some amplification mechanisms responsible for past temperature changes. On the other hand, if the number of proxy records is too low, it could not be used to provide reconstructions of past changes at a regional scale.     

Consistency of observed winter precipitation trends in northern Europe with regional climate change projections

Often it is claimed that the recent changes in northern European climate are at least partly anthropogenic even though a human influence has not yet been successfully detected. Hence we investigate whether the recent changes are consistent with regional climate change projections. Therefore, trends in winter (DJF) mean precipitation in northern Europe are compared to human induced changes as predicted by a set of four regional climate model simulations. The patterns of recent trends and predicted changes match reasonably well as indicated by pattern correlation and the similarity is very likely not random. However, the model projections generally underestimate the recent change in winter precipitation. That is, the signal-to-noise ratio of the anthropogenic precipitation change is either rather low or the presently used simulations are significantly flawed in their ability to project changes into the future. European trends contain large signals related to the North Atlantic Oscillation (NAO), of which a major unknown part may be unrelated to the anthropogenic signal. Therefore, we also examine the consistency of recent and projected changes after subtracting the NAO signal in both the observations and in the projections. It turns out that even after the removal of the NAO signal, the pattern of trends in the observations is similar to those projected by the models. At the same time, the magnitude of the trends is considerably reduced and closer to the magnitude of the change in the projections.     

Modeled seasonality of glacial abrupt climate events

Greenland ice cores, as well as many other paleo-archives from the northern hemisphere, recorded a series of 25 warm interstadial events, the so-called Dansgaard-Oeschger (D-O) events, during the last glacial period. We use the three-dimensional coupled global ocean–atmosphere–sea ice model ECBILT-CLIO and force it with freshwater input into the North Atlantic to simulate abrupt glacial climate events, which we use as analogues for D-O events. We focus our analysis on the Northern Hemisphere. The simulated events show large differences in the regional and seasonal distribution of the temperature and precipitation changes. While the temperature changes in high northern latitudes and in the North Atlantic region are dominated by winter changes, the largest temperature increases in most other land regions are seen in spring. Smallest changes over land are found during the summer months. Our model simulations also demonstrate that the temperature and precipitation change patterns for different intensifications of the Atlantic meridional overturning circulation are not linear. The extent of the transitions varies, and local non-linearities influence the amplitude of the annual mean response as well as the response in different seasons. Implications for the interpretation of paleo-records are discussed.     

Ensemble standardization constraints on the influence of the tree growth trends in dendroclimatology

Climate Dynamics - Tree growth trends can affect the interpretation of the response of tree-ring proxies (especially tree-ring width) to climate in the low-frequency band, which in turn may limit...     

Methodology for construction, calibration and validation of a national hydrological model for Denmark

An integrated groundwater/surface water hydrological model with a 1 km² grid has been constructed for Denmark covering 43,000 km². The model is composed of a relatively simple root zone component for estimating the net precipitation, a comprehensive three-dimensional groundwater component for estimating recharge to and hydraulic heads in different geological layers, and a river component for streamflow routing and calculating stream–aquifer interaction. The model was constructed on the basis of the MIKE SHE code and by utilising comprehensive national databases on geology, soil, topography, river systems, climate and hydrology. The present paper describes the modelling process for the 7330 km² island of Sjælland with emphasis on the problems experienced in combining the classical paradigms of groundwater modelling, such as inverse modelling of steady-state conditions, and catchment modelling, focussing on dynamic conditions and discharge simulation. Three model versions with different assumptions on input data and parameter values were required until the performance of the final, according to pre-defined accuracy criteria, model was evaluated as being satisfactory. The paper highlights the methodological issues related to establishment of performance criteria, parameterisation and assessment of parameter values from field data, calibration and validation test schemes. Most of the parameter values were assessed directly from field data, while about 10 ‘free’ parameters were subject to calibration using a combination of inverse steady-state groundwater modelling and manual trial-and-error dynamic groundwater/surface water modelling. Emphasising the importance of tests against independent data, the validation schemes included combinations of split-sample tests (another period) and proxy-basin tests (another area).     

Modeling undertow due to random waves

A numerical model of undertow due to random waves is developed. The model includes three sub-models: (i) a model for multi-directional and multi-frequency random wave transformation, (ii) a surface roller evolution model, and (iii) a model for calculating the vertical distribution and the mean value of the undertow velocity. The calculation of wave trough level is performed based on a theory for the wave asymmetry. The model was successfully validated against small- and large-scale laboratory experiments. Thus, the model is expected to provide reliable input for the modeling of sediment transport and morphological change due to waves and currents.     

Spectra of poorly resolved ocean data

Standard temperature and pressure sensors on Aanderaa RCM8 current meters have a resolution of 0.024 °C and 0.6 bar, each equal to 1 digital number (value) over a range of 1024. It is shown that an 11-month deep-ocean temperature record using only four values can contain useful spectral information on internal wave motions. This is partially due to the modulation of high-frequency data by non-zero low-frequency (subinertial) variations. This result follows from the comparison of this record with artificial three- and four-value data constructed from temperature records observed in stronger stratified waters nearby. These artificial records show main features of the internal wave band similar to those observed in the original data spectra. Peaks at tidal harmonic frequencies and enhancements at sum-tidal-inertial interaction frequencies are preserved in the artificial data, but overall noise level (and thus the continuum spectral slope) is enhanced with respect to the properly resolved records (using 15 and 100 values). As a demonstration of the stable accuracy of the temperature sensors, the poorly resolved records provided an estimate of mean stratification to within 5% of the estimate using Seabird CTD data.Responsible Editor: Hans Burchard     

Detection of PAHs in seawater using surface-enhanced Raman scattering (SERS)

The laboratory characterization of a field-operable surface-enhanced Raman scattering sensor (SERS optode) is presented for the detection of aromatic hydrocarbons in seawater. The sensor has been developed for deployment with a robust underwater spectrograph. To meet the demands of the harsh seawater application, sol-gel derived SERS substrates were used. The calibration curves of six PAHs were determined to be of Langmuir adsorption isotherm type with limits of detection ranging from the microg l(-1) to ng l(-1) level. The experimentally determined adsorption constants varied strongly with the molecular weight of the analytes and correlated with their solubility. A mixture of five PAHs dissolved in seawater was investigated to demonstrate the utility of this method for screening. Emphasis was put on the interference from suspended particulate matter (SPM). The Raman measurement with backscattering configuration was shown to be immune against turbidities up to 1000 NTU. The physico-chemical interference arising from adsorption by the sediment was measured on-line by adding sediment to a PAH-spiked solution. According to the calibration curve, the PAH concentration decrease corresponded to more than 98% of the analyte being scavenged by the sediment.     

Modelling transverse dunes

Transverse dunes appear in regions of mainly unidirectional wind and high sand availability. A dune model is extended to two‐dimensional calculation of the shear stress. It is applied to simulate dynamics and morphology of three‐dimensional transverse dunes. In the simulations they seem to reach translational invariance and do not stop growing. Hence, simulations of two‐dimensional dune ?elds have been performed. Characteristic laws were found for the time evolution of transverse dunes. Bagnold's law of the dune velocity is modi?ed and reproduced. The interaction between transverse dunes led to the interesting conclusion that small dunes can travel over bigger ones. Copyright © 2004 John Wiley & Sons, Ltd.     

Comparison of abstraction scenarios simulated by SWAT and SWAT-MODFLOW

In hydrological modelling of catchments, wherein streams are groundwater-fed, an accurate representation of groundwater processes and their interaction with surface water is crucial. With this purpose, a coupled model was recently developed linking SWAT (Soil and Water Assessment Tool) with the fully-distributed groundwater model MODFLOW (Modular Groundwater Flow). In this study, SWAT and SWAT-MODFLOW were applied to a Danish groundwater-dominant catchment, simulating groundwater abstraction scenarios and assessing the benefits and drawbacks of SWAT-MODFLOW. Both models demonstrated good performance. However, SWAT-MODFLOW provided more realistic outputs when simulating abstraction: the decrease in streamflow was similar to the volume of water abstracted, while in SWAT the impact was negligible. SWAT also showed impacts on streamflow only when abstractions were taken from the shallow aquifer, not from the deep aquifer. Overall, SWAT-MODFLOW demonstrated wider possibilities for groundwater analysis, providing more insights than SWAT in supporting decision making in relation to environmental assessment.