An inter-comparison of tidal solutions computed with a range of unstructured grid models of the Irish and Celtic Sea Regions

Three finite element codes, namely TELEMAC, ADCIRC and QUODDY, are used to compute the spatial distributions of the M₂, M₄ and M₆ components of the tide in the sea region off the west coast of Britain. This region is chosen because there is an accurate topographic dataset in the area and detailed open boundary M₂ tidal forcing for driving the model. In addition, accurate solutions (based upon comparisons with extensive observations) using uniform grid finite difference models forced with these open boundary data exist for comparison purposes. By using boundary forcing, bottom topography and bottom drag coefficients identical to those used in an earlier finite difference model, there is no danger of comparing finite element solutions for “untuned unoptimised solutions” with those from a “tuned optimised solution”. In addition, by placing the open boundary in all finite element calculations at the same location as that used in a previous finite difference model and using the same M₂ tidal boundary forcing and water depths, a like with like comparison of solutions derived with the various finite element models was possible. In addition, this open boundary was well removed from the shallow water region, namely the eastern Irish Sea where the higher harmonics were generated. Since these are not included in the open boundary, forcing their generation was determined by physical processes within the models. Consequently, an inter-comparison of these higher harmonics generated by the various finite element codes gives some indication of the degree of variability in the solution particularly in coastal regions from one finite element model to another. Initial calculations using high-resolution near-shore topography in the eastern Irish Sea and including “wetting and drying” showed that M₂ tidal amplitudes and phases in the region computed with TELEMAC were in good agreement with observations. The ADCIRC code gave amplitudes about 30 cm lower and phases about 8° higher. For the M₄ tide, in the eastern Irish Sea amplitudes computed with TELEMAC were about 4 cm higher than ADCIRC on average, with phase differences of order 5°. For the M₆ component, amplitudes and phases showed significant small-scale variability in the eastern Irish Sea, and no clear bias between the models could be found. Although setting a minimum water depth of 5 m in the near-shore region, hence removing wetting and drying, reduced the small-scale variability in the models, the differences in M₂ and M₄ tide between models remained. For M₆, a significant reduction in variability occurred in the eastern Irish Sea when a minimum 5-m water depth was specified. In this case, TELEMAC gave amplitudes that were 1 cm higher and phases 30° lower than ADCIRC on average. For QUODDY in the eastern Irish Sea, average M₂ tidal amplitudes were about 10 cm higher and phase 8° higher than those computed with TELEMAC. For M₄, amplitudes were approximately 2 cm higher with phases of order 15° higher in the northern part of the region and 15° lower in the southern part. For M₆ in the north of the region, amplitudes were 2 cm higher and about 2 cm lower in the south. Very rapid M₆ tidal-phase changes occurred in the near-shore regions. The lessons learned from this model inter-comparison study are summarised in the final section of the paper. In addition, the problems of performing a detailed model–model inter-comparison are discussed, as are the enormous difficulties of conducting a true model skill assessment that would require detailed measurements of tidal boundary forcing, near-shore topography and precise knowledge of bed types and bed forms. Such data are at present not available.     

Data reporting norms for 40Ar/39Ar geochronology

Data reported in ⁴⁰Ar/³⁹Ar geochronology studies are commonly insufficient to allow computation of ages. This deficiency renders it difficult to compare ages based on different standards or constants, and often hinders critical evaluation of the results. Herein are presented an enumeration of the data that should be reported in all ⁴⁰Ar/³⁹Ar studies, including a discussion in support of these requirements. The minimum required data are identified and distinguished from parameters that are useful but may be derived from them by calculation. Finally, recommendations are made for metadata needed to document age calculations (e.g., from age spectrum or isochron analyses).     

Large-scale structure, the cosmic microwave background and primordial non-Gaussianity

Cosmic microwave background and large-scale structure data will shortly improve dramatically with the Microwave Anisotropy Probe and Planck Surveyor , and the Anglo-Australian 2-Degree Field and Sloan Digital Sky Survey. It is therefore timely to ask which of the microwave background and large-scale structure will provide a better probe of primordial non-Gaussianity. In this paper we consider this question, using the bispectrum as a discriminating statistic. We consider several non-Gaussian models and find that in each case the microwave background will provide a better probe of primordial non-Gaussianity. Our results suggest that if microwave background maps appear Gaussian, then apparent deviations from Gaussian initial conditions in galaxy surveys can be attributed with confidence to the effects of biasing. We demonstrate this precisely for the spatial bispectrum induced by local non-linear biasing.     

Significance of geomorphological and subsurface drainage controls on failures of peat‐covered hillslopes triggered by extreme rainfall

