ALBERT

Notes: This paper is concerned with a new technique (the terrestrial/photogrammetric (TP) technique) for the detection and compensation of systematic height errors in block aerial triangulation. The technique improves the height accuracy of such triangulation by reducing the original bridging distance, using new bands of corrected photogrammetric tie points as height control together with the original bands of the terrestrial (or ground) height control. A number of photogrammetric blocks with different characteristics and configurations have been tested in a comprehensive manner. The tests show that the TP technique gives consistent results over a wide range of circumstances and that it is a most effective, simple and inexpensive method for the compensation of systematic errors.〈section xml:id="abs1-2"〉〈title type="main"〉RésuméCet article relate une technique nouvelle permettant de détecter et de compenser les erreurs altimétriques systématiques dans l'aérotriangulation par blocs: la technique terrestre/photogrammétrique (TP). Cette technique améliore la précision altimétrique de l'aérotriangulation en réduisant la distance de l'aérocheminement original, grâce à l'utilisation d'un canevas de points d'appuis altimétriques composés de points de liaison photogrammétriques corrigés, en supplément du canevas original de points d'appui altimétriques déterminés sur le terrain.On a essayé de manière approfondie un grand nombre de blocs dotés de caractéristiques et de configurations diverses. Les essais montrent que la technique TP fournit des résultats fiables dans un large éventail de circonstances et que c'est une méthode de compensation des erreurs systématiques efficace, simple et peu coûteuse.〈section xml:id="abs1-3"〉〈title type="main"〉ZusammenfassungIm Artikel wird ein neues Verfahren behandelt (das terrestrisch-photogrammetrische Verfahren (TP)), das zur Aufdeckung und Kompensation system-atischer Höhenfehler in der Block triangulation dient. Es wird die Höhengenauigkeit durch Verkürzung der ursprünglichen Überbrückungsdistanz verbessert, indem neue Streifen korrigierter Verbindungspunkte als Höhenpasspunkte zusammen mit den ursprünglichen Reihen der terrestrischen Höhenpasspunkte verwendet werden. Eine Anzahl photogrammetrischer Blöcke unterschiedlicher Figur und Charakteristik wurde umfassend getestet. Die Tests zeigen, dass die TP-Methode konstante Ergebnisse für einen grossen Bereich van Möglichkeiten ergibt und dass es ein höchst effektives, einfaches und billiges Verfahren zur Eliminierung systematischer Fehler ist.

Notes: During June and July 1981, synthetic aperture radar imagery was acquired of selected sites in ten European countries as part of the European SAR 580 campaign. Evaluation of this imagery in the UK will be carried out by 18 approved experiments (out of 56 for Europe) when the processed data becomes available to investigators. To assist the evaluation process, a review is given of some recent research on applications of imaging radar in agricultural and urban areas. Research is described which considered the effects of terrain parameters and radar parameters on crop classification accuracy. The relative merits of a quantitative or a qualitative approach in classification of radar data are considered. Combined use of radar and Landsat data achieves a marked improvement in classification accuracy. The development of spatial filtering techniques may assist in the problem in urban areas resulting from the dependence of image tone on the relative alignment of ground features and radar look direction.〈section xml:id="abs1-2"〉〈title type="main"〉RésuméPendant les mois dejuin et juillet 1981, on a, dans le cadre de la campagne européenne SAR 580, recueilli des données de radar imageur à synthèse d'ouverture sur des sites choisis dans une dizaine de pays d'Europe. Ces données seront évaluées au Royaume-Uni au moyen de 18 expériences ayant été acceptées (sur un total de 56 menées en Europe); ce travail commencera dès que les données une fois traitées seront distribuées aux chercheurs. Pour accompagner cette procédure d'evaluation, on fournit une analyse de quelques recherches récentes sur les applications de l'imagerie radar aux zones agricoles et urbaines.L'une de ces recherches prenait en compte les effets des paramètres liés au terrain sur l'exactitude de la classification des cultures. Les mérites réciproques des approches quantitatives et qualitatives dans la classification des données radar sont également examinés.Le recours à une combinaison de données radar et Landsat permet d'améliorer de manière significative l'exactitude de la classification. La mise en oeuvre des techniques de filtrage spatial peut aider à réduire, dans les zones urbaines, les conséquences néfastes sur la tonalité de l'image, d'objets au sol en alignement relatif avec la direction d'observation du radar.〈section xml:id="abs1-3"〉〈title type="main"〉ZusammenfassungIm Juni und Juli 1981 erfolgten Radaraufnahmen mittels synthetischer Apertur von ausgewählten Gebieten in 10 europäischen Ländern als Teil der Europäischen SAR 580 Kampagne. Die Analyse dieser Aufzeichnungen wird im Vereinigten Königreich durch 18 (aus insgesamt 56 für Europa) bestätigte Experimente erfolgen, nachdem die Daten zur Verfügung stehen. Um den Analysenprozess zu unterstützen, wird ein Überblick über kürzliche Forschungen zu Radaraufzeichnungen in landwirtschaftlichen und städtischen Gebieten gegeben. Es werden auch Forschungen beschrieben, bei denen die Einflüsse von Geländeparametern auf die Genauigkeit der Ernteklassifizierung untersucht wurden. Der relative Wert eines quantitativen oder qualitativen Versuchs zur Klassifizierung von Radardaten wird abgeschätzt. Eine Kombination von Radarund Landsat-Daten ergibt eine merkliche Verbesserung bei Klassifizierungsgenauigkeit. Die Entwicklung von Verfahren der räumlichen Filterung kann von Wert bei städtischen Gebieten sein, weil es eine Abhängigkeit des Grauwertes von der relativen Ausdehnung von Geländecharakteristika und der Richtung der Radaraufnahme gibt.

Notes: Some aspects of the XYZBLC suite of programs for aerial triangulation are described. The Ordnance Survey Bute block of photography at a scale of 1: 24000 is computed by independent models using XYZBLC. Three of the strips in the block are treated as a test block and are computed with different configurations of control. Conclusions are drawn with regard to the usefulness of height control doublets and other dense control configurations.〈section xml:id="abs1-2"〉〈title type="main"〉Résumé Quelques aspects de l'ensemble XYZBLC de logiciels d'aérotriangulation sont décrits. Le bloc Bute de l'Ordnance Survey, composé de photographies à l'échelle du 1: 24000, a été calculé par modèles indépendants avec les logiciels XYZBLC. On a considéré trois des bandes du bloc comme formant un bloc d'essai et on les a calculées avec diverses combinaisons de canevas de points d'appui. On en tire des conclusions quant à l'utilité des doublets de points d'appui altimétriques et de diverses autres configurations denses de canevas d'appui. 〈section xml:id="abs1-3"〉〈title type="main"〉Zusammenfassung Beschreibung einiger Aspekte des XYZBLC-Programmsystems zur Aero-triangulation. Ein Block des Ordnance Survey im Bildmassstab 1: 24000 wurde damit nach der Methode der unabhängigen Modelle berechnet. Drei der Streifen des Blocks wurden als Testblock behandelt und mit verschiedenen Passpunktanordnungen berechnet. Es werden Schlussfolgerungen in bezug auf die Nützlichkeit von Zwillingspunkten zur Höhenorientierung und andere dichte Passpunktanordnungen gezogen.

Notes: The wide availability of reconnaissance frame photography, together with the advantage of high image resolution that such photography provides and the extreme flight conditions under which it can be taken, means that it is of great interest to photogrammetrists. However, the geometry of the image is invariably distorted due to the mode of operation of the camera. In this paper, the effects of the focal plane shutter and apparent image motion compensation on the geometry of the photographs are considered.〈section xml:id="abs1-2"〉〈title type="main"〉RésuméLe fait que la photographic de reconnaissance soil très répandue, conjointe-ment aux avantages que présente la grande résolution des images fournies par ces systèmes et les conditions extrêmes de vol dans lesquelles on peut les prendre, milite pour une prise en compte de lew intérêt par les photogrammètres. Toutefois la géométrie de ces images est systématiquement déformée par suite du mode de fonctionnement de la chambre. On traite dans cet article des effets de l'obturateur focal et de la compensation du filé de l'image sur la géométrie des photographies.〈section xml:id="abs1-3"〉〈title type="main"〉ZusammenfassungDie grosse Verbreitung von Erkundungsfotos zusammen mit der grossen Auflösung solcher Fotos und den extremen Bedingungen unter denen sie auf-genommen wurden, lassen diese auch für den Photogrammeter interessant erscheinen. Die Geometrie dieser Bilder ist jedoch infolge der speziellen Kamera-funktion verzerrt. Im Artikel werden die Einflüsse des Schlitzverschlusses und des Bildwanderungsausgleichs auf die Bildgeometrie untersucht.

