Medieval Warming, Little Ice Age, and European impact on the environment during the last millennium in the lower Hudson Valley, New York, USA

Establishing natural climate variability becomes particularly important in large urban areas in anticipation of droughts. We present a well-dated bi-decadal record of vegetation, climate, land use, and fire frequency from a tidal marsh in the Hudson River Estuary. The classic Medieval Warm Period is evident through striking increases in charcoal and Pinus dominance from 800–1300 A.D., paralleling paleorecords southward along the Atlantic seaboard. Higher inputs of inorganic sediment during this interval suggest increased watershed erosion during drought conditions. The presence of the Little Ice Age ensues with increases in Picea and Tsuga, coupled with increasing organic percentages due to cooler, moister conditions. European impact is manifested by a decline in arboreal pollen due to land clearance, increased weedy plant cover (i.e., Ambrosia, Plantago, and Rumex), and an increase in inorganic particles to the watershed.     

The pressure dependence of the zirconium-in-rutile thermometer

The solubility of ZrO₂ in rutile is strongly temperature-dependent and has been identified as a potentially powerful thermometer when the rutile coexists with an appropriate buffer assemblage, e.g. zircon + quartz. In combination with experimental data at 10 kbar, previous consideration of data on natural rutile has not identified a pressure dependence for the thermometer. However, the expected volume change as a result of substitution of the larger Zr⁴⁺ cation for Ti⁴⁺ suggests that the Zr content of rutile should decrease with increasing pressure. To investigate the pressure dependence of the thermometer, piston cylinder (at 10, 20 & 30 kbar) and 1 atm furnace experiments were performed in the system ZrO₂-TiO₂-SiO₂. The solubility of ZrO₂ in rutile, in the presence of zircon and quartz was reversed at each pressure value. From these experiments, the thermodynamics of the end-member reaction ZrSiO₄ = SiO₂ + ZrO₂ (in rutile) have been determined. There is a secondary pressure effect accompanying the primary temperature dependence of the Zr content of rutile. New thermometer equations are, in the α -quartz field: in the β -quartz field and in the coesite field in which φ is ppm Zr, P is in kbar and R is the gas constant, 0.0083144 kJ K⁻¹. Thermometric results using these equations are shown for a range of geological settings.     

The effects of sub-blocking temperature metamorphism on the K/Ar systematics of hornblendes: 40Ar/39Ar dating of polymetamorphic garnet amphibolite from the Franciscan Complex,California

The effects of sub-blocking temperature metamorphism on the K/Ar system in hornblende, as revealed by ⁴⁰Ar/³⁹Ar release spectra, have been studied in a polymetamorphic knocker from the Franciscan Complex (FC), California. A primary amphibolite assemblage of horn-blende +rutile+epidote+apatite±garnet±sphene is variably overprinted by a blueschist facies assemblage of blue amphibole +lawsonite+chlorite+white mica+pumpellyite±sphene. The secondary assemblage formed at a temperature of 370° C, below that at which rapid Ar diffusion is expected in hornblende. Hornblendes from three, variably-altered samples of garnet amphibolite yield total gas ages of 147 to 161 Ma, but the corresponding plateau ages of 163.0±2.8, 160.6±2.2, and 161.8±2.2 Ma are identical within error. Hornblende separates with lower total gas ages come from more highly overprinted rocks, have excess K compared to that expected on the basis of electron probe analyses, and exhibit anomalously high K/Ca ratios in the low-temperature fractions of their incremental heating spectra. The reduced total gas ages result from the presence of thin (2 m) sheets of younger white mica in hornblendes from the moderately and highly altered amphibolites. The secondary micas are difficult to detect because of their small size and low abundance (2%), but because their K content is 50 to 100 times that of the host hornblende, they contribute significantly to the K and Ar budgets of the sample. The mica intergrowths are not removed by normal sample preparation, but because the mica inclusions degas at lower temperatures than hornblende during vacuum extraction, incremental heating analyses can provide precise cooling ages for the hornblendes as well as useful estimates of the age of the mica inclusions. The hornblende separate from the most altered sample contained 20±10% younger blue amphibole replacing hornblende, but its plateau age was not significantly affected. This is consistent with a replacement process in which K and Ar loss from the hornblende are coupled, leaving the K/Ar system undisturbed in relict primary grains. The K and Ar budgets of the sample are not strongly affected by the blue amphibole because of its very low K content. Because partial replacement of primary amphibole by high-K phyllosilicates occurs in many geological environments, effects like those described here could be widespread.     

Analysis of the spatial distribution of gamma-ray bursts in their host galaxies

We compare the radial distributions of known localized gamma-ray bursts (GRBs) relative to the centers of their host galaxies with the distributions of known objects in nearby galaxies (supernovae of various types, X-ray binaries), the hypothetical dark-matter profiles, and the distribution of luminous matter in galaxies in the model of an exponential disk. By comparing the moments of empirical distributions, we show that the radial distribution of GRBs in galaxies differs significantly from that of other sources. We suggest a new statistical method for comparing empirical samples that is based on estimating the number of objects within a given radius. The exponential disk profile was found to be in best agreement with the radial distribution of GRBs. The distribution of GRBs relative to the centers of their host galaxies also agrees with the dark matter profile at certain model parameters.     

A trajectory-oriented carriageway-based road network data model, part 1: Background

Introduction Withthedramaticdevelopmentofpositioning andtelecommunicationtechnologies,lotsoftraj ectorydataofmovingvehiclescanbecalculated,collected,andtransferred.Inthemeanwhile,thisdevelopmentandtheavailabilityoftrajectory datahavealsomotivatedresearche…     

Tritium in the Dead Sea

