Efficient GOCE satellite gravity field recovery based on least-squares using QR decomposition

We develop and apply an efficient strategy for Earth gravity field recovery from satellite gravity gradiometry data. Our approach is based upon the Paige-Saunders iterative least-squares method using QR decomposition (LSQR). We modify the original algorithm for space-geodetic applications: firstly, we investigate how convergence can be accelerated by means of both subspace and block-diagonal preconditioning. The efficiency of the latter dominates if the design matrix exhibits block-dominant structure. Secondly, we address Tikhonov-Phillips regularization in general. Thirdly, we demonstrate an effective implementation of the algorithm in a high-performance computing environment. In this context, an important issue is to avoid the twofold computation of the design matrix in each iteration. The computational platform is a 64-processor shared-memory supercomputer. The runtime results prove the successful parallelization of the LSQR solver. The numerical examples are chosen in view of the forthcoming satellite mission GOCE (Gravity field and steady-state Ocean Circulation Explorer). The closed-loop scenario covers 1 month of simulated data with 5 s sampling. We focus exclusively on the analysis of radial components of satellite accelerations and gravity gradients. Our extensions to the basic algorithm enable the method to be competitive with well-established inversion strategies in satellite geodesy, such as conjugate gradient methods or the brute-force approach. In its current development stage, the LSQR method appears ready to deal with real-data applications.     

Field scale heterogeneity of redox conditions in till-upscaling to a catchment nitrate model

Point scale studies in different settings of glacial geology show a large local variation of redox conditions. There is a need to develop an upscaling methodology for catchment scale models. This paper describes a study of field-scale heterogeneity of redox-interfaces in a till aquitard within an area of 600?×?600 m. The results showed significant variation of the depths to the redox-interface and thickness of the aquitard. Nitrate was present above the redox-interface but reduced to non-detectable levels a few metres below the interface. An upscaling approach for an area of 92 km² is proposed. Two models are proposed to predict the depth to the redox-interface in the aquitard and the resulting nitrate recharge concentrations to an underlying aquifer. The first model assumes that the depth to the redox-interface reflects the hydraulic head in the aquitard, and the second model assumes that the depth of the redox-interface is randomly distributed according to a log-normal probability distribution function. The upscaling approach using the random redox model estimated recharge concentrations comparable to the observed concentration in the underlying aquifer. The presented upscaling approach is applicable in distributed catchment models where sub-grid variability cannot be represented by the large grids.     

Interstellar neutral hydrogen filaments at high galactic latitudes and the bennett pinch

Observed properties of interstellar neutral hydrogen filaments suggest the presence of the Bennett pinch as described by the Carlqvist relationship with rotation around the filament axes included. A brief summary is first given of three ways in which a filament model for interstellar cloud structure was tested. Preliminary results from highresolution HI mapping of gas and dust in an apparent HI cloud indicate that the neutral gas and dust within and around its boundary is itself highly filamentary. An attempt to detect magnetic fields in this and similar features using the Zeeman effect technique at the 21-cm wavelength of interstellar neutral hydrogen set upper limits of a fewµG. In contrast, the strength of the toroidal magnetic field expected from the examination of the Carlqvist relationship is of order 5µG, which would be produced by a current of 1.4 · 10¹³ A. Zeeman effect technology is at present not able to detect toroidal magnetic fields of this order at the edge of barely resolved HI filaments. Nevertheless, currently available high-resolution HI data suggest that interstellar filament physics should take into account the role of currents and pinches for creating and stabilizing the structures.     

Infrared albedos and rotation curves of the Galilean satellites

Infrared (1.5–5 μm) albedos and rotation curves of the Galilean satellites have been obtained. The data suggest that the rotational variation in the infrared is less than ±10% for all four satellites. While no conclusion about rotational variation could be reached for Io, the 1.57 μm data for the outer three satellites marginally suggest phase correlation with the visual variation. The geometric albedos obtained are in general agreement with earlier results. For Io, the absorption feature near 1.5 μm found by Pilcher et al. (1972) is confirmed, thus contradicting the flat spectrum measured by Fink et al. (1973). Io and Ganymede were observed in the 1.57 μm bandpass as they reappeared from eclipse. The curve for Io shows a slight (<10%) overshoot similar to those sometimes reported for visual measurements. This result is based on a single reappearance, and is extremely tentative.     

A distinctive type of ascending prominence — ‘fountain’

Cinematographic observations of solar prominences made at Mauna Loa during the past couple of years suggest that there is a well-defined sub-class of ascending prominences characterized by closed-system transference of chromospheric material along an arch or loop (up one leg and down the other); meanwhile the entire prominence envelope steadily rises upward and expands through the corona. We denote these prominences as fountains. Several examples are described. Fountains appear to be well contained by coronal magnetic fields. Their total kinetic energy is in the order of 10³⁰ erg but dissipation is typically quite slow (over time periods like 100 min) so that the correlative disturbances (radio bursts, coronal transients, chromopsheric brightenings, etc.) are generally unspectacular or non-existent.This research was started when E. T.-H. was still on the staff of the High Altitude Observatory.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

