Erosion rates of alpine bedrock summit surfaces deduced from in situ Be and Al

We have measured the concentration of in situ produced cosmogenic ¹⁰Be and ²⁶Al from bare bedrock surfaces on summit flats in four western U.S. mountain ranges. The maximum mean bare-bedrock erosion rate from these alpine environments is 7.6 ± 3.9 m My⁻¹. Individual measurements vary between 2 and 19 m My⁻¹. These erosion rates are similar to previous cosmogenic radionuclide (CRN) erosion rates measured in other environments, except for those from extremely arid regions. This indicates that bare bedrock is not weathered into transportable material more rapidly in alpine environments than in other environments, even though frost weathering should be intense in these areas. Our CRN-deduced point measurements of bedrock erosion are slower than typical basin-averaged denudation rates ( 50 m My⁻¹). If our measured CRN erosion rates are accurate indicators of the rate at which summit flats are lowered by erosion, then relief in the mountain ranges examined here is probably increasing.

We develop a model of outcrop erosion to investigate the magnitude of errors associated with applying the steady-state erosion model to episodically eroding outcrops. Our simulations show that interpreting measurements with the steady-state erosion model can yield erosion rates which are either greater or less than the actual long-term mean erosion rate. While errors resulting from episodic erosion are potentially greater than both measurement and production rate errors for single samples, the mean value of many steady-state erosion rate measurements provides a much better estimate of the long-term erosion rate.     

Numerical Studies with a Regional Atmospheric Climate Model Based on Changes in the Roughness Length for Momentum and Heat Over Antarctica

A regional atmospheric climate model is used toexamine the effect of changes in the roughnesslengths of momentum (_z0m) and heat (_z0h)on the structure of the lower atmosphere and on thesurface energy fluxes over Antarctica. Fourexperiments were carried out in which _z0mand/or _z0h were altered with respect to acontrol experiment. The changes consisted of (1) alowering of _z0m from a field aggregated froma vegetation map with an orographic correction basedon the European Centre for Medium-Range WeatherForecasts _z0m field, to a constant value of10^-3 m; and (2) a lowering of _z0h from a valueequal to _z0m to a constant value of 10^-3 mor a value dependent on the wind speed via a surfacerenewal model. A reduction of _z0m results in theexpected increase in near-surface wind speed. It alsoresults in an increase in the depth of the layer in whichsouth-easterly near-surface winds prevail, and in adecrease in the strength of the large-scale flow overthe continent, in particular in summer. In theescarpment region a decrease of _z0m is foundto result in too high wind speeds. Surface temperatureson average decrease while atmospheric temperaturesincrease, resulting in an increase of near-surfacestatic stability. Changes in roughness lengths donot significantly change the temperature profiles.The surface fluxes, on average found reduced, aremodelled best by using the _z0h based on thesurface renewal method.     

Use of microtremor measurement for assessing site effects in Northern Belgium – interpretation of the observed intensity during the MS = 5.0 June 11 1938 earthquake

To provide quantitative information on site effects in the northern part of Belgium, forty-seven H/V microtremor measurements were performed with 5 second seismometers over an area of about 15.000 km². Most of the results show a northward regular increase of the fundamental period in agreement with the augmentation of the Mesozoic and Cainozoic soft sedimentsthickness from a few meters 40 km south of Brussels to 900 m at the Netherlands-Belgium border. The measured resonance frequency values wereconsistent with theoretical computations performed at different sites onthe basis of existing information and shallow seismic experiments. At oneparticular site (Uccle) where borehole data were available, microtremor measurements using an array of four seismological stations with differentapertures allowed to obtain the low frequency part of the Rayleigh wave dispersion curve, extending the range covered by the analysis of surfacewaves artificially generated. The Vs profile derived from the surface waveinversion corroborates the 1 Hz natural frequency of the site. Comparison of these results with the macroseismic information concerning the M_S = 5.0 1938 earthquake which occurred 50 km west of Brussels,confirmed the hypothesis that the geological structure of the Brabant massifis likely to control damage distribution during such an earthquake. Comparisonbetween the intensity map of the 1938 earthquake and the resonance period ofsediments obtained by our microtremor study shows a clear relation betweenthe two parameters. During the 1938 earthquake, site effects played a prominent role due to the dimension of the source whose corner frequency wasabout 1 Hz.     

Petrogenesis of orthopyroxene- and amphibole-bearing andesites, Mustique, Grenadine Islands, Lesser Antilles Arc: Isotope, trace element and physical constraints

