Drainage, sediment transport, and denudation rates on the Nanga Parbat Himalaya, Pakistan

The Nanga Parbat Himalaya presents some of the greatest relief on Earth, yet sediment production and denudation rates have only been sporadically addressed. We utilized field measurements and computer models to estimate bank full discharge, sediment transport, and denudation rates for the Raikot and Buldar drainage basins (north slope of Nanga Parbat) and the upper reach of the Rupal drainage basin (south slope).The overall tasks of determining stream flow conditions in such a dynamic geomorphic setting is challenging. No gage data exist for these drainage basins, and the overall character of the drainage basins (high relief, steep flow gradients, and turbulent flow conditions) does not lend itself to either ready access or complete profiling.Cross-sectional profiles were surveyed through selected reaches of these drainage basins. These data were then incorporated into software (WinXSPRO) that aids in the characterization (stage, discharge, velocity, and shear stress) of high altitude, steep mountain stream conditions.Complete field measurements of channel depths were rarely possible (except at several bridges where the middle of the channel could actually be straddled and probed) and, when coupled with velocity measurements, provided discrete points of field-measured discharge calculations. These points were then used to calibrate WinXSPRO results for the same reach and provided a confidence level for computer-generated results.Flow calculations suggest that under near bank full conditions, the upper Raikot drainage basin produces discharges of 61 cm and moves about 11,000 tons day⁻¹ (9980 tons day⁻¹) of sediment through its channel. Bank full conditions on the upper portion of the Rupal drainage basin generate discharges of 84 cm and moves only about 3800 tons day⁻¹ (3450 tons day⁻¹) of sediment. Although the upper Rupal drainage basin moves more water, the lower slope of the drainage basin (0.03) generates a much smaller shear stress (461 Pa) than does the higher slope (0.12) of the upper Raikot drainage basin (1925 Pa).Dissolved and suspended sediment loads were measured from water/sediment samples collected throughout the day and night over a period of 10 days at the height of the summer melt season but proved to be a minor variable in transport flux. Channel bed loads were measured using a pebble count method of bank material and then used to generate ratings curves of bed loads relative to discharge volumes. When coupled with discharge data and basin area, mean annual sediment yield and denudation rates for Nanga Parbat are produced. Denudation rates calculated in this fashion range from 0.2 mm year⁻¹ in the slower, more sluggish Rupal drainage basin to almost 6 mm year⁻¹ in the steeper, faster flowing Raikot and Buldar drainage basins.     

Tidal dwarfs in the M81 group: the second generation?

We present a quantitative star formation history derivation of the four suspected tidal dwarf galaxies in the M 81 group: Holmberg IX, BK3N,Arp-loop (A0952+69) and Garland using HST/WFPC2 images of these galaxies. We construct a library of synthetic Colour-Magnitude Diagrams(CMDs) based on theoretical isochrones and data-derived determinations of photometric errors. These synthetic CMDs were combined linearly andχ²-compared to observed photometry. All the galaxies show continuous star formation between about 20 and 200 Myr ago with star formation rates between 7.5⋅10^-3 M_⊙/yr and 7.67⋅10^-4 M_⊙/yr. The metallicity of the detected stars is spanning rather a wide range, being lower than solar abundance. We suppose, that all the galaxies were formed out of material from metal-poor outer part of the giant spiral galaxy M81after tidal interaction about 200 Myr ago. However, this suggestion requires significantly more deep color-magnitude diagrams to be sure with the scenario of the galaxy evolution. This revised version was published online in July 2006 with corrections to the Cover Date.     

Volcanic and Scientific Activity at Kick ’em Jenny Submarine Volcano 2001–2002: Implications for Volcanic Hazard in the Southern Grenadines,Lesser Antilles

