Fluid infiltration during contact metamorphism of interbedded marble and calc-silicate hornfels, Twin Lakes area, central Sierra Nevada, California

Abstract In the Twin Lakes area, central Sierra Nevada, California, most contact metamorphosed marbles contain calcite + dolomite + forsterite ± diopside ± phlogopite ± tremolite, and most calc-silicate hornfelses contain calcite + diopside + wollastonite + quartz ± anorthite ± K-feldspar ± grossular ± titanite. Mineral-fluid equilibria involving calcite + dolomite + tremolite + diopside + forsterite in two marble samples and wollastonite + anorthite + quartz + grossular in three hornfels samples record P± 3 kbar and T± 630° C. Various isobaric univariant assemblages record CO₂-H₂O fluid compositions of χCO₂= 0.61–0.74 in the marbles and χCO₂= 0.11 in the hornfelses. Assuming a siliceous dolomitic limestone protolith consisting of dolomite + quartz ° Calcite ± K-feldspar ± muscovite ± rutile, all plausible prograde reaction pathways were deduced for marble and hornfels on isobaric T-XCO₂ diagrams in the model system K₂O-CaO-MgO-Al₂O₃-SiO₂-H₂O-CO₂. Progress of the prograde reactions was estimated from measured modes and mass-balance calculations. Time-integrated fluxes of reactive fluid which infiltrated samples were computed for a temperature gradient of 150 °C/km along the fluid flow path, calculated fluid compositions, and estimated reaction progress using the mass-continuity equation. Marbles and hornfelses record values in the range 0.1–3.6 × 10⁴ cm³/cm² and 4.8–12.9 × 10⁴ cm³/cm², respectively. For an estimated duration of metamorphism of 10⁵ years, average in situ metamorphic rock permeabilities, calculated from Darcy's Law, are 0.1–8 × 10^?6 D in the marbles and 10–27 × 10^?6 D in the hornfelses. Reactive metamorphic fluids flowed up-temperature, and were preferentially channellized in hornfelses relative to the marbles. These results appear to give a general characterization of hydrothermal activity during contact metamorphism of small pendants and screens (dimensions ± 1 km or less) associated with emplacement of the Sierra Nevada batholith.     

A laboratory test of the soil chemical submodels of two models of catchment acidification

Parameters for ion exchange selectivity and aluminium hydroxide dissolution in the soil chemical submodels used in applications of the Birkenes model and of MAGIC are compared and several discrepancies identified for organic soils. A laboratory column simulation of the soil chemical submodels is proposed and applied to soils from the Loch Dee area in Galloway. Experimental results were well predicted by a simplified version of MAGIC, with ion exchange selectivity parameters similar to those used in a previous simulation of one subcatchment of Loch Dee. The aluminium hydroxide dissolution parameter used previously was found to be too low for the organic soil materials, where a value of 106 predicted the experimental results more closely. The model developed also included a simple silicate weathering reaction to release base cations into the system. It is concluded that such simple laboratory simulations are useful for independent calibration of the soil chemical submodel of catchment models.     

Palaeotemperature reconstructions of the European permafrost zone during marine oxygen isotope Stage 3 compared with climate model results

A palaeotemperature reconstruction based on periglacial phenomena in Europe north of approximately 51 °N, is compared with high‐resolution regional climate model simulations of the marine oxygen isotope Stage 3 (Stage 3) palaeoclimate. The experiments represent Stage 3 warm (interstadial), Stage 3 cold (stadial) and Last Glacial Maximum climatic conditions. The palaeotemperature reconstruction deviates considerably for the Stage 3 cold climate experiments, with mismatches up to 11 °C for the mean annual air temperature and up to 15 °C for the winter temperature. However, in this reconstruction various factors linking climate and permafrost have not been taken into account. In particular a relatively thin snow cover and high climatic variability of the glacial climate could have influenced temperature limits for ice‐wedge growth. Based on modelling the 0 °C mean annual ground temperature proves to be an appropriate upper temperature limit. Using this limit, mismatches with the Stage 3 cold climate experiments have been reduced but still remain. We therefore assume that the Stage 3 ice wedges were generated during short (decadal time‐scale) intervals of extreme cold climate, below the mean temperatures indicated by the Stage 3 cold climate model simulations. Copyright © 2003 John Wiley & Sons, Ltd.     

A THREE-DIMENSIONAL ELASTIC NESTED-GRID MESO-(β-γ)SCALE ATMOSPHERIC MODEL——PARTI:MODEL DESCRIPTION

A three-dimensional elastic nonhydrostatic mesoscale(β-γ)model with nested-grid is presented.It uses a set of fullequations in terrain-following coordinates as its basic dynamic frame,which is solved with a time-splitting algorithmfor acoustic and gravity waves.The model physical parameterization includes a K-theory subgrid eddy mixing for cloudand free atmosphere,a bulk planetary boundary layer parameterization,and three types of sofisticated cloudmicrophysics schemes with double-parameters for hail-bearing clouds,warm clouds and snowing clouds respectively.The model is designed to be used flexibly for simulations of a variety of meso-and small-scale atmospheric processes,and can be improved as a regional and local operational NWP system in future.     

