Two‐dimensional physically based finite element runoff model for small agricultural watersheds: I. Model development

Soil erosion by water is the root cause of ecological degradation in the Shiwalik foothills of Northern India. Simulation of runoff and its component processes is a pre‐requisite to develop the management strategies to tackle the problem, successfully. A two‐dimensional physically based distributed numerical model, ROMO2D has been developed to simulate runoff from small agricultural watersheds on an event basis. The model employs the 2‐D Richards equation with sink term to simulate infiltration and soil moisture dynamics in the vadoze zone under variable rainfall conditions, and 2‐D Saint‐Venant equations under the kinematic wave approximation along with Manning's equation as the stage‐discharge equation for runoff routing. The various flow‐governing equations have been solved numerically by employing a Galerkin finite element method for spatial discretization using quadrilateral elements and finite difference techniques for temporal solutions. The ROMO2D computer program has been developed as a class‐based program, coded in C + + in such a way that with minor modifications, the model can be used to simulate runoff on a continuous basis. The model writes output for a runoff hydrograph of each storm. Model development is described in this paper and the results of model testing and field application are to be presented in a subsequent paper. Copyright © 2008 John Wiley & Sons, Ltd.     

Influence of subdaily polar motion model on nutation offsets estimated by very long baseline interferometry

This paper studies the connection between the subdaily model for polar motion used in the processing of very long baseline interferometry (VLBI) observations and the estimated nutation offsets. By convention accepted by the International Earth Rotation Service, the subdaily model for polar motion recommended for routine processing of geodetic observations does not contain any daily retrograde terms due to their one-to-one correlation with the nutation. Nevertheless, for a 24-h VLBI solution a part of the signal contained in the polar motion given by the used subdaily model is numerically mistaken for a retrograde daily sidereal signal. This fictitious retrograde daily signal contributes to the estimated nutation, leading to systematic differences between the nutation offsets from VLBI solutions computed with different subdaily polar motion models. We demonstrate this effect using solutions for all suitable 24-h VLBI sessions over a time span of 11 years (2000–2011). By changing the amplitudes of one tidal term in the underlying subdaily model for polar motion and comparing the estimated parameters to the solutions computed with the unchanged subdaily model, the paper shows and explains theoretically the effects produced by the individual subdaily terms on the VLBI nutation estimates.     

Ra/Ra and Ra/Ba ratios to track barite formation and transport in the water column

We measured ²²⁸Ra_ex/²²⁶Ra_ex and ²²⁶Ra_ex/Ba_ex ratios in suspended and sinking particles collected at the Oceanic Flux Program (OFP) time-series site in the western Sargasso Sea and compared them to seawater ratios to provide information on the origin and transport of barite (BaSO₄) in the water column. The ²²⁸Ra_ex/²²⁶Ra_ex ratios of the suspended particles down to 2000 m are nearly identical to those of seawater at the same water depth. These ratios are much lower than expected if suspended barite was produced in surface waters and indicate that barite is produced throughout the mesopelagic layer. The ²²⁸Ra_ex/²²⁶Ra_ex activity ratios of sinking particles collected at 1500 and 3200 m varied mostly between 0.1 and 0.2, which is intermediate between the seawater ratio at these depths (<0.03) and the seawater ratios found in the upper 250 m (0.31-0.42). This suggests that excess Ba (i.e., Ba_ex = Ba_total − Ba_lithogenic), considered to be mainly barite, present in the sinking flux is a mixture of crystals formed recently in the upper water column, formed several years earlier in the upper water column, or formed recently in deeper waters. We observe a sizeable temporal variability in the ²²⁸Ra_ex/²²⁶Ra_ex ratios of sinking particles, which indicates temporal variability in the relative proportion of barite crystals originating from surface (with a high ²²⁸Ra_ex/²²⁶Ra_ex ratio) and mesopelagic (with a low ²²⁸Ra_ex/²²⁶Ra_ex ratio) sources. However, we could not discern a clear pattern that would elucidate the factors that control this variability. The ²²⁶Ra/Ba ratios measured in seawater are consistent with the value reported from the GEOSECS expeditions (2.3 dpm μmol⁻¹) below 500 m depth, but are significantly lower in the upper 500 m. High ²²⁶Ra_ex/Ba_ex ratios and elevated Sr concentrations in suspended particles from the upper water column suggest preferential uptake of ²²⁶Ra over Ba during formation of SrSO₄ skeletons by acantharians, which must contribute to barite formation in shallow waters. Deeper in the water column the ²²⁶Ra_ex/Ba_ex ratios of suspended particles are lower than those of seawater. Since ²²⁸Ra_ex/²²⁶Ra_ex ratios demonstrate that suspended barite at these depths has been produced recently and in situ, their low ²²⁶Ra_ex/Ba_ex ratios indicate preferential uptake of Ba over Ra in barite formed in mesopelagic water.     

Transmission electron microscope texture and crystal chemistry of coexisting ortho- and clinopyroxene in the Antarctic ureilite Frontier Mountain 90054: Implications for thermal history

