Wiener optimal filtering of GRACE data

We present a spatial averaging method for Gravity Recovery and Climate Experiment (GRACE) gravity-field solutions based on the Wiener optimal filtering. The optimal filter is designed from the least-square minimization of the difference between the desired and filtered signals. It requires information about the power spectra of the desired gravitational signal and the contaminating noise, which is inferred from the average GRACE degree-power spectrum. We show that the signal decreases with increasing spherical harmonic degree j with approximately j^−b, where b = 1.5 for GRACE data investigations. This is termed the Second Kaula rule of thumb for temporal variations of the Earth’s gravity field. The degree power of the noise increases, in the logarithmic scale, linearly with increasing j. The Wiener optimal filter obtained for the signal model with b = 1.5 closely corresponds to a Gaussian filter with a spatial half width of 4° (∼440 km). We find that the filtered GRACE gravity signal is relatively insensitive to the exponent b of the signal model, which indicates the robustness of Wiener optimal filtering. This is demonstrated using the GFZ-GRACE gravity-field solution for April 2004.     

Modeling of circulation and salinity in Hooghly estuary

The circulation and salinity distribution in the Hooghly Estuary have been studied by developing a two‐dimensional depth‐averaged numerical model for the lower estuary, where the flow is vertically well mixed. This has been coupled with a one‐dimensional model for the upper estuary, where the flow is assumed to be unidirectional and well mixed over the depth and breadth. The Hooghly River receives high freshwater discharge during the monsoon season (June to September), which has significant effect on the salinity distribution in the estuary. The model‐simulated currents, elevations, and salinities are in good agreement with observations during the dry season. However, during the wet season the computed salinities seem to deviate slightly from the observed values.     

An inter-comparison of tidal solutions computed with a range of unstructured grid models of the Irish and Celtic Sea Regions

Three finite element codes, namely TELEMAC, ADCIRC and QUODDY, are used to compute the spatial distributions of the M₂, M₄ and M₆ components of the tide in the sea region off the west coast of Britain. This region is chosen because there is an accurate topographic dataset in the area and detailed open boundary M₂ tidal forcing for driving the model. In addition, accurate solutions (based upon comparisons with extensive observations) using uniform grid finite difference models forced with these open boundary data exist for comparison purposes. By using boundary forcing, bottom topography and bottom drag coefficients identical to those used in an earlier finite difference model, there is no danger of comparing finite element solutions for “untuned unoptimised solutions” with those from a “tuned optimised solution”. In addition, by placing the open boundary in all finite element calculations at the same location as that used in a previous finite difference model and using the same M₂ tidal boundary forcing and water depths, a like with like comparison of solutions derived with the various finite element models was possible. In addition, this open boundary was well removed from the shallow water region, namely the eastern Irish Sea where the higher harmonics were generated. Since these are not included in the open boundary, forcing their generation was determined by physical processes within the models. Consequently, an inter-comparison of these higher harmonics generated by the various finite element codes gives some indication of the degree of variability in the solution particularly in coastal regions from one finite element model to another. Initial calculations using high-resolution near-shore topography in the eastern Irish Sea and including “wetting and drying” showed that M₂ tidal amplitudes and phases in the region computed with TELEMAC were in good agreement with observations. The ADCIRC code gave amplitudes about 30 cm lower and phases about 8° higher. For the M₄ tide, in the eastern Irish Sea amplitudes computed with TELEMAC were about 4 cm higher than ADCIRC on average, with phase differences of order 5°. For the M₆ component, amplitudes and phases showed significant small-scale variability in the eastern Irish Sea, and no clear bias between the models could be found. Although setting a minimum water depth of 5 m in the near-shore region, hence removing wetting and drying, reduced the small-scale variability in the models, the differences in M₂ and M₄ tide between models remained. For M₆, a significant reduction in variability occurred in the eastern Irish Sea when a minimum 5-m water depth was specified. In this case, TELEMAC gave amplitudes that were 1 cm higher and phases 30° lower than ADCIRC on average. For QUODDY in the eastern Irish Sea, average M₂ tidal amplitudes were about 10 cm higher and phase 8° higher than those computed with TELEMAC. For M₄, amplitudes were approximately 2 cm higher with phases of order 15° higher in the northern part of the region and 15° lower in the southern part. For M₆ in the north of the region, amplitudes were 2 cm higher and about 2 cm lower in the south. Very rapid M₆ tidal-phase changes occurred in the near-shore regions. The lessons learned from this model inter-comparison study are summarised in the final section of the paper. In addition, the problems of performing a detailed model–model inter-comparison are discussed, as are the enormous difficulties of conducting a true model skill assessment that would require detailed measurements of tidal boundary forcing, near-shore topography and precise knowledge of bed types and bed forms. Such data are at present not available.     

Amplitudes of long-period PcP and the core-mantle boundary

The amplitudes of the core reflection PcP are sensitive to the wave velocities and densities in the neighborhood of the core-mantle boundary (CMB). We study the amplitude ratio of the long-period phases PcP and P from two South American deep-focus earthquakes with favorable fault-plane solution, depth and magnitude, as recorded by WWNSS and CSN stations in North America.Comparison is made with long-period PcP/P amplitude ratios, derived from theoretical seismograms for a variety of CMB models. Models from previous studies, which were mainly derived from short-period PcP observations and which are characterized by discrete layers above the CMB, are almost all inconsistent with the long-period data. The data also discriminate against low nonzero S velocities below the CMB. Simple first-order-discontinuity models of the CMB, for instance according to the Jeffreys-Bullen earth model or according to recent models based mainly on free oscillations, explain the data reasonably well.Model improvements are attempted by varying the P-velocity gradient above the CMB. The best amplitude fit is obtained for a rather strong decrease in P velocity with depth in this zone which, however, gives no acceptable traveltime fit for PcP. The scatter in body-wave amplitudes is considerable even for long-period waves and may prevent the correct assessment of that part of the amplitude variation of a phase with distance that is due to the variation of velocities and densities with depth alone.     

Evolution and observational signature of diffused antiworld

The existence of macroscopic regions with antibaryon excess in the baryon asymmetric Universe with general baryon excess is the possible consequence of practically all models of baryosynthesis. Diffusion of matter and antimatter to the border of antimatter domains defines the minimal scale of the antimatter domains surviving to the present time. A model of diffused antiworld is considered, in which the density within the surviving antimatter domains is too low to form gravitationally bound objects. The possibility to test this model by measurements of cosmic gamma ray fluxes is discussed. The expected gamma ray flux is found to be acceptable for modern cosmic gamma ray detectors and for those planned for the near future.     

Problems associated with direct displacement-based design of concrete bridges with single-column piers, and some suggested improvements

Currently available displacement-based design (DBD) procedures for bridges are critically evaluated with a view to identifying extensions and/or modifications of the procedure, for it to be applicable to final design of a fairly broad class of bridges. An improved direct DBD procedure is presented, including a suite of comprehensive design criteria and proper consideration of the degree of fixity of the pier top. The design of an overpass bridge (originally designed to a current European Code), applying the improved ‘direct’ displacement-based design (DDBD) procedure is presented and both ‘conventional’ and displacement-based designs are assessed using non-linear response-history analysis (NLRHA); comparisons are made in terms of both economy and seismic performance of the different designs. It is seen that DDBD provided a more rational base shear distribution among piers and abutments when compared to the force-based design procedure and adequately captured the displacement pattern, closely matching the results of the more rigorous NLRHA.