On 19 September 2003, 40 landslides of 140–18 000 m³ volume occurred within 2·5 km² on the slopes of Dooncarton Mountain (Republic of Ireland) during a storm that may have exceeded 90 mm within 90 minutes. The landslides were investigated to determine the reasons for such a high density of slope failures. All of the landslides were surveyed within four months, and nine of them were investigated in detail. The six largest landslides, all peat failures, accounted for 57% of the more than 100 000 m³ of material displaced during the event. A consistent sequence of superficial materials was found on the failed hillslopes, including an extensive iron pan at the base of a buried soil horizon 0·3 m below the base of the peat. Morphologically, almost all of the landslides occurred on steep planar slopes or around sharp convexities, with the latter failures developing retrogressively upslope. The only significant relationship found from analysis of 371 subsurface pipes and 142 seepage cracks (defined here as contiguous fissures conducting concentrated subsurface flow) across all the failures was that the thinner the peat cover, the deeper the pipes and seepage cracks occurred below the base of peat. It is concluded that most of the landslides were probably caused by a combination of excess water pressures in the buried soil horizon and the thinner overburden of peat or peaty soil associated with the steeper slope segments. Pipes and seepage cracks formed on the iron pan probably existed prior to the failure event and may have contributed to the high water pressures as rainwater inputs exceeded their discharge capacities. One large peat slide was probably triggered by excess water pressures developed within and between artificial tine cuts. The properties of the blanket peat were generally of little consequence in the occurrence of the landslides, but relict desiccation cracks and other structural weaknesses through the peat mass were probably highly significant. Although several aspects of the peat failures correspond to previously published examples, the context of these failures in terms of the topography and upland catena is distinctive. Copyright © 2007 John Wiley & Sons, Ltd.     

Impact of irrigated agriculture on groundwater-recharge salinity: a major sustainability concern in semi-arid regions

Intensive irrigated agriculture substantially modifies the hydrological cycle and often has major environmental impacts. The article focuses upon a specific concern—the tendency for progressive long-term increases in the salinity of groundwater recharge derived from irrigated permeable soils and replenishment of unconfined aquifers in more arid regions. This process has received only scant attention in the water-resource literature and has not been considered by agricultural science. This work makes an original contribution by analysing, from scientific principles, how the salinisation of groundwater recharge arises and identifies the factors affecting its severity. If not proactively managed, the process eventually will impact irrigation waterwell salinity, the productivity of agriculture itself, and can even lead to land abandonment. The types of management measure required for mitigation are discussed through three detailed case histories of areas with high-value groundwater-irrigated agriculture (in Spain, Argentina and Pakistan), which provide a long-term perspective on the evolution of the problem over various decades.     

Effects of variable winds on biological productivity on continental shelves in coastal upwelling systems

The production and distribution of biological material in wind-driven coastal upwelling systems are of global importance, yet they remain poorly understood. Production is frequently presumed to be proportional to upwelling rate, yet high winds can lead to advective losses from continental shelves, where many species at higher trophic levels reside. An idealized mixed-layer conveyor (MLC) model of biological production from constant upwelling winds demonstrated previously that the amount of new production available to shelf species increased with upwelling at low winds, but declined at high winds [Botsford, L.W., Lawrence, C.A., Dever, E.P., Hastings, A., Largier, J., 2003. Wind strength and biological productivity in upwelling systems: an idealized study. Fisheries Oceanography 12, 245–259]. Here we analyze the response of this model to time-varying winds for parameter values and observed winds from the Wind Events and Shelf Transport (WEST) study region. We compare this response to the conventional view that the results of upwelling are proportional to upwelled volume. Most new production per volume upwelled available to shelf species occurs following rapid increases in shelf transit time due to decreases in wind (i.e. relaxations). However, on synoptic, event time-scales shelf production is positively correlated with upwelling rate. This is primarily due to the effect of synchronous periods of low values in these time series, paradoxically due to wind relaxations. On inter-annual time-scales, computing model production from wind forcing from 20 previous years shows that these synchronous periods of low values have little effect on correlations between upwelling and production. Comparison of model production from 20 years of wind data over a range of shelf widths shows that upwelling rate will predict biological production well only in locations where cross-shelf transit times are greater than the time required for phytoplankton or zooplankton production. For stronger mean winds (narrower shelves), annual production falls below the peak of constant wind prediction [Botsford et al., 2003. Wind strength and biological productivity in upwelling systems: an idealized study. Fisheries Oceanography 12, 245–259], then as winds increase further (shelves become narrower) production does not decline as steeply as the constant wind prediction.     

Oxidation of pyrite during extraction of carbonate associated sulfate

