High-performance direct gravitational N-body simulations on graphics processing units

We present the results of gravitational direct N-body simulations using the commercial graphics processing units (GPU) NVIDIA Quadro FX1400 and GeForce 8800GTX, and compare the results with GRAPE-6Af special purpose hardware. The force evaluation of the N-body problem was implemented in Cg using the GPU directly to speed-up the calculations. The integration of the equations of motions were, running on the host computer, implemented in C using the 4th order predictor–corrector Hermite integrator with block time steps. We find that for a large number of particles (N  10⁴) modern graphics processing units offer an attractive low cost alternative to GRAPE special purpose hardware. A modern GPU continues to give a relatively flat scaling with the number of particles, comparable to that of the GRAPE. The GRAPE is designed to reach double precision, whereas the GPU is intrinsically single-precision. For relatively large time steps, the total energy of the N-body system was conserved better than to one in 10⁶ on the GPU, which is impressive given the single-precision nature of the GPU. For the same time steps, the GRAPE gave somewhat more accurate results, by about an order of magnitude. However, smaller time steps allowed more energy accuracy on the grape, around 10⁻¹¹, whereas for the GPU machine precision saturates around 10⁻⁶ For N  10⁶ the GeForce 8800GTX was about 20 times faster than the host computer. Though still about a factor of a few slower than GRAPE, modern GPUs outperform GRAPE in their low cost, long mean time between failure and the much larger onboard memory; the GRAPE-6Af holds at most 256k particles whereas the GeForce 8800GTX can hold 9 million particles in memory.     

Biostratigraphic and aminostratigraphic constraints on the age of the Middle Pleistocene glacial succession in north Norfolk,UK

