Ice flow physical processes derived from the ERS-1 high-resolution map of the Antarctica and Greenland ice sheets

The ERS-1 satellite, launched in 1991, has provided altimetric observations of the Greenland Ice Sheet and 80 per cent of the Antarctica Ice Sheet north of 82°S. It was placed in a geodetic (168-day repeat) orbit between April 1994 and March 1995, yielding a 1.5  km across-track spacing at latitude 70° with a higher along-track sampling of 350  m. We have analysed the waveform altimetric data from this period to compute maps with a 1/30° grid size. Data processing consists of correcting for environmental factors and editing and retracking the waveforms. A further step consists of reducing the radial orbit error through crossover analysis and correcting the slope error to second order. The high-resolution topography of both ice sheets reveals numerous details. A kilometre-scale surface roughness running at 45° from the flow direction is the dominant topographic characteristic of both continents. Antarctica also exhibits many scars due to local flow anomalies. Several physical processes can be identified: abrupt transitions from deformation to sliding and vice versa, and impressive strike-slip phenomena, inducing en echelon folds.     

Improved ocean-geoid resolution from retracked ERS-1 satellite altimeter waveforms

The ocean geoid can be inferred from the topography of the mean sea surface. Satellite altimeters transmit radar pulses and determine the return traveltime to measure sea-surface height. The ERS-1 altimeter stacks 51 consecutive radar reflections on board the satellite to a single waveform. Tracking the time shift of the waveform gives an estimate of the distance to the sea surface. We retrack the ERS-1 radar traveltimes using a model that is focused on the leading edge of the waveforms. While earlier methods regarded adjacent waveforms as independent statistical events, we invert a whole sequence of waveforms simultaneously for a spline geoid solution. Smoothness is controlled by spectral constraints on the spline coefficients. Our geoid solutions have an average spectral density equal to the expected power spectrum of the true geoid. The coherence of repeat track solutions indicates a spatial resolution of 31  km, as compared to 41  km resolution for the ERS-1 Ocean Product. While the resolution of the latter deteriorates to 47  km for wave heights above 2  m, our geoid solution maintains its resolution of 31  km for rough sea. Retracking altimeter waveform data and constraining the solution by a spectral model leads to a realistic geoid solution with significantly improved along-track resolution.     

Combined use of field observations and SAR interferometry to study ice dynamics and mass balance in Dronning Maud Land, Antarctica

In the region of the Schirmacheroase (71 °S, 12°E) various geodetic and glaciological research activities have been carried out in the last decade. Several times three geodetic-glaciological traverses were undertaken to study ice velocity, accumulation and ablation, and ice surface height changes. Repeated ground surveys show a significant decrease in surface heights by about 15 cm/y for a large blue-ice area. This paper presents the first interferometrically derived ice velocity field of the inland ice close to the Schirmacheroase. The interferometric analysis of the synthetic aperture radar (SAR) data is performed in combination with ground-based information. Since only ERS-1&2 tandem mission image couples are available for this region a digital elevation model (DEM) is used to remove the effect of topography. Ice velocities up to 100 m/y are proved interferometrically for this part of the inland ice.     

Fully normalized spherical cap harmonics: application to the analysis of sea-level data from TOPEX/POSEIDON and ERS-1

We introduce two sets of fully normalized harmonics for the spectral analysis of functions defined on a spherical cap. The harmonics are the products of Fourier functions and the fully normalized associated Legendre functions of non-integer degree. Using Sturm-Liouville theory for boundary-value problems, we present two convenient and stable formulae for computing the zeros of the associated Legendre functions that form two sets of orthogonal functions. Formulae for the stable numerical evaluation of the fully normalized associated Legendre functions of non-integer degree that avoid the gamma function are also derived. The result from the expansions of sea-level anomaly from altimetry into Set 2 fully normalized cap harmonics shows fast convergence of the series, and the degree variances decay rapidly without aliasing effects. The zero-degree coefficients (Set 2) of sea-level anomaly from TOPEX/ POSEIDON (T/P) and ERS-1 indicate an El Niño event during 1993 January-1993 July, and a La Niño event during 1993 November-1994 July, although the ERS-1 result is less obvious. Ocean circulations over the South China Sea and the Kuroshio area are clearly identified with the low-degree expansions of sea-surface topography (SST) from T/P and ERS-1. A cold-core eddy of 4° in diameter centred at 17.5°N, 118 E was detected with the expansion of SST from T/P cycle 47, and a property of the cap harmonics is used to compute this eddy's kinetic energy. The kinetic energy is at a low in winter and high in summer, and its variation seems to be periodic with an amplitude of 0.4 m² S^-2.     

