Ensemble-Based Observation Targeting for Improving Ozone Prediction in Houston and the Surrounding Area

This study examines the effectiveness of targeted meteorological observations for improving ozone prediction in Houston and the surrounding area based on perfect-model simulation experiments. Supplementary observations are targeted for the location that has the highest impact factor (maximum Kalman gain) estimated from an ensemble and is expected to minimize ozone forecast uncertainty at the verification time. It is found that the observational impact factor field varies with time and is sensitive to ensemble resolutions and physics parameterizations. The efficiency of observation targeting is further examined through assimilating observations in areas with different impact factors using an ensemble Kalman filter. It is found that the ensemble sensitivity analysis is capable of locating supplementary observations that may reduce meteorological and ozone forecast error, but not as effectively as expected.     

Dead Sea Basin Imaged by Ambient Seismic Noise Tomography

In the framework of the Dead Sea Integrated Research project (DESIRE), 59 seismological stations were deployed in the region of the Dead Sea Basin. Twenty of these stations recorded data of sufficiently high quality between May and September 2007 to be used for ambient seismic noise analysis. Empirical Green’s functions are extracted from cross-correlations of long term recordings. These functions are dominated by Rayleigh waves, whose group velocities can be measured in the frequency range from 0.1 to 0.5 Hz. Analysis of positive and negative correlation lags of the Green’s functions makes it possible to identify the direction of the source of the incoming energy. Signals with frequencies higher than 0.2 Hz originate from the Mediterranean Sea, while low frequencies arrive from the direction of the Red Sea. Travel times of the extracted Rayleigh waves were measured between station pairs for different frequencies, and tomographically inverted to provide independent velocity models. Four such 2D models were computed for a set of frequencies, all corresponding to different sampling depths, and thus together giving an indication of the velocity variations in 3D extending to a depth of 10 km. The results show low velocities in the Dead Sea Basin, consistent with previous studies suggesting up to 8 km of recent sedimentary infill in the Basin. The complex structure of the western margin of the Basin is also observed, with sedimentary infill present to depths not exceeding 5 km west of the southern part of the Dead Sea. The high velocities associated with the Lisan salt diapir are also observed down to a depth of ~5 km. The reliability of the results is confirmed by checkerboard recovery tests.     

The Effects of the Earth’s Curvature on Gravity and Geoid Calculations

While it is obvious that large-scale gravity studies should account for the sphericity of the Earth, each case should be examined. If a geometry model is very large for the 3D-gravity calculation, it cannot be correctly defined in Cartesian coordinates. Because of the Earth’s curvature it is necessary to use spherical coordinates, the importance of which is shown in this paper. The calculation of the gravity for a cylinder reveals, 1 m above the center of the cylinder, a relative difference of 13% between the models with Cartesian and spherical coordinates.     

Homogenization of Earthquake Catalog for Northeast India and Adjoining Region

A catalog for northeast India and the adjoining region for the period 1897–2009 with 4,497 earthquakes events is compiled for homogenization to moment magnitude M _w,GCMT in the magnitude range 3–8.7. Relations for conversion of m _b and M _s magnitudes to M _w,GCMT are derived using three different methods, namely, linear standard regression, inverted standard regression (ISR) and orthogonal standard regression (OSR), for different magnitude ranges based on events data for the catalog period 1976–2006. The OSR relations for M _s to M _w,GCMT conversion derived in this paper have significantly lower errors in regression parameters compared to the relations reported in other studies. Since the number of events with magnitude ≥7 for this region is scanty, we, therefore, considered whole India region to obtain the regression relationships between M _w,GCMT and M _s,ISC. A relationship between M _w,GCMT and M _w,NEIC is also obtained based on 17 events for the range 5.2 ≤ magnitude ≤ 6.6. A unified homogeneous catalog prepared using the conversion relations derived in this paper can serve as a reference catalog for seismic hazard assessment studies in northeast India and the adjoining region.     

Joint inversion of teleseismic P waveforms and surface-wave group velocities from ambient seismic noise in the Bohemian Massif

Joint inversion of teleseismic P-waveforms and local group velocities of surface waves retrieved from ambient seismic noise has been performed to model velocity structure of the crust and uppermost mantle of the Bohemian Massif. We analysed P-waveforms of 381 teleseismic earthquakes recorded at 54 broadband seismic stations located on the territory of the Czech Republic and in its close surroundings. Group velocities of Rayleigh and Love surface waves were obtained by cross-correlating long-term recordings of seismic noise. The basis for waveform inversion is the well-known methodology of P-to-S receiver functions constructed from converted phases. Due to instabilities in direct inversion of receiver functions caused by the necessity of applying deconvolution, we propose an alternative formulation to fit observed and calculated radial components of P waveforms. The joint inversion is transformed into a search for the minimum of the cost function defined as a weighted sum of waveform and group velocity misfits. With the use of the robust stochastic optimizer (Differential Evolution Algorithm), neither derivatives nor a starting model are needed. The task was solved for 1D layered isotropic models of the crust and the uppermost mantle. We have performed a sequence of inversions with models containing one, two, three and four layers above a half-space. By using statistical criteria (F-test) we were able to select the simplest velocity models satisfying data and representing local geological structures. Complex crustal models are typical for stations located close to boundaries of major tectonic units. The relatively low average P to S wave-velocity ratio is in agreement with the generally accepted view that the BM crust is predominantly felsic.     

