Continuous zircon growth during long-lived granulite facies metamorphism: a microtextural,U–Pb,Lu–Hf and trace element study of Caledonian rocks from the Arctic

Microtextural, U–Pb, trace element and Lu–Hf analyses of zircons from gneisses dredged from the Chukchi Borderland indicate a long-lived, Cambrian–Ordovician, granulite facies metamorphism. These results reveal a complete prograde, peak and cooling history of zircon growth during anatexis. Early increasing temperatures caused modification and Pb-loss of Precambrian zircons by recrystallization and dissolution/re-precipitation of existing grains. Small variations in initial ¹⁷⁶Hf/¹⁷⁷Hf results (0.282325–0.282042) and flat HREE patterns of these zircons indicate that they grew by dissolution/re-precipitation in the presence of garnet. Zircons subsequently crystallized from a partial melt during peak to post-peak metamorphism from 530 to 485 Ma. A broad range of initial ¹⁷⁶Hf/¹⁷⁷Hf ratios (0.282693–0.282050) and mineral inclusions within zircons suggest that this phase of growth incorporated Zr and Hf obtained from the breakdown of Zr-enriched phases. Microtextural evidence along with trace element and isotopic data suggests that final growth of metamorphic rims on zircon occurred during slow cooling and crystallization of residual partial melts during the early Ordovician (485–470 Ma). Younger, late Ordovician–Silurian (420–450 Ma) euhedral, oscillatory-zoned, trace element-enriched zircons crystallized within leucocratic veins that intrude the gneisses. Their age corresponds to granitoids dated from this same dredge. The intrusives and veins provide evidence that the Chukchi Borderland rifted from a position near Pearya and northwest Svalbard, which represent the northern continuation of the Caledonian orogen. Evidence for earlier Cambrian metamorphism has not been reported from this region. The age of granulite facies metamorphism reported here represents the earliest phase of deformation in the Arctic Caledonides.     

Calculating the two-dimensional magnetotelluric Jacobian in finite elements using reciprocity

To speed up the calculation of the field Jacobian for 2-D magnetoteliuric inversion using finite elements, the principle of electromagnetic reciprocity is applied. The governing relationship for the Jacobian of the field along strike is obtained by differentiating the Helmholtz equation with respect to the resistivity of each region in the finite-element mesh. The result is a similar Helmholtz equation for the Jacobian, with new sources distributed over all nodes within the parameter medium. However, according to the principle of electromagnetic reciprocity, the roles of sources and receivers are interchangeable. Utilizing reciprocity, the field values obtained from the original forward problem and for new unit sources imposed at the receivers are then utilized in the calculation of the Jacobian by simple multiplication and summation with finite-element terms at each rectangle in the mesh. For the auxiliary (across-strike) fields, the Jacobian terms are obtained by solving source vectors loaded with parabola coefficients used in the approximation to Maxwell's equations. Jacobian terms for the apparent resistivity ( p _a), the impedance phase (φ) and the vertical magnetic field ( K _zy) are then calculated utilizing the parallel- and auxiliary-field Jacobians. Comparison of Jacobian values obtained from reciprocity calculations and by differencing two forward solutions show that the reciprocity method is accurate and can be used to decrease the number of calculations required to obtain sensitivities by one to two orders of magnitude.     

Crustal structure of the Ruby Mountains and southern Carlin Trend region, Nevada, from magnetotelluric data