Spectral analysis using orthonormal functions with a case study on the sea surface topography

b
Spherical harmonics are orthonormalized using the Gram-Schmidt process in a function space. The problem of linear dependence of spherical harmonics over the oceans is studied using the Gram matrices and consequently three sets of orthonormal (ON) functions have been constructed. For the process an efficient formula for computing inner products of spherical harmonics has been developed. Important spectral properties of the ON functions are addressed. The ON functions may be used for representing the sea surface topography (SST) in the analysis of satellite altimeter data. The geoid error can be transformed to a representation by the ON functions and hence the comparison of powers of the geoid error and the SST signal only over the oceans is possible, leading to a better way of determining the cut-off frequency of the SST in the simultaneous solution using satellite altimeter data. As a case study, the modified Levitus SST is expanded into the ON functions. The results show that 99.90 per cent of that signal's energy is contained within degree 24 of the orthonormal functions. Such expansions also render better spectral behaviour of oceanic signals as compared to that from spherical harmonic expansions. The study shows that these generalized Fourier functions are suitable for spectral analyses of oceanic signals and they can be applied to future altimetric mission where the geoid and the SST are to be recovered.     

Detailed 1×1° gravimetric Indian Ocean geoid and comparison with GEOS-3 radar altimeter geoid profiles

Summary. A new set of 1×1° mean free-air anomalies in the Indian Ocean is determined on the basis of previously published free-air anomaly maps (Talwani & Kahle) and the most recent Lamont surface ship gravity measurements. The data are then used to compute a (total) 1×1° gravimetric Indian Ocean geoid. The computation is carried out by combining the Goddard Space Flight Center (GSFC) GEM-6 geoid and a difference geoid that corresponds to the differences between the set of 1×1° surface gravity values and the GEM-6 gravity anomalies. The difference geoid is highest over the Madagascar Ridge (+ 20 m) and lowest over the Timor Trough (-30 m). The total geoid is compared with GEOS-3 radar altimeter derived geoid profiles and geophysical implications are discussed.     

Combination of TOPEX/POSEIDON data with a hydrographic inversion for determination of the oceanic general circulation and its relation to geoid accuracy

A global estimate of the absolute oceanic general circulation from a geostrophic inversion of in situ hydrographic data is tested against and then combined with an estimate obtained from TOPEX/POSEIDON altimetric data and a geoid model computed using the JGM-3 gravity-field solution. Within the quantitative uncertainties of both the hydrographic inversion and the geoid estimate, the two estimates derived by very different methods are consistent. When the in situ inversion is combined with the altimetry/geoid scheme using a recursive inverse procedure, a new solution, fully consistent with both hydrography and altimetry, is found. There is, however, little reduction in the uncertainties of the calculated ocean circulation and its mass and heat fluxes because the best available geoid estimate remains noisy relative to the purely oceano-graphic inferences. The conclusion drawn from this is that the comparatively large errors present in the existing geoid models now limit the ability of satellite altimeter data to improve directly the general ocean circulation models derived from in situ measurements. Because improvements in the geoid could be realized through a dedicated spaceborne gravity recovery mission, the impact of hypothetical much better, future geoid estimates on the circulation uncertainty is also quantified, showing significant hypothetical reductions in the uncertainties of oceanic transport calculations, Full ocean general circulation models could better exploit both existing oceanographic data and future gravity-mission data, but their present use is severely limited by the inability to quantify their error budgets.     

On the quality of marine magnetic anomaly sources and sea-floor topography

