Laplace-transform analytic element solution of transient flow in porous media

A Laplace-transform analytic element method (LT-AEM) is described for the solution of transient flow problems in porous media. Following Laplace transformation of the original flow problem, the analytic element method (AEM) is used to solve the resultant time-independent modified Helmholtz equation, and the solution is inverted numerically back into the time domain. The solution is entirely general, retaining the mathematical elegance and computational efficiency of the AEM while being amenable to parallel computation. It is especially well suited for problems in which a solution is required at a limited number of points in space–time, and for problems involving materials with sharply contrasting hydraulic properties. We illustrate the LT-AEM on transient flow through a uniform confined aquifer with a circular inclusion of contrasting hydraulic conductivity and specific storage. Our results compare well with published analytical solutions in the special case of radial flow.     

Feldspathic lunar meteorites and their implications for compositional remote sensing of the lunar surface and the composition of the lunar crust

We present new compositional data for six feldspathic lunar meteorites, two from cold deserts (Yamato 791197 and 82192) and four from hot deserts (Dhofar 025, Northwest Africa 482, and Dar al Gani 262 and 400). The concentrations of FeO (or Al₂O₃) and Th (or any other incompatible element) together provide first-order compositional information about lunar polymict samples (breccias and regoliths) and regions of the lunar surface observed from orbit. Concentrations of both elements on the lunar surface have been determined from data acquired by orbiting spacecraft, although the derived concentrations have large uncertainties and some systematic errors compared to sample data. Within the uncertainties and errors in the concentrations derived from orbital data, the distribution of FeO and Th concentrations among lunar meteorites, which represent ∼18 source regions on the lunar surface, is consistent with that of 18 random samples from the surface. Approximately 11 of the lunar meteorites are low-FeO and low-Th breccias, consistent with large regions of the lunar surface, particularly the northern farside highlands. Almost all regoliths from Apollo sites, on the other hand, have larger concentrations of both elements because they contain Fe-rich volcanic lithologies from the nearside maria and Th-rich lithologies from the high-Th anomaly in the northwestern nearside. The feldspathic lunar meteorites thus offer our best estimate of the composition of the surface of the feldspathic highlands, and we provide such an estimate based on the eight most well-characterized feldspathic lunar meteorites. The variable but high (on average) Mg/Fe ratio of the feldspathic lunar meteorites compared to ferroan anorthosites confirms a hypothesis that much of the plagioclase at the surface of the feldspathic highlands is associated with high-Mg/Fe feldspathic rocks such as magnesian granulitic breccia, not ferroan anorthosite. Geochemically, the high-Mg/Fe breccias appear to be unrelated to the mafic magnesian-suite rocks of the Apollo collection. Models for the formation of the upper lunar crust as a simple flotation cumulate composed mainly of ferroan anorthosite do not account for the complexity of the crust as inferred from the feldspathic lunar meteorites.     

Increased sediment accretion rates following invasion byPhragmites australis: The role of litter

Negative connotations of invasive plants worldwide have implicated them as the bearers of unfavorable ecosystem change. We contrasted 5-yr-old and 20-yr-oldPhragmites populations with pre-invasion areas occupied byTypha spp. andPanicum virgatum in an oligohaline tidal marsh of Chesapeake Bay. Peak live biomass was 3 times greater, while standing dead and litter was twice as great in the 20-yr-oldPhragmites. It is this abundance of concentrated litter on the marsh surface of maturePhragmites populations that we implicate as encouraging the trapping of organic and mineral matter. The rate of vertical accretion in 20-yr-oldPhragmites populations is 3–4 mm yr⁻¹ above the adjacent populations. By integrating the constant initial concentration and constant rate of supply models on individual²¹⁰Pb cores, we estimate thatPhragmites populations require a minimum of 7-yr post-colonization to enhance rates of accretion in this system. In ligh of the considerable loss of marsh habitat from relative sea-level rise, this finding contests the view that invasion creates strictly undesirable change at the ecosystem level.     

A method for precise measurement of argon 40/36 and krypton/argon ratios in trapped air in polar ice with applications to past firn thickness and abrupt climate change in Greenland and at Siple Dome, Antarctica