Comparison of the FORT approximation of the coupling ray theory with the Fourier pseudospectral method

The standard ray theory (RT) for inhomogeneous anisotropic media does not work properly or even fails when applied to S-wave propagation in inhomogeneous weakly anisotropic media or in the vicinity of shear-wave singularities. In both cases, the two shear waves propagate with similar phase velocities. The coupling ray theory was proposed to avoid this problem. In it, amplitudes of the two S waves are computed by solving two coupled, frequency-dependent differential equations along a common S-wave ray. In this paper, we test the recently developed approximation of coupling ray theory (CRT) based on the common S-wave rays obtained by first-order ray tracing (FORT). As a reference, we use the Fourier pseudospectral method (FM), which does not suffer from the limitations of the ray method and yields very accurate results. We study the behaviour of shear waves in weakly anisotropic media as well as in the vicinity of intersection, kiss or conical singularities. By comparing CRT and RT results with results of the FM, we demonstrate the clear superiority of CRT over RT in the mentioned regions as well as the dangers of using RT there.     

Sea level variability at Antalya and Mente? tide gauges in Turkey: atmospheric, steric and land motion contributions

Sea level trends and interannual variability at Antalya and Menteş tide gauges are investigated during the 1985–2001 period, quantifying the roles of atmospheric, steric and local land motion contributions. Tide gauge sea level measurements, temperature/salinity climatologies and GPS data are used in the analyses and the results are compared with the output of a barotropic model forced by atmospheric pressure and wind. The overall sea level trends at two tide gauges collocated with GPS are in the range of 5.5 to 7.9 mm/yr during the study period, but showing different behaviour in the sub-periods 1985–1993 and 1993–2001 due to variations in the contributing factors both in space and time. After the removal of the atmospheric forcing and steric contribution from sea level records, the resulting trends vary between 1.9 to 4.5 mm/yr in Antalya and −1.2 to −11.6 mm/yr in Menteş depending on the period considered. Vertical land movement estimated from GPS data seems to explain the high positive residual trend in Antalya during the whole period. On the other hand, the source of the highly negative sea level trend of about −14 mm/yr in Menteş during 1985–1993 could not be resolved with the available datasets. Interannual variability of wind and atmospheric pressure appear to dominate the sea level at both tide gauges during the study period. Atmospheric and steric contributions together account for ∼50% of the total sea level variance at interannual time scales. Mass induced sea level variations which were not considered in this study may help to close the sea level trend budgets as well as to better explain the interannual sea level variance.     

Emergence patterns and microhabitat preference of aquatic dance flies (Empididae; Clinocerinae and Hemerodromiinae) on a longitudinal gradient of barrage lake system

Aquatic dance flies (Empididae; Clinocerinae and Hemerodromiinae) are important components of freshwater assemblages, especially in running waters. They are predators as larvae and adults and thus essential for understanding aquatic food webs. This study was conducted in Plitvice lakes National Park (Croatia) representing a wide variety of freshwater habitats (springs, streams, lakes and tufa barriers). Adults were collected monthly from March 2007 until March 2009 using pyramid-type emergence traps at 13 locations. A total of 3865 specimens comprising 18 species were collected. The dominant genus was Chelifera, while the most abundant species was Hemerodromia unilineata. All species were univoltine except Chelifera precabunda, Chelifera pyrenaica and Chelifera stigmatica that were bivoltine. Considerable differences in composition and structure of aquatic dance flies assemblages were recorded along a longitudinal gradient of studied sites, primarily related to differences in physical and chemical parameters of water. Water temperature was the main factor influencing the timing of emergence. Hemerodromia species preferred variable water temperature throughout the year while the majority of the Chelifera species preferred stable water temperature characteristic of spring sites. Furthermore, discharge affected assemblage composition of aquatic dance flies. The highest abundance of aquatic dance flies was recorded in lotic habitats with fast water current over substrates of moss, gravel and particulate tufa with detritus. These results give a new insight on microhabitat preference and their distribution on unique karstic habitats.     

Zero-order ray-theory Green tensor in a heterogeneous anisotropic elastic medium

Although synthetic seismograms can, in many cases, be generated by direct numerical methods, approximate analytical ray-theory solutions are often very useful in forward and especially inverse problems of seismic wave propagation in complex 3-D media. In this paper, the frequency-domain zero-order ray-theory Green tensor in a heterogeneous anisotropic elastic medium is derived from the zero-order ray-theory approximation using the representation theorem applied in ray-centred coordinates.