Summary. The quality of marine magnetic anomaly sources is described with the power-density representation of a stochastic model of random temporal and spatial emplacement of the marine magnetic anomaly source in the oceanic crust. Typical values of sea-floor spreading and emplacement parameters define a high-fidelity process of recording and sea-surface detection of the palaeomagnetic field reversals for spreading rates over 20 mm yr^-1. An analogous stochastic model is developed for the formation of sea-floor topography by normal faulting. It is shown that the random process of normal faulting observed in the inner walls of the FAMOUS rift valley can account for the quality of the adjacent West Rift Mountains topography.     

CryoSat-2卫星测高计划及其应用

作为欧空局地球探测计划的一项重要任务,Cryo Sat-2卫星于2010年携带着Ku波段SIRAL(干涉/合成孔径雷达高度计)发射升空。Cryo Sat-2任务的优势在于高精度地测量两极海冰厚度以及监测南极大陆和格陵兰冰盖的变化,同时92°的近极地轨道一定程度上填补了先前遥测卫星的高纬数据缺口。本文对Cryo Sat-2卫星的科学目标、科学需求、任务概况、仪器载荷以及数据产品等方面做了详细的介绍,重点强调了Cryo Sat-2的首要载荷SIRAL和数据校准,最后探讨了Cryo Sat-2在极地领域的应用。     

哨兵卫星Sentinel-1A数据特性及应用潜力分析

Sentinel-1A是欧空局"哥白尼计划"发射的首颗对地观测卫星,其搭载的C波段SAR传感器有效地延续了ERS-1/2和ENVISAT ASAR对地观测任务。简要介绍了Sentinel-1A卫星的特点和轨道参数,并与ERS-1/2和ENVISAT ASAR进行了对比。详细介绍了Sentinel-1A的成像模式(SM,IW,EW,WV)以及每种模式下数据产品(Level-0,Level-1,Level-2)的特征。最后分析了Sentinel-1A数据在不同领域的实际应用,为其在对地观测中的广泛应用提供了参考。     

Evidence for underthrusting beneath the Queen Charlotte Margin, British Columbia, from teleseismic receiver function analysis

The Queen Charlotte Fault zone is the transpressive boundary between the North America and Pacific Plates along the northwestern margin of British Columbia. Two models have been suggested for the accommodation of the ∼20 mm yr⁻¹ of convergence along the fault boundary: (1) underthrusting; (2) internal crustal deformation. Strong evidence supporting an underthrusting model is provided by a detailed teleseismic receiver function analysis that defines the underthrusting slab. Forward and inverse modelling techniques were applied to receiver function data calculated at two permanent and four temporary seismic stations within the Queen Charlotte Islands. The modelling reveals a ∼10 km thick low-velocity zone dipping eastward at 28° interpreted to be underthrusting oceanic crust. The oceanic crust is located beneath a thin (28 km) eastward thickening (10°) continental crust.     

Crustal thickness estimation beneath the southern central Andes at 30°S and 36°S from S wave receiver function analysis

We have used the S wave receiver function (SRF) technique to investigate the crustal thickness beneath two seismic profiles from the CHARGE project in the southern central Andes. A previous study employing the P wave receiver function method has observed the Moho interface beneath much of the profiles. They found, however, that the amplitude of the P to S conversion was diminished in the western part of the profiles and have attributed it to a reduction of the impedance contrast at the Moho due to lower crustal ecologitization. With SRF, we have successfully detected S to P converted waves from the Moho as well as possible conversions from other lithospheric boundaries. The continental South American crust reaches its maximum thickness of ∼70 km (along 30°S between 70°W and 68.5°W) beneath the Principal Cordillera and the Famatina system and becomes thinner towards the Sierras Pampeanas with a thickness of ∼40 km. Negative phases, possibly related to the base of the continental and oceanic lithosphere, can be recognized in the summation traces at different depths. By comparing our results with data obtained from previous investigations, we are able to further constrain the thickness of the crust and lithosphere beneath the central Andes.     

Pn anisotropy studies in northern Britain and the eastern and western United States

