The role of prior model calibration on predictions with ensemble Kalman filter

This paper, based on a real world case study (Limmat aquifer, Switzerland), compares inverse groundwater flow models calibrated with specified numbers of monitoring head locations. These models are updated in real time with the ensemble Kalman filter (EnKF) and the prediction improvement is assessed in relation to the amount of monitoring locations used for calibration and updating. The prediction errors of the models calibrated in transient state are smaller if the amount of monitoring locations used for the calibration is larger. For highly dynamic groundwater flow systems a transient calibration is recommended as a model calibrated in steady state can lead to worse results than a noncalibrated model with a well-chosen uniform conductivity. The model predictions can be improved further with the assimilation of new measurement data from on-line sensors with the EnKF. Within all the studied models the reduction of 1-day hydraulic head prediction error (in terms of mean absolute error [MAE]) with EnKF lies between 31% (assimilation of head data from 5 locations) and 72% (assimilation of head data from 85 locations). The largest prediction improvements are expected for models that were calibrated with only a limited amount of historical information. It is worthwhile to update the model even with few monitoring locations as it seems that the error reduction with EnKF decreases exponentially with the amount of monitoring locations used. These results prove the feasibility of data assimilation with EnKF also for a real world case and show that improved predictions of groundwater levels can be obtained.     

Rock salt mass in the Paleoproterozoic sequence of the Onega trough in Karelia (from the Onega parametric well data)

A thick (200 m) rock salt mass covering Late Archean granitoids was exposed for the first time in the Early Proterozoic volcanogenic-sedimentary formations in the Onega trough of the east Baltic Shield by the Onega parametric well. The mineral composition of salts, their geochemical features, and the isotopic composition of carbonate carbon and oxygen have been studied. After fluid inclusions present in salts, their metamorphism temperature and isotopic composition of helium and argon were determined. The obtained results give evidence of the fact that rock salts and magnesites associated with them were formed in an evaporate basin with participation of deep crustal processes. The age of the underlying granitoids (2.716 ± 9 Ma) is determined using the Pb—Pb method.     

Developing a modern pollen–climate calibration data set for Norway

Bjune, A. E., Birks, H. J. B., Peglar, S. M. & Odland, A. 2010: Developing a modern pollen–climate calibration data set for Norway. Boreas, Vol. 39, pp. 674–688. 10.1111/j.1502‐3885.2010.00158.x. ISSN 0300‐9483. Modern pollen–climate data sets consisting of modern pollen assemblages and modern climate data (mean July temperature and mean annual precipitation) have been developed for Norway based on 191 lakes and 321 lakes. The original 191‐lake data set was designed to optimize the distribution of the lakes sampled along the mean July temperature gradient, thereby fulfilling one of the most critical assumptions of weighted‐averaging regression and calibration and its relative, weighted‐averaging partial least‐squares regression. A further 130 surface samples of comparable taphonomy, taxonomic detail and analyst became available as a result of other projects. These 130 samples, all from new lakes, were added to the 191‐lake data set to create the 321‐lake data set. The collection and construction of these data sets are outlined. Numerical analyses involving generalized linear modelling, constrained ordination techniques, weighted‐averaging partial least‐squares regression, and two different cross‐validation procedures are used to asses the effects of increasing the size of the calibration data set from 191 to 321 lakes. The two data sets are used to reconstruct mean July temperature and mean annual precipitation for a Holocene site in northwest Norway and a Lateglacial site in west‐central Norway. Overall, little is to be gained by increasing the modern data set beyond about 200 lakes in terms of modern model performance statistics, but the down‐core reconstructions show less between‐sample variability and are thus potentially more plausible and realistic when based on the 321‐lake data set.     

A robust algorithm for ellipse-based discrete element modelling of granular materials

With recent research indicating the importance of the rolling mechanism of deformation in granular systems consisting of perfectly round particles, it has become popular to use ellipse-shaped particles in the Discrete Element Method (DEM) numerical model. Inherent in this technique is the need for accurately computing ellipse to ellipse intersection, in order to properly detect contact formation and compute relative contact velocities. However, the commonly used algorithms for computing ellipse-ellipse intersection are generally poorly conditioned and can be inaccurate. An alternate method for computing ellipse-ellipse intersection is developed and presented which results in a well-conditioned, stable and accurate contact detection method. These modification are incorporated into the general DEM algorithm.     

X-ray transient sources

The galactic distribution and physical nature of X-ray transient sources is investigated. Two types of transients are considered. The observational data on 41 X-ray transient sources are given, and the average parameters of hard and soft X-ray transients are estimated.     

Observation of Precipitating Systems over Complex Orography with Meteorological Doppler Radars: A Feasibility Study

Summary  This paper concerns the use of airborne or ground-based Doppler radars to observe precipitating systems over complex orography. As nearly all of the previous experiments involving Doppler radars were conducted over flat surfaces over the continents or the oceans, new techniques are needed firstly to separate ground clutter from meteorological signal and, in the case of airborne Doppler observations, to deduce navigational errors. Secondly, it is necessary to take the atmospheric circulation induced by orography into account in the three-dimensional wind field analysis. Variational techniques are presented to solve these problems. The proposed methods are tested with simulated ground-based and airborne Doppler radar observations for analytic flows over analytic terrains and for numerically simulated wind and reflectivity fields for the Brig event (22 September 1993) of heavy precipitation over the southern flank of the Alps (Cosma and Richard, 1998), and with actual airborne Doppler data relative to weak snow showers over the Rocky Mountains on 12 March 1995. Received March 22, 1999/Revised June 1, 1999     

Linear mechanism of generation and intencification of internal gravity waves in the ionosphere at their interaction with a nonuniform zonal wind: 1. Model of the medium and initial dynamic equations

The specific features of the generation and intensification of internal gravity wave structures in different atmospheric-ionospheric regions, caused by zonal local nonuniform winds (shear flows), are studied. The model of the medium has been explained and an initial closed system of equations has been obtained in order to study the linear and nonlinear dynamics of internal gravity waves (IGWs) when they interact with the geomagnetic field in a dissipative ionosphere (for the D, E, and F regions).     

A ten-to-twelve year variation in the stratosphere of the Northern Hemisphere

This review paper for STIB (Stratosphere-Troposphere Interaction and the Biosphere; a proposed core project for IGBP) summarizes several features of a recently discovered 10–12 year oscillation in the atmosphere on the Northern Hemisphere. The oscillation is especially strong in the stratosphere during the warmer half of the year, but it is evident in the stratosphere and troposphere also in winter if the data are grouped according to the phase of the Quasi-Biennial Oscillation of the wind in the equatorial stratosphere. During the 40 years with data available to describe the oscillation it was phase locked with the 11-year solar cycle.Affiliate Scientist, NCAR.The National Center for Atmospheric Research is sponsored by the National Science Foundation.