地理信息系统数据库设计

首先分析了地理信息系统的组成、地理信息系统数据的特征、地理信息系统数据库的功能，然后进行了中心式导航系统数据库的设计，给出了数据库模型的E-R图，最后描述了空间地理信息的查询功能。     

Sensitivity analysis of a deterministic water temperature model to forest canopy and soil temperature in Catamaran Brook (New Brunswick,Canada)

A coupled deterministic hydrological and water temperature model, CEQUEAU, was modified to include soil temperature and crown closure in its calculation of local advective terms in the heat budget. The modified model was than tested to verify its sensitivity to these modifications. An analysis of the heat budget of a small forested catchment in eastern Canada revealed that the advective term related to interflow plays a significant role in the daily water heat budget, providing on average 28% of the local advective budget (which also includes advective heat terms from surface runoff and groundwater) and nearly 14% of the total heat budget (which includes all radiative terms at the water surface, convection and evaporation, as well as the local advective terms). Relative sensitivity indices (RSIs) were used to verify the impact of the newly introduced parameters and variables. Among them, parameters related to the forest cover (crown closure and leaf area index) have a maximum RSI of ?0·6; i.e. a 100% increase in value produces a 60% decrease in the local advective term. Parameters with the greatest influence are the volume of water contributing to interflow and the amplitude of the net radiative flux at the soil surface, which, if doubled, would double the contribution of the local interflow advective term to the heat budget. Copyright © 2003 John Wiley & Sons, Ltd.     

Modelling the effect of low soil temperatures on transpiration by Scots pine

For ecosystem modelling of the Boreal forest it is important to include processes associated with low soil temperature during spring‐early summer, as these affect the tree water uptake. The COUP model, a physically based SVAT model, was tested with 2 years of soil and snow physical measurements and sap flow measurements in a 70‐year‐old Scots pine stand in the boreal zone of northern Sweden. During the first year the extent and duration of soil frost was manipulated in the field. The model was successful in reproducing the timing of the soil warming after the snowmelt and frost thaw. A delayed soil warming, into the growing season, severely reduced the transpiration. We demonstrated the potential for considerable overestimation of transpiration by the model if the reduction of the trees' capacity to transpire due to low soil temperatures is not taken into account. We also demonstrated that the accumulated effect of aboveground conditions could be included when simulating the relationship between soil temperature and tree water uptake. This improved the estimated transpiration for the control plot and when soil warming was delayed into the growing season. The study illustrates the need of including antecedent conditions on root growth in the model in order to catch these effects on transpiration. The COUP model is a promising tool for predicting transpiration in high‐latitude stands. Copyright © 2006 John Wiley & Sons, Ltd.     

Screening of polymers on selected Hawaii soils for erosion reduction and particle settling

In recent years, high‐molecular‐weight anionic polyacrylamides (PAMs) have been tested on a variety of soils, primarily in temperate climates. However, little information is available regarding the effectiveness of PAM for preventing soil loss through runoff in tropical settings. Screening tests were performed using three negatively charged PAMs and one positively charged PAM on five Hawaii soils (two Oxisols, one Vertisol, and two Aridisols) to determine erosion loss, sediment settling, and aggregate stability. A laboratory‐scale rainfall simulator was used to apply erosive rainfall at intensities from 5 to 8·5 cm h^?1 at various PAM doses applied in both dry and solution forms. Soil detachment due to splash and runoff, as well as the runoff and percolate water volumes, were measured for initial and successive storms. The impact of PAM on particle settling and aggregate stability was also evaluated for selected soil‐treatment combinations. Among the PAMs, Superfloc A‐836 was most effective, and significantly reduced runoff and splash sediment loss for the Wahiawa Oxisol and Pakini Andisol at rates varying between 10 and 50 kg ha^?1. Reduced runoff and splash sediment loss were also noted for PAM Aerotil‐D when applied in solution form to the Wahiawa Oxisol. Significant reductions in soil loss were not noted for either the Lualualei Vertisol or the Holomua Oxisol. It is believed that the high montmorillonite content of the Lualualei Vertisol and the low cation‐exchange capacity of the Holomua Oxisol diminished the effectiveness of the various PAMs tested. The polymers were also found to enhance sediment settling of all soils and helped improve their aggregate stability. This screening study shows the potential use of PAM for tropical soils for applications such as infiltration enhancement, runoff reduction, and enhanced sedimentation of detention ponds. Copyright © 2005 John Wiley & Sons, Ltd.     

Analytical solutions of one‐dimensional infiltration before and after ponding

This paper is a continuation of previous research, which obtained a convenient solution for arbitrary surface fluxes before ponding. By means of Fourier Transformation this has been extended to derive analytical solutions of a linearized Richards' equation for arbitrary input fluxes after surface saturation. Copyright © 2003 John Wiley & Sons, Ltd.