Ice‐cover variability on shallow lakes at high latitudes: model simulations and observations

A one‐dimensional thermodynamic model for simulating lake‐ice phenology is presented and evaluated. The model can be driven with observed daily or hourly atmospheric forcing of air temperature, relative humidity, wind speed, cloud amount and snowfall. In addition to computing the energy balance components, key model output includes the temperature profile at an arbitrary number of levels within the ice/snow (or the water temperature if there is no ice) and ice thickness (clear ice and snow‐ice) on a daily basis, as well as freeze‐up and break‐up dates. The lake‐ice model is used to simulate ice‐growth processes on shallow lakes in arctic, sub‐arctic, and high‐boreal forest environments. Model output is compared with field and remote sensing observations gathered over several ice seasons. Simulated ice thickness, including snow‐ice formation, compares favourably with field measurements. Ice‐on and ice‐off dates are also well simulated when compared with field and satellite observations, with a mean absolute difference of 2 days. Model simulations and observations illustrate the key role that snow cover plays on the seasonal evolution of ice thickness and the timing of spring break‐up. It is also shown that lake morphometry, depth in particular, is a determinant of ice‐off dates for shallow lakes at high latitudes. Copyright © 2003 John Wiley & Sons, Ltd.     

Phase Relations and Melting of Anhydrous K-bearing Eclogite from 1200 to 1600{degrees}C and 3 to 5 GPa

To investigate eclogite melting under mantle conditions, wehave performed a series of piston-cylinder experiments usinga homogeneous synthetic starting material (GA2) that is representativeof altered mid-ocean ridge basalt. Experiments were conductedat pressures of 3·0, 4·0 and 5·0 GPa andover a temperature range of 1200–1600°C. The subsolidusmineralogy of GA2 consists of garnet and clinopyroxene withminor quartz–coesite, rutile and feldspar. Solidus temperaturesare located at 1230°C at 3·0 GPa and 1300°C at5·0 GPa, giving a steep solidus slope of 30–40°C/GPa.Melting intervals are in excess of 200°C and increase withpressure up to 5·0 GPa. At 3·0 GPa feldspar, rutileand quartz are residual phases up to 40°C above the solidus,whereas at higher pressures feldspar and rutile are rapidlymelted out above the solidus. Garnet and clinopyroxene are theonly residual phases once melt fractions exceed 20% and garnetis the sole liquidus phase over the investigated pressure range.With increasing melt fraction garnet and clinopyroxene becomeprogressively more Mg-rich, whereas coexisting melts vary fromK-rich dacites at low degrees of melting to basaltic andesitesat high melt fractions. Increasing pressure tends to increasethe jadeite and Ca-eskolaite components in clinopyroxene andenhance the modal proportion of garnet at low melt fractions,which effects a marked reduction in the Al₂O₃ and Na₂O contentof the melt with pressure. In contrast, the TiO₂ and K₂O contentsof the low-degree melts increase with increasing pressure; thusNa₂O and K₂O behave in a contrasted manner as a function ofpressure. Altered oceanic basalt is an important component ofcrust returned to the mantle via plate subduction, so GA2 maybe representative of one of many different mafic lithologiespresent in the upper mantle. During upwelling of heterogeneousmantle domains, these mafic rock-types may undergo extensivemelting at great depths, because of their low solidus temperaturescompared with mantle peridotite. Melt batches may be highlyvariable in composition depending on the composition and degreeof melting of the source, the depth of melting, and the degreeof magma mixing. Some of the eclogite-derived melts may alsoreact with and refertilize surrounding peridotite, which itselfmay partially melt with further upwelling. Such complex magma-genesisconditions may partly explain the wide spectrum of primitivemagma compositions found within oceanic basalt suites. KEY WORDS: eclogite; experimental petrology; mafic magmatism; mantle melting; oceanic basalts     

Comparison of oscillatory and stationary flow through porous media

Experiments were conducted in an oscillatory water tunnel to investigate what effects temporal inertia has on the resistance of a granular medium. The flow law governing stationary porous media flow is reasonably well accepted and understood but the effects of unsteady flow have commonly been neglected. The present research was designed to assess the magnitude of the acceleration effects on media of uniformly packed spheres of equal diameter and on one sample of randomly placed stone.Oscillatory flow tests were made in a large oscillatory flume tunnel with periods varying from 3 to 12 seconds. The influence of properties of the medium (grain size and porosity) were tested by using spheres with two different diameters and packing each size sphere in different geometric arrangements. Tests made on a stone sample provided a qualitative assessment of the effects of more random material properties.For the experiments described in this paper, the Forchheimer unsteady-stationary flow law described the oscillatory measurements well when velocities and energy losses were maximum. Empirical coefficients determined from steady-stationary flow were generally found to apply to the unsteady flow, however some evidence of dependency on the period of oscillation was noted.     

Ultra-rapid landscape response and sediment yield following glacier retreat, Icy Bay, southern Alaska

