Analytic Element Modeling of Steady Interface Flow in Multilayer Aquifers Using AnAqSim

This paper presents the analytic element modeling approach implemented in the software AnAqSim for simulating steady groundwater flow with a sharp fresh‐salt interface in multilayer (three‐dimensional) aquifer systems. Compared with numerical methods for variable‐density interface modeling, this approach allows quick model construction and can yield useful guidance about the three‐dimensional configuration of an interface even at a large scale. The approach employs subdomains and multiple layers as outlined by Fitts (2010) with the addition of discharge potentials for shallow interface flow (Strack 1989). The following simplifying assumptions are made: steady flow, a sharp interface between fresh‐ and salt water, static salt water, and no resistance to vertical flow and hydrostatic heads within each fresh water layer. A key component of this approach is a transition to a thin fixed minimum fresh water thickness mode when the fresh water thickness approaches zero. This allows the solution to converge and determine the steady interface position without a long transient simulation. The approach is checked against the widely used numerical codes SEAWAT and SWI/MODFLOW and a hypothetical application of the method to a coastal wellfield is presented.     

Mineralogy and thermal properties of V-type Asteroid 956 Elisa: Evidence for diogenitic material from the Spitzer IRS (5-35 μm) spectrum

We present the thermal infrared (5-35 μm) spectrum of 956 Elisa as measured by the Spitzer Infrared Spectrograph (“IRS”; Houck, J.R. et al. [2004]. Astrophys. J. Suppl. 154, 18-24) together with new groundbased lightcurve data and near-IR spectra. From the visible lightcurve photometry, we determine a rotation period of 16.494 ± 0.001 h, identify the rotational phase of the Spitzer observations, and estimate the visible absolute magnitude (H_V) at that rotational phase to be 12.58 ± 0.04. From radiometric analysis of the thermal flux spectrum, we find that at the time of observation 956 Elisa had a projected radius of 5.3 ± 0.4 km with a visible albedo p_V = 0.142 ± 0.022, significantly lower than that of the prototype V-type asteroid, 4 Vesta. (This corresponds to a radius of 5.2 ± 0.4 km at lightcurve mean.) Analysis with the standard thermal model (STM) results in a sub-solar temperature of 292.3 ±  2.8 K and beaming parameter η = 1.16 ± 0.05. Thermophysical modeling places a lower limit of on the thermal inertia of the asteroid’s surface layer (if the surface is very smooth) but more likely values fall between 30 and depending on the sense of rotation.The emissivity spectrum, calculated by dividing the measured thermal flux spectrum by the modeled thermal continuum, exhibits mineralogically interpretable spectral features within the 9-12 μm reststrahlen band, the 15-16.5 μm Si-O-Si stretching region, and the 16-25 μm reststrahlen region that are consistent with pyroxene of diogenitic composition: extant diogenitic pyroxenes fall within the narrow compositional range Wo_2±1En_74±2Fs_24±1. Spectral deconvolution of the 9-12 μm reststrahlen features indicates that up to ≈20% olivine may also be present, suggesting an olivine-diogenite-like mineralogy. The mid-IR spectrum is inconsistent with non-cumulate eucrite as the major component on the surface of 956 Elisa, although cumulate eucrite material may be present at abundances lower than that of the diogenite component.Analysis of new near-IR spectra of 956 Elisa with the Modified Gaussian Model (MGM; Sunshine, J.M., Pieters, C.M., Pratt, S.F. [1990]. J. Geophys. Res. 95 (May), 6955-6966) results in two pyroxene compositions: 75% magnesian low-Ca pyroxene and 25% high-Ca pyroxene. High-Ca pyroxene is not evident in the mid-IR data, but may belong to a component that is underrepresented in the mid-IR spectrum either because of its spatial distribution on the asteroid or because of its particle size. High-Ca pyroxenes that occur as exsolution lamellae may also be more evident spectrally in the NIR than in the mid-IR. In any case, we find that the mid-IR spectrum of 956 Elisa is dominated by emission from material of diogenite-like composition, which has very rarely been observed among asteroids.     

Regolith depth growth on an icy body orbiting Saturn and evolution of bidirectional reflectance due to surface composition changes