Notes: Book Reviewed in this article:AERIAL PHOTOGRAPHY. By G. S. Kumar.MANUAL OF PHOTOGRAMMETRY (Fourth Edition). Editor-in-Chief: C. C. Slama.

Notes: The algorithm presented assumes a preliminary determination of the elements of inner orientation and bases its mathematical model on the collinearity condition, from which projective equations are derived. Linearisation of condition equations is based on Taylor's expansion, whilst the subsequent solution follows iteration cycles. The method aims to derive composite corrections to image coordinates, with the aid of test field data.〈section xml:id="abs1-2"〉〈title type="main"〉RésuméL'auteur présente un algorithme oú l'on suppose déterminés au préalable les éléments d'orientation interne. La modélisation mathématique s'appuie sur les conditions de co-linéarité et sur les équations qui en résultent en projection. On linéarise les équations de conditions en utilisant un développement en série de Taylor, et l'on obtient la solution correspondante par itération. La méthode vise à apporter des corrections aux coordonnées images en utilisant les données d'un polygone d'essai.〈section xml:id="abs1-3"〉〈title type="main"〉ZusammenfassungBei dem dargestellten Algorithmus wird die vorläufige Bestimmung der Elemente der inneren Orientierung angenommen, und das mathematische Modell basiert aufder Kollinearitätsgleichung, woraus perspektive Gleichungen abgeleitet werden. Die Linearisierung der Bedingungsgleichungen erfolgt durch Anwendung der Taylor-Reihe, was für die folgenden Berechnungen einen iterativen Prozess bedingt. Mit der Methode wird angestrebt, Bildkoordinaten mit Hilfe von Testfelddaten zu verbessern.

Notes: There is a growing interest in the potential of synthetic aperture radar (SAR) for mapping and for gathering intelligence for map revision. This paper first gives a simple guide to the characteristics of SAR and then describes a study carried out to investigate the geometric properties of the imagery and the information content which is of interest for topographic mapping and the production of large scale plans. The conclusions cautiously indicate a potential although better quality images are required and more information must be known about the optimum conditions which should obtain when the data are recorded.〈section xml:id="abs1-2"〉〈title type="main"〉RésuméOn assiste à un intérét croissant envers le radar à ouverture synthétique (SAR) et ses possibilityés pour la cartographie et pour la collecte de renseignements utiles à la révision. Ce texte commence par donner un simple rappel des caractéristiques de ces radars pour décrire ensuite une étude effectuée sur le contenu informatif et les proprietyés géométriques de l'imagerie radar, qui intèresse la cartographie topographique et la confection de plans à grande échelle.Les conclusions font prudemment état d'une potentialité mais soulignent qu'il faudrait disposer d'images de meilleure qualityé et avoir davantage d'informations sur les conditions optimales à respecter lors de l'enregistrement des données.〈section xml:id="abs1-3"〉〈title type="main"〉ZusammenfassungEs besteht zunehmendes Interesse an den Möglichkeiten von Radaraufnahmen mittels synthetischer Apertur (SAR) zur Kartenherstellung und Erfassung von Daten zur Laufendhaltung. Zunächst wird im Artikel eine einfache Einführung in die SAR-Charakteristika gegeben. Dann erfolgt eine Beschreibung eines Versuches zur Untersuchung der geometrischen Eigenschaften und des Informationsgehalts, der von Interesse für die Herstellung topographischer und grossmassstäbiger Karten ist. Vorsichtige Schlussfolgerungen ergeben Möglichkeiten zur Verwendung solcher Informationen, wobeijedoch eine bessere Bildqualität zu fordern ist. Auch müssen noch mehr Kenntnisse über die optimalen Bedingungen bezüglich der Datenregistrierung bestehen.

Notes: Landsat multispectral scanner (MSS) and return beam vidicon (RBV) image data in both analogue and digital formats were used to assess the performance of the MSS and RBV sensor systems, and to determine their potential for recording information associated with topographic or land use maps of selected Chinese study areas (Beijing, Tianjin and the Huang He elbow). MSS images permitted the detection of 40 per cent to 50 per cent of the features identified on maps of 1: 250000 scale and larger, whereas the higher resolution Landsat 3 RBV images provided 10 per cent to 20 per cent additional information. Completeness of information extracted from a colour cathode ray tube display exceeded that obtained from analogue single band black and white MSS images by only 5 per cent to 10 per cent. Only through major improvements in spatial resolution will it be possible to compile maps at scales of 1: 250000 and larger from satellite data. Significantly, it appears that instantaneous field of view values of 5 m to 10 m will be required to produce image data of sufficient resolution for topographic mapping or detailed thematic studies of the Asian environment.〈section xml:id="abs1-2"〉〈title type="main"〉Résumé On a utilisé les données de l'imagerie spatiale provenant du scanneur multispectral (MSS) et du tube vidicon à retour de faisceau (RBV) de Landsat, à la fois sous forme analogique et numérique, pour évaluer la performance des systèmes de capteurs MSS et RBV, et pour déterminer leur aptitude à contenir l'information figurant dans les cartes topographiques et les cartes d'occupation du sol, sur des polygones d'essais sélectionnés en Chine (Pekin, Tianjin, et le coude du Huan He). Les images MSS ont permis de détecter 40 à 50 pour cent des détails figurant sur les cartes à une échelle supérieure ou égale à 1: 250000, tandis que les images RBV à meilleure résolution de Landsat C fournissaient 10 à 20 pour cent d'informations supplémentaires. Le volume des informations extraites sur un écran cathodique de visualisation en couleurs n'a dépassé le volume obtenu à partir d'images d'une seule bande MSS en noir et blanc que de 5 à 10 pour cent seulement. Il semble que ce ne soit que par des améliorations importantes dans la résolution spatiale que l'on pourra tirer des données de satellites des cartes à des échelles supérieures ou égales à 1: 250000. Il apparait de manière significative qu'un champ instantané de prise de vues de 5 à 10 m (pixel) sera nécessaire pour obtenir des images ayant une résolution suffisante pour la cartographie topographique ou pour des études thématiques détaillées sur l'environnement asiatique.〈section xml:id="abs1-3"〉〈title type="main"〉Zusammenfassung Verwendung von Landsat-aufzeichnungen, die mit Hilfe des multispektralen Scanners (MSS) und des Return Beam Vidicon (RBV) aufgenommen wurden und sowohl in analoger als auch digitaler Form vorlagen. um die Leistungsfähigkeit dieser Aufnahmesysteme zu ermitteln und die Möglichkeiten zur Gewinnung von Informationen zur Herstellung von topographischen oder Landnutzungskarten in ausgewählten chinesischen Versuchsgebieten (Beijing, Tianjin undder Huang He Bogen). MSS-Bilder gestatteten die Erkennung von 40 bis 50% der auf Karten im Massstab 1: 250 000 oder grösser identifizierten Objekte, während die RBV-Bilder von Landsat 3, die eine grössere Auflösung hatten, 10 bis 20% zusätzliche Informationen lieferten.Die Vollständigkeit der Information, die von einem Farbdisplay erhalten wurde, überstieg die von analogen einkanaligen schwarzweissen MSS-Bildern erhaltenen nur um 5 bis 10%. Nur durch wesentliche Verbesserungen der Objektauflösung wird es möglich sein, Karten 1: 250000 und grösser nach Satellitendaten herzustellen. Es wird offensichtlich eine Objecktauflösung von 5 bis 10m gefordert, um Bilddaten genügender Auflösung für die Herstellung topographischer Karten oder thematischer Studien des asiatischen Raums zu erzeugen.