The páramo is a neotropical alpine ecosystem that covers more than 75,000 km² of the northern Andes of Colombia, Ecuador, Venezuela, and Peru. It provides important environmental services: more than 10 million people in the Andean highlands benefit from the water supply and regulation function, which is attributed to the volcanic soils that underlie the ecosystem. The soils are also major carbon sinks of global significance. Severe land use changes and soil degradation threaten both the hydrology and carbon sink function. Nevertheless, soil genesis and properties in the páramo is rather poorly understood, nor are their ecological functions well documented. The impact of the geomorphology of the páramo on soil genesis was studied in the rio Paute basin, south Ecuador. Two toposequences were described and analysed. In each toposequence, four pedons were selected representing summit, backslope, undrained plain situation, and valley bottom positions in the landscape. The soils are classified as Hydric Andosols in the World Reference Base for Soil Resources and Epiaquands or Hydrudands in Soil Taxonomy. They are very acidic and have a high organic matter content, high P deficiency, and Al toxicity. Their water content ranges from 2.64 g g^− 1 at saturation, down to 1.24 g g^− 1 at wilting point, resulting in a large water storage capacity. Two major soil forming processes are identified: (1) volcanic ash deposition and (2) accumulation of organic carbon. Volcanic ash deposits may vary in depth as a result of regional geomorphological factors such as parent material, orientation, slope, and altitude. Organic carbon accumulation is an interaction of both waterlogging, which depends on the position in the landscape, and the formation of organometallic complexes with Al and Fe released during volcanic ash breakdown. Despite the high variability in parent material and topography, the soil is characterised by a notable homogeneity in physico-chemical properties. Statistical analysis reveals that only topographic location has a slight but significant influence on soil pH as well as the organic matter content, saturated conductivity and water retention at high pressure. Finally, the exceptional properties of these soils provide useful insights to improve classification of the Andosols reference group of the FAO World reference Base for Soil Resources.     

Parameter identification of inelastic structures under dynamic loads

The generalized model of differential hysteresis contains 13 control parameters with which it can curve‐fit practically any hysteretic trace. Three identification algorithms are developed to estimate the control parameters of hysteresis for different classes of inelastic structures. These algorithms are based upon the simplex, extended Kalman filter, and generalized reduced gradient methods. Novel techniques such as global search and internal constraints are incorporated to facilitate convergence and stability. Effectiveness of the devised algorithms is demonstrated through simulations of two inelastic systems with both pinching and degradation characteristics in their hysteretic traces. Owing to very modest computing requirements, these identification algorithms may become acceptable as a design tool for mapping the hysteretic traces of inelastic structures. Copyright © 2002 John Wiley & Sons, Ltd.     

Synclinal-horst basins: examples from the southern Rio Grande rift and southern transition zone of southwestern New Mexico, USA

In areas of broadly distributed extensional strain, the back‐tilted edges of a wider than normal horst block may create a synclinal‐horst basin. Three Neogene synclinal‐horst basins are described from the southern Rio Grande rift and southern Transition Zone of southwestern New Mexico, USA. The late Miocene–Quaternary Uvas Valley basin developed between two fault blocks that dip 6–8° toward one another. Containing a maximum of 200 m of sediment, the Uvas Valley basin has a nearly symmetrical distribution of sediment thickness and appears to have been hydrologically closed throughout its history. The Miocene Gila Wilderness synclinal‐horst basin is bordered on three sides by gently tilted (10°, 15°, 20°) fault blocks. Despite evidence of an axial drainage that may have exited the northern edge of the basin, 200–300 m of sediment accumulated in the basin, probably as a result of high sediment yields from the large, high‐relief catchments. The Jornada del Muerto synclinal‐horst basin is positioned between the east‐tilted Caballo and west‐tilted San Andres fault blocks. Despite uplift and probable tilting of the adjacent fault blocks in the latest Oligocene and Miocene time, sediment was transported off the horst and deposited in an adjacent basin to the south. Sediment only began to accumulate in the Jornada del Muerto basin in Pliocene and Quaternary time, when an east‐dipping normal fault along the axis of the syncline created a small half graben. Overall, synclinal‐horst basins are rare, because horsts wide enough to develop broad synclines are uncommon in extensional terrains. Synclinal‐horst basins may be most common along the margins of extensional terrains, where thicker, colder crust results in wider fault spacing.     

Petrography and uplift history of the Quaternary Takidani Granodiorite: could it have hosted a supercritical (HDR) geothermal reservoir?