We describe a method for measuring the ⁴⁰Ar/³⁶Ar ratio and the ⁸⁴Kr/³⁶Ar ratio in air from bubbles trapped in ice cores. These ratios can provide constraints on the past thickness of the firn layer at the ice core site and on the magnitude of past rapid temperature variations when combined with measured ¹⁵N/¹⁴N. Both variables contribute to paleoclimatic studies and ultimately to the understanding of the controls on Earth’s climate. The overall precision of the ⁴⁰Ar/³⁶Ar method (1 standard error of the mean) is 0.012‰ for a sample analyzed in duplicate, corresponding to ±0.6 m in reconstructed firn thickness. We use conventional dynamic isotope ratio mass spectrometry with minor modifications and special gas handling techniques designed to avoid fractionation. About 100 g of ice is used for a duplicate pair of analyses. An example of the technique applied to the GISP2 ice core yields an estimate of 11 ± 3K of abrupt warming at the end of the last glacial period 15,000 years ago. The krypton/argon ratio can provide a diagnostic of argon leakage out of the bubbles, which may happen (naturally) during bubble close-off or (artifactually) if samples are warmed near the freezing point during core retrieval or storage. Argon leakage may fractionate the remaining ⁴⁰Ar/³⁶Ar ratio by +0.007‰ per ‰ change in ⁸⁴Kr/³⁶Ar, introducing a possible bias in reconstructed firn thickness of about +2 m if thermal diffusion is not accounted for or +6 m if thermal diffusion effects are quantified with measured ¹⁵N/¹⁴N. Reproducibility of ⁸⁴Kr/³⁶Ar measured in air is about ±0.2‰ (1 standard error of the mean) but is about ±1‰ for ice core samples. Ice core samples are systematically enriched in ⁸⁴Kr/³⁶Ar relative to atmosphere by ∼5‰, probably reflecting preferential size-dependent exclusion of the smaller argon atom during bubble entrapment. Recent results from the Siple Dome ice core reveal two climate events during the last deglaciation, including an 18-m reduction in firn thickness associated with an abrupt warming at sometime between 18 and 22 kyr BP and a partial or total removal of the firn during an ablation event at 15.3 kyr BP.     

Serpentinization of abyssal peridotites from the MARK area, Mid-Atlantic Ridge: sulfur geochemistry and reaction modeling

The opaque mineralogy and the contents and isotope compositions of sulfur in serpentinized peridotites from the MARK (Mid-Atlantic Ridge, Kane Fracture Zone) area were examined to understand the conditions of serpentinization and evaluate this process as a sink for seawater sulfur. The serpentinites contain a sulfur-rich secondary mineral assemblage and have high sulfur contents (up to 1 wt.%) and elevated δ³⁴S_sulfide (3.7 to 12.7‰). Geochemical reaction modeling indicates that seawater-peridotite interaction at 300 to 400°C alone cannot account for both the high sulfur contents and high δ³⁴S_sulfide. These require a multistage reaction with leaching of sulfide from subjacent gabbro during higher temperature (∼400°C) reactions with seawater and subsequent deposition of sulfide during serpentinization of peridotite at ∼300°C. Serpentinization produces highly reducing conditions and significant amounts of H₂ and results in the partial reduction of seawater carbonate to methane. The latter is documented by formation of carbonate veins enriched in ¹³C (up to 4.5‰) at temperatures above 250°C. Although different processes produce variable sulfur isotope effects in other oceanic serpentinites, sulfur is consistently added to abyssal peridotites during serpentinization. Data for serpentinites drilled and dredged from oceanic crust and from ophiolites indicate that oceanic peridotites are a sink for up to 0.4 to 6.0 × 10¹² g seawater S yr⁻¹. This is comparable to sulfur exchange that occurs in hydrothermal systems in mafic oceanic crust at midocean ridges and on ridge flanks and amounts to 2 to 30% of the riverine sulfate source and sedimentary sulfide sink in the oceans. The high concentrations and modified isotope compositions of sulfur in serpentinites could be important for mantle metasomatism during subduction of crust generated at slow spreading rates.     

Estimation of recent glacial variations in baspa basin using remote sensing technique

The Himalayas has one of the largest concentrations of glaciers outside the Polar Regions. Various reports suggest that significant number of mountain glaciers is shrinking due to climatic variations. Monitoring of these glaciers is important to assess future availability of water resources in the Himalayan region. However, Himalayan glaciers are normally difficult to monitor due to the rugged, mountainous terrain. Therefore, images of Indian Remote Sensing Satellite were used to monitor glaciers in the Baspa basin. Investigations have shown the presence of 30 glaciers in the basin, with areal extent of 167 km2. Out of these, 19 glaciers, with areal extent of 140 km² were selected to estimate retreat. Investigation suggests that almost all glaciers are retreating in the study basin and overall 19% deglaciation has been observed from 1962 to 2001. In general, altitude distribution appears to have significant influence on glacial retreat. Glaciers located around 5000 m altitude range are showing 24% loss as compared to 14% by glaciers located in altitude range higher than 5400 m. In addition, mean altitude of glacier terminus is shifted upward by 88 m, i.e. from 4482 to 4570 m in last 39 years. The glacial volumes were estimated using regression relationship between area and depth. The investigations have suggested that 19.10 km³ of glacial water stored in the 19 glaciers in 1962, has been reduced to 14.71 km³ in 2001, respectively, an overall loss of 23 percent in a period between 1962 and 2001. These investigations suggest that all glaciers in the Baspa Basin are reducing and in long term, such reducing trend can create scarcity of water in the region.     