Summary. Available seismic refraction data from three different continental areas, northern Britain and the eastern and western United States, has been studied for possible P_n , velocity anisotropy using the methods described by Bamford. There are various deficiencies in the time—distance data used in each case but, while the uppermost mantle beneath northern Britain and the eastern United States seems to be isotropic within the limits of measurement error, there is a small but significant anisotropy beneath the western United States.
Both the amount (up to 3 per cent) and the direction (70–80° east of north) of this anisotropy are very similar to the results obtained in the Pacific Ocean off California. We tentatively conclude that this anisotropy is present as a consequence of the subduction of oceanic lithosphere beneath the western United States.     

Lithoprobe seismic reflection profiling across Vancouver Island: results from reprocessing

Summary. Multichannel seismic reflection sections recorded across Vancouver Island have revealed two extensive zones of deep seismic reflections that dip gently to the northeast, and a number of moderate northeasterly dipping reflections that can be traced to the surface where major faults are exposed. Based on an integrated interpretation of these data with information from gravity, heat flow, seismicity, seismic refraction, magnetotelluric and geological studies it is concluded that the lower zone of gently dipping reflections is due to underplated oceanic sediments and igneous rocks associated with the current subduction of the Juan de Fuca plate, and that the upper zone represents a similar sequence of accreted rocks associated with an earlier episode of subduction. The high density/high velocity material between the two reflection zones is either an underplated slab of oceanic lithosphere or an imbricated package of mafic rocks. Reprocessing of data from two of the seismic lines has produced a remarkable image of the terrane bounding Leech River fault, with its dip undulating from >60° near the surface to 20° at 3 km depth and ∼38° at 6 km depth.     

Gravimetric geoid in the Northwest Pacific Ocean

Summary. A total of 3708 1 × 1° free-air gravity anomaly averages have been used to construct a new 1 × 1° gravimetric geoid of the Northwest Pacific Ocean. The 1 × 1° averages are based on a compilation of 147000 surface ship and pendulum gravity measurements. The gravimetric geoid reveals information in the geoid of the Northwest Pacific not present in currently used satellite derived models. The RMS difference between the 1 × 1° geoid and satellite derived models is about ±6 m. Difference geoid undulations range from a maximum of +19 m over the Hawaiian ridge to a minimum of −31 m over the junction of the Kuril and Aleutian trenches. The Hawaiian swell is associated with a geoidal high of up to +15 m with wavelengths of about 2200 km and the topographic rises seaward of deep-sea trenches are associated with geoidal highs of up to 4m with wavelengths of about 220–900 km. The main difference between the gravimetric geoid and the satellite derived models occurs over the Pacific basin where discrepancies reach +10 m with wavelengths of 4000 km. The agreement between the gravi-metric geoid and Skylab-4 and Geos-3 altimeter data is close for wavelengths greater than about 300 km but poor for shorter wavelengths.     

Oceanic basalt continuous thermal demagnetization curves

The palaeomagnetic standard technique of stepwise thermal demagnetization (STD), long regarded as unreliable for oceanic basalts that have undergone low temperature alteration, has recently been applied in a number of studies to characterize the natural remanent magnetization (NRM) of such rocks. In order to better understand STD data of oceanic basalts, and to possibly identify the magnetominerals that are carrying the NRM, we have carried out a number of continuous and STD experiments on seven oceanic basalt samples. During continuous thermal demagnetization (CTD), a sample is heated to a certain temperature and its NRM is measured during heating and subsequent cooling. Even when CTD reveals only titanomaghemite unblocking at 400°C as the remanence carrier, STD behaviour can be very complex and unblocking is observed at temperatures of up to 500°C and higher. CTD also allowed to identify a partial or full self-reversal of NRM due to interaction between two types of magnetominerals in one sample. The higher degree of maghemitization of smaller titanomaghemite grains with respect to larger ones, which are less efficient in carrying the remanence, was seen for three samples by a shift of 80°C between the strong field thermomagnetic curve and the NRM measured at elevated temperature. In several cases, the identification of the NRM-carrying magnetomineral was not possible from CTD data due to the ambiguity of Curie temperatures in the titanomagnetite/titanomaghemite system.     

The late Maastrichtian palaeomagnetic pole of the Pacific Plate