Abstract On the island of Mustique, fresh and propylitized olivine–plagioclase–clinopyroxene basalt, plagioclase–clinopyroxene–orthopyroxene and plagioclase–clinopyroxene–amphibole andesite lavas and minor intrusions are interbedded with Oligocene pyroclastic and epiclastic rocks. Chemical data show that two isotopically identical, but chemically different, suites of lava are present: (i) the OPXS (⁸⁷Sr/⁸⁶Sr 0.70403–0.70454; ¹⁴³Nd/¹⁴⁴Nd 0.512952–0.512986; δ¹⁸O_cpx 5.49 and 5.61), comprising basalts and orthopyroxene‐bearing andesites; and (ii) the AMPHS (⁸⁷Sr/⁸⁶Sr 0.70401–0.70457; ¹⁴³Nd/¹⁴⁴Nd 0.512981–0.513037; δ¹⁸O_cpx 5.54), made up of basalts and amphibole‐bearing andesites. The OPXS has higher contents of TiO₂, P₂O₅, light rare earth elements, Sm, Pb, Th, U, Zr, Y and Nb, and higher La/Yb ratios than the AMPHS. The isotopic data suggest that both suites formed from melts derived from the same subduction‐modified depleted mantle source as the volcanic rocks of nearby St Vincent and Bequia, and the northern islands of the Lesser Antilles Arc. The immobile trace element contents, and La/Yb ratios, of the OPXS are indicative of ～10% partial melting of the source, whereas those of the AMPHS are indicative of ～25% partial melting. The within‐suite chemical variation of the OPXS is consistent with ～45% fractional crystallization of its intratelluric mineral assemblages, and that of the AMPHS is consistent with the removal of ～65% of its intratelluric assemblages. Experimental evidence suggests that both suites of basalt crystallized at pressures <8 kbar from melts containing 1–2 wt% water. After extensive fractional crystallization, the andesites crystallized at pressures between approximately 5 and 2 kbar. The OPXS magmas appear to have lost more of their water content than the AMPHS magmas. Thus, the OPXS andesites formed from melts with an estimated water content of 2–3 wt%, whereas the AMPHS andesites formed from melts containing at least 4.5 wt% water.     

The International Laser Ranging Service and its support for IGGOS

The International Laser Ranging Service (ILRS) was established in September 1998 as a service within the IAG to support programs in geodetic, geophysical, and lunar research activities and to provide data products to the International Earth Rotation Service (IERS) in support of its prime objectives. Now in operation for 5 years, the ILRS develops: (1) the standards and specifications necessary for product consistency and (2) the priorities and tracking strategies required to maximize network efficiency. The service collects, merges, analyzes, archives and distributes satellite and lunar laser ranging data to satisfy a variety of scientific, engineering, and operational needs and encourages the application of new technologies to enhance the quality, quantity, and cost effectiveness of its data products. The ILRS works with: (1) the global network to improve station performance; (2) new satellite missions in the design and building of retroreflector targets to maximize data quality and quantity and (3) science programs to optimize scientific data yield. The ILRS Central Bureau maintains a comprehensive web site as the primary vehicle for the distribution of information within the ILRS community. The site, which can be accessed at: http://ilrs.gsfc.nasa.gov is also available at mirrored sites at the Communications Research Laboratory (CRL) in Tokyo and the European Data Center (EDC) in Munich.During the last 2 years, the ILRS has addressed very important challenges: (1) data from the field stations are now submitted hourly and made available immediately through the data centers for access by the user community; (2) tracking on low satellites has been significantly improved through the sub-daily issue of predictions, drag functions, and the real-time exchange of time biases; (3) analysis products are now submitted in SINEX format for compatibility with the other space geodesy techniques; (4) the Analysis Working Group is heavily engaged in Pilot Projects as it works toward an ILRS “standard” global solution and (5) SLR has significantly increased its participation in the International Terrestrial Reference Frame (ITRF) activity, which is important to the success of IGGOS.     

Data assimilation (4D-VAR) to forecast flood in shallow-waters with sediment erosion

In this paper, the four-dimensional variational data assimilation technique (4D-VAR) is presented as a tool to forecast floods. Our study is limited to purely hydrological flows and supposes that the weather, here a big rain, has been already forecasted by meteorological services. The technique consists in minimizing, in the sense of Lagrange, the cost function: a measure of the difference between calculated data and available observations, here the water level. This is done under constraints that are the equations of the physical model. In our case, we modified the shallow-water equations to include a simplified sediment transport model. The steepest descent algorithm is then used to find the minimum. This is made possible because we can compute analytically the gradient of the cost function by using the adjoint equations of the model. As an application of the 4D-VAR technique, the overflowing of the Chicoutimi River at the Chute-Garneau dam, during the 1996 flood, is investigated. It is found that the 4D-VAR method reduces the error in the water height forecast even when the erosion model is not activated. In terms of Lyapunov exponents, we estimate the predictability horizon of such an event to be about half-an-hour after a big rain. However, this limit of predictability can be increased by using more observations or by using a finer computational grid.     

