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.     

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.     

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.     

Runoff quality prediction from small urban catchments using SWMM

The RUNOFF block of EPA's storm water management model (SWMM) was used to simulate the quantity and quality of urban storm water runoff from four relatively small sites (i.e. 5·97–23·56 ha) in South Florida, each with a specific predominant land use (i.e. low density residential, high density residential, highway and commercial). The objectives of the study were to test the applicability of this model in small subtropical urban catchments and provide modellers with a way to select appropriate input parameters to be used in planning studies. A total of 58 storm events, measured by the US Geological Survey (USGS), provided hyetographs, hydrographs and pollutant loadings for biological oxygen demand (BOD₅), total suspended solids (TSS), total Kjeldahl nitrogen (TKN) and lead (Pb), and were used for calibration of the model. Several other catchment characteristics, also measured or estimated by USGS, were used in model input preparation. Application of the model was done using the Green–Ampt equation for infiltration loss computation, a pollutant accumulation equation using a power build-up equation dependent on the number of dry days, and a power wash-off equation dependent on the predicted runoff rate. Calibrated quantity input parameters are presented and compared with suggested values in the literature. The impervious depression storage was generally found to be the most sensitive calibration parameter, followed by the Manning's roughness coefficients of conduit and overland flow, the Green–Ampt infiltration parameters and, finally, the pervious depression storage. Calibrated quality input parameters are presented in the form of regression equations, as a function of rainfall depth and the number of antecedent dry days. A total of 16 independent rainfall events were used for verification of the model, which showed a good comparison with observed data for both hydrographs and pollutant loadings. Average model predictions for the four constituent concentrations from the verification runs also showed good agreement with NURP published values in Florida and US sites. © 1998 John Wiley & Sons, Ltd.     

Tectonics of southwestern Mexico,isotopic evidence,nuclear Central America,Late Cretaceous break up

In the paleogeographic reconstruction of Mexico and northern Central America, an ever-increasing amount of evidence shows that the entire region is a collage of suspect terranes transported from abroad, whose timing and sense of motion are now beginning to be understood. Among these, the Chortis block (nuclear Central America) and the Baja California Peninsula have been proposed as pieces of continent separated from the Pacific coast of southwestern Mexico, that have moved either southeastward by the Farallon plate or northwestward by the Kula plate. Previous studies mainly confined to the northern margin of the Chortis block, confirmed a left-lateral displacement of 130 km in Neogene time. Further studies made northwestward along the Mexican coast provided a better understanding of magmatic and metamorphic processes in the area, and suggested times of detachment increased to 30 Ma, 40 Ma, and 66 Ma. The pre-detachment westernmost position of the block has changed, depending on the model chosen, from Puerto Vallarta and beyond, to the current position. Here we show that the isotopic mineral ages from coastal granites along the coast from Puerto Vallarta, Jalisco (80 Ma) to Puerto Angel, Oaxaca (11 Ma) record systematic decrease of cooling ages from NW to SE. This pattern is interpreted to result from the progressive uplift of rocks exposed at the present-day coast in that direction, such uplift occurred in response to the development of the Middle America Trench at the newly formed continental margin when the Chortis block was sliding at an average rate of 1.5 cm/year in a sinistral sense to its present position. Our results also constrain the position of the Kula-Farallon spreading axis north of Puerto Vallarta. These observations led us to conclude that several indicators point to this time and region for the onset of strike-slip drifting of the Chortis block toward its current position. Here, we also present several view points in terms of other possilble interpretations to different tectonic, geologic and isotopic data sets published recently by different authors.     

Decay of a highly eccentric satellite

The re-entry phase of a highly eccentric satellite is discussed. Numerical simulations allowing the prediction of the exact date of re-entry of a highly eccentric satellite are exposed.It is shown that under very particular circumstances the life of the satellite can be extended by a few days. The number of final revolutions of the rapidly contracting orbit depends critically on the air density between 70 km and 100 km.Re-entry of the European scientific satellite HEOS-1 predicted for 28 October, 1975 is near such a situation.     

Paleoceanography of surface waters in the Gulf of Mexico during the late Quaternary

Late Quaternary paleotemperatures and paleosalinities of surface waters of the Gulf of Mexico were estimated using a multivariate statistical analysis of census data of planktonic foraminifera. Two climatic extremes were selected for detailed basinwide study, the climatic optimum 125,000 yr ago and the glacial maximum 18,000 yr ago. In addition, patterns of climatic change were examined in seven piston cores from 127,000 yr ago to the present day. During the climatic optimum 125,000 yr ago temperature distributions in surface waters were similar to those of the present. The 22°C winter isotherm trended northeastward across the central basin and paleotemperatures decreased northward. Summer distributions were nearly homogeneous and ranged between 28° and 29°C. Winter salinities were 1‰ fresher than present values in the northmost Gulf and 0.4‰ fresher in the central basin. Summer salinities were similar during both times. In contrast, during the last glacial maximum temperatures were 1° to 2°C cooler in winter and 1°C cooler in summer, and isotherms formed a circular pattern in the Gulf during both seasons. Salinity was 0.3‰ fresher in winter than at present but 0.6‰ saltier in summer. Conditions deteriorated from the climatic optimum to the glacial maximum. In the Mexico Basin, winter temperatures were 2°C cooler from 75,000 to 45,000 yr ago (Y6 to Y3 Subzones), summer temperatures reached a minimum (3°C cooler) 32,000 yr ago (Y2–Y3 boundary), and seasonality reached minimal values (5°C) from 45,000 to 15,000 yr ago. All three parameters became similar in value to those in the Straits of Florida from 45,000 to 15,000 yr ago, suggesting that the exchange of surface waters was enhanced at this time between the two regions. Summer salinities remained similar to present conditions in the Mexico Basin, whereas, winter salinities increased 2‰ by 32,000 yr ago and then fell 0.5‰ until the glacial maximum ended. The Westerlies may have migrated southward over the Mexico Basin in winter from 32,000 to 15,000 yr ago.     

Facies Structure and Quantitative Parameters of Pleistocene Sediments of the Bering Sea

Lithofacies zoning is described for the first time for the Neo- and Eopleistocene of the Bering Sea. Four lithofacies sedimentation zones are distinguished: (I) terrigenous; (II) siliceous–terrigenous; (III) siliceous, and (IV) volcanoterrigenous ones. Corresponding maps were treated using Ronov volumetric method to quantify sedimentation parameters for distinguished lithofacies zones (subzones) and types of Pleistocene sediments. It was revealed that terrigenous sediments predominate over other sediments. Accumulation of the terrigenous sediments was more intense (by 1.4 times) in the Neopleistocene than in the Eopleistocene. The sedimentation rate of siliceous sediments of the Bowers Ridge in the Eopleistocene was two times higher than in the Neopleistocene.     

Evidence for short term changes in the intensity and spectrum of hard X-rays from Cyg X-1

The paper presents experimental evidence for the existence of fast intensity fluctuations with time scales of the order of a minute in the X-ray emission from Cyg X-1 at energies greater than 29 keV. Spectral variations over time intervals of 20–25 min are also observed in the same energy range. Whereas, similar intensity and spectral fluctuations have been reported earlier at lower energies the observations presented here is the first evidence for the existence of similar fluctuations at high energies.