Kick em Jenny submarine volcano, ~8 km north of Grenada, has erupted at least 12 times since it was first discovered in 1939, making it the most frequently active volcano in the Lesser Antilles arc. The volcano lies in shallow water close to significant population centres and directly beneath a major shipping route, and as a consequence an understanding of the eruptive behaviour and potential hazards at the volcano is critical. The most recent eruption at Kick em Jenny occurred on December 4 2001, and differed significantly from past eruptions in that it was preceded by an intensive volcanic earthquake swarm. In March 2002 a multi-beam bathymetric survey of the volcano and its surroundings was carried out by the NOAA ship Ronald H Brown. This survey provided detailed three-dimensional images of the volcano, revealing the detailed morphology of the summit area. The volcano is capped by a summit crater which is breached to the northeast and which varies in diameter from 300 to 370 m. The depth to the summit (highest point on the crater rim) is 185 m and the depth to the lowest point inside the crater is 264 m. No dome is present within the crater. The crater and summit region of Kick em Jenny are located at the top of an asymmetrical cone which is about 1300 m from top to bottom on its western side. It lies within what appear to be the remnants of a much larger arcuate collapse structure. An evaluation of the morphology, bathymetry and eruptive history of the volcano indicates that the threat of eruption-generated tsunamis is considerably lower than previously thought, mainly because the volcano is no longer thought to be growing towards the surface. Of more major and immediate concern are the direct hazards associated with the volcano, such as ballistic ejecta, water disturbances and lowered water density due to degassing.     

Ages and metallicities of globular clusters in NGC 4472

We have derived ages and metallicities from co-added spectra of 131 globular clusters associated with the giant elliptical galaxy NGC 4472. Based upon a calibration with Galactic globular clusters, we find that our sample of globular clusters in NGC 4472 span a metallicity range of approximately −1.6≤[Fe/H]≤0 dex. There is evidence of a radial metallicity gradient in the globular cluster system which is steeper than that seen in the underlying starlight. Determination of the absolute ages of the globular clusters is uncertain, but formally, the metal-poor population of globular clusters has an age of 14.5±4 Gyr and the metal-rich population is 13.8±6 Gyr old. Monte Carlo simulations indicate that the globular cluster populations present in these data are older than 6 Gyr at the 95 per cent confidence level. We find that within the uncertainties, the globular clusters are old and coeval, implying that the bimodality seen in the broadband colours primarily reflects metallicity and not age differences.     

Structural investigation of the synthetic CaAn(PO4)2 (An = Th and Np) cheralite-like phosphates

The crystallographic structures of the synthetic cheralite, CaTh(PO₄)₂, and its homolog CaNp(PO₄)₂ have been investigated by X-ray diffraction at room temperature. Rietveld analyses showed that both compounds crystallize in the monoclinic system and are isostructural to monazite LnPO₄ (Ln = La to Gd). The space group is P2₁/n (I.T. = 14) with Z = 2. The refined lattice parameters of CaTh(PO₄)₂ are a = 6.7085(8) Å, b = 6.9160(6) Å, c = 6.4152(6) Å, and β = 103.71(1)° with best fit parameters R _wp = 4.87%, R _p = 3.69% and R _B = 3.99%. For CaNp(PO₄)₂, we obtained a = 6.6509(5) Å, b = 6.8390(3) Å, c = 6.3537(8) Å, and β = 104.12(6)° and R _wp = 6.74%, R _p = 5.23%, and R _B = 6.05%. The results indicate significant distortions of bond length and angles of the PO₄ tetrahedra in CaTh(PO₄)₂ and to a lesser extent in CaNp(PO₄)₂. The structural distortions were confirmed by Raman spectroscopy of CaTh(PO₄)₂. A comparison with the isostructural compounds LnPO₄ (Ln = Ce and Sm) confirmed that the substitution of the large rare earth trivalent cations with Ca²⁺ and Th⁴⁺ introduces a distortion of the PO₄ tetrahedra.     

Modelling converted seismic waveforms in isotropic and anisotropic 1-D gradients: discontinuous versus continuous gradient representations

Over the past decade, there have been numerous receiver function studies directed at imaging the lithosphere-asthenosphere boundary (LAB). Although it is generally accepted that receiver function phases observed in these studies are derived from physical mode conversions at depth within the lithosphere-asthenosphere transition, it is still debatable as to whether these phases are directly indicative of the LAB. This is because interpretation of receiver function LAB signals relies on understanding the elastic characteristics of the Earth??s outer thermal boundary layer. The main issues for receiver function imaging are the sharpness of the elastic material property transition and, more importantly, what specifically are the material gradients. To test the various transition models, a forward modelling approach is required that allows accurate waveform synthetics for a range of discontinuous and continuous gradients in anisotropic, elastic media. We present a derivation of the reflection and transmission response for continuous one-dimensional (1-D) gradients in generally anisotropic elastic media. We evaluate the influence of 1-D isotropic and anisotropic elastic gradients on the seismic waveform by comparing numerical results of models for discontinuous and continuous transitions. The results indicate that discontinuous representations using layers each with uniform parameters and with thicknesses on the order of approximately 1/3 to 1/8 of the dominant seismic wavelength can be used to accurately model P-to-S and S-to-P mode conversions due to continuous transitions of both isotropic and anisotropic elastic properties. From a practical point of view, when comparing synthetic modelling with observation, this constraint can be relaxed further. The presence of signal noise and/or the result of receiver function stacking techniques will likely obscure these subtle waveform e ff ects. Hence this study suggests that accurate synthetic waveforms for LAB transitions can be modelled with discontinuous gradient representations using a reasonable number of discrete transition layers with layer thicknesses no greater than 1/2 to 1/3 the dominant seismic wavelength.     