Numerical modelling of interface scattering of seismic wavefield from a random rough interface in an acoustic medium: comparison between 2D and 3D cases

Seismic wavefield scattering from a statistically randomly rough interface in a multilayered piecewise homogeneous medium is studied in 3D. The influence of the surface roughness on the scattered wavefield is analysed numerically by using a finite‐difference operator in the acoustic domain. Since interface scattering in the real practical sense is a 3D physical phenomenon, we show in this work that the scattering response of a randomly rough interface is not the same in 3D situations as in the 2D cases described in some earlier works. For a given interface roughness height in 3D, an interface roughness height at least three times greater is required to produce an equivalent phase scattering effect in 2D situations, for a given correlation length of the interface roughness scale. Based on observations from spectral analysis, we show that scattering results principally in de‐phasing and frequency band‐limiting of the incident wavefront, the frequency band‐limiting properties being comparable to cases reported in the literature for absorption and thin‐layer filtering. The interface scattering phenomenon should be critically considered when using amplitude and phase information from seismic signal during inversion processes.     

From water to tillage erosion dominated landform evolution

While water and wind erosion are still considered to be the dominant soil erosion processes on agricultural land, there is growing recognition that tillage erosion plays an important role in the redistribution of soil on agricultural land. In this study, we examined soil redistribution rates and patterns for an agricultural field in the Belgian loess belt. ¹³⁷Cs derived soil erosion rates have been confronted with historical patterns of soil erosion based on soil profile truncation. This allowed an assessment of historical and contemporary landform evolution on agricultural land and its interpretation in relation to the dominant geomorphic process. The results clearly show that an important shift in the relative contribution of tillage and water erosion to total soil redistribution on agricultural land has occurred during recent decades. Historical soil redistribution is dominated by high losses on steep midslope positions and concavities as a result of water erosion, leading to landscape incision and steepening of the topography. In contrast, contemporary soil redistribution is dominated by high losses on convex upperslopes and infilling of slope and valley concavities as a result of tillage, resulting in topographic flattening. This shift must be attributed to the increased mechanization of agriculture during recent decades. This study shows that the typical topographical dependency of soil redistribution processes and their spatial interactions must be accounted for when assessing landform and soil profile evolution.     

Preliminary Characterisation of New Glass Reference Materials (GSA-1G, GSC-1G, GSD-1G and GSE-1G) by Laser Ablation-Inductively Coupled Plasma-Mass Spectrometry Using 193 nm, 213 nm and 266 nm Wavelengths

New glass reference materials GSA-1G, GSC-1G, GSD-1G and GSE-1G have been characterised using a prototype solid state laser ablation system capable of producing wavelengths of 193 nm, 213 nm and 266 nm. This system allowed comparison of the effects of different laser wavelengths under nearly identical ablation and ICP operating conditions. The wavelengths 213 nm and 266 nm were also used at higher energy densities to evaluate the influence of energy density on quantitative analysis. In addition, the glass reference materials were analysed using commercially available 266 nm Nd:YAG and 193 nm ArF excimer lasers. Laser ablation analysis was carried out using both single spot and scanning mode ablation. Using laser ablation ICP-MS, concentrations of fifty-eight elements were determined with external calibration to the NIST SRM 610 glass reference material. Instead of applying the more common internal standardisation procedure, the total concentration of all element oxide concentrations was normalised to 100%. Major element concentrations were compared with those determined by electron microprobe. In addition to NIST SRM 610 for external calibration, USGS BCR-2G was used as a more closely matrix-matched reference material in order to compare the effect of matrix-matched and non matrix-matched calibration on quantitative analysis. The results show that the various laser wavelengths and energy densities applied produced similar results, with the exception of scanning mode ablation at 266 nm without matrix-matched calibration where deviations up to 60% from the average were found. However, results acquired using a scanning mode with a matrix-matched calibration agreed with results obtained by spot analysis. The increased abundance of large particles produced when using a scanning ablation mode with NIST SRM 610, is responsible for elemental fractionation effects caused by incomplete vaporisation of large particles in the ICP.