The Quaternary Takidani Granodiorite (Japan Alps) is analogous to the type of deep-seated (3–5 km deep) intrusive-hosted fracture network system that might support (supercritical) hot dry/wet rock (HDR/HWR) energy extraction. The I-type Takidani Granodiorite comprises: porphyritic granodiorite, porphyritic granite, biotite-hornblende granodiorite, hornblende-biotite granodiorite, biotite-hornblende granite and biotite granite facies; the intrusion has a reverse chemical zonation, characterized by >70 wt% SiO₂ at its inferred margin and <67 wt% SiO₂ at the core. Fluid inclusion evidence indicates that fractured Takidani Granodiorite at one time hosted a liquid-dominated, convective hydrothermal system, with <380°C, low-salinity reservoir fluids at hydrostatic (mesothermal) pressure conditions. ‘Healed’ microfractures also trapped >600°C, hypersaline (35 wt% NaCl_eq) fluids of magmatic origin, with inferred minimum pressures of formation being 600–750 bar, which corresponds to fluid entrapment at 2.4–3.0 km depth. Al-in-hornblende geobarometry indicates that hornblende crystallization occurred at about 1.45 Ma (7.7–9.4 km depth) in the (marginal) eastern Takidani Granodiorite, but later (at 1.25 Ma) and shallower (6.5–7.0 km) near the core of the intrusion. The average rate of uplift across the Takidani Granodiorite from the time of hornblende crystallization has been 5.1–5.9 mm/yr (although uplift was about 7.5 mm/yr prior to 1.2 Ma), which is faster than average uplift rates in the Japan Alps (3 mm/yr during the last 2 million years). A temperature–depth–time window, when the Takidani Granodiorite had potential to host an HDR system, would have been when the internal temperature of the intrusive was cooling from 500°C to 400°C. Taking into account the initial (7.5 mm/yr) rate of uplift and effects of erosion, an optimal temperature–time–depth window is proposed: for 500°C at 1.54–1.57 Ma and 5.2±0.9 km (drilling) depth; and 400°C at 1.36–1.38 Ma and 3.3±0.8 km (drilling) depth, which is within the capabilities of modern drilling technologies, and similar to measured temperature–depth profiles in other active hydrothermal systems (e.g. at Kakkonda, Japan).     

Long term atmospheric deposition as the source of nitrate and other salts in the Atacama Desert, Chile: New evidence from mass-independent oxygen isotopic compositions

Isotopic analysis of nitrate and sulfate minerals from the nitrate ore fields of the Atacama Desert in northern Chile has shown anomalous ¹⁷O enrichments in both minerals. Δ¹⁷O values of 14-21 ‰ in nitrate and 0.4 to 4 ‰ in sulfate are the most positive found in terrestrial minerals to date. Modeling of atmospheric processes indicates that the Δ¹⁷O signatures are the result of photochemical reactions in the troposphere and stratosphere. We conclude that the bulk of the nitrate, sulfate and other soluble salts in some parts of the Atacama Desert must be the result of atmospheric deposition of particles produced by gas to particle conversion, with minor but varying amounts from sea spray and local terrestrial sources. Flux calculations indicate that the major salt deposits could have accumulated from atmospheric deposition in a period of 200,000 to 2.0 M years during hyper-arid conditions similar to those currently found in the Atacama Desert. Correlations between Δ¹⁷O and δ¹⁸O in nitrate salts from the Atacama Desert and Mojave Desert, California, indicate varying fractions of microbial and photochemical end-member sources. The photochemical nitrate isotope signature is well preserved in the driest surficial environments that are almost lifeless, whereas the microbial nitrate isotope signature becomes dominant rapidly with increasing moisture, biologic activity, and nitrogen cycling. These isotopic signatures have important implications for paleoclimate, astrobiology, and N cycling studies.     

Hydrothermal Quartz Vein Formation, Revealed by Coupled SEM-CL Imaging and Fluid Inclusion Microthermometry: Shuteen Complex, South Gobi, Mongolia

Abstract. Scanning electron microscopy-cathodoluminescence (SEM-CL) imaging of vein quartz in the Cu-mineralised, Shuteen Complex (South Gobi, Mongolia) has revealed a complex history of crystal growth, dissolution and microfracture healing, associated with several hydrothermal events that could not be detected using other observational techniques (e.g. transmitted/reflected light microscopy, back-scattered electron imaging, or secondary electron imaging).
The quartz initially grew as CL-bright/grey crystals in a 345±30C liquid reservoir, as inferred by the analysis of primary liquid fluid inclusions (average Th of 343C; 6.6∼7.7 wt% NaCl_eq). Quartz precipitation occurred at the edge of the crystals as reservoir fluids cooled to 260±25C, as indicated by micron-scale CL-dark/CL-bright quartz growth bands containing abundant fluid inclusions (with an average T_h values of 261C). Pressure fluctuations were the likely cause of dissolution, as SEM-CL imaging reveals the quartz have corroded or rounded crystal edges, and precipitation of later quartz into open space. SEM-CL imaging shows the quartz contains healed microfractures that trapped low salinity fluids (3.9 wt% NaC1_eq) with Th values of 173±15C.
SEM-CL imaging provides a means of deciphering the thermal and chemical evolution of the fossil Shuteen hydrothermal system, and the nature of hydrothermal quartz vein-forming processes, by facilitating the correlation of distinct fluid inclusion populations and their relative chronology, with specific hydrothermal events.