On the computation of the spherical harmonic terms needed during the numerical integration of the orbital motion of an artificial satellite

A method is presented for the accurate and efficient computation of the forces and their first derivatives arising from any number of zonal and tesseral terms in the Earth's gravitational potential. The basic formulae are recurrence relations between some solid spherical harmonics,V _n,m, associated with the standard polynomial ones.     

Visible and infrared radiometer on Seasat-1

The visible and infrared radiometer (VIRR) is a scanner for the collection of digital data in the visible and thermal-infrared portions of the spectrum over broad swaths of the earth's surface [1]. Some of the essential electronic and mechanical details of the system are described; and the results of an engineering assessment of its operation in orbit are given. The set of algorithms used to transform the raw VIRR data to calibrated geophysical measurements is discussed. Some of the results of a preliminary geophysical evaluation are presented.     

Molecular line mapping of the giant molecular cloud associated with RCW 106 – II. Column density and dynamical state of the clumps

We present a fully sampled C¹⁸O (1–0) map towards the southern giant molecular cloud (GMC) associated with the H  ii region RCW 106, and use it in combination with previous ¹³CO (1–0) mapping to estimate the gas column density as a function of position and velocity. We find localized regions of significant ¹³CO optical depth in the northern part of the cloud, with several of the high-opacity clouds in this region likely associated with a limb-brightened shell around the H  ii region G333.6−0.2. Optical depth corrections broaden the distribution of column densities in the cloud, yielding a lognormal distribution as predicted by simulations of turbulence. Decomposing the ¹³CO and C¹⁸O data cubes into clumps, we find relatively weak correlations between size and linewidth, and a more sensitive dependence of luminosity on size than would be predicted by a constant average column density. The clump mass spectrum has a slope near −1.7, consistent with previous studies. The most massive clumps appear to have gravitational binding energies well in excess of virial equilibrium; we discuss possible explanations, which include magnetic support and neglect of time-varying surface terms in the virial theorem. Unlike molecular clouds as a whole, the clumps within the RCW 106 GMC, while elongated, appear to show random orientations with respect to the Galactic plane.     

Analytical model for screening potential CO<Subscript>2</Subscript> repositories

Assessing potential repositories for geologic sequestration of carbon dioxide using numerical models can be complicated, costly, and time-consuming, especially when faced with the challenge of selecting a repository from a multitude of potential repositories. This paper presents a set of simple analytical equations (model), based on the work of previous researchers, that could be used to evaluate the suitability of candidate repositories for subsurface sequestration of carbon dioxide. We considered the injection of carbon dioxide at a constant rate into a confined saline aquifer via a fully perforated vertical injection well. The validity of the analytical model was assessed via comparison with the TOUGH2 numerical model. The metrics used in comparing the two models include (1) spatial variations in formation pressure and (2) vertically integrated brine saturation profile. The analytical model and TOUGH2 show excellent agreement in their results when similar input conditions and assumptions are applied in both. The analytical model neglects capillary pressure and the pressure dependence of fluid properties. However, simulations in TOUGH2 indicate that little error is introduced by these simplifications. Sensitivity studies indicate that the agreement between the analytical model and TOUGH2 depends strongly on (1) the residual brine saturation, (2) the difference in density between carbon dioxide and resident brine (buoyancy), and (3) the relationship between relative permeability and brine saturation. The results achieved suggest that the analytical model is valid when the relationship between relative permeability and brine saturation is linear or quasi-linear and when the irreducible saturation of brine is zero or very small.     

A two-dimensional analytical solution of groundwater responses to tidal loading in an estuary and ocean

Previous studies on tidal dynamics of coastal aquifers have focussed on the inland propagation of oceanic tides in the cross-shore direction, a configuration that is essentially one-dimensional. Aquifers at natural coasts can also be influenced by tidal waves in nearby estuaries, resulting in a more complex behaviour of head fluctuations in the aquifers. We present an analytical solution to the two-dimensional depth-averaged groundwater flow equation for a semi-infinite aquifer subject to oscillating head conditions at the boundaries. The solution describes the tidal dynamics of a coastal aquifer that is adjacent to a cross-shore estuary. Both the effects of oceanic and estuarine tides on the aquifer are included in the solution. The analytical prediction of the head fluctuations is verified by comparison with numerical solutions computed using a standard finite-difference method. An essential feature of the present analytical solution is the interaction between the cross- and along-shore tidal waves in the aquifer area near the estuary’s entry. As the distance from the estuary or coastline increases, the wave interaction is weakened and the aquifer response is reduced, respectively, to the one-dimensional solution for oceanic tides or the solution of Sun (Sun H. A two-dimensional analytical solution of groundwater response to tidal loading in an estuary, Water Resour Res 1997;33:1429–35) for two-dimensional non-interacting tidal waves.     

Rapid magmatic processes accompany arc–continent collision: the Western Bismarck arc, Papua New Guinea