Linear theory of the response of a two layer ocean to a moving hurricane‡

Abstract

This paper describes the linear response of an inviscid two‐layer model of a deep ocean on an f‐plane to a hurricane translating across the surface at constant speed. The forcing is a localized, radially‐symmetric pattern of positive wind stress curl and negative pressure anomaly. Only the steady state response is considered. The principal result is the identification of an internal wake in the lee of the storm, present when the translation speed of the storm exceeds the baroclinic long wave speed. The amplitude of the wake depends on the length of time over which the stress is experienced at a given point. The angle of the wedge filled by the wake is small, an effect due to the fact that the scale of a hurricane is typically larger than the baroclinic radius of deformation. After the wake disperses, a geostrophically balanced baroclinic ridge remains along the storm track.     

Stability analysis of the 2007 Chehalis lake landslide based on long-range terrestrial photogrammetry and airborne LiDAR data

On December 4th 2007, a 3-Mm³ landslide occurred along the northwestern shore of Chehalis Lake. The initiation zone is located at the intersection of the main valley slope and the northern sidewall of a prominent gully. The slope failure caused a displacement wave that ran up to 38 m on the opposite shore of the lake. The landslide is temporally associated with a rain-on-snow meteorological event which is thought to have triggered it. This paper describes the Chehalis Lake landslide and presents a comparison of discontinuity orientation datasets obtained using three techniques: field measurements, terrestrial photogrammetric 3D models and an airborne LiDAR digital elevation model to describe the orientation and characteristics of the five discontinuity sets present. The discontinuity orientation data are used to perform kinematic, surface wedge limit equilibrium and three-dimensional distinct element analyses. The kinematic and surface wedge analyses suggest that the location of the slope failure (intersection of the valley slope and a gully wall) has facilitated the development of the unstable rock mass which initiated as a planar sliding failure. Results from the three-dimensional distinct element analyses suggest that the presence, orientation and high persistence of a discontinuity set dipping obliquely to the slope were critical to the development of the landslide and led to a failure mechanism dominated by planar sliding. The three-dimensional distinct element modelling also suggests that the presence of a steeply dipping discontinuity set striking perpendicular to the slope and associated with a fault exerted a significant control on the volume and extent of the failed rock mass but not on the overall stability of the slope.     

The GPS Toolbox

The GPS Toolbox is dedicated to highlighting algorithms utilized by GPS engineers and scientists. If you have an interesting algorithm you would like to share with our readers or if you have a topic you would like to see covered in a future column, contact us at gps-toolbox@ngs.noaa.gov. To comment on the algorithms presented here, or to leave a request for an algorithm you may be looking for, visit our Web site (http://www.ngs.noaa.gov/gps-toolbox). ? 2000 John Wiley & Sons, Inc.     

Inversion of electromagnetic data: An overview of new techniques

This paper explores some of the newer techniques for acquiring and inverting electromagnetic data. Attention is confined primarily to the 2d magnetotelluric (MT) problem but the inverse methods are applicable to all areas of EM induction. The basis of the EMAP technique of Bostick is presented along with examples to illustrate the efficacy of that method in structural imaging and in overcoming the deleterious effects of near-surface distortions of the electric field. Reflectivity imaging methods and the application of seismic migration techniques to EM problems are also explored as imaging tools. Two new approaches to the solution of the inverse problem are presented. The AIM (Approximate Inverse Mapping) inversion of Oldenburg and Ellis uses a new way to estimate a perturbation in an iterative solution which does not involve linearization of the equations. The RRI (Rapid Relaxation Inverse) of Smith and Booker shows how approximate Fréchet derivatives and sequences of 1d inversions can be used to develop a practical inversion algorithm. The overview is structured to provide insight about the latest inversion techniques and also to touch upon most areas of the inverse problem that must be considered to carry out a practical inversion. These include model parameterization, methods of calculating first order sensitivities, and methods for setting up a linearized inversion.