R-script modified VS-Lite model

Drought is a key factor affecting forest ecosystem processes at different spatio-temporal scales. For accurately modeling tree functioning ? and thus for producing reliable simulations of forest dynamics ? the consideration of the variability in the timing and extent of drought effects on tree growth is essential, particularly in strongly seasonal climates such as in the Mediterranean area. Yet, most dynamic vegetation models (DVMs) do not include this intra-annual variability of drought effects on tree growth. We present a novel approach for linking tree-ring data to drought simulations in DVMs. A modified forward model of tree-ring width (VS-Lite) was used to estimate seasonal- and site-specific growth responses to drought of Scots pine (Pinus sylvestris L.), which were subsequently implemented in the DVM ForClim. Ring-width data from sixteen sites along a moisture gradient from Central Spain to the Swiss Alps, including the dry inner Alpine valleys, were used to calibrate the forward ring-width model, and inventory data from managed Scots pine stands were used to evaluate ForClim performance. The modified VS-Lite accurately estimated the year-to-year variability in ring-width indices and produced realistic intra-annual growth responses to soil drought, showing a stronger relationship between growth and drought in spring than in the other seasons and thus capturing the strategy of Scots pine to cope with drought. The ForClim version including seasonal variability in growth responses to drought showed improved predictions of stand basal area and stem number, indicating the need to consider intra-annual differences in climate-growth relationships in DVMs when simulating forest dynamics. Forward modeling of ring-width growth may be a powerful tool to calibrate growth functions in DVMs that aim to simulate forest properties in across multiple environments at large spatial scales.

(Table 3) Air temperatures, and snow and ice characteristics of fresh and brackish water lakes in the Northwest Territories, Canada

The algorithms designed to estimate snow water equivalent (SWE) using passive microwave measurements falter in lake-rich high-latitude environments due to the emission properties of ice covered lakes on low frequency measurements. Microwave emission models have been used to simulate brightness temperatures (Tbs) for snowpack characteristics in terrestrial environments but cannot be applied to snow on lakes because of the differing subsurface emissivities and scattering matrices present in ice. This paper examines the performance of a modified version of the Helsinki University of Technology (HUT) snow emission model that incorporates microwave emission from lake ice and sub-ice water. Inputs to the HUT model include measurements collected over brackish and freshwater lakes north of Inuvik, Northwest Territories, Canada in April 2008, consisting of snowpack (depth, density, and snow water equivalent) and lake ice (thickness and ice type). Coincident airborne radiometer measurements at a resolution of 80x100 m were used as ground-truth to evaluate the simulations. The results indicate that subsurface media are simulated best when utilizing a modeled effective grain size and a 1 mm RMS surface roughness at the ice/water interface compared to using measured grain size and a flat Fresnel reflective surface as input. Simulations at 37 GHz (vertical polarization) produce the best results compared to airborne Tbs, with a Root Mean Square Error (RMSE) of 6.2 K and 7.9 K, as well as Mean Bias Errors (MBEs) of -8.4 K and -8.8 K for brackish and freshwater sites respectively. Freshwater simulations at 6.9 and 19 GHz H exhibited low RMSE (10.53 and 6.15 K respectively) and MBE (-5.37 and 8.36 K respectively) but did not accurately simulate Tb variability (R= -0.15 and 0.01 respectively). Over brackish water, 6.9 GHz simulations had poor agreement with airborne Tbs, while 19 GHz V exhibited a low RMSE (6.15 K), MBE (-4.52 K) and improved relative agreement to airborne measurements (R = 0.47). Salinity considerations reduced 6.9 GHz errors substantially, with a drop in RMSE from 51.48 K and 57.18 K for H and V polarizations respectively, to 26.2 K and 31.6 K, although Tb variability was not well simulated. With best results at 37 GHz, HUT simulations exhibit the potential to track Tb evolution, and therefore SWE through the winter season.

Pore water methane characteristics and sulphate reduction rates of ODP Leg 164 sediments

Anaerobic methane oxidation (AMO) was characterized in sediment cores from the Blake Ridge collected during Ocean Drilling Program (ODP) Leg 164. Three independent lines of evidence support the occurrence and scale of AMO at Sites 994 and 995. First, concentration depth profiles of methane from Hole 995B exhibit a region of upward concavity suggestive of methane consumption. Diagenetic modeling of the concentration profile indicates a 1.85-m-thick zone of AMO centered at 21.22 mbsf, with a peak rate of 12.4 nM/d. Second, subsurface maxima in tracer-based sulfate reduction rates from Holes 994B and 995B were observed at depths that coincide with the model-predicted AMO zone. The subsurface zone of sulfate reduction was 2 m thick and had a depth integrated rate that compared favorably to that of AMO (1.3 vs. 1.1 nmol/cm**2/d, respectively). These features suggest close coupling of AMO and sulfate reduction in the Blake Ridge sediments. Third, measured d13CH4 values are lightest at the point of peak model-predicted methane oxidation and become increasingly 13C-enriched with decreasing sediment depth, consistent with kinetic isotope fractionation during bacterially mediated methane oxidation. The isotopic data predict a somewhat (60 cm) shallower maximum depth of methane oxidation than do the model and sulfate reduction data.

Petrology and geochemistry at DSDP Site 89-462

The 16 samples of Deep Sea Drilling Project (DSDP) Leg 89 basalts that we analyzed for whole rock major and trace elements and for mineralogic compositions are identical to some of the basalts recovered during Leg 61. Leg 89 samples are mostly olivine-plagioclase-clinopyroxene sparsely phyric basalts and exhibit a wide variety of textures. These basalts have lower TiO2 at a given Mg/(Mg+Fe2+)*100 than MORB (midocean ridge basalt). We recognize three major chemical types of basalts in the Nauru Basin. We believe that different degrees of partial melting, modified by fractional crystallization and possibly by magma mixing at shallow depths, can explain the chemical differences among the three groups. This petrogenetic model is consistent with the observed downhole chemical-chronostratigraphic relations of the samples. New 87Sr/86Sr and U3Nd/144Nd analyses of basalt samples from DSDP Site 462 indicate that the Nauru Basin igneous complex is within the Sr-Nd isotopic range of ocean island basalt. Thus the Nauru Basin igneous complex resembles MORB in many aspects of its chemistry, morphology, and secondary alteration patterns (Larson, Schlanger, et al., 1981), but not in its isotopic characteristics. If it were not for the unambiguous evidence that the Nauru Basin complex was erupted off-ridge, the complex could easily be interpreted as normal oceanic layer 2. For this reason, we speculate that the Nauru Basin igneous complex was produced in an oceanic riftlike environment when multiple, fast-propagating rifts were formed during the fast seafloor spreading episode in the Cretaceous.