Tritium data in the Dead Sea for the period 1960–1979 are given. Tritium levels have increased until 1965 in the upper layers of the Dead Sea reaching a level of 170 TU, in response to the atmospheric buildup of tritium from thermonuclear testing. The levels have been decreasing ever since, both because of rapidly declining atmospheric concentrations of tritium and because of mixing of the surface layers with tritium deficient, deeper water masses. The depth of penetration of the tracer delineated the depth of meromictic stratification and successfully monitored the deepening of the pycnocline, until the overturn in 1979 homogenised the entire tritium profile. Modelling the changing tritium inventory over this period showed the predominance of the direct exchange across the air/sea interface, both in the buildup of tritium in the lake and also in its subsequent removal from it. The good fit between calculated and measured tritium inventories confirmed the evaporation estimate of 1.46 m/yr (the mean value for the period) with a precision unattained by other methods.     

Tourmaline in Proterozoic Massive Sulfide Deposits from Rajasthan, India

We have analyzed the chemical composition and boron isotope composition of tourmaline from tourmalinites, granite and a quartz-tourmaline vein from the Deri ore zone and from a pegmatitic band in the Rampura-Agucha ore body. These two Proterozoic massive sulfide deposits occur in the Aravalli-Delhi orogenic belt, Rajasthan, northwest India. Tourmaline from stratiform tourmalinites closely associated with the massive sulfides in the Deri deposit have preserved their original chemical compositions despite regional and thermal metamorphism in the area. These tourmalines have low Fe/(Fe + Mg) ratios (0.19–0.30; mean 0.26) that suggest formation close to the sediment-sea water interface. The δ¹¹B values (−15.5 and −16.4‰) are compatible with boron derived from leaching of argillaceous sediments and/or felsic volcanics underlying the original massive sulfide deposit during its formation. Boron isotope compositions measured in tourmaline from a post-ore granite and quartz-tourmaline vein in the Deri deposit indicate that boron in these tourmalines was derived from the tourmalinites produced during ore formation. The boron isotope systematics of a coarse brown tourmaline crystal from a pegmatitic band on the hanging wall contact of the Rampura-Agucha deposit indicate that 45 ± 25% of the boron within the original tourmaline was lost during upper amphibolite facies regional metamorphism. Received: 3 April 1996 / Accepted: 11 April 1996     

Efficient methods for solving water flow in variably saturated soils under prescribed flux infiltration

The non-linear solvers in numerical solutions of water flow in variably saturated soils are prone to convergence difficulties. Many aspects can give rise to such difficulties and in this paper we address the gravity term and the prescribed-flux boundary in the Picard iteration. The problem of the gravity term in the Picard iteration is iteration-to-iteration oscillation as the gravity term is treated, by analogy with the time-step advance technique, ‘explicitly’ in the iteration. The proposed method for the gravity term is an improvement of the ‘implicit’ approach of Zhang and Ewen [Water Resour. Res. 36 (2000) 2777] by extending it to heterogeneous soil and approximating the inter-nodal hydraulic conductivity in the diffusive term and the gravity term with the same scheme. The prescribed-flux boundary in traditional methods also gives rise to iteration-to-iteration oscillation because there is no feedback to the flux in the solution at the new iteration. To reduce such oscillation, a new method is proposed to provide such a feedback to the flux. Comparison with traditional Picard and Newton iteration methods for a wide range of problems show that a combination of these two proposed methods greatly improves the stability and consequently the computational efficiency, making the use of small time step and/or under-relaxation solely for convergence unnecessary.     

Lunar gravity field determination using SELENE same-beam differential VLBI tracking data

A lunar gravity field model up to degree and order 100 in spherical harmonics, named SGM100i, has been determined from SELENE and historical tracking data, with an emphasis on using same-beam S-band differential VLBI data obtained in the SELENE mission between January 2008 and February 2009. Orbit consistency throughout the entire mission period of SELENE as determined from orbit overlaps for the two sub-satellites of SELENE involved in the VLBI tracking improved consistently from several hundreds of metres to several tens of metres by including differential VLBI data. Through orbits that are better determined, the gravity field model is also improved by including these data. Orbit determination performance for the new model shows improvements over earlier 100th degree and order models, especially for edge-on orbits over the deep far side. Lunar Prospector orbit determination shows an improvement of orbit consistency from 1-day predictions for 2-day arcs of 6 m in a total sense, with most improvement in the along and cross-track directions. Data fit for the types and satellites involved is also improved. Formal errors for the lower degrees are smaller, and the new model also shows increased correlations with topography over the far side. The estimated value for the lunar GM for this model equals 4902.80080±0.0009 km³/s² (10 sigma). The lunar degree 2 potential Love number k ₂ was also estimated, and has a value of 0.0255 ± 0.0016 (10 sigma as well).     

Seeding agent dispersion within convective cloud as simulated by a 3-D numerical model

Summary The three-dimensional mesoscale cloud-resolving model ARPS (Advanced Regional Prediction System) was used to investigate the dispersal of an inert seeding agent within a cumulonimbus (Cb) cloud developing from two different initial states. In this paper, we stress the influence vortices in the cloud have on seeding agent dispersion. If a strong directional ambient wind shear is present in the lowest layer, a vortex pair formed at the flanks of the simulated cloud. Following the velocity field, a considerable amount of the injected seeding agent would be thrown out to the rear of the cloud, where both updrafts associated with vortices and downdrafts occurred. After a short time the agent was present only in the cloud periphery. If the Cb cloud developed under conditions where directional ambient wind shear did not exist, seeding agent dispersion would be quite different. In this case, almost all the seeding agent was transported into the main updraft, while the residence time of the agent within the cloud was longer due to the weaker cloud dynamics. Therefore, we must pay attention to whether or not the cloud contains vortices when we make the decision where to seed. This is necessary in order to minimize the loss of seeding material.