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.     

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.     

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     

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.     

A mathematical programming approach for scheduling equipment in a surface coal mining operation

Summary A mathematical programming model for scheduling open pit mining was developed and validated using data from a surface mining operation. A two-phase solution procedure was used involving repeated evaluations of an integer scheduling model and a simple transportation model.     

Distinct modes of bidecadal and multidecadal variability in a climate reconstruction of the last centuries from a South Pacific coral

The Rarotonga coral Sr/Ca time series (Linsley et al. in Science 290:1145–1148, 2000) provides a near-monthly resolved proxy record of South Pacific climate variability over the last ~300 years. Here we show that two distinct interdecadal, quasi-periodic time components with periods of ~80 and ~25 years can be identified in this time series by Singular Spectrum Analysis. Their associated spatial patterns in the global sea surface temperature (SST) field show notable differences. Whereas the multidecadal component is associated with a global SST pattern that was recently associated with solar forcing on multidecadal timescales, the bidecadal component is associated with a well known pattern of Pacific decadal to interdecadal SST variability.     

On the robustness of ENSO teleconnections

Using observations covering the last 128 years we show that apparent changes in El Niño-Southern Oscillation (ENSO) teleconnections can be explained by chance and stem from sampling variability. This result is backed by experiments in which an atmosphere model is driven by 123 years of observed sea surface temperature. The possibility of ENSO teleconnection changes in a warming climate is further investigated using coupled GCMs driven by past and projected future greenhouse gas concentrations. These runs do not exclude physical changes in the teleconnection strength but do not agree on their magnitude and location. If existing at all, changes in the strength of ENSO teleconnection, other than obtained by chance, are small and will only be detectable on centennial time scales.     

Response of precipitation over Greenland and the adjacent ocean to North Pacific warm spells during Dansgaard–Oeschger stadials

ABSTRACT Palaeoceanographic reconstructions from the North Atlantic indicate massive ice breakouts from East Greenland near the onset of cold Dansgaard–Oeschger (D–O) stadials. In contrast to these coolings in the North Atlantic area, a new sea-surface temperature record reveals concomitant warm spells in the northern North Pacific. A sensitivity experiment with an atmospheric general circulation model is used to test the potential impact of sea-surface warmings by 3.5 °C in the North Pacific, on top of otherwise cold stadial climate conditions, on the precipitation regime over the Northern Hemisphere ice sheets. The model predicts a maximum response over East Greenland and the Greenland Sea, where a 40% increase in net annual snow accumulation occurs. This remote effect of North Pacific warm spells on the East Greenland snow-accumulation rate may play an important role in generating D–O cycles by rebuilding the ice lost during ice breakouts. In addition, the increased precipitation over the Greenland Sea may help to sustain the D–O stadial climate state.     

The challenge to detect and attribute effects of climate change on human and natural systems

Future scenarios of the energy system under greenhouse gas emission constraints depict dramatic growth in a range of energy technologies. Technological growth dynamics observed historically provide a useful comparator for these future trajectories. We find that historical time series data reveal a consistent relationship between how much a technology’s cumulative installed capacity grows, and how long this growth takes. This relationship between extent (how much) and duration (for how long) is consistent across both energy supply and end-use technologies, and both established and emerging technologies. We then develop and test an approach for using this historical relationship to assess technological trajectories in future scenarios. Our approach for “learning from the past” contributes to the assessment and verification of integrated assessment and energy-economic models used to generate quantitative scenarios. Using data on power generation technologies from two such models, we also find a consistent extent - duration relationship across both technologies and scenarios. This relationship describes future low carbon technological growth in the power sector which appears to be conservative relative to what has been evidenced historically. Specifically, future extents of capacity growth are comparatively low given the lengthy time duration of that growth. We treat this finding with caution due to the low number of data points. Yet it remains counter-intuitive given the extremely rapid growth rates of certain low carbon technologies under stringent emission constraints. We explore possible reasons for the apparent scenario conservatism, and find parametric or structural conservatism in the underlying models to be one possible explanation.     

Eemian tropical and subtropical African moisture transport: an isotope modelling study

During the last interglacial insolation maximum (Eemian, MIS 5e) the tropical and subtropical African hydrological cycle was enhanced during boreal summer months. The climate anomalies are examined with a General Circulation Model (ECHAM4) that is equipped with a module for the direct simulation of ¹⁸O and deuterium (H ₂ ¹⁸ O and HDO, respectively) in all components of the hydrological cycle. A mechanism is proposed to explain the physical processes that lead to the modelled anomalies. Differential surface heating due to anomalies in orbital insolation forcing induce a zonal flow which results in enhanced moisture advection and precipitation. Increased cloud cover reduces incoming short wave radiation and induces a cooling between 10°N and 20°N. The isotopic composition of rainfall at these latitudes is therefore significantly altered. Increased amount of precipitation and stronger advection of moisture from the Atlantic result in isotopically more depleted rainfall in the Eemian East African subtropics compared to pre-industrial climate. The East–West gradient of the isotopic rainfall composition reverses in the Eemian simulation towards depleted values in the east, compared to more depleted western African rainfall in the pre-industrial simulation. The modelled re-distribution of δ¹⁸O and δD is the result of a change in the forcing of the zonal flow anomaly. We conclude that the orbitally induced forcing for African monsoon maxima extends further eastward over the continent and leaves a distinct isotopic signal that can be tested against proxy archives, such as lake sediment cores from the Ethiopian region.