We have collected about 150 magnetotelluric (MT) soundings in northeastern Nevada in the region of the Ruby Mountains metamorphic core complex uplift and southern Carlin mineral trend, in an effort to illuminate controls on core complex evolution and deposition of world-class gold deposits. The region has experienced a broad range of tectonic events including several periods of compressional and extensional deformation, which have contributed to the total expression of electrical resistivity. Most of the soundings reside in three east–west profiles across increasing degrees of core uplift to the north (Bald Mountain, Harrison Pass, and Secret Pass latitudes). One short cross-line was also taken to assess an east–west structure to the north of the northern profile. Model resistivity cross-sections were derived from the MT data using a 2-D inversion algorithm, which damps departures of model parameters from an a priori structure. Geological interpretation of the resistivity combines previous seismic, potential field and isotope models, structural and petrological models for regional compression and extension, and detailed structural/stratigraphic interpretations incorporating drilling for petroleum and mineral exploration. To first order, the resistivity structure is one of a moderately conductive, Phanerozoic sedimentary section fundamentally disrupted by intrusion and uplift of resistive crystalline rocks. Late Devonian and early Mississippian shales of the Pilot and Chainman Formations together form an important conductive marker sequence in the stratigraphy and show pronounced increases in conductance (conductivity–thickness product) from east to west. These increases are attributed to graphitization caused by Elko–Sevier era compressional shear deformation and possibly by intrusive heating. The resistive crystalline central massifs adjoin the host stratigraphy across crustal-scale, steeply dipping fault zones. The zones provide pathways to the lower crust for heterogeneous, upper crustal induced, electric current flow. Resistive core complex crust appears steeply bounded under the middle of the neighboring grabens and not to deepen at a shallow angle to arbitrary distances to the west. The numerous crustal breaks imaged with MT may contribute to the low effective elastic thickness (T_e) estimated regionally for the Great Basin and exemplify the mid-crustal, steeply dipping slip zones in which major earthquakes nucleate. An east–west oriented conductor in the crystalline upper crust spans the East Humboldt Range and northern Ruby Mountains. The conductor may be related to nearby graphitic metasediments, with possible alteration by middle Tertiary magmatism. Lower crustal resistivity everywhere under the profiles is low and appears quasi one-dimensional. It is consistent with a low rock porosity (<1 vol.%) containing hypersaline brines and possible water-undersaturated crustal melts, residual to the mostly Miocene regional extension. The resistivity expression of the southern Carlin Trend (CT) in the Pinon Range is not a simple lineament but rather a family of structures attributed to Eocene intrusion, stratal deformation, and alteration/graphitization. Substantial reactivation or overprinting by core complex uplift or Basin–Range extensional events seems likely. We concur with others that the Carlin Trend may result in part from overlap of the large Eocene Northeast Nevada Volcanic Field with Precambrian–Paleozoic deep-water clastic source rocks thickening abruptly to the west of the Pinon Range, and projecting to the north–northwest.     

Testing a model for predicting the timing and location of shallow landslide initiation in soil‐mantled landscapes

The growing availability of digital topographic data and the increased reliability of precipitation forecasts invite modelling efforts to predict the timing and location of shallow landslides in hilly and mountainous areas in order to reduce risk to an ever‐expanding human population. Here, we exploit a rare data set to develop and test such a model. In a 1·7 km² catchment a near‐annual aerial photographic coverage records just three single storm events over a 45 year period that produced multiple landslides. Such data enable us to test model performance by running the entire rainfall time series and determine whether just those three storms are correctly detected. To do this, we link a dynamic and spatially distributed shallow subsurface runoff model (similar to TOPMODEL) to an in?nite slope model to predict the spatial distribution of shallow landsliding. The spatial distribution of soil depth, a strong control on local landsliding, is predicted from a process‐based model. Because of its common availability, daily rainfall data were used to drive the model. Topographic data were derived from digitized 1 : 24 000 US Geological Survey contour maps. Analysis of the landslides shows that 97 occurred in 1955, 37 in 1982 and ?ve in 1998, although the heaviest rainfall was in 1982. Furthermore, intensity–duration analysis of available daily and hourly rainfall from the closest raingauges does not discriminate those three storms from others that did not generate failures. We explore the question of whether a mechanistic modelling approach is better able to identify landslide‐producing storms. Landslide and soil production parameters were ?xed from studies elsewhere. Four hydrologic parameters characterizing the saturated hydraulic conductivity of the soil and underlying bedrock and its decline with depth were ?rst calibrated on the 1955 landslide record. Success was characterized as the most number of actual landslides predicted with the least amount of total area predicted to be unstable. Because landslide area was consistently overpredicted, a threshold catchment area of predicted slope instability was used to de?ne whether a rainstorm was a signi?cant landslide producer. Many combinations of the four hydrological parameters performed equally well for the 1955 event, but only one combination successfully identi?ed the 1982 storm as the only landslide‐producing storm during the period 1980–86. Application of this parameter combination to the entire 45 year record successfully identi?ed the three events, but also predicted that two other landslide‐producing events should have occurred. This performance is signi?cantly better than the empirical intensity–duration threshold approach, but requires considerable calibration effort. Overprediction of instability, both for storms that produced landslides and for non‐producing storms, appears to arise from at least four causes: (1) coarse rainfall data time scale and inability to document short rainfall bursts and predict pressure wave response; (2) absence of local rainfall data; (3) legacy effect of previous landslides; and (4) inaccurate topographic and soil property data. Greater resolution of spatial and rainfall data, as well as topographic data, coupled with systematic documentation of landslides to create time series to test models, should lead to signi?cant improvements in shallow landslides forecasting. Copyright © 2003 John Wiley & Sons, Ltd.     

One-dimensional checkerboards and blending heights

