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.     

Meteorological excitation of the annual polar motion

Summary. Numerous studies have indicated that the annual term in the polar motion cannot be explained in any detail by meteorological/hydrological excitation and no reasonable alternative excitations have been put forward. Part of the problem has been that the hydrostatic adjustment of the oceans to the atmospheric pressure changes has traditionally been computed using the inverse barometer approach. This approach does not properly model the gravitational interaction between the atmosphere and oceans, and the inverse barometer theory is modified in this paper to account for this properly. The information necessary to compute the ocean tide and polar excitation caused by any change in the atmospheric pressure pattern is presented. The results of the application of this theory to two global atmospheric pressure data sets are examined and compared to results of other workers.
It is concluded that the atmosphere is observed well enough to answer the question of the annual excitation of polar motion and it is argued that the ground water excitation is the component with the largest error and remains the chief obstacle to the successful solution of this problem.     

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     

Refinement of the modified time-to-failure method for intermediate-term earthquake prediction

The modified time-to-failure method for intermediate-term earthquake prediction utilizes empirical relationships to reduce the number of unknown parameters providing a stable and unique solution set. The only unknown parameters in the modified time-to-failure method are the time and size of the impending main shock. The modified time-to-failure equation is used to model the precursory events and a prediction contour diagram is constructed with the magnitude and time-of-failure as the axes of the diagram. The root-mean-square (rms) is calculated for each set of time and magnitude on the prediction diagram representing the difference between the model (calculated) acceleration and the actual accelerated energy release of the precursory events. A small region, corresponding to the low rms region on the diagram, defines the prediction. The prediction has been shown to consistently under-estimate the magnitude and over-estimate the time-of-failure. These shortcomings are caused by an underestimation in energy release of the modified time-to-failure equation at the very end of the sequence. An empirical correction can be applied to the predicted results to minimize this problem. A main shock location search technique has been developed for use with the modified time-to-failure method. The location technique is used to systematically search an earthquake catalog and identify locations corresponding to precursory sequences that display accelerated energy releases. It has shown good results when applied in retrospective predictions, and is essential for the practical application of the modified time-to-failure method. In addition, an observed linear characteristic in long-term energy release can be used to minimize false predictions. The refined empirical relationships that eliminate or constrain unknown constants used in the modified time-to-failure method and the main shock location search technique are used in a practical application in the New Madrid Seismic Zone (NMSZ). The NMSZ, which is over due for a magnitude 6 event according to recurrence rates (Johnston and Nava, 1985), makes this region ideal for testing the method. One location was identified in the NMSZ as a high risk area for an event in the magnitude 4.5 range. The prediction, if accurate, is of scientific interest only because of the relatively small size of the main shock.     

Faunas from the upper silurian (upper ludlovian) in the ludlow-much wenlock district,England

Faunal lists for 17 sections in the Overton Formation are presented. Faunal diversity appears to increase towards the formations base while faunal dominance appears to increase towards its top.     

Brachiopods from China seas

A collection of brachiopods by the Institute of Oceanology, Academia Sinica (Qingdao), contains eight species from seven genera. Six of the species have been recorded previously from China seas—Lingula adamsi, L. anatina, Discinisca stella, Pelagodiscus atlanticus, Campages mariae, Terebratalia coreanica. Two species (Terebratulina hataiana andFrenulina sanguinolenta) have been described from other parts of the Pacific area. The apparent absence of any endemic species is a noteworthy feature of Chinese in comparison with Japanese faunas.     

Assessment of atmosphere-ocean general circulation model simulations of winter northern hemisphere atmospheric blocking

An assessment of six coupled Atmosphere-Ocean General Circulation Models (AOGCMs) is undertaken in order to evaluate their ability in simulating winter atmospheric blocking highs in the northern hemisphere. The poor representation of atmospheric blocking in climate models is a long-standing problem (e.g. D’Andrea et?al. in Clim Dyn 4:385–407, 1998), and despite considerable effort in model development, there is only a moderate improvement in blocking simulation. A modified version of the Tibaldi and Molteni (in Tellus A 42:343–365, 1990) blocking index is applied to daily averaged 500?hPa geopotential fields, from the ERA-40 reanalysis and as simulated by the climate models, during the winter periods from 1957 to 1999. The two preferred regions of blocking development, in the Euro-Atlantic and North Pacific, are relatively well captured by most of the models. However, the prominent error in blocking simulations consists of an underestimation of the total frequency of blocking episodes over both regions. A more detailed analysis revealed that this error was due to an insufficient number of medium spells and long-lasting episodes, and a shift in blocking lifetime distributions towards shorter blocks in the Euro-Atlantic sector. In the Pacific, results are more diverse; the models are equally likely to overestimate or underestimate the frequency at different spell lengths. Blocking spatial signatures are relatively well simulated in the Euro-Atlantic sector, while errors in the intensity and geographical location of the blocks emerge in the Pacific. The impact of models’ systematic errors on blocking simulation has also been analysed. The time-mean atmospheric circulation biases affect the frequency of blocking episodes, and the maximum event duration in the Euro-Atlantic region, while they sometimes cause geographical mislocations in the Pacific sector. The analysis of the systematic error in time-variability has revealed a negative relationship between the high-frequency variability of the transient eddies in the areas affected by blocking and blocking frequency. The blocking responses to errors in the low-frequency variability are different according to the region considered; the amplitude of the low-frequency variability is positively related to the blocking frequency and persistence in the Euro-Atlantic sector, while no such consistency is observed in the Pacific.