Abstract— Frontier Mountain (FRO) 90054, from Antarctica, is a rare clino- and orthopyroxene-bearing ureilite with a coarse equigranular oriented texture (grains up to 3 mm); it is classified as a low-shock Ca-rich type. The crystal chemistry of its clinopyroxene (Wo_39.3En_54.6Fs_6.1), orthopyroxene (En_{84 2}Fs₁₁Wo_4.8) and olivine (Fa_12.6Fo_86.9) was investigated by single-crystal x-ray structural refinements and transmission electron microscope (TEM) observations to obtain data on the evolutionary history of the parent body. The M1 octahedron and unit cell volumes of the orthopyroxene and clinopyroxene are consistent with low-pressure crystallization. The closure temperatures for intracrystalline Mg-Fe²⁺ ordering yielded values of 674 °C and 804 °C for opx and 596 °C for cpx, which indicate high-temperature equilibration and fast cooling. Trasmission electron microscope investigations were performed on clinopyroxene, orthopyroxene and pigeonite. The (100) twin lamellae in the clinopyroxene and intergrowth of clino- and orthoenstatite lamellae in orthopyroxene most probably originated by deformation. Exsolution was not observed in any of the phases, which suggests rapid cooling. Analysis by TEM also revealed interstitial Na-rich glass and pigeonite with sharp h + k odd reflections and rare stacking faults parallel to (100). Textural and crystal chemical data, obtained by TEM, indicated rapid cooling that was probably due to fast radiative heat loss as a result of the disintegration of the parent body into small fragments, which subsequently reassembled into a larger body. One or more collisional events caused fine-scale stacking faults and partial melting.     

Soil-water dynamics and tree water uptake in the Sacramento Mountains of New Mexico (USA): a stable isotope study

In the southwestern United States, precipitation in the high mountains is a primary source of groundwater recharge. Precipitation patterns, soil properties and vegetation largely control the rate and timing of groundwater recharge. The interactions between climate, soil and mountain vegetation thus have important implications for the groundwater supply. This study took place in the Sacramento Mountains, which is the recharge area for multiple regional aquifers in southern New Mexico. The stable isotopes of oxygen and hydrogen were used to determine whether infiltration of precipitation is homogeneously distributed in the soil or whether it is partitioned among soil-water ‘compartments’, from which trees extract water for transpiration as a function of the season. The results indicate that “immobile” or “slow” soil water, which is derived primarily from snowmelt, infiltrates soils in a relatively uniform fashion, filling small pores in the shallow soils. “Mobile” or “fast” soil water, which is mostly associated with summer thunderstorms, infiltrates very quickly through macropores and along preferential flow paths, evading evaporative loss. It was found that throughout the entire year, trees principally use immobile water derived from snowmelt mixed to differing degrees with seasonally available mobile-water sources. The replenishment of these different water pools in soils appears to depend on initial soil-water content, the manner in which the water was introduced to the soil (snowmelt versus intense thunderstorms), and the seasonal variability of the precipitation and evapotranspiration. These results have important implications for the effect of climate change on recharge mechanisms in the Sacramento Mountains.     

Micro-sheeting of granite and its relationship with landsliding specifically after the heavy rainstorm in June 1999, Hiroshima Prefecture, Japan

Certain types of granite in mountainous areas are microscopically sheeted to a depth of 50 m due to unloading under the stress field that reflects slope morphology. Micro-sheets generally strike parallel to major slope surfaces and gently dip downslope, forming cataclinal overdip slopes. The cataclinal overdip slope accelerates creep movement of micro-sheeted granite, which in turn loosens and disintegrates granite via the widening or neoformation of cracks, probably in combination with stress release, temperature change, and changes in water content near the ground surface. The surface portion of micro-sheeted granite is thus loosened with a well-defined basal front, which finally slides in response to heavy rain. Innumerable landslides of this type occurred in Hiroshima Prefecture, western Japan, following the heavy rainstorm of 29 June 1999. Following such landslides, the weathering of micro-sheeted granite exposed on the landslide scar recommences, setting the stage for future landslide.     

Estimates of source parameters of <Emphasis Type=Italic>M</Emphasis>4.9 Kharsali earthquake using waveform modelling

This paper presents the computation of time series of the 22 July 2007 M 4.9 Kharsali earthquake. It occurred close to the Main Central Thrust (MCT) where seismic gap exists. The main shock and 17 aftershocks were located by closely spaced eleven seismograph stations in a network that involved VSAT based real-time seismic monitoring. The largest aftershock of M 3.5 and other aftershocks occurred within a small volume of 4 × 4 km horizontal extent and between depths of 10 and 14 km. The values of seismic moment (M _∘) determined using P-wave spectra and Brune’s model based on f ² spectral shape ranges from 10¹⁸ to 10²³ dyne-cm. The initial aftershocks occurred at greater depth compared to the later aftershocks. The time series of ground motion have been computed for recording sites using geometric ray theory and Green’s function approach. The method for computing time series consists in integrating the far-field contributions of Green’s function for a number of distributed point source. The generated waveforms have been compared with the observed ones. It has been inferred that the Kharsali earthquake occurred due to a northerly dipping low angle thrust fault at a depth of 14 km taking strike N279°E, dip 14° and rake 117°. There are two regions on the fault surface which have larger slip amplitudes (asperities) and the rupture which has been considered as circular in nature initiated from the asperity at a greater depth shifting gradually upwards. The two asperities cover only 10% of the total area of the causative fault plane. However, detailed seismic imaging of these two asperities can be corroborated with structural heterogeneities associated with causative fault to understand how seismogenesis is influenced by strong or weak structural barriers in the region.