We analyze the checkerboard problem of many alternating surfaces with different properties, on scales up to (say) 3,000 m. Power-law representations of the vertical profiles of mean wind speed and eddy diffusivity lead to solutions in terms of Kelvin and trigonometric functions.These solutions are used to determine blending heights (_*), where deviations from the mean of concentration, or of vertical flux density, fall to some small fraction, , of their value at the surface. Values of _*are important for regional and larger-scale meteorological models. In smaller scale micrometeorological studies, they may serve also as the top levels of surface boundary layers.An important result for both theoretical and experimental contexts is that deviations of flux persist with elevation much more strongly than those of concentration, so that, in general, _* should be based on flux rather than concentration. Representative values of _*, for = 0.05, are of order 5 and 30 m for surface pattern wavelengths of 10² and 10³ m, respectively. Values of _* are robust to changes in adopted power-law indices, and are independent of wind speed. Surface roughness has a mild but calculable effect.     

The Venus ionosphere at grazing incidence of solar radiation: Transport of plasma to the night ionosphere

Using a quasi-two-dimensional model of the Venus ionosphere, we calculated the ion number densities and horizontal ion bulk velocities expected for a range of solar zenith angles near the terminator (80 to 100°), and compared them with data obtained from the Pioneer Venus Orbiter retarding potential analyzer. The calculated ion bulk velocity arises entirely from the solar EUV-induced plasma pressure gradient and has a magnitude consistent with observations; ionization by suprathermal electrons is neglected in those computations. We find that while photoionization is the dominant source of ionospheric plasma for solar zenith angles less than 92°, plasma transport from the dayside is the dominant plasma source for solar zenith angles greater than 95°. We also show that the main nightside plasma peak at approximately 140 km altitude is of the F₂ type (i.e., is diffusion controlled). Its altitude and shape are thus quite insensitive to the altitude of the ion source.     

Climatic and volcanic forcing revealed in a 50,000-year diatom record from Lake Massoko, Tanzania

The interactions between climatic and volcanic forcing on diatom communities contained in a 50,000-year sedimentary sequence from Lake Massoko, Tanzania, were examined. At the century scale, 19 discrete tephra inputs to the lake isolated the sedimentary nutrient supply and shifted the diatom communities to those tolerant of low phosphorus levels, whereas at the millennial scale, diatom-inferred shifts in precipitation–evaporation based on conductivity optima and diatom life-form ratios were broadly similar to lake-level reconstructions from Lake Rukwa, Lake Malawi, and others in the region. Some fluctuations of Lake Massoko are consistent with the precession-driven changes in insolation, but the major climate shifts do not relate directly to orbital forcing of summer insolation south of the equator and show more consistency with records from the equatorial and northern tropics that receive rainfall from the passing of the intertropical convergence zone. Sea surface temperatures are strongly correlated to multimillennial-scale climate patterns over this region of Africa.     

A note on the effects of a thin visco-elastic mud layer on small amplitude water-wave propagation

In this note we investigated the effects of a thin visco-elastic mud layer on wave propagation. Within the framework of linear water-wave theory, analytical solutions are obtained for damping rate, dispersion relation between wave frequency and wave number, and velocity components in the water column and mud layer. The wave attenuation rate reaches a maximum value when the mud layer thickness is about the same as the mud boundary layer thickness. Heavier mud has a weaker effect on the wave damping. However, the wave attenuation rate does not always decrease as the elastic shear modulus increases. In the range of small values for elastic shear modulus, the wave attenuation can be amplified quite significantly. The current solutions are compared with experimental data with different wave conditions and mud properties. In general, good agreements are observed.     

Evaluating a modified point-based method to downscale cell-based climate variable data to high-resolution grids

To address the demand for high spatial resolution gridded climate data, we have advanced the Daymet point-based interpolation algorithm for downscaling global, coarsely gridded data with additional output variables. The updated algorithm, High-Resolution Climate Downscaler (HRCD), performs very good downscaling of daily, global, historical reanalysis data from 1° input resolution to 2.5 arcmin output resolution for day length, downward longwave radiation, pressure, maximum and minimum temperature, and vapor pressure deficit. It gives good results for monthly and yearly cumulative precipitation and fair results for wind speed distributions and modeled downward shortwave radiation. Over complex terrain, 2.5 arcmin resolution is likely too low and aggregating it up to 15 arcmin preserves accuracy. HRCD performs comparably to existing daily and monthly US datasets but with a global extent for nine daily climate variables spanning 1948–2006. Furthermore, HRCD can readily be applied to other gridded climate datasets.