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.     

The role of boundary organizations in climate change adaptation from the perspective of municipal practitioners

In climate change impact research it is crucial to carefully select the meteorological input for impact models. We present a method for model selection that enables the user to shrink the ensemble to a few representative members, conserving the model spread and accounting for model similarity. This is done in three steps: First, using principal component analysis for a multitude of meteorological parameters, to find common patterns of climate change within the multi-model ensemble. Second, detecting model similarities with regard to these multivariate patterns using cluster analysis. And third, sampling models from each cluster, to generate a subset of representative simulations. We present an application based on the ENSEMBLES regional multi-model ensemble with the aim to provide input for a variety of climate impact studies. We find that the two most dominant patterns of climate change relate to temperature and humidity patterns. The ensemble can be reduced from 25 to 5 simulations while still maintaining its essential characteristics. Having such a representative subset of simulations reduces computational costs for climate impact modeling and enhances the quality of the ensemble at the same time, as it prevents double-counting of dependent simulations that would lead to biased statistics.     

Determination of mercury biogeochemical fluxes in the remote Mackenzie River Basin,northwest Canada,using speciation of sulfur and organic carbon

The Mackenzie River Basin (MRB) in NW Canada is one of the least human-impacted large watersheds in the world. The western and eastern sub-basins of the MRB are also marked by contrasting geology, geomorphology, hydrology, and biology. These remarkable differences in a remote river basin provide a unique opportunity to probe the biogeochemical processes governing the sources, transport, and bioavailability of Hg at the terrestrial-marine interface and ultimately in the Arctic Ocean. Based on a large dataset of the concentration and speciation of Hg, S and organic matter in surface sediments across the MRB, a source-apportioned budget was constructed for Hg in the MRB. The results indicate that the Hg flux in the basin originates primarily from the weathering of sulfide minerals in the western mountainous sub-basin (∼78% of the total), followed by the erosion of coal deposits along the mainstream of the Mackenzie River (∼10%), with the remainder split between atmospheric inputs and organic matter-bound Hg (6% and 5%, respectively). Although the Hg flux from the eastern peatland sub-basin only accounts for ∼10% of the total riverine Hg flux, Hg in this region correlates strongly with labile organic matter, and will likely have a much stronger influence on local biota.     

Melting of metasomatized subcontinental lithosphere: undersaturated mafic lavas from Rungwe,Tanzania

 This paper uses the geochemistry of primitive mafic lavas from the Rungwe volcanic province (southwestern Tanzania) to infer the source mineralogy and melting history. Post-Miocene mafic lavas from Rungwe include alkali basalts, basanites, nephelinites and picrites with up to 18.9 wt% MgO; nephelinites (>13.5% normative nepheline) are restricted to Kiejo volcano in the southern portion of the province. Rungwe lavas differ from most Western Rift volcanics in that they are not unusually potassic (K₂O/Na₂O ca. 0.40). Sparsely phyric mafic lavas contain phenocrysts and xenocrysts of plagioclase (An_82–90), clinopyroxene (4.5–9.5 wt% Al₂O₃), and olivine (Fo_79–88); one basanite contains a 1 mm xenocryst of apatite included in magnesian clinopyroxene. All samples have high abundances of incompatible elements (e.g., 0.7–2.2 wt% P₂O₅) and are enriched in REE relative to HFSE (Hf, Zr, Ti, Y), Cs, Ba, and K. Some incompatible element ratios are constant throughout the Rungwe suite (e.g., Zr/Nb, Sr/Ce, K/Rb), but other ratios are extremely variable and exceed the range measured in global Ocean Island Basalts (OIB) (e.g., Ba/Nb, Sm/Zr, La/Nb, Pb/Ce, Nb/U). The range in degree of silica saturation, and its excellent correlation with P₂O₅/Al₂O₃, indicate that the Rungwe suite records variable degrees of melting. Variations of individual incompatible trace element abundances in nephelinite and basanite samples suggest that the source contains metasomatic amphibole, ilmenite, apatite, and zircon. The Rungwe suite is interpreted as a series of low-percentage melts of CO₂-rich peridotite at pressures that span the garnet-spinel transition. A geochemical comparison of Rungwe samples to lavas from other Western Rift volcanic centers requires that the source mineralogy varies along the rift axis, although each province is underlain by metasomatized peridotite. The incompatible trace element signatures of Western Rift lavas indicate that the source area is typically homogeneous on the scale of individual volcanoes, although lavas from each volcano reflect a range in degree of melting. Significantly, volcanoes with distinct geochemistry are always separated by major rift faults, suggesting that volcanic and tectonic surface features may correspond to metasomatic provinces within the subcontinental lithospheric mantle. Received: 30 May 1994 / Accepted: 5 April 1995     

A review of the assessments of the major fisheries resources in South Africa

The waters off South Africa's coastline boast a rich mix of commercially fished species. Quantitative assessments of these marine resources have developed from simple methods first applied in the 1970s, to models that encompass a wide range of methodologies. The more valuable resources have undergone regular assessments in recent decades, with frequencies closely related to the management approach employed for each fishery. Many of these assessments form the operating models used to simulation-test candidate management procedures. This paper provides a comprehensive review of the assessments of 11 of the most important fisheries resources in South Africa. Some assessments use simple biomass dynamics models, whereas others are a hybrid of age- and length-based models, each designed to model the specific characteristics of the resource and fishery concerned. Many of the assessments have been disaggregated by species/stock and/or area as related multispecies/stock/ distribution hypotheses have arisen. This paper explores the similarities and differences in the data available and the methods applied. The review indicates that, whereas the status of three of these resources cannot be estimated reliably at present, the status of six resources is considered to be reasonable to good, whereas that of abalone Haliotis midae and West Coast rock lobster Jasus lalandii remains poor.     