P-SH conversions in a flat-layered medium with anisotropy of arbitrary orientation

P-SH conversion is commonly observed in teleseismic P waves, and is often attributed to dipping interfaces beneath the receiver. Our modelling suggests an alternative explanation in terms of flat-layered anisotropy. We use reflectivity techniques to compute three-component synthetic seismograms in a 1-D anisotropic layered medium. For each layer of the medium, we prescribe values of seismic velocities and hexagonally symmetric anisotropy about a common symmetry axis of arbitrary orientation. A compressional wave in an anisotropic velocity structure suffers conversion to both SV -and SH -polarized shear waves, unless the axis of symmetry is everywhere vertical or the wave travels parallel to all symmetry axes. The P-SV conversion forms the basis of the widely used 'receiver function' technique. The P-SH conversion occurs at interfaces where one or both layers are anisotropic. A tilted axis of symmetry and a dipping interface in isotropic media produce similar amplitudes of both direct ( P ) and converted ( Ps ) phases, leaving the backazimuth variation of the P-Ps delay as the main discriminant. Seismic anisotropy with a tilted symmetry axis leads to complex synthetic seismograms in velocity models composed of just a few flat homogeneous layers. It is possible therefore to model observations of P coda with prominent transverse components with relatively simple 1-D velocity structures. Successful retrieval of salient model characteristics appears possible using multiple realizations of a genetic-algorithm (GA) inversion of P coda from several backazimuths. Using GA inversion, we determine that six P coda recorded at station ARU in central Russia are consistent with models that possess strong (> 10 per cent) anisotropy in the top 5 km and between 30 and 43 km depth. The symmetry axes are tilted, and appear aligned with the seismic anisotropy orientation in the mantle under ARU suggested by SKS splitting.     

Mitigating the Effects of Climate Change on the Water Resources of the Columbia River Basin

The potential effects of climate change on the hydrology and water resources of the Columbia River Basin (CRB) were evaluated using simulations from the U.S. Department of Energy and National Center for Atmospheric Research Parallel Climate Model (DOE/NCAR PCM). This study focuses on three climate projections for the 21st century based on a `business as usual' (BAU) global emissions scenario, evaluated with respect to a control climate scenario based on static 1995 emissions. Time-varying monthly PCM temperature and precipitation changes were statistically downscaled and temporally disaggregated to produce daily forcings that drove a macro-scale hydrologic simulation model of the Columbia River basin at 1/4-degree spatial resolution. For comparison with the direct statistical downscaling approach, a dynamical downscaling approach using a regional climate model (RCM) was also used to derive hydrologic model forcings for 20-year subsets from the PCM control climate (1995–2015) scenario and from the three BAU climate(2040–2060) projections. The statistically downscaled PCM scenario results were assessed for three analysis periods (denoted Periods 1–3: 2010–2039,2040–2069, 2070–2098) in which changes in annual average temperature were +0.5,+1.3 and +2.1 °C, respectively, while critical winter season precipitation changes were –3, +5 and +1 percent. For RCM, the predicted temperature change for the 2040–2060 period was +1.2 °C and the average winter precipitation change was –3 percent, relative to the RCM controlclimate. Due to the modest changes in winter precipitation, temperature changes dominated the simulated hydrologic effects by reducing winter snow accumulation, thus shifting summer streamflow to the winter. The hydrologic changes caused increased competition for reservoir storage between firm hydropower and instream flow targets developed pursuant to the Endangered Species Act listing of Columbia River salmonids. We examined several alternative reservoir operating policies designed to mitigate reservoir system performance losses. In general, the combination of earlier reservoir refill with greater storage allocations for instream flow targets mitigated some of the negative impacts to flow, but only with significant losses in firm hydropower production (ranging from –9 percent in Period1 to –35 percent for RCM). Simulated hydropower revenue changes were lessthan 5 percent for all scenarios, however, primarily due to small changes inannual runoff.     

Dislocating modernity: Identity, space and representations of street trade in Durban, South Africa

South Africa’s cities have experienced dramatic changes over the past decade. Cities are now home to a multiracial population, and have been transformed by new forms of economic and social interaction. For some, these changes have become a significant source of fear and anxiety. In this paper, we examine reactions to urban spatial change in the city of Durban, as expressed in local newspapers and interviews with suburban residents. We describe how the discourses of urban change in Durban have centred on the increased presence of street traders within the city’s public spaces, and the various ways in which the activity of street trade has disrupted long-established modernist norms governing the occupation and use of the urban space. Specifically, we offer a detailed reading of three prominent narratives within the discussion of street traders in Durban--chaos, congestion and pollution. We argue that street traders have come to embody a wide range of more deeply seated cultural anxieties, which have been brought to the fore in the context of South Africa’s transition. These anxieties arise from the ways in which modern understandings of order, agency and subjectivity have been called into question by material changes in the city, and have implications for the nature of citizenship and civic engagement in post-apartheid South Africa.