Third-order Uranus-Neptune planetary theory

In this part we calculate the secular and critical terms arising from the indirect part of the classical planetary Hamiltonian for Uranus and Neptune. We neglect in our expansions powers higher than the second in the eccentricity-inclination. Our required results, are expressed in terms of Poincaré variables.     

Approaches to environmental restoration of a polluted harbour with submerged archaeology: the Alexandria case study

Many invaluable underwater buildings of archaeological interest in Alexandria were discovered in 1996 at different sites in the Eastern Harbour of Alexandria. There is a belief that the best way to protect these invaluable heritages is to transfer them to an underwater park or museum. Obviously, the execution of such a project depends essentially upon the water quality (including water transparency) improving in the future. The harbour is presently polluted by discharge of wastewater effluents from different sources. It has recently been decided to restore this important coastal area through: (1) stopping the direct discharge of wastewater effluents into this semi-enclosed harbour in 1993 and (2) gradually reducing the discharge of the municipal wastewater through marine outfalls at two sites lying at the outer sides of the harbour. Zero discharge is expected to be effective by the end of the year 2001. The present work, therefore, is a follow up of the study of water quality in the harbour after 1993: in 1996 and 1999-2000. The water quality of an open sea reference station was also studied for comparison. The results reveal occurrence of an improvement of the environmental conditions in the harbour. The water has turned from being eutrophic to mesotrophic. The harbour is expected to become meso/oligotrophic as soon as the complete cessation of the discharge from the two outside sources is attained.     

Forward Modeling with Forced Neural Networks for Gravity Anomaly Profıle

In this paper, we introduce a new method called Forced Neural Network (FNN) to find the parameters of the object in geophysical section respect to gravity anomaly assuming the prismatic model. The aim of the geological modeling is to find the shape and location of underground structure, which cause the anomalies, in 2D cross section. At the first stage, we use one neuron to model the system and apply back propagation algorithm to find out the density difference. At the second level, quantization is applied to the density differences and mean square error of the system is computed. This process goes on until the mean square error of the system is small enough. First, we use FNN to two synthetic data, and then the Sivas–Gürün basin map in Turkey is chosen as a real data application. Anomaly values of the cross section, which is taken from the gravity anomaly map of Sivas–Gürün basin, are very close to those obtained from the proposed method.     

Geochemistry of the Pliocene basalts erupted along the Malatya-Ovacik fault zone (MOFZ), eastern Anatolia, Turkey: Implications for source characteristics and partial melting processes

The study area is located on the middle sector of the Malatya-Ovacık Fault Zone (MOFZ) in the eastern Anatolia. Four basaltic flows from bottom to top, which are tholeiitic in character and intercalated with Pliocene sedimentary rocks, were erupted along this fault zone. Chemical compositions of these flows reveal some differences between the first flow and others in terms of high-field strength elements (HFSEs) (e.g. Ti, Zr, Nb). Limited variations in compositions within the first flow and upper flows suggest a limited fractionation range. Trace-element patterns exhibit that all the flows have similar and OIB-like patterns without positive peak at Pb and a trough at Nb—Ta, indicating minimal or no crustal contribution. Rare-earth element (REE) patterns indicate that the first flow has flat patterns with negative Eu anomaly, whereas the upper flows have variable enrichments in LREE and depletions in HREE. La/Sm_N, Dy/Yb_N and Zr/Y ratios exhibit that the degree of partial melting decreases from the first flow to upper flows. Higher values of La/Yb_N ratio for the upper flows and depletions at Y and Yb on the trace-element patterns suggest the presence of garnet as a residual phase, which imply that the depth of partial melting took place solely in the garnet-stability field. OIB-like trace-element patterns and trace-element ratios (e.g. La/Nb, Ce/Y and Zr/Nb) emphasize that the melts forming the Arguvan basalt were originated from the asthenospheric mantle rather than the lithospheric mantle.     

Expansion of the first order planetary disturbing function by Smart's method

In this part we expand the indirect part of the planetary perturbing function by Smart's method, via Taylor's theorem. We neglect, in our expansion, terms of degree higher than the fourth with regard to the eccentricities and tangents of the inclinations.     

Mapping and attenuation of surface waves side scattered by near-surface diffractors

Near-surface diffractors are one of the problems in land seismic exploration. They can scatter the surface wave energy emanating from the seismic source and contaminate the signal received by seismic receivers. The scattered energy from the near-surface diffractors manifests itself on seismic shot gathers as strong hyperbolic events, called diffractions, masking the weakly reflected body waves. Diffractions present complications to most of surface-wave suppression schemes, especially when they have been scattered by scatterers outside the receiver line. Different methods have been used to eliminate diffractions from seismic data, including geophone arrays, filtering, and inverse scattering. Each of those methods has its own limitations. In this study, we present processing algorithms to map and attenuate near-surface diffractors of surface waves in seismic shot gathers. The mapping algorithm is based on semblance measurements and time–offset relations, while the attenuation algorithm is based on the least-square fitting of a source wavelet. The algorithms are applied on synthetic data from two different models. The first model has three near-surface diffractors, while the second model has three clusters of near-surface diffractors. Each cluster consists of three near-surface diffractors with a different geometry for each cluster. The results show that the proposed algorithms are successful in locating and attenuating most near-surface diffractors, except when the separation between individual diffractors is below the wavelength of the diffracted surface wave.     

Stellar Population Analysis of Galaxies based on Genetic Algorithms

We present a new method for determining the age and relative contribution of different stellar populations in galaxies based on the genetic algorithm. We apply this method to the barred spiral galaxy NGC 3384, using CCD images in U, B, V, R and I bands. This analysis indicates that the galaxy NGC 3384 is mainly inhabited by old stellar population (age >109yr). Some problems were encountered when numerical simulations are used for determining the contribution of different stellar populations in the integrated color of a galaxy. The results show that the proposed genetic algorithm can search efficiently through the very large space of the possible ages.     

Response of a tunnel deeply embedded in a viscoelastic medium

This paper presents a three‐dimensional energy‐based solution for the time‐dependent response of a deeply embedded and unsupported semi‐infinite tunnel of circular cross‐section. The tunnel is taken to be excavated quasi‐instantaneously from an infinite rock body that initially exhibits an isotropic stress state and that is made up of a homogeneous, isotropic and viscoelastic material. The viscoelastic behaviour is modelled by means of Burger's model, and the rock is taken to behave volumetrically linear elastic and to exhibit exclusively deviatoric creep. This viscoelastic problem is transformed into the Laplace domain, where it represents a quasi‐elastic problem. The displacement fields in the new solution are taken to be the products of independent functions that vary in the radial and longitudinal directions. The differential equations governing the displacements of the system and appropriate boundary conditions are obtained using the principle of minimum potential energy. The solutions for these governing equations in the Laplace domain are then obtained analytically and numerically using a one‐dimensional finite difference technique. The results are then transformed back into the time domain using an efficient numerical scheme. The accuracy of the new solution is comparable with that of a finite element analysis but requires much less computation effort. Copyright © 2011 John Wiley & Sons, Ltd.