New U–Th–Ra, major and trace element, and Sr–Nd–Pb isotope data are presented for young lavas from the New Britain and Western Bismarck arcs in Papua New Guinea. New Britain is an oceanic arc, whereas the latter is the site of an arc–continent collision. Building on a recent study of the Manus Basin, contrasts between the two arcs are used to evaluate the processes and timescales of magma generation accompanying arc–continent collision and possible slab detachment. All three suites share many attributes characteristic of arc lavas that can be ascribed to the addition of a regionally uniform subduction component derived from the subducting altered oceanic crust and sediment followed by dynamic melting of the modified mantle. However, the Western Bismarck arc lavas diverge from the Pb isotope mixing array formed by the New Britain and the Manus Basin lavas toward elevated ²⁰⁸Pb/²⁰⁴Pb. We interpret this to reflect a second and subsequent addition of sediment melt at crustal depth during collision. ²³⁸U and ²²⁶Ra excesses are preserved in all of the lavas and are greatest in the Western Bismarck arc. High-Mg andesites with high Sr/Y ratios in the westernmost arc are attributed to recent shallow mantle flux melting at the slab edge. Data for two historical rhyolites are also presented. Although these rhyolites formed in quite different tectonic settings and display different geochemical and isotopic compositions, both formed from mafic parents within millennia.     

Application of airborne LiDAR data and airborne multispectral imagery to structural mapping of the upper section of the Troodos ophiolite, Cyprus

Structural maps are traditionally produced by mapping features such as faults, folds, fabrics, fractures and joints in the field. However, large map areas and the spatially limited ground perspective of the field geologist can potentially increase the likelihood that not all structural features will be identified within a given area. The ability to recognise and map both local and regional structural features using high-resolution remote sensing data provides an opportunity to complement field-based mapping to help generate more comprehensive structural maps. Nonetheless, vegetation cover can adversely affect the extraction of structural information from remotely sensed data as it can mask the appearance of subtle spectral and geomorphological features that correspond to geological structures. This study investigates the utility of airborne Light Detection And Ranging (LiDAR) data and airborne multispectral imagery for detailed structural mapping in vegetated ophiolitic rocks and sedimentary cover of a section of the northern Troodos ophiolite, Cyprus. Visual enhancement techniques were applied to a 4-m airborne LiDAR digital terrain model and 4-m airborne multispectral imagery to assist the generation of structural lineament maps. Despite widespread vegetation cover, dykes and faults were recognisable as lineaments in both data sets, and the predominant strike trends of lineaments in all resulting maps were found to be in agreement with field-based structural data. Interestingly, prior to fieldwork, most lineaments were assumed to be faults, but were ground-verified as dykes instead, emphasising the importance of ground-truthing. Dyke and fault trends documented in this study define a pervasive structural fabric in the upper Troodos ophiolite that reflects the original sea-floor spreading history in the Larnaca graben. This structural fabric has not previously been observed in such detail and is likely to be continuous in adjacent regions under sedimentary cover. This information may be useful to future exploration efforts in the region focused on identification of structurally controlled mineral and groundwater resources. Overall, our case study highlights the efficacy of airborne LiDAR data and airborne multispectral imagery for extracting detailed and accurate structural information in hard-rock terrain to help complement field-based mapping.     

Surface geostrophic currents across the Antarctic circumpolar current in Drake Passage from 1992 to 2004

The Southern Ocean plays an important role in the global overturning circulation as a significant proportion of deep water is converted into intermediate and deeper water masses in this region. Recently, a secular trend has been reported in wind stress around the Southern Ocean and it is thought theoretically that the strength of the ACC is closely related to wind stress, so one consequence should be a corresponding increase in ACC transport and hence changes in the rate of the global overturning. There are no long-term data sets of ACC transport and so we must examine other data that may also respond to changing wind stress. Here we calculate surface currents in Drake Passage every seven days over 11.25 years from 1992 to 2004. We combine surface velocity anomalies calculated from satellite altimeter sea surface heights with measured surface currents. Since 1992, the UK has regularly occupied WOCE hydrographic section SR1b across the ACC in Drake Passage. From seven hydrographic sections surface currents are estimated by referencing relative geostrophic velocities from CTD sections with current measurements made by shipboard and lowered acoustic Doppler current profilers. Combining the seven estimates of surface currents with the altimeter data reduces bias in the estimates of average currents over time through Drake Passage and we show that surface current anomalies estimated by satellite and in situ observations are in good agreement. The strongest surface currents are found in the Subantarctic and Polar Fronts with average speeds of 50 cm/s and 35 cm/s, respectively and are inversely correlated, so that maximum westward flow in one corresponds to minimum westward flow in the other. The average cross-sectional weighted surface velocity from 1992 to 2004 is 16.7 ± 0.2 cm/s. A spectral analysis of the average surface current has only weakly increasing energy at higher frequencies and there is no dominant mode of variability. The standard deviation of the seven day currents is 0.68 cm/s and a running 12 month average has only a slightly smaller standard deviation of 0.52 ± 0.16 cm/s. The southern annular mode (SAM) measures the circumpolar average of wind stress and like the surface currents its spectrum has slightly increased energy at frequencies greater than 1 cpy. A cospectral analysis of these, averaging cospectra of five slightly overlapping 36 month segments improve statistical reliability, suggests that there is coherence between them at 1 cpy with the currents leading changes in the Southern annular mode. We conclude that the SAM and average Drake Passage surface currents are weakly correlated with no dominant co-varying modes, and hence predicting Southern Ocean transport variability from the SAM is not likely to give significant results and that secular trends in surface currents are likely to be masked by weekly and interannual variability.