Notes: There are two types of masked layers in seismic refraction work: the velocity reversal (low-velocity layer) and the hidden layer (insufficient velocity contrast or layer thickness). On the basis of an analytical formulation of the general case of a masked layer under an overburden of plane and parallel multiple refractors the two limiting cases are discussed: the solution resulting from an uncritical interpretation of the measured time-distance curve and the blind zone solution. Between these two limiting cases there is a variety of possible masked layer solutions. These no-blind zone solutions—as well as the blind zone solution itself—are formulated separately for the velocity inversion and the hidden layer case.For the evaluation of some no-blind zone solution a diagram is presented which can be used for any case of multiple refractors in the overburden of the masked layer. However, it is only for the three- and the four-layer case that a blind zone interpretation by use of diagrams is advisable. Such diagrams are presented together with the basic sets of formulae which contain as parameters only ratios of velocities and layer thicknesses. As the velocity of the masked layer is usually unknown the diagrams are principally constructed to show the dependence on the masked layer velocity. This is useful for estimation of the largest possible error.

Notes: Wave equation migration techniques have shown the limits of traditional stacking methods with data from tectonically complicated areas. An improved stack can be obtained utilizing the dip-moveout correction technique based on offset continuation. The properties and the limits of the algorithms used are summarized briefly.Several synthetic and real data examples are shown and compared with the results obtained using conventional processing in order to show the focusing effects and the strong improvement in signal-to-noise ratios, both at the stacked and migrated section level. The possibility of exploiting this technique to transform multiple coverage into increased spatial resolution is illustrated with examples.

Notes: The accuracy of short length digital linear filter operators can be substantially increased if the sampling interval as well as the abscissa shift are properly adjusted. This may be done by a trial and error process of adjustment of these parameters until the error made by the filter operator, applied to a suitably chosen test function, is smallest.As an illustration of the application of this method, 7-, 11- and 19-point filters for the calculation of Schlumberger apparent resistivity from a known resistivity transform are designed. Errors with the new 7-point filter are seen to be less than those with a 19-point filter of conventional design. The errors with the new 19-point filter are two to three orders of magnitude smaller than those made by the conventional 19-point filter.The new method should provide digital linear operators that allow significant improvements in accuracy for comparable computation efforts, or substantial reduction in computation for comparable accuracy of results, or something of both.

Notes: A well-known technique for the migration of normal-incidence two-way travel-time maps is extended to common-source-point travel-time data. The travel time and the travel-time gradient are used to compute the parameters defining the tangent plane of the reflecting interface. It is also shown how the curvature matrix of the received wavefront can be used to compute the curvature of the reflecting interface. The method is initially derived for common-source-point data and then extended to common-midpoint data.In a three-dimensional medium the wavefront curvature matrix is computed by solving a 2 × 2 symmetric matrix Riccati equation. In a two-dimensional medium and in a medium with constant velocity gradient, the wavefront curvature matrix is computed by solving a scalar Riccati equation and two linear equations. The migration procedures are also simplified.When the velocity function is unknown, the migration procedures cannot be used. An inverse modeling algorithm which simultaneously performs the migration and estimates the velocity function must then be applied. Two different inversion schemes are discussed briefly.

Notes: Time-domain-induced polarization (IP) laboratory measurements were performed on about 200 fine sediment samples with varying water content. The results permitted an analysis of IP properties of clays, loams, silts, and sands.Particular emphasis has been given to the analysis of the chargeability m as a function of lithotype and the water content.By analyzing decay curves, a new parameter was identified. It is a statistically specific characteristic of the lithotype and is independent of the water content. Therefore, it provides a diagnostic parameter for lithotype identification. In association with the values of chargeability and electrical resistivity, this parameter permits a reliable evaluation of water content and yields useful information about the porosity and permeability of the lithotype.

Notes: Approaches to the reduction of bias in the computation of the elements of the magnetotelluric impedance tensor have been proposed in the past by several authors. In this paper a clear distinction is made between random errors and bias errors. No effort is made to reduce either, but the emphasis is on their estimation. Both types of errors depend critically upon the polarization of the magnetic field. The random error increases with increasing noise-to-signal ratio in the electrical field, and it is rather insensitive to noise in the magnetic field. The bias error increases with increasing noise-to-signal ratio in the magnetic field. Expressions for random errors and maximum bias errors are developed and discussed using a single station set-up. Random errors with a reference station set-up are also calculated.

Notes: A numerical technique to compute the resistivity transform directly from the observed Wenner sounding data has been developed. In principle, the procedure is based on a decomposition method and consists of two steps: the first step determines a function that approximates the apparent resistivity data and the second step transforms this function into the corresponding kernel by an analytical operation. The proposed method is tested on some theoretical master curves. A high degree of precision is achieved with very little computer time. The applicability is shown on two field examples.

Notes: A suite of electrical, radiation, and mechanical borehole probes were run in a 76-mm-diameter borehole drilled to a slant depth of 380 m in leptite and granite. The hole is located in Precambrian bedrock in central Sweden where a site is dedicated to in-situ experiments pertaining to the disposal of radioactive wastes. The challenge to borehole logging methods for such site investigations is to resolve geological features and fluid flow parameters in geological sites which are initially chosen for their homogeneity, low porosity, and minimal fracturing. The Stripa borehole is characterized by high electrical resistivity values in the 20–100 kΩm range, by acoustic velocities around 5800 m s-1 (which is close to laboratory values on intact specimens), and by total porosity of around one volume percent. In this context, probe resolution was adequate to produce interpretable information on almost all of the logs.Two principal rock types were encountered in the hole: granite, of quartz monzonitic composition, and leptite. The granite and leptite intercepts are subdivided into units characterized by mafic mineral content, sulfide mineral content, and electrical and radiation properties. Iron-rich zones in the leptite are highly anomalous on the gamma-gamma and neutron logs; thin mafic zones in the granite can also be distinguished. Occurrences of a few percent pyrite are detected by the electrical, gamma-gamma, and neutron logs. Although overall porosity is quite low throughout the hole, analysis of the resistivity and neutron logs indicates the porosity increases by a few volume percent at fracture zones. The differential resistance and caliper probes detect borehole diameter roughness of less than 1 mm, helping to confirm acoustic waveform anomalies which are indicative of fracture zones. Compres-sional wave transit time and shear-wave interference patterns usually occur coincident with open fractures observed in core, the correlation being especially good at major fracture zones.

Notes: The use of conventional CDP-stacking in the processing of reflection data imposes restrictions on the horizontal and vertical resolution. Ideally, the final seismic section should consist only of short offset or, in practice, near-trace primary energy. Through the use of the iterative stacking algorithm, the signal-to-noise ratio on a single trace in the CDP-gather may be improved to an extent comparable to what occurs on a conventional stacked trace. By using this approach and treating the near-trace section after iterative stacking as the final section, the seismic resolution can be improved.

Notes: A Fredholm integral equation of the first kind with respect to the surface charge density is presented via a weighted residual formulation in the standard isotropic problem. The surface charge densities are numerically obtained by solving the well-known Fredholm integral equation of the second kind, where the new equation can be regarded as a constraint. The accuracy of boundary element solutions is examined in connection with the violation of such a constraint and the “modified reciprocal averaging techniques” are proposed.

Notes: The roots of the so-called shaly-sand problem in hydrocarbon evaluation lie in the effect of relatively fine-grained minerals upon measured electrical parameters of granular reservoirs. This influence manifests itself as an excess conductivity, over and above that due to the purely geometric effects of electrolyte distribution within the pore space. For formations with low shaliness, this excess conductivity is usually insignificant in typical oilfield situations.The influence of shaliness upon observed values of formation resistivity has been appraised by collating core-sample data from four different reservoirs. It has been demonstrated that during the course of electrical measurement under conditions of full electrolyte saturation, any given lithology can exhibit both negligible and highly significant shale effects depending upon the resistivity of the interstitial aqueous electrolyte. The effects of shaliness are also governed by the degree of water saturation.Because the manifestation of shaliness in electrical data is not a function of lithology alone, recourse is made to a more realistic concept of shale effects whereby a formation, or section of a formation, is classified as “effectively clean” or “effectively shaly” according to whether it obeys or defies, respectively, the fundamental empirical laws of Archie (1942). In particular, since an intrinsic formation factor can be obtained directly in fully-saturated effectively clean reservoirs, whereas only an apparent quantity may be recorded directly in fully-saturated effectively shaly reservoirs, the ratio of apparent to intrinsic formation factor serves as a useful conceptual indicator of shale effects, attaining the limiting value of unity only under effectively clean conditions.In the context of electrical measurement the terms “shaliness” and “shale effects” are evidently not synonymous and it is the latter which should be considered when selecting equations for the computation of water saturation. The implications for well-log analysis follow through formulated guidelines that describe the relative levels of shale effects in different zones of lithologically uniform reservoirs.