Analysis of Cassini/CIRS limb spectra of Titan acquired during the nominal mission II: Aerosol extinction profiles in the 600–1420 cm spectral range

We have analyzed the continuum emission of limb spectra acquired by the Cassini/CIRS infrared spectrometer in order to derive information on haze extinction in the 3–0.02 mbar range (∼150–350 km). We focused on the 600–1420 cm⁻¹ spectral range and studied nine different limb observations acquired during the Cassini nominal mission at 55°S, 20°S, 5°N, 30°N, 40°N, 45°N, 55°N, 70°N and 80°N. By means of an inversion algorithm solving the radiative transfer equation, we derived the vertical profiles of haze extinction coefficients from 17 spectral ranges of 20-cm⁻¹ wide at each of the nine latitudes. At a given latitude, all extinction vertical profiles retrieved from various spectral intervals between 600 and 1120 cm⁻¹ display similar vertical slopes implying similar spectral characteristics of the material at all altitudes. We calculated a mean vertical extinction profile for each latitude and derived the ratio of the haze scale height (H_haze) to the pressure scale height (H_gas) as a function of altitude. We inferred H_haze/H_gas values varying from 0.8 to 2.4. The aerosol scale height varies with altitude and also with latitude. Overall, the haze extinction does not show strong latitudinal variations but, at 1 mbar, an increase by a factor of 1.5 is observed at the north pole compared to high southern latitudes. The vertical optical depths at 0.5 and 1.7 mbar increase from 55°S to 5°N, remain constant between 5°N and 30°N and display little variation at higher latitudes, except the presence of a slight local maximum at 45°N. The spectral dependence of the haze vertical optical depth is uniform with latitude and displays three main spectral features centered at 630 cm⁻¹, 745 cm⁻¹ and 1390 cm⁻¹, the latter showing a wide tail extending down to ∼1000 cm⁻¹. From 600 to 750 cm⁻¹, the optical depth increases by a factor of 3 in contrast with the absorbance of laboratory tholins, which is generally constant. We derived the mass mixing ratio profiles of haze at the nine latitudes. Below the 0.4-mbar level all mass mixing ratio profiles increase with height. Above this pressure level, the profiles at 40°N, 45°N, 55°N, at the edge of the polar vortex, display a decrease-with-height whereas the other profiles increase. The global increase with height of the haze mass mixing ratio suggest a source at high altitudes and a sink at low altitudes. An enrichment of haze is observed at 0.1 mbar around the equator, which could be due to a more efficient photochemistry because of the strongest insolation there or an accumulation of haze due to a balance between sedimentation and upward vertical drag.     

Effects of catastrophic flooding on stream biogeochemistry in a headwater stream in Shenandoah National Park,USA

This study examined the stream chemistry changes in Staunton River (a second‐order headwater stream with an average annual discharge 704 m³ ha^?1 yr^?1, Shenandoah National Park, Virginia) resulting from a catastrophic flood in June 1995. This flood, which followed after 800 mm of rain in a 4‐day period, caused large‐scale debris flows and complete scouring of riparian soils down to bedrock in the lower 2 km of the stream, and has been estimated to be a 1000‐year flood. The flood affected stream chemistry on both short‐ and long‐term time scales. The primary short‐term response was elevations in stream concentration of Ca²⁺, Mg²⁺, and K⁺ by 59%, 87%, and 49%, respectively, for 6 months immediately following the flood. The long‐term impact of decreased concentration of all base cations and SiO₂ during summer months (8% average) lasted about 2 years. At the episodic time scale, Ca²⁺, Mg²⁺, and K⁺ flushed from soil sources during pre‐flood storms while Na⁺ and SiO₂ diluted; these trends generally reversed during post‐flood storms for 2 years. Short‐term effects are attributed to the leaching of unconsolidated soil and upturned organic matter that clogged the streambed after the flood. The long‐term and superimposed episodic impacts may have resulted from the loss of riparian soils and vegetation in the flood. Copyright © 2008 John Wiley & Sons, Ltd.     

The problem of auroral NO+ concentration and the intensity of 1·27μ bands of O2

The problem of the theoretical computation of the emission intensities and ion composition in a weak aurora which has been preceded by a stronger event is examined. For this purpose a model auroral precipitation consisting of biexponential primaries is considered. The softer of the two components is brighter, and begins to decay after remaining steady for ten to fifteen minutes. The other, harder component starts to build up at that instant. Our results suggest that at least a part of the high $\text{n(}\text{NO}{}^{\mathrm{+}}\text{)}\text{n(}\text{O}{}_2{}^{\mathrm{+}}\text{)}\text{or}\text{I}\text{(1·27 μ)}\text{I}\text{(3914}\text{A}\text{?}\text{)}$ ratios could be attributed to the retention, by the atmosphere, of the memory of previous auroral precipitations. Thus, the serious energy paradox in the context of 1·27 μ intensity need not arise, and, in the context of the large NO⁺ density, it may perhaps be unnecessary to invoke any major conversion of O₂ to NO thus avoiding the associated energy problem.