The sources of material comprising a mass ejection coronal transient

The origin of the material which is ejected during a white light coronal transient has not been determined heretofore. Study of a disturbance on 26 and 27 August 1973, during which a slowly ascending prominence and a more rapid accompanying coronal transient were simultaneously observed, helps to resolve this question. Prominence images obtained in Hα 6563 Å and in He II 304 Å are nearly identical. The mass ejection transient observed in white light (3700–7000 Å) appeared to be a loop about 1 R_⊙ higher than the top of the ascending prominence; it accelerated away from the prominence below it. These observations imply: (1) the bulk of the ejected material did not originate in the ascending prominence; (2) therefore, most of the material must have come from the low corona above the prominence, (and was at coronal temperatures during its outward passage); and (3) the total event - ascending prominence accompanied by coronal mass ejection - was far larger, more energetic, and longer lasting than would be inferred from the prominence observations alone. The transient of 26–27 August was slow and of atypical shape compared to other mass ejection transients, but we believe that these three conclusions apply to most, if not all, of the more than 60 loop-shaped coronal transients observed by the High Altitude Observatory's coronagraph during the nine-month flight of Skylab.     

Interpretation of single-well tracer tests using fractional-flow dimensions. Part 1: Theory and mathematical models

Generalized equations using fractional-flow dimensions were derived to estimate the Darcy and seepage velocities obtained from the point-dilution and the single-well injection-withdrawal field tests conducted in fractured-rock aquifers. Seepage velocities can only be estimated from single-well tests if the hydraulic conductivity and the hydraulic gradient are known a priori. However, if a radial-convergent test is also performed between two boreholes, the kinematic porosity can be estimated and be used to estimate the seepage velocity from the single-well test results. To apply the generalized equations, the flow dimension and the extent of the flow region must be known. Therefore, the generalized radial flow (GRF) model of Barker (1988; a generalized radial flow model for hydraulic tests in fractured rock. Water Resour Res 24(10):1796–1804) is used to estimate the flow dimension because of its wide range of applications. A pumping test performed on the boreholes will yield an estimate of the fractional-flow dimension by applying the GRF model. Electronic Publication     

The fall and rise of metropolitan government in Copenhagen

In this paper it is argued that the re-establishment of a metropolitan government in the Greater Copenhagen Area in 2000 (Greater Copenhagen Authority) is not to be seen as a return to the kind of regional, coordinative authority (Greater Copenhagen Council) that was abolished in the Thatcherite climate of the 1980s. Throughout the 1980s and 1990s new forms of urban governance, entrepreneurialism and a higher national priority to the capital region were introduced, which changed the conditions under which the new metropolitan government is operating. This `contextualization' of the analysis of the two generations of metropolitan government points to the fact that typologies of metropolitan government, such as the one put forward by Sharpe (1995), often lacks analytical cogency. This revised version was published online in August 2006 with corrections to the Cover Date.     

Atmospheric nitrogen inputs into the North Sea: effect on productivity

The ANICE (Atmospheric Nitrogen Inputs into the Coastal Ecosystem) project addressed the atmospheric deposition of nitrogen to the North Sea, with emphasis on coastal effects. ANICE focused on quantifying the deposition of inorganic nitrogen compounds to the North Sea and the governing processes. An overview of the results from modelling and experimental efforts is presented. They serve to identify the role of the atmosphere as a source of biologically essential chemical species to the marine biota. Data from the Weybourne Atmospheric Observatory (UK) are used to evaluate the effect of short episodes with very high atmospheric nitrogen concentrations. One such episode resulted in an average deposition of 0.8 mmol N m⁻² day⁻¹, which has the potential to promote primary productivity of 5.3 mmol C m⁻² day⁻¹. This value is compared to long-term effects determined from model results. The total calculated atmospheric deposition to the North Sea in 1999 is 948 kg N km⁻¹, i.e. 0.19 mmol N m⁻² day⁻¹ which has the potential to promote primary productivity of 1.2 mmol C m⁻² day⁻¹. Detailed results for August 1999 show strong gradients across the North Sea due to adjacent areas where emissions of NO_x and NH₃ are among the highest in Europe. The average atmospheric deposition to the southern part of the North Sea in August 1999 could potentially promote primary production of 2.0 mmol C m⁻² day⁻¹, i.e. 5.5% of the total production at this time of the year in this area of the North Sea. For the entire study area the atmospheric contribution to the primary production per m² is about two-third of this value. Most of the deposition occurs during short periods with high atmospheric concentrations. This atmospheric nitrogen is almost entirely anthropogenic in origin and thus represents a human-induced perturbation of the ecosystem.