Notes: Shear waves can today be generated and observed, though not with the flexibility and the technical standard of compressional waves, and they can be identified in seismograms by various means. Their potential lies not so much in their lower velocity (corresponding—for the same frequency—to shorter wavelength and higher resolution) but in the fact that they probe the earth with stresses and strains that differ from those of compressional waves. Full utilization of the information potential of shear waves, therefore, requires combined use of P-and S-waves.Complications in the combined use of different wave types should be regarded as opportunities to obtain additional information. A typical example is the observation that the depth of one and the same interface estimated on the bases of P- and SH-reflections, respectively, can differ significantly. This discrepancy may be due to the anisotropy of a finely layered medium. Under favorable circumstances some of the parameters describing this anisotropy can be deduced from the different depth estimates and the curvature of the squared-offset/squared-time representation of the different reflections. Since in anisotropic media vertically polarized shear waves are significantly different from horizontally polarized ones, the combined observation of all three waves opens up additional possibilities.

Notes: An interpretation scheme for horizontal-loop EM measurements is presented for a permeable sphere model. The induced multipole moments are found to contribute significantly even at very low frequencies for a permeable conductor. The anomaly profiles are computed considering multipole excitation (up to 20) to study the effect of depth of burial and permeability of conductor. The anomaly half-width along with the inphase and quadrature anomaly amplitudes allow direct interpretation of the parameters of the sphere. The above scheme is suitable for results of Dighem II (coplanar configuration), Slingram and Max Min II measurements.

Notes: Wavenumber aliasing is the main limitation of conventional optimum least-squares linear moveout filters: it prevents adequate reject domain weighting for efficient coherent noise rejection. A general frequency domain multichannel filter design technique based on a one-to-one mapping method between two-dimensional (2D) space and one-dimensional (1D) space is presented. The 2D desired response is mapped to the 1D frequency axis after a suitable sorting of the coefficients. A min-max or Tchebycheff approximation to the desired response is obtained in the 1D frequency domain and mapped back to the 2D frequency domain. The algorithm is suitable for multiband 2D filter design. No aliasing damage is inherent in the linear moveout filters designed using this technique because the approximation is done in the frequency-wavenumber (f, k)-domain.Linear moveout filters designed by using the present coefficient mapping technique achieve better pass domain approximations than the corresponding conventional least-squares filters. Compatible reject domain approximations can be obtained from suitable mappings of the origin coefficient of the desired (f k)-response to the 1D frequency axis. The (fk)-responses of linear moveout filters designed by using the new technique show equi-ripple behavior. Synthetic and real data applications show that the present technique is superior to the optimum least-squares filters and straight stacking in recovering and enhancing the signal events with relatively high residual statics. Their outputs also show higher resolution than those of the optimum least-squares filters.

Notes: The inversion of seismic traces allows the estimation of reservoir porosity from an analysis of transit times derived from the pseudo-velocity logs. A four-step computational procedure is illustrated consisting of (i) inversion of seismic traces and calculation of interval velocities; (ii) accurate stratigraphic interpretation; (iii) determination of the petrophysical parameters for the porosity evaluation; (iv) analysis of the reliability of the results and final corrections.Both the possibilities and the limits of the method are discussed. One of the causes of error is the fact that impedances—and not velocities—are readily obtainable from seismic data. Moreover, the porosity due to fracturation contributes only slightly to velocity, while it often contributes most importantly to the permeability.Results are shown for two of the most significant reservoir types, i.e. carbonatic and clastic. Two cases belonging to the latter type will be examined. In the first case the primary porosity is dominant. The second case is very complex and both primary and secondary porosity are present.

Notes: To investigate the effect of permeability on the propagation of seismo-acoustic waves through marine sediments, a theoretical model based on Biot's equations is established which relates the compressional wave velocity measured at a fixed frequency to computed velocities at zero and infinite frequencies in terms of sediment porosity and permeability. The model is examined experimentally in a standard soil mechanics consolidation test (itself dependent, among other things, on sediment porosity and permeability) which has been modified to include measurements of compressional wave velocity at 1 MHz and shear-wave velocity at 5 kHz. This test allows the elastic modulus of the sediment frame to be assessed under different load conditions simultaneous with the velocity determinations.From a number of tests on different samples, five samples are chosen to typify the range of sediment sizes. The results show that the difference between the measured velocity at 1 MHz and the model-derived velocity at zero frequency increases with increasing particle size (from clays to fine sand), with decreasing porosity, and with increasing permeability. For sediments coarser than fine sand the simple model breaks down, possibly because of the dominance of scattering/diffraction effects at the high frequency of the experiment. Within this limitation the model seems satisfactory to offer a capability of predicting the permeability of a sea floor sediment to an order of magnitude by the in situ measurement of seismic velocities over a wide range of frequencies; the prediction process requires a good in situ determination of sediment porosity such as that offered by electrical formation factor measurements.

Notes: A methodology is outlined for using general complex compliance functions—which contain the constant-Q case and any kind of fluid saturant as special cases—to determine the nature of that saturant from seismic reflections observed on seismic records. The important attribute of such a formulation based upon complex compliance functions is that the presence of any kind of saturant is indicated by the values of the parameters contained within the complex compliance function, the latter having the same mathematical form for all rocks under any conditions as long as it remains a crystalline solid.

Notes: Gravity measurements were made in the Mailaram copper mines, Andhra Pradesh. The observations were distributed between the two shafts situated about 220 m apart and in the three levels up to a maximum depth of 100 m. Assuming a normal free-air gradient, average densities for the three layers were determined as 2.631, 2.604, and 2.823 g cm-3, respectively. Upon incorporating the weighted mean density values from measurements on samples, the free-air gradients were found to be 0.315 mGal m-1 for the second layer (i.e. between the first and second levels) and 0.2978 mGal m-1 for the third layer (i.e. between the second and third levels). The density variation map obtained from the gravity data, the deduced anomalies, the weighted mean density values from measurements on rock samples, and the varying free-air coefficients all suggest correspondence with the concentration of ore lodes.

Notes: Common midpoint data are now being collected with ever increasing source-receiver offsets. For wide aperture seismic data classical methods of interpretation fail, since velocity analyses and signal-to-noise enhancement methods based on hyperbolic traveltime curves are no longer appropriate. Therefore, the goals of increased velocity resolution and signal enhancement, which motivate the increase in offset, are not achieved. Approximate methods, involving higher order traveltime curves or extrapolations, have been developed for velocity analysis but these are ineffective in the presence of refracted arrivals, and lack a physical basis. These problems can be minimized by transforming the observational data to the domain of intercept or vertical delay time τ and horizontal ray parameter p. In this domain headwave refractions are collapsed into points and both near vertical and wide angle reflections can be analyzed simultaneously to derive velocity-depth information, even in the presence of velocity gradients or low velocity zones.

Notes: Using the concept of reciprocity and the known solution for the electric field of a vertical oscillating magnetic dipole source placed over a two-layered halfspace, an integral expression for the vertical magnetic field produced by a horizontal rectangular loop, carrying an oscillating current and placed on the surface of the same halfspace, is deduced. This expression is such that it can be evaluated by a combination of straightforward numerical integration and digital linear filter techniques. Displacement currents everywhere in space are neglected.Fields both inside and outside the loop are presented as sounding curves. A little over a decade of frequencies (100–3000 Hz) appears to be the minimum requirement and two decades (100–10000 Hz) are sufficient for most situations. Parametric and geometric sounding curves both show good resolution of subsurface layering. Phase shows better resolution than amplitude. Similarly, resolution is higher when conductivity increases with depth than when it decreases.Besides being useful for the electromagnetic depth sounding for layered earth structures, such computations can be applied to determine normal corrections to Turam observations, whenever a rectangular loop of finite size is used as a transmitter.