The sulfur isotopic composition of carbonate associated sulfate (CAS) has been used to investigate the geochemistry of ancient seawater sulfate. However, few studies have quantified the reliability of δ³⁴S of CAS as a seawater sulfate proxy, especially with respect to later diagenetic overprinting. Pyrite, which typically has depleted δ³⁴S values due to authigenic fractionation associated with bacterial sulfate reduction, is a common constituent of marine sedimentary rocks. The oxidation of pyrite, whether during diagenesis or sample preparation, could thus adversely influence the sulfur isotopic composition of CAS. Here, we report the results of CAS extractions using HCl and acetic acid with samples spiked with varying amounts of pyrite. The results show a very strong linear relationship between the abundance of fine-grained pyrite added to the sample and the resultant abundance and δ³⁴S value of CAS. This data represents the first unequivocal evidence that pyrite is oxidized during the CAS extraction process. Our mixing models indicate that in samples with much less than 1 wt.% pyrite and relatively high δ³⁴S_pyrite values, the isotopic offset imparted by oxidation of pyrite should be much less than ? 4‰. A wealth of literature exists on the oxidation of pyrite by Fe³⁺ and we believe this mechanism drives the oxidation of pyrite during CAS extraction, during which the oxygen used to form sulfate is taken from H₂O, not O₂. Consequently, extracting CAS under anaerobic conditions would only slow, but not halt, the oxidation of pyrite. Future studies of CAS should attempt to quantify pyrite abundance and isotopic composition.     

西藏高原南部雅鲁藏布江缝合带地区地壳电性结构研究

为了探测西藏高原南部雅鲁藏布江缝合带地区地壳浅部和深部构造沿东西和南北方向的变化特征,在雅鲁藏布江缝合带地区布设了三条南北向剖面(错那—墨竹工卡、亚东—雪古拉、吉隆—措勤),采用超宽频带大地电磁测深方法进行了地壳、上地幔电性结构探测研究,发现该区主要电性结构特征为：1. 雅鲁藏布江缝合带附近表层发育大规模的高阻体,岩体延伸最深达30km以上,是冈底斯花岗岩体的反映. 2. 雅鲁藏布江缝合带的南部有小规模的良导体,在其下方和北侧发育有大规模良导体. 3. 沿剖面从南往北壳内普遍发育良导体,各良导体主体间是不连续的,规模逐渐增大,总体北倾,在缝合带附近产状较陡. 4. 在雅鲁藏布江缝合带附近良导体由西往东规模逐渐增大,导电性逐渐变好,相对雅鲁藏布江在剖面上的位置逐渐南移. 这些重要的电性特征可能是印度板块向北俯冲所形成的,深部大规模的良导体特征沿东西向的差异可能是板块碰撞引起物质沿东西向运移作用的结果.     

FRAGMENTAL BRECCIAS AND THE COLLISIONAL EVOLUTION OF ORDINARY CHONDRITE PARENT BODIES

Our survey of type 4–6 ordinary chondrites indicates that gas-poor, melt-rock and/or exotic clast-bearing fragmental breccias constitute 5%, 22% and 23%, respectively, of H, L and LL chondrites. These abundances contrast with the percentages of solar-gas-rich regolith breccias among ordinary chondrites: H (14%), L (3%) and LL (8%) (Crabb and Schultz, 1981). Petrologic study of several melt-rock-clast-bearing fragmental breccias indicates that some acquired their clasts prior to breccia metamorphism and others acquired them after metamorphism of host material. In general, the melt-rock clasts in gas-poor H chondrite fragmental breccias were acquired after breccia metamorphism and were probably formed by impacts into boulders or exposed outcrops of H4-6 material in the H chondrite parent body regolith. In contrast, most of the melt-rock clasts in gas-poor L and LL fragmental breccias were acquired prior to breccia metamorphism. The low abundance of regolith breccias among L chondrites and evidence that at least two-thirds of the L chondrites suffered a major shock event 0.5 Gyr ago, suggest that the L parent body may have been disrupted by a major collision at that time and that the remaining parent body fragments were too small to develop substantial regoliths (e.g., Heymann, 1967; Crabb and Schultz, 1981). Such a disruption would have exposed a large amount of L chondrite bedrock, some of which would consist of fragmental breccias that acquired melt-rock clasts very early in solar system history, prior to metamorphism. The exposed bedrock would serve as a potential target for sporadic meteoroid impacts to produce a few fragmental breccias with unmetamorphosed melt-rock clasts. The high proportion of genomict brecciated LL chondrites reflects a complex collisional history, probably including several episodes of parent body disruption and gravitational reassembly. Differences in the abundances of different kinds of breccias among the ordinary chondrite groups are probably due to the stochastic nature of major asteroidal collisions.