Considerable debate surrounds the age of the Middle Pleistocene glacial succession in East Anglia following some recent stratigraphical reinterpretations. Resolution of the stratigraphy here is important since it not only concerns the glacial history of the region but also has a bearing on our understanding of the earliest human occupation of north‐western Europe. The orthodox consensus that all the tills were emplaced during the Anglian (Marine Isotope Stage (MIS) 12) has recently been challenged by a view assigning each major till to a different glacial stage, before, during and after MIS 12. Between Trimingham and Sidestrand on the north Norfolk coast, datable organic sediments occur immediately below and above the glacial succession. The oldest glacial deposit (Happisburgh Till) directly overlies the ‘Sidestrand Unio‐bed’, here defined as the Sidestrand Hall Member of the Cromer Forest‐bed Formation. Dating of these sediments therefore has a bearing on the maximum age of the glacial sequence. This paper reviews the palaeobotany and describes the faunal assemblages recovered from the Sidestrand Unio‐bed, which accumulated in a fluvial environment in a fully temperate climate with regional deciduous woodland. There are indications from the ostracods for weakly brackish conditions. Significant differences are apparent between the Sidestrand assemblages and those from West Runton, the type site of the Cromerian Stage. These differences do not result from contrasting facies or taphonomy but reflect warmer palaeotemperatures at Sidestrand and a much younger age. This conclusion is suggested by the higher proportion of thermophiles at Sidestrand and the occurrence of a water vole with unrooted molars (Arvicola) rather than its ancestor Mimomys savini with rooted molars. Amino acid racemisation data also indicate that Sidestrand is significantly younger than West Runton. These data further highlight the stratigraphical complexity of the ‘Cromerian Complex’ and support the conventional view that the Happisburgh Till was emplaced during the Anglian rather than the recently advanced view that it dates from MIS 16. Moreover, new evidence from the Trimingham lake bed (Sidestrand Cliff Formation) above the youngest glacial outwash sediments (Briton's Lane Formation) indicates that they also accumulated during a Middle Pleistocene interglacial – probably MIS 11. All of this evidence is consistent with a short chronology placing the glacial deposits within MIS 12, rather than invoking multiple episodes of glaciation envisaged in the ‘new glacial stratigraphy’ during MIS 16, 12, 10 and 6. Copyright © 2009 John Wiley & Sons, Ltd.     

粤港澳大湾区发展的理论框架与发展战略探究

粤港澳大湾区是在“一带一路”倡议下,由珠三角经济圈的强化合作而产生的新地理概念。作为新兴的全球化港湾,粤港澳大湾区如何突破新自由主义框架理论,探索出符合社会主义市场经济体制的发展路径,是当前决策者和规划者需要思考的重要问题。本文简要梳理了现有的城市群、区域增长极理论,新经济地理集群理论,全球产业链网络理论并结合大湾区现状优势,对湾区发展战略定位与发展路径进行探究,认为可通过粤港澳三方通力合作,创造一个内生型的经济与产业本土增长模式为主要方向。具体表现为：首先,要以先进制造业为立足点、实现自我创新的产业升级,形成完善的制造业产业链,并成为全球生产网络的重要节点与区域性枢纽;其次,大湾区还应利用自身的科研与教育、金融与创新资源优势,推进“一国两制三关税区”、尤其是与香港在大湾区建设中的全方位参与联动,完善产权制度,加强合规和规则经济的市场经济体系建设,将大湾区打造为中国的科创中心及全球性金融中心。     

Vulnerability of population and transportation infrastructure at the east bank of Delaware Bay due to coastal flooding in sea-level rise conditions

Catastrophic flooding associated with sea-level rise and change of hurricane patterns has put the northeastern coastal regions of the United States at a greater risk. In this paper, we predict coastal flooding at the east bank of Delaware Bay and analyze the resulting impact on residents and transportation infrastructure. The three-dimensional coastal ocean model FVCOM coupled with a two-dimensional shallow water model is used to simulate hydrodynamic flooding from coastal ocean water with fine-resolution meshes, and a topography-based hydrologic method is applied to estimate inland flooding due to precipitation. The entire flooded areas with a range of storm intensity (i.e., no storm, 10-, and 50-year storm) and sea-level rise (i.e., current, 10-, and 50-year sea level) are thus determined. The populations in the study region in 10 and 50 years are predicted using an economic-demographic model. With the aid of ArcGIS, detailed analysis of affected population and transportation systems including highway networks, railroads, and bridges is presented for all of the flood scenarios. It is concluded that sea-level rise will lead to a substantial increase in vulnerability of residents and transportation infrastructure to storm floods, and such a flood tends to affect more population in Cape May County but more transportation facilities in Cumberland County, New Jersey.     

Using Spatial Data to Explore Landholder Awareness and Concern about Dryland Salinity

This paper draws upon research in Victoria's Goulburn Broken Catchment exploring landholder responses to dryland salinity. It has been assumed that part of the explanation for limited adoption of recommended practices was that landholders were unaware of the extent of dryland salinity. Socio-economic data from a mail survey was combined in a Geographic Information System (GIS) with other layers, including the location of discharge sites and depth to saline ground water. Comparisons of expert maps and landholder identified salinity sites suggested that landholders in the upper catchment had excellent knowledge of the current extent of salinity on their property. At the same time, the expert maps failed to predict half of the saline-affected sites identified by landholders. It seems that the extensive community education effort undertaken in this region had been successful in raising salinity awareness. Our research also highlighted that most landholders were not concerned about the impacts of dryland salinity and appear to believe they can 'live with salt'. This is an important issue because the small, diffused amounts of salt exported from these properties in the upper catchment are/will have a substantial impact downstream. The research methodology and findings have important policy and management implications, and these are discussed.     

Fine‐scale diet of the Australian sea lion (Neophoca cinerea) using DNA‐based analysis of faeces

We applied DNA‐based faecal analysis to determine the diet of female Australian sea lions (n = 12) from two breeding colonies in South Australia. DNA dietary components of fish and cephalopods were amplified using the polymerase chain reaction and mitochondrial DNA primers targeting the short (~100 base pair) section of the 16S gene region. Prey diversity was determined by sequencing ~50 amplicons generated from clone libraries developed for each individual. Faecal DNA was also combined and cloned from multiple individuals at each colony and fish diversity determined. Diets varied between individuals and sites. Overall, DNA analysis identified a broad diversity of prey comprising 23 fish and five cephalopod taxa, including many species not previously described as prey of the Australian sea lion. Labridae (wrasse), Monacanthidae (leatherjackets) and Mullidae (goat fish) were important fish prey taxa. Commonly identified cephalopods were Octopodidae (octopus), Loliginidae (calamary squid) and Sepiidae (cuttlefish). Comparisons of fish prey diversity determined by pooling faecal DNA from several samples provided a reasonable but incomplete resemblance (55–71%) to the total fish diversity identified across individual diets at each site. Interpretation of diet based on the recovery of prey hard‐parts identified one cephalopod beak (Octopus sp.) and one fish otolith (Parapriacanthus elongatus). The present study highlights the value of DNA‐based analyses and their capabilities to enhance information of trophic interactions.     

Co-diagenesis of iron and phosphorus in hydrothermal sediments from the southern East Pacific Rise: Implications for the evaluation of paleoseawater phosphate concentrations

We present a detailed study of the co-diagenesis of Fe and P in hydrothermal plume fallout sediments from ∼19°S on the southern East Pacific Rise. Three distal sediment cores from 340-1130 km from the ridge crest, collected during DSDP Leg 92, were analysed for solid phase Fe and P associations using sequential chemical extraction techniques. The sediments at all sites are enriched in hydrothermal Fe (oxyhydr)oxides, but during diagenesis a large proportion of the primary ferrihydrite precipitates are transformed to the more stable mineral form of goethite and to a lesser extent to clay minerals, resulting in the release to solution of scavenged P. However, a significant proportion of this P is retained within the sediment, by incorporation into secondary goethite, by precipitation as authigenic apatite, and by readsorption to Fe (oxyhydr)oxides. Molar P/Fe ratios for these sediments are significantly lower than those measured in plume particles from more northern localities along the southern East Pacific Rise, and show a distinct downcore decrease to a depth of ∼12 m. Molar P/Fe ratios are then relatively constant to a depth of ∼35 m. The Fe and P speciation data indicate that diagenetic modification of the sediments is largely complete by a depth of 2.5 m, and thus depth trends in molar P/Fe ratios can not solely be explained by losses of P from the sediment by diffusion to the overlying water column during early diagenesis. Instead, these sediments are likely recording changes in dissolved P concentrations off the SEPR, possibly as a result of redistribution of nutrients in response to changes in oceanic circulation over the last 10 million years. Furthermore, the relatively low molar P/Fe ratios observed throughout these sediments are not necessarily solely due to losses of scavenged P by diffusion to the overlying water column during diagenesis, but may also reflect post-depositional oxidation of pyrite originating from the volatile-rich vents of the southern East Pacific Rise. This study suggests that the molar P/Fe ratio of oxic Fe-rich sediments may serve as a proxy of relative changes in paleoseawater phosphate concentrations, particularly if Fe sulfide minerals are not an important component during transport and deposition.     

Geomorphologic and Palaeoenvironmental Development of Holocene Perched Coastal Dune Systems in Brittany, France

Two perched dune systems are investigated in the Baie d'Audierne (Finistère, Brittany, northwest France). A maximum age of 4545–4160 cal. yr bp for the onset of sand accretion is provided through radiocarbon (¹⁴C) analysis of an organic-rich basal palaeosol, corresponding closely to deceleration, but before stabilisation, of regional Holocene relative sea levels (RSLs). Molluscan analysis through the dune sand indicates an initial bare sand or sparsely vegetated surface, subsequently maturing and colonised by scrub vegetation, and finally returning to more open dune grassland conditions. These data suggest a three-stage Perched Dune Development Model (PDDM), which may be applicable to perched dunes on at least a regional scale: (1) initial sand sheet inundation as RSL rises; (2) stabilisation of RSL allowing dune to mature; and (3) sea-level stabilisation stimulates cliff formation, perching and progressive landward retreat of the dune shore. This study recognises that perched dunes hold greater potential for establishing chronologies for regional Holocene dune development than the more widely studied lowland dune systems, where evidence for early sand inundation may be obscured or lost by modest rises in sea level.     

Statistical landslide susceptibility assessment of the Mansehra and Torghar districts,Khyber Pakhtunkhwa Province,Pakistan

This paper presents laboratory experiments and numerical simulations of effects of submerged obstacles on tsunami-like solitary wave and its run-up. This study was carried out for the breaking and non-breaking solitary waves on 1:19.85 uniform slope which contains a submerged obstacle. New laboratory experiments are performed to describe the mitigation of tsunami amplitude and run-up under the effect of submerged obstacles. We are based on experimental results obtained to validate the numerical model. The numerical modeling using COULWAVE aims essentially to show the effect of the obstacle on the shape of solitary wave and the limit of this effect. Using a multiple nonlinear regression, we have determined a model to estimate height of run-up according to the amplitude of the wave and the obstacle peak depth.