Notes: The surface resistivity method has been used to study a few exposed coal seams located in the northwestern part of the Raniganj Coalfield (belonging to the Damodar Valley Gondwana basin), India. Different electrode configurations, viz., Wenner, two-electrode and half-Schlumberger, have been used with different electrode spacings in horizontal profiling. The vertical distribution of resistivity has been studied using Schlumberger configuration at different locations along the profile.Laboratory studies of resistivity (at partial to full water saturation conditions) and porosity of different coal, shaly sandstone and sandstone samples of the survey area had already been carried out before the field survey. The results indicated the presence of a good resistivity contrast between the coal seams and the surrounding formations.The field results did not bear this out: the resistivity responses of the coal formations on the profiles and sounding curves are not as clear as one would expect for such high contents. The Wenner profiles show a broad resistivity anomaly over the coal seams. Two-electrode profiles are less noisy than Wenner profiles. Sharp peaks have been observed over the coal seams. The half-Schlumberger configuration seems the best: all the coal seams and their edges can be accurately outlined on the resistivity curves.Borehole data close to the profile have been used for correlating the field results. There is good agreement with vertical electrical soundings.

Notes: Source strength of an airgun array may be increased by:〈list style="custom"〉—utilizing higher pressure,—increasing total array volume,—employing more guns,—improving gun efficiency.One measure of gun efficiency is “specific source strength”, Pa*, defined as source strength per unit quantity of air used. Typical units are MPa m/l. Most developments are directed toward increasing gun pressure and/or gun volume to increase source strength of the array. These efforts require that more air compressors be installed onboard the ship. Consequently, a larger ship may be needed for the additional compressors, guns, and auxiliary equipment.A development program was initiated in 1976 to increase source strength of the array without using a larger ship. New guns were designed and built—one for 41.4 MPa and 7.37 liter (6000 p.s.i./450 in3) operation and another with 13.8 MPa and 4.92 liter (2000 p.s.i./300 in3) capability. Experiments were conducted with these new guns (and existing guns) over a range of pressures from 13.8 to 41.4 MPa (2000 to 6000 p.s.i.).Design of the new guns was aided by a mathematical model. The model relates physical dimensions of the airgun to acoustic pressure in the water. It consists of four nonlinear differential equations relating〈list style="custom"〉—shuttle motion,—bubble pressure,—chamber pressure,—bubble radius.The last equation is the “free-bubble-oscillation equation” and represents the ideal case of a pressurized bubble released instantaneously in water. The three other equations modify this ideal case; the four equations together model an airgun of the type manufactured by Bolt Associates, Inc.

Notes: Geometrical acoustic and wave theory lead to a second-order partial differential equation that links seismic sections with different offsets. In this equation a time-shift term appears that corresponds to normal moveout; a second term, dependent on offset and time only, corrects the moveout of dipping events.The zero-offset stacked section can thus be obtained by continuing the section with maximum offset towards zero, and stacking along the way the other common-offset sections.Without the correction for dip moveout, the spatial resolution of the section is noticeably impaired, thus limiting the advantages that could be obtained with expensive migration procedures. Trade-offs exist between multiplicity of coverage, spatial resolution, and signal-to-noise; in some cases the spatial resolution on the surface can be doubled and the aliasing noise averaged out.Velocity analyses carried out on data continued to zero offset show a better resolution and improved discrimination against multiples. For instance, sea-floor multiples always appear at water velocity, so that their removal is simplified.This offset continuation can be carried out either in the time-space domain or in the time-wave number domain. The methods are applied both to synthetic and real data.

Notes: Geoelectric soundings were used to investigate the younger fault activity (Middle and Early Pleistocene) of the southwestern fault margin area of the Central Graben (and extension of the Rhine Graben) in the southern Netherlands and northern Belgium. Some effects of this fault activity can still be observed in the present geomorphology and hydrography. Investigations have been concentrated on the uppermost 20–30 m, consisting mainly of Middle and Early Pleistocene fluviatile deposits. Marshy clays and loams are most conductive. They alternate with large sand beds and gravelly gully infillings, characterized by high specific resistivities. Problems in the interpretation of the geoelectric data, caused by the variability of the deposits, are partially solved by in-situ resistivity measurements (“mini-electric” measurements). In spite of the lithologic inhomogeneity, a few marker horizons allow the geologic structure to be determined. Some new tectonic boundaries have been traced. The movements along the faults during the Quaternary are very small (less than 10 m). Therefore, the network of soundings has to be very dense and a very intensive analysis of the soundings is necessary. The results also have a hydrogeologic significance. This project illustrates the possibility of locating small tectonic structures in relatively inhomogeneous deposits by a detailed and carefully designed geoelectric survey.

Notes: This paper is directed at modeling layered media. We extend the plane-wave normal-incidence state-space model developed by Mendel, Nahi and Chan in 1979, to the non-normal incidence case. To do this we introduce a shifting principle, a zero-offset wavefront, and zero-offset travel times for different layers. We also develop an algorithm for obtaining a synthetic line source reflection seismogram. In this algorithm non-normal incidence plane-wave seismograms are summed over a range of incident angles. The algorithm is based on a modified version of Sommerfield's (1896) theorem. Simulations of acoustic and elastic media are included which illustrate the applicability of our plane-wave and line source seismograms for both elastic and acoustic cases.

Notes: The principles for ray-tracing and wavefront curvature calculations in a three-dimensional medium are reviewed. A new derivation of the transformation of the wavefront curvature matrix at an interface between two inhomogeneous media is given. The derivation is based on a Taylor series expansion of the ray refraction equation at the interface between two inhomogeneous media, and only elementary geometric arguments are used. The wave-front curvature transformation at the interface is obtained by neglecting all terms in the direction of the surface normal.With proper definition of the variables, the derivation is also valid for a reflected wave-front. A simplified transformation rule is derived for a reflected wave of the same type as the incident wave.

Notes: A brief history of the development of the inverse problem in resistivity sounding is presented with the development of the equations governing the least-squares inverse. Five algorithms for finding the minimum of the least-square problem are described and their speed of convergence is compared on data from two planar earth models. Of the five algorithms studied, the ridge-regression algorithm required the fewest numbers of forward problem evaluations to reach a desired minimum.Solution space statistics, including (1) parameter-standard errors, (2) parameter correlation coefficients, (3) model parameter eigenvectors, and (4) data eigenvectors are discussed. The type of weighting applied to the data affects these statistical parameters. Weighting the data by taking log10 of the observed and calculated values is comparable to weighting by the inverse of a constant data error. The most reliable parameter standard errors are obtained by weighting by the inverse of observed data errors. All other solution statistics, such as dataparameter eigenvector pairs, have more physical significance when inverse data error weighting is used.

Notes: An analytical relationship for the autocorrelation function of an upsweep with high-frequency attenuation is used in the construction of synthetic seismograms.Field experiments were conducted in two areas to investigate the attenuation of upsweep where the near-surface materials were different. The results showed that the attenuation of high frequencies occurs at the source point depending on the near-surface lithology.The attenuation effect is usually neglected in the construction of the input wavelet of synthetic seismograms for Vibroseis data. In this study, the high-frequency attenuation of upsweep was considered in the construction of the input wavelet for the synthetic seismogram in an area where the Vibroseis technique was used. The synthetic seismogram generated in this manner had a better correlation with the Vibroseis section than that of corresponding synthetics using minimum-phase and the unattenuated autocorrelation wavelet of the upsweep.

Notes: The frequency domain response of a grounded linear circuit is deduced. The model of the ground is a homogeneous conducting half-space whose horizontal and vertical conductivities are specified. The source itself is formulated as an insulated current-carrying wire of finite length lying on the surface that is grounded at its end points. The tangential electric field expressions are used to calculate the induced voltage between two additional electrodes that are connected by an insulated wire. Particular attention is paid to the transient coupling response in the potential circuit when the source current is a step-function of time. Among other things, it is shown that the voltage response is a step-function only if the circuits are at right-angles and if the ground is effectively isotropic and nonpolarizable.