Using a Markov chain model, we consider the regolith growth on a small body in orbit around Saturn, subject to meteoritic bombardment, and assuming all impact ejecta are re-collected. We calculate the growth of regolith and the fractional pollution, assuming an initial pure ice body and amorphous carbon as a pollutant. We extend the meteorite flux of Cuzzi and Estrada (Cuzzi, J., Estrada, P. [1998]. Icarus 132, 1-35) to larger sizes to consider the effect of disruption of the moonlet on other moonlets in the ensemble. This is a relatively small effect, completely negligible for moonlets of 1 m radius. For the given impact model, fractional pollution reaches 22% for 1 m bodies, but only 3% for 10 m bodies, 1.7% for 20 m bodies, and 1% for 30 m bodies after 4 byr. By considering an ensemble of moonlets, which have identical cross-sections for releasing and capturing ejecta, this analysis can be extended to a model of particles in Saturn’s rings, where the calculated spectra can be compared to observed ring spectra. The measured spectral reflectance of Saturn’s rings from Cassini observations therefore constrains the size and age of the ring particles. The comparison between 1 m, 10 m, 20 m, and 30 m particles confirms that for larger ring mass, the current rings would be less polluted; for the largest particles, we expect negligible changes in the UV spectrum after 4 byr of meteoritic bombardment. We consider two end members for mixing of the meteoritic material: areal and intimate. Given the uncertainties in the actual mixing of the meteoritic infall and in its composition (as a worst case, we assume the meteoritic material is 100% amorphous carbon, intimately mixed) initially pure ice 30 m ring particles would darken after 4 byr of exposure by 15%.     

Lethal concentrations of suspended solids for common native fish species that are rare in New Zealand rivers with high suspended solids loads

Abstract

Common native fish species that are rare in New Zealand rivers with high suspended solids (SS) concentrations were identified to determine which species were likely to be most affected by high SS concentrations. Mortality rates for the three most sensitive species (banded kokopu Galaxias fasciatus, redfin bully Gobiomorphus huttoni, common smelt Retropinna retropinnd) were measured over a range of SS concentrations to identify lethal concentrations. The survival of banded kokopu and redfin bully was not reduced by SS concentrations up to 43 000 g m^?3; however, the 24‐h survival of common smelt was reduced by SS concentrations over 1000 g m^?3, and its median lethal concentration (LC50) was close to 3000 g m^?3. SS concentrations over 3000 g m^?3 for more than 24 h will therefore reduce smelt but not affect the mortality of most other common native fish species. Although mortality from high concentrations of SS can account for the scarcity of smelt in rivers with high SS concentrations, mortality from high SS does not explain the scarcity of banded kokopu or redfin bully in such rivers. High turbidity (>20 nephelometric turbidity units) is known to limit the upstream migration and recruitment of juvenile banded kokopu in rivers with increased SS concentrations. Banded kokopu is therefore affected by increased turbidity from SS rather than lethal SS concentrations. As redfin bully were not affected by either high turbidity levels or high SS concentrations, the low abundance of this species in rivers with high SS concentrations is likely to be related to the siltation of benthic habitats.     

Setting new thermal standards for transitional and coastal (TraC) waters

The construction of a new generation of coastal power stations in the UK and other western European countries has highlighted the absence of robust standards for thermal discharges to transitional and coastal waters. The proposed discharge volumes are greater than hitherto, yet there has been little independent critical examination of their potential impact, whilst much of the existing guidance has been adapted from freshwater practice. This review considers the available knowledge on the tolerance and behaviour of fish and other marine biota to heated effluents. Appropriate ways are suggested of grouping fish species to reflect their sensitivity to thermal effects. The plethora of existing standards are considered and their validity assessed in a framework of predicted seawater temperature rise. Those species or groups of organisms most likely to be affected are identified and finally specific recommendations for thermal standards consistent with long term sustainability are proposed.     

Influence of anthropogenic land‐use change on hillslope erosion in the Waipaoa River Basin,New Zealand

European settlement of the Poverty Bay Region resulted in deforestation and conversion of > 90% of the landscape to pastureland. The resulting loss of vegetation triggered a rapid increase in hillslope erosion as widespread landslide complexes and gully systems developed on weak lithologic units in the Waipaoa Basin. To quantify the rate and volume of historic hillslope degradation, we used a 1956–2010 sequence of aerial photographs for a ~16 km² catchment to map temporal changes in the spatial extent of active landslides. Then we created a ‘turf index’ based on the extent and style of pastoral ground disruption, which correlates with downslope velocity. Based on the movement of trees and other features, we assigned average velocities to the turf classes as follows: (1) minimal disrupted ground: 0.6 m/yr, (2) a mix of disrupted ground and intact blocks: 3.4 m/yr, and (3) no intact blocks or vegetation: > 6 m/yr. We then calculated the average annual sediment flux using these turf‐derived velocities, the width of the landslide‐channel intersection, and an average toe depth of 4.4 ± 1.3 m (mean ± standard deviation [SD]) from 37 field measurements. The resulting catchment averaged erosion rates are (mean ± SD): 29.9 ± 12.9 mm/yr (1956), 28.8 ± 13.7 mm/yr (1969), 13.4 ± 4.9 mm/yr (1979), 17.0 ± 6.2 mm/yr (1988), and 9.9 ± 3.6 mm/yr (2010). Compared with long‐term (post‐18 ka) erosion rates (1.6 mm/yr) and the long‐term uplift rate (~1 mm/yr) for this site, the 50‐year anthropogenically‐driven rate is an order of magnitude larger (~20 mm/yr). Previously, we measured an increase in erosion over the past 3.4 kyr (2.2 mm/yr), and here, we demonstrate this increase could be primarily due to human land‐use change – showing that a century of rapid erosion superimposed on the background geologic rate can profoundly skew the interpretation of erosion rates. Copyright © 2016 John Wiley & Sons, Ltd.     