Summary. The Pacific plate's late Maastrichtian (∼ 69 Ma) palaeomagnetic pole, which constrains the northward motion of the Pacific plate during the Cainozoic and latest Cretaceous, was studied. A recently proposed method for obtaining oceanic plate palaeomagnetic poles by combining dissimilar data was extended to accept, as input, the relative amplitudes of magnetic lineations with different azimuths or widely separated sites or both. Combining late Maastrichtian palaeomagnetic data-the relative amplitudes and skewness of magnetic lineations, palaeolatitudes from a palaeomagnetic study of basalt and sediment in vertical cores, a pole from the inversion of the magnetic anomaly over a seamount, and present locations of equatorial sediment facies—yielded a best fit pole of 71°N, 9°E and a 95 per cent confidence ellipse with the major semiaxis of 6° striking 91° clockwise from north and the minor semiaxis of 2° striking 1° clockwise from north. This best fit pole, when compared to the pole expected if the hotspots have been fixed with respect to the spin axis, demonstrates that the hotspots in the Pacific Ocean have shifted ∼ 10° south with respect to the spin axis during the Cainozoic. This best fit pole, when compared to the best fit Campanian pole of the Pacific plate, demonstrates that the pole wandered rapidly, 1.1° Ma^-1, with respect to the Pacific plate during the latest Cretaceous.     

ICESat与ICESat-2应用进展与展望

近十年来,格陵兰和南极冰盖消融增速,两极地区显著变化.美国国家航空航天管理局(NASA)的冰、云和陆地高程卫星(ICESat)自2003年发射至2009年的在轨期间,用于定量分析冰盖和海冰的变化速率,并对这些变化的驱动机制进行检测与监测.ICESat在极地冰雪遥感应用等方面展示出了无可比拟的优势,对精确估算冰盖对海平面...     

Interpretation of aeromagnetic data across the central crystalline shield area of Nigeria

Summary. A residual map of the total magnetic field (above 25 000 nT base) is presented for a portion of the central crystalline shield area of Nigeria and overlapping small portions of the Chad basin and the Benue rift (8°30'−12° 00'lat, and 7°−10°30' long). The map (based on a dataset digitized from recently released aeromagnetic sheets of Nigeria) leads to four results. (1) A magnetic boundary, evident on the map, separates the Younger Granite complexes into two groups. The groups are petrologically different, and the boundary may be a fault line with uplift to the south. (2) South of the boundary the map is dominated by a system of sub-parallel anomalies striking NE–SW, possibly representing major tectonic trends, and a set of fractures through which the Younger Granite complexes were intruded. The trend of the system parallels the Benue rift and lineaments in the oceanic crust off West Africa. (3) Negative magnetic anomalies lie over most of the known ring complexes, and over some suspected buried ring complexes and other intrusions. (4) 2½-and 3-D modelling shows that the larger complexes extend to 12 km depth, and the smaller ones to 6 km. They have nearly vertical sides, and magnetization contrasts range from 0.3 to 0.5 A m⁻¹.     

Anisotropy in the oceanic lithosphere of the North-western Pacific Basin

Summary. A series of long-range explosion seismological experiments has been conducted by the use of specially designed ocean bottom seismographs (OBSs) in the Western Pacific. OBS studies of apparent velocity measurements by the use of natural earthquakes have also been made. The experiments have made clear that large-scale P -wave anisotropy exists in the entire thickness of the oceanic lithosphere. The existence of the large-scale anisotropy in the oceanic lithosphere has been demonstrated for the first time by seismic body-wave studies. Previously, anisotropy had been found only in the uppermost oceanic mantle in the Eastern Pacific.
The azimuth of the maximum velocity, 8.6 km s^-1, is about 155° clock-wise from north. The direction is perpendicular to the magnetic lineation of the region, however, the direction differs from the direction of the present plate motion by about 30°. So it appears that the anisotropy has been 'frozen' at least since the change of the plate motion that occurred 40 Myr ago. The frozen anisotropy should set important constraints on the mechanical properties of the lithosphere such as the viscosity and temperature of the lower lithosphere.