Notes: Application of dip-moveout attempts to correct pre-stack data in such a way that they stack correctly, however far the reflector geometry departs from an ideal plane horizontal interface. Even for plane dipping reflectors under a constant velocity overburden a common mid-point gather contains reflections distributed over a finite segment of the reflector. It is shown that under these conditions application of dip-moveout moves dipping energy on common offset gathers in such a manner that common mid-point gathers become true common depth point gathers.

Notes: Two methods for estimating the true formation temperature from well logs are compared. One method requires knowledge of the circulation time, whereas the other requires an estimate of the thermal diffusivity of the contents of the well. Both methods require three or more successive bottom-hole temperature measurements. Data from 157 wells have been analyzed, and the calculated formation temperatures from the two methods agree well. By an analysis of best matching of the two methods, it is found that 0.35 × 10-6 m2 s-1 is a good estimate for the thermal diffusivity of the well contents.

Notes: A new presentation for magnetotelluric data is proposed to be used together with the traditional representation. It is easily computed by applying the Bostick transform to all field data points. The result is an approximate depth-resistivity diagram which contains an indication of the data quality as well as the homogeneity of the subsurface. This resistivity section should be used in the interpretation of surveys in which accuracy is either not of prime concern or cannot be achieved. Since this mode of presentation is easily understood it can be used to familiarize earth scientists with the use of the MT method in exploration.

Notes: Seismic dynamic deconvolution is the mathematical basis on which a degree of unification in different prospecting methods is possible, relative to the parameter identification in horizontally stratified media. There is a basic structure which has some immediate applications to the inversion of resistivity data and possibly to other problems. For resistivity soundings there exists a key equation which is parallel to the energy conservation law in the theory of synthetic seismograms.

Notes: Magnetotelluric response is studied for an inhomogeneous medium having conductivity varying linearly with depth as σ(z) =σ1+αz. For a medium having conductivity increasing linearly with depth, the phase of the impedance approaches 60° at long periods and the apparent resistivity becomes log (ρa) = 2 log (1.36/α1/3) — 1/3 log (T'). The asymptote of log (ρa, T'→∞) when plotted against log (T') has a constant gradient —1/3 and has an intercept on the log (T') axis, which equals 6 log (1.36/α1/3). When a homogeneous layer with a moderate thickness overlies an inhomogeneous half-space, this layer does not affect the asymptote, but it affects the cut-off period and pushes this toward the long period direction. For a medium having conductivity decreasing linearly with depth, the impedance is equivalent to that of a Cagniard two-layer model; the intercept period related to the thickness is T'0=σ1(h2/2)2. Homogeneous multilayer approximations to an inhomogeneous layer are also investigated, and it is shown that the fit to the model variation depends on the number of layers and the layer parameters chosen.

Notes: Modern seismic processing techniques developed in recent years have provided the explorationist with more meaningful data than would have been predicted even by optimists. Correct migration of seismic data, relative amplitude preservation of reflections, and seismic trace inversion represent the necessary efforts to ensure that the best possible picture of in situ physical properties of the subsurface section is revealed. Moreover, compacted and over-pressured zones can be predicted from surface data prior to drilling a well through them.The basic tool for predicting overpressured zones from the surface is still the velocity analysis derived from good reflection data with few erratic multiples. The extraction of regional normal compaction trends from the seismic velocities allows one—where velocities deviate from the trend—to locate the top of overpressure. Moreover, the statistical behavior of the ratios of the sonic log vs pore pressure in existing boreholes enables one to convert the deviation from the trend of the seismic velocities into overpressure rates expected at the seismic reflection horizon.This paper presents a field case study showing how the knowledge of well site lithology together with the more detailed information extracted from inverted seismic data enables the prediction to match well conditions with high reliability.

Notes: Complete analytical expressions are developed for the first and second derivatives of the Newtonian gravitational potential in arbitrary directions due to the homogeneous revolutional compartment with a polygonal vertical section by applying the Gaussian divergence theorem in the cylindrical coordinate system. Elementary solutions presented can easily be translated into magnetic anomalies caused by a uniformly magnetized body.The divergence approach in the polar coordinate system is also described, and gravity attractions in the radial direction are presented in closed form associated with a homogeneous cap compartment. The explicit solutions are tested against well-known formulae for a cylinder, cone, infinite plate, and sphere.

Notes: The notion of nonvertical incidence plane wave seismograms and nonzero effect point source seismograms are sometimes confused with each other in the literature. This short note attempts to clarify this confusion which is observed in Douze (1979).

Notes: In the Gulf of Cadiz, a Tertiary basin became filled by clastic series during Miocene and Pliocene times. This terrigenous influx, derived from the Iberic Meseta in the north, is characterized by a sandy episode during the Tortonian and Messinian. The sand deposits were probably connected with uplift and major erosion of the Meseta during the sliding of the olistostrome, which occupied the south of the basin from late Helvetian to middle Tortonian.High resolution seismic techniques produced a good picture of the stratigraphy and of the depositional environment of the sands.A further study, using the amplitude of the reflections, inversion of seismic traces into acoustic impedance traces, and modeling, provides a remarkable example of the possibilities of seismic stratigraphy for depicting the lateral evolution of facies and localizing hydrocarbon occurrences. Out of seven exploratory wells based upon seismic information, six encountered gas-bearing sands with economic potential.

Notes: Problems of the resistivity sounding method may be reparametrized by means of Dar Zarrouk—or generalized D.Z.—transformations.A complete set of equivalent models is subdivided into three classes: initial models, intermediate models, and all other models. Equivalent reparametrized models are subdivided into three analogous classes. There is no one-to-one mapping of intermediate primary and intermediate reparametrized models. Non-trivial equivalent models are generated from initial models via intermediate reparametrized models. A method is developed to find equivalent models relative to the uncertainty of resistivity soundings.There are different forms of approximate equivalence laws, apart from those commonly accepted. The equivalence depends considerably on the “flatness” of the sounding curve. For mean flatness, the electric models exhibit the (hp±2)-equivalence rather than the (hp±1)-equivalence.

Notes: AEM (airborne electromagnetic) anomalies caused by massive sulphide conductors and superficial conductors can be recognized with a statistical method, as shown by an analysis of Input AEM data from Karnataka State. The weathering in the survey area is of tropical type. Parameters, such as various amplitude ratios and time parameters (inverse of decay rate) for exponential and power-law decay were analyzed for sulphide bodies, conducting soil, superficial conductors, and cultural conductors. Time parameters τ1 (exponential decay) is defined as ratio of time differences between the third and fourth channel to the logarithmic value of the relative amplitude of the two channels. Time parameter K1 (power-law decay) is defined as ratio of the difference of the logarithmic values of the delay times of the third and fourth channels to the logarithmic value of the relative amplitude of the two channels. The two parameters have been useful in recognizing sulphide conductors. Also the first channel Input amplitude and logarithmic plot of the transients appear to be helpful in conductor identification. Channel ratios seem to be the least effective parameters of conductor identification. In the area studied both power-law and the conventional exponential decay were found equally suitable for approximating Input AEM transients.

Notes: A seismic source consisting of a 700 kg weight that could be dropped vertically or projected down a ramp inclined at 45° to the vertical was tested as a source of P, SV and SH waves within a crystalline rock body at Chalk River, Ontario. The seismic energy was recorded by arrays of both horizontal and vertical-component geophones at distances between 30 and 600 m from the source, which was operated over glacial overburden varying in thickness from less than a meter to a few tens of meters. Seismic energy was more efficiently generated when the overburden thickness was at least several meters. The signals identified visually as S are generally true S, though some may be the converted wave PS. The SV amplitudes are generally larger than those of P, regardless of the type of shot, while the signal frequencies are roughly 60 Hz and 90 Hz, respectively. The horizontal-component seismograms for the inclined shots showed no evidence of SH polarization, and the SH amplitudes were only rarely enhanced relative to P and SV amplitudes on changing from vertical to inclined shots. These unexpected results are attributed to the combined effect of the high velocity and density contrasts and the irregularity of the boundary between the glacial overburden and the rock body.