Tsunami damage to coastal defences and buildings in the March 11th 2011 M w 9.0 Great East Japan earthquake and tsunami

On March 11th 2011 a M _w 9.0 mega-thrust interface subduction earthquake, the Great East Japan Earthquake, occurred 130 km off the northeast coast of Japan in the Pacific Ocean at the Japan Trench, triggering tsunami which caused damage along 600 km of coastline. Observations of damage to buildings (including vertical evacuation facilities) and coastal defences in Tōhoku are presented following investigation by the Earthquake Engineering Field Investigation Team (EEFIT) at 10 locations in Iwate and Miyagi Prefectures. Observations are presented in the context of the coastal setting and tsunami characteristics experienced at each location. Damage surveys were carried out in Kamaishi City and Kesennuma City using a damage scale for reinforced concrete (RC), timber and steel frame buildings adapted from an earlier EEFIT tsunami damage scale. Observations show that many sea walls and breakwaters were overtopped, overturned, or broken up, but provided some degree of protection. We show the extreme variability of damage in a local area due to inundation depth, flow direction, velocity variations and sheltering. Survival of many RC shear wall structures shows their high potential to withstand local earthquake and significant tsunami inundation but further research is required into mitigation of scour, liquefaction, debris impact, and the prevention of overturning failure. Damage to steel and timber buildings are also discussed. These observations are intended to contribute to mitigation of future earthquake and tsunami damage by highlighting the key features which influence damage level and local variability of damage sustained by urban coastal infrastructure when subjected to extreme tsunami inundation depths.     

Seasonal and Spatial Effects of Wastewater Effluent on Growth,Survival, and Accumulation of Microbial Contaminants by Oysters in Mobile Bay,Alabama

We measured seasonal effects of wastewater treatment plant (WTP) effluent on growth, survival, and accumulation of microbes in oysters near a major WTP in Mobile Bay, AL. Despite higher nutrients near the WTP, seasonal conditions rather than distance affected chlorophyll a concentration and oyster growth. In summer and fall, when oyster growth was higher, δ¹⁵N‰ in oysters near the WTP changed through time to reflect δ¹⁵N‰ in effluent (approx. −4‰). Microbial indicators (male-specific coliphage, fecal coliforms) were highest in oysters near the WTP in all seasons and correlated with δ¹⁵N‰ in fall and summer. Increased riverine discharge and slower acquisition of δ¹⁵N‰ likely confounded correlations in winter/spring. Although we did not detect gross ecological effects of wastewater exposure for oysters, data indicated wastewater-derived particles entered the local food web and accumulated in oysters. These data highlight the importance of using multiple indicators of wastewater exposure and considering both seasonal and spatial effects when defining wastewater influence on a system or species.     

Dental microwear texture analysis of late Pliocene Procynocephalus subhimalayanus (Primates: Cercopithecidae) of the Upper Siwaliks, India

Late Pliocene Procynocephalus subhimalayanus from the Upper Siwaliks, India is known from only three specimens. The dietary proclivities of this taxon have implications for reconstructing the paleoecology of the Upper Siwaliks. The dental microwear texture properties of Procynocephalus are compared to those from extant tropical forest primates including Alouatta palliata (n = 11), Cebus apella (n = 13), Gorilla gorilla (n = 9), Lophocebus albigena (n = 15) and Trachypithecus cristatus (n = 12). Dental microwear textures are generated by scanning the surface enamel of Facet 9 using white-light confocal microscopy at 100x. Four variables were extracted from scale-sensitive fractal analysis, and the data were ranked before ANOVA with post-hoc tests of significance and multivariate analyses were performed. Procynocephalus clusters closest to Lophocebus, Cebus and some Gorilla specimens suggesting hard-object feeding characterized a portion of its diet. The dental microwear texture of Procynocephalus supports interpretations of widespread grasslands of the Late Pliocene Kansal Formation (Pinjor zone). The extreme enamel complexity characterizing Procynocephalus may derive from consumption of underground storage organs, or other foods with high grit loads. Foods consumed near ground level carry a heavy load of abrasive minerals possibly contributing to greater enamel surface complexity and textural fill volume.