Notes: Entrenchment and subsequent filling of a glacial valley have led to localization of iron-enriched ground water in the Silurian carbonate aquifer at the Vistron plant, Lima, Ohio. All production wells are open exclusively to the carbonate aquifer. Water from the saturated glacial deposits of relatively low permeability and high ironcontent is drawn laterally into the more permeable carbonate aquifer containing water with relatively low concentrations of iron. Wells closest to the saturated glacial deposits yield water with iron concentrations greater than 1.4 milligrams per liter. These wells, in the southwestern part of the site, apparently form a sink preventing iron-enriched water from migrating into other pumping wells in the central and northern parts of the area.Iron concentrations have not increased in individual wells between 1971 and 1981. They are not expected to increase in the central and northern wells as long as about 20% of the plant's ground water is withdrawn from at least two of the southwestern wells.

Notes: Nitrate-N concentrations during the Summers of 1980 and 1981 exceeded 10 mg/1 in 68 of the 164 ground-water samples collected from a shallow water-table aquifer underlying a sand plain near Alliston, Ontario. Three extensive zones of nitrate contamination were associated with major potato-growing areas on the sand plain. Nitrate concentrations were positively correlated with both the percentage area of heavily fertilized crops (potatoes, corn, sod and asparagus) and nitrogen fertilizer application rates in the vicinity of ground-water sampling sites. Chloride levels in ground water exhibited a positive association with KCl fertilizer application rates. Ground water under potato fields had relatively consistent C1/NO3-N ratios despite considerable variations in nitrate-N and chloride concentrations. Ground water with a low nitrate content was found beneath forest and permanent pasture. These data suggest that the use of commercial nitrogen fertilizer is a major source of nitrate in the aquifer.

Notes: Abstract. A simple and convenient method for calculating two-dimensional dispersion on a TI 58 or 59 programmable calculator utilizing the normal distribution program available in the Solid State Software Master Library module as a subroutine is presented with an example.

Notes: Three simplistic models are developed for evaluating the transport of organic pollutants through soil to ground water. The models consider mobility and first-order degradation. The first calculates linear sorption/desorption of the pollutant and first-order degradation without considering dispersion. The second is similar to the first but also considers dispersion. The third considers nonlinear sorption following a Freundlich equation and first-order degradation but does not consider dispersion. The models are compared to field data for the pesticides aldicarb and DDT. The models projected a lower mobility for DDT than was observed in the field.

Notes: Darcy's equation and unsaturated-flow theory are used to calculate flow through earth linings that are placed in ponds or channels to reduce seepage. The procedure utilizes the relation between unsaturated hydraulic conductivity and pressure head of the underlying material as calculated from the measured saturated hydraulic conductivity and the relation between water content and pressure head. The method enables the selection of the liner (thickness and hydraulic conductivity) that will keep seepage below a certain maximum limit. Since the hydraulic conductivity of clays is affected by the cationic composition and the salt concentration of the soil solution, the chemical composition of the liquid moving through the liner must be taken into account. This applies also to other chemicals, including solvents and other organic compounds that may be in the water. Travel times of water from the surface impoundment to the underlying ground water are calculated from the seepage rate and the corresponding water content in the vadose zone. Accumulation of solids (mine tailings, for example) can further reduce the seepage from the pond. Proper design of waste-water ponds also requires analysis of the response of the underlying ground water (mound buildup) and the movement of pollutants in the vadose zone and aquifer. In view of the high costs of earth liners, prediction of the seepage is necessary to make sure that the selected lining material, the thickness of the liner itself, and the method of construction will produce the desired results.

Notes: A multidimensional, finite-difference model for ground-water flow and heat transport is used to analyze the thermal energy storage experiment conducted by Auburn University in Mobile, Alabama. The experiment consisted of three stages–namely, injection, storage and recovery occurring for 80, 51 and 41 days, respectively. This application demonstrates the validation evidence that the model adequately and accurately simulates the field experiment. The numerical model includes the effects of: hydraulic anisotropy, thermal convection and conduction, and heat loss to the adjacent confining strata. Observed aquifer isotherms at the end of each stage are compared with predicted values on a cylindrical grid situated about the well. The degree of vertical discretization used in the model is shown to impact the predicted temperature profiles at each stage, but has minimal effect on the recovery water temperature.

Notes: The author has used the concept of depression cone volume to derive type curve equations for large-diameter wells in aquifers of finite extent. Comparisons with the corresponding type curves in an infinitely extensive aquifer have been given. These solutions are based on joint exploitation of the ground-water movement equation (Darcy's law) and the continuity equation for large-diameter wells.

Notes: The majority of well fields in the river region in the Netherlands, where anaerobic water is withdrawn from the shallow aquifer, have problems with well clogging. In order to test the supposition that sulfate-reducing bacteria play a role in this clogging process, sulfate-reducing bacteria in water from wells on well fields with and without the occurrence of clogging were enumerated. In water withdrawn from nonclogging wells, the Most Probable Number of sulfate-reducing bacteria averaged 5 per 100 ml, whereas in wells subject to clogging, the number averaged 25 per 100 ml. A statistical analysis by Wilcoxon's order test confirmed that a significant difference exists between the numbers of sulfate-reducing bacteria in clogging and nonclogging wells. The mechanism whereby sulfate-reducing bacteria contribute towards this type of well clogging is as yet unclear.

Notes: The Kansas Geological Survey is pursuing an effort to automate some of the more common methods of aquifer pumping-test analysis. This paper discusses the results of work done on the leaky artesian aquifer as defined by Hantush and Jacob (1955). The paper covers the basic theory of the aquifer type, the numerical solution of the leaky artesian-well function, and the methodology of achieving the “best fit” parameters in the least squares' sense. Several data sets are used to demonstrate the applicability of the proposed technique. These examples indicate the generally satisfactory results produced by the automated analysis documented here.The algorithm has good convergence properties. Initial estimates for the aquifer parameters may vary by about three orders of magnitude above or below the correct values. For typical data sets the rms fitting error should be less than a few tenths of a foot. If this is not the case, one is probably not dealing with a simple leaky aquifer. This method of pumping-test analysis does not eliminate the role of an experienced hydrologist to define the local hydrogeology and aquifer type. However, once the decision is made as to which aquifer configuration is being observed, this program will, in a quick and unbiased fashion, give an accurate assessment of the leaky-aquifer parameters within the limits of the theoretical approximations and the data quality.

Notes: The software for simulation of the three-dimensional ground-water flow, developed by Trescott (1975), is executable only on bigger machines with large memory capacities. To popularize the sophisticated software for easily accessible small machines, the organization is restructured and various programming facilities are availed. The new design is tested with an earlier quoted example. The program requires about 11K bytes of memory as against 72K bytes of the original program on IBM/370. The methodology borne out of the present work to implement large programs for small computers is presented along with the modified code.

Notes: In the design of ground-water heat pump systems the quantitative analysis of heat transport and heat storage is of great interest. A unified finite element approach to the transient nonlinear heat transport and heat storage problem is presented. The theory presented includes nonlinear physical properties and boundary conditions, coupled conductive and convective heat flow, freezing (phase change), and time-dependent heat input and output. The given theory is applied to two problems. The first problem is a transient study during a few years of heat storage in saturated clay. A pure conductive theory is used due to the small amount of convective water flow in clay. The second problem deals with transient heat transport in an aquifer according to coupled conduction-convection theory on a one-dimensional model problem.

Notes: The principal water-bearing units (in descending stratigraphie order) in Pipeline Canyon are the alluvium, the Dilco Coal Member of the Crevasse Canyon Formation, the Torrivio Sandstone Member of the Gallup Formation, and the Upper and Lower Gallup Sandstones. Presently, the alluvium is recharged by a perennial, southward-flowing stream sustained by mine dewatering discharge in addition to natural precipitation and runoff. Localized infiltration has created ground-water mounds in the alluvium. Artesian conditions may exist in sandstone units which are adjacent and hydraulically connected to the alluvium in these areas. Basement faults have produced the Fort Wingate and Pipeline Canyon lineaments which intersect in the area of investigation. Fault-related fracturing of the sandstones allows for significant ground-water recharge via the overlying alluvium. Flow in the sedimentary formations is generally to the northeast, while the alluvial system flows to the southwest. The natural water quality for the aquifers in the area is dominated by sodium and sulfate ions. However, complete characterization of the natural water quality is complex because of ground-water contamination by acidic tailings fluids from a local uranium mill.

Notes: This paper describes a systematic and straightforward method for the estimation of velocity components in three dimensions from hydraulic head data. Groups of four measurement points are connected to form tetrahedrons, and a linear interpolation scheme is used to obtain a head gradient estimate for each tetrahedron. Application of Darcy's law then yields the desired velocity component values. A sample calculation and comparison between this method and a two-dimensional approach are also included.

Notes: Two earthfill sections of Wallace Dam on the Oconee River near Eatonton, Georgia were constructed with vertical drainage filters (chimney drains) in a clay fill zone. In order to evaluate the performance of the filter in the west dike of Wallace Dam, the finite element Galerkin method was utilized in formulating a numerical model to study the steady-state saturated-unsaturated seepage characteristics through the earth dam. The resulting model is applied to Station 58+00 of the west dike of Wallace Dam. Numerical results for the four cases analyzed in this study describe the location of the zero pressure isobar and total hydraulic head values ranging from 425 feet (130 m) to 365 feet (111 m). Model results indicate a maximum seepage velocity of 2.62 feet per day (0.80 m/day) using a saturated horizontal hydraulic conductivity of 0.283 foot per day (8.64 cm/day). A maximum seepage rate is calculated to be 0.266 cubic foot per second (7.52 × 10−3 m3/sec). Analysis of the hydrostatic uplift forces along the base of the dam indicates an average pressure head reduction of 51 feet (16 m) from the upstream to the downstream side of the dam. Comparison of local seepage velocities to the critical seepage velocity upstream of the filter and inside the filter indicate a factor of safety against piping (a factor of safety against a quick condition arising in the soil) ranging from 0.3 to 3.7.

Notes: The ground-water observation well network in many parts of Kansas has been developed and expanded through the years without serious attempt to determine the adequacy of the network for any specified purpose or to assess its cost effectiveness. This study was undertaken to examine the existing well network in northwest Kansas and to determine the arrangement that offers the most satisfactory accuracy for the purpose of monitoring it. To achieve this goal, we have employed the theory of regionalized variables to estimate the amount of spatial variability of the water table. The error analysis produced by universal kriging indicates that a significant reduction in the number of wells could be achieved by employing a regular 4-mile (6.4-km) network, without affecting the present level of accuracy. It also indicates that it is not practical to reduce the estimation error in the water-table surface uniformly throughout the region because to do so would increase the cost of monitoring wells drastically. For example, reducing the presently existing error by 50 percent throughout the area would require 16 times more wells than the currently existing well network.

Notes: In river basins where aquifers are closely interrelated with a stream, large-scale ground-water development can lower water tables near the stream and diminish stream flows. Junior surface right holders are adversely affected. A digital computer simulation of the hydrologic-economic system on the lower South Platte River in Colorado is employed to study economic impacts of two water management policies. Open access management is found to yield high income benefits but imposes substantial costs on surface-water users when water supplies are limited. Incorporation of ground water into the appropriation system can help avoid the losses to surface-water users, but greatly reduces the income of ground-water users.

Notes: Basalts are a major source of ground water throughout the Columbia River Plateau region of the Pacific Northwest. Development and management of ground water in these basalts are complicated by the spatial variability of the hydrologic characteristics of the deep, stratified lava flows, but new irrigation developments and municipal and industrial water needs are placing increasingly larger demands on the ground-water resource. Water management decisions are aided by individual basin studies that contribute to greater understanding of the regional ground-water system. A distributed-system multiple-storage model for the Deschutes River Basin, Oregon, reveals the magnitude of spatial differences in ground-water recharge, storage, and discharge for this watershed. Input-output analysis elucidates the functional characteristics of the basin groundwater system, and it identifies the presence and magnitude of interbasin linkages in the ground-water system. Implementation of ground-water development strategies based on storage and transmission characteristics simulated by the model demonstrates that basin or regional perspectives are necessary to fully utilize ground-water storage in basalts.

Notes: Bioherms found in the Lower Permian Tyrrellfjellet Member of Central Spitsbergen contain bindstones produced by the ?hydrozoan taxon Palaeoaplysina. We briefly describe the bioherms’stratigraphical and environmental setting and note their development on hardground surfaces. The association of these structures with bituminous limestones is interesting in view of the petroleum resources of similar buildups in the USSR.

Notes: Aeolian deposits are described in terms of their areal distribution, sedimentary and pedological characteristics. These deposits, which have accumulated at the valley margins, result from deflation of fluvian sediments deposited in the valley bottom. Fine horizontal lamination in the upper horizons is tentatively attributed to winnowing of primary depositional units following partial cementation by salt precipitates. This stratification breaks down at depth and is replaced by a gleyed horizon, resulting from drainage impedance by permafrost. Proximal aeolian accumulations on flu-vially inactive areas of the valley bottom may typify many other valleys in Spitsbergen.

Notes: Occurrence of Salterella rugosa Billings on Nordaustlandet, Svalbard, is reported for the first time. The rocks that produced these fossils belong to the upper part of the Kapp Sparre Formation. This unit can be correlated with the Oslobreen Dolomite in Ny Friesland (Spitsbergen, Svalbard) from which the same species has been reported earlier; these fossils have not been illustrated previously. In both areas many specimens are poorly preserved, and it is becoming increasingly evident that the difference between Salterella and Volborthella simply may be the result of loss of the outer shell in the Volborthella.

Notes: Foraminiferal faunas were analysed from 22 cores (length 30-165 cm), three grab samples and one dredge sample taken from Spitsbergenbanken and the trenches Storfjordrenna and Bjørnøyrenna. The distribution pattern of selected taxa and eight assemblages are discussed in relation to water depth, currents and sedimentation. The assemblages found on Spitsbergenbanken where currents are particularly active, are clearly different from the assemblages occurring in the more quiet environments of the trenches.In the Quaternary evolution of the area three successive stages are distinguished: glacial shelf environment of the Weichselian; semiglacial, or transitional environment of the Late Weichselian and Early Holocene; open, or deglaciated shelf environment of the Holocene. Minor faunal changes in Storfjordrenna are explained by changes in the current pattern.In some assemblages marked postmortem changes have taken place which include solution of calcareous shells, breakage and abrasion by currents, and transport of shells out of the high energy bank area with assumed redeposition in deeper water.

Notes: The pre-Carboniferous miogeosynclinal sedimentary and volcanic sequences of south-central Prins Karls Forland belong to a suite of rocks known as the Forland Complex. During early Palaeozoic diastrophism the Complex of Prins Karls Forland experienced predeformational, low grade metamorphism. A knowledge of the stratigraphy and structure of the region, together with an analysis of metamorphic mineral assemblages and reactions enable the conditions of metamorphism to be determined within certain limits. It is found that temperatures during metamorphism were between 400°C and 500°C, and pressures were about 6 kb. These conditions provide an estimate of palaeogeothermal gradient, consistent with observed mineralogies, of 18 to 21°C.km−1. In addition, compositional variations within mineral phases are found to be largely dependent upon bulk-rock chemistry, and are of limited use in determining grade changes.

Notes: The Orustdalen Formation of Brøggerhalvaya is well exposed on two coastal sections and is 200–250 m thick. It consists dominantly of siliceous conglomerates and sandstones with subordinate shales and impure coal. One horizon low in the sequence yielded spores of uppermost Dinantian/lowermost Namurian age. Three sedimentary fades have been recognised. The fluvial channel facies (interbedded conglomerates and cross-stratified sandstones) are interpreted as braided stream deposits with flow directions to south and west. A shale facies, sometimes with drifted plant remains, is interpreted as overbank in origin. Highly quartzose, medium-to-coarse-grained, cross-stratified sandstones make up the reworked facies where palaeocurrents are bimodal or indicate movement to the NE. The Formation is thought to represent sediments shed from an early Carboniferous fault scarp eroding’a siliceous Lower Palaeozoic source terrain. The proximity of a nearby coastline led to the construction of fan deltas where fluvial deposits were reworked by waves and/or tides. Early diagenetic events include the local development of pyrite, kaolinite and calcite cements. Reddening’probably occurred during Middle Carboniferous times beneath an exposed land surface. Following burial, pressure solution and quartz cementation eliminated porosity at depths greater than 1000 m.