Effect of spatial variation characteristics on contouring of design storm depth

The ordinary kriging method, a geostatistical interpolation technique, was applied for developing contour maps of design storm depth in northern Taiwan using intensity–duration–frequency (IDF) data. Results of variogram modelling on design storm depths indicate that the design storms can be categorized into two distinct storm types: (i) storms of short duration and high spatial variation and (ii) storms of long duration and less spatial variation. For storms of the first category, the influence range of rainfall depth decreases when the recurrence interval increases, owing to the increasing degree of their spatial independence. However, for storms of the second category, the influence range of rainfall depth does not change significantly and has an average of approximately 72 km. For very extreme events, such as events of short duration and long recurrence interval, we do not recommend usage of the established design storm contours, because most of the interstation distances exceed the influence ranges. Our study concludes that the influence range of the design storm depth is dependent on the design duration and recurrence interval and is a key factor in developing design storm contours. Copyright © 2003 John Wiley & Sons, Ltd.     

Boundary curves for families of planar orbits

For monoparametric familiesf(x,y)=c of planar orbits, created by a planar potentialV(x,y), we introduce the notion of the family boundary curves (FBC). All members of the familyf(x,y)=c are traced in an allowable region of thexy plane, defined by the corresponding FBC, with total energyE=E(c) varying along the family. Family boundary curves are also found for two-parametric familiesf(x,y,b)=c. The relation of equilibrium points and asymptotic orbits, possibly possessed by the potentialV(x,y), to be FBC is studied.     

The origin of the correlation between the spin parameter and the baryon fraction of galactic disks

The puzzling correlation between the spin parameter λ of galactic disks and the disk-to-halo mass fraction f _disk is investigated. It is shown that such a correlation arises naturally from large uncertainties in determining the virial masses of dark matter halos. This result indicates that halo properties derived from fits to observed rotation curves are still very uncertain. An analysis of λ versus f _disk as function of the adopted halo virial mass shows that for reasonable halo concentrations f _disk ≈ 0.01 - 0.07 which is significantly smaller than the universal baryon fraction. Most of the available gas either never settled into the galactic disks or was ejected subsequently. This revised version was published online in August 2006 with corrections to the Cover Date.     

The Homoclinic Tangle of A 1:2 Resonance in a 2-D Hamiltonian System

We study in great detail the geometry of the homoclinic tangle, with respect to the energy, corresponding to an unstable periodic orbit of type 1:2, on a surface of section representing a 2-D Hamiltonian system. The tangle consists of two resonance areas, in contrast with the tangles of type-l or -{l, m, k, x = 0} considered in previous studies, that consist of only one resonance area. We study the intersections of the inner and outer lobes of the same resonance area and of the two resonance areas. The intersections of the lobes follow certain rules. The detailed study of these rules allows us to derive quantitative relations about the number of intersections and to understand the complex behavior of the higher order lobes by studying the lower order lobes. We find 1st, 2nd, 3rd, etc. order intersections formed by lobes making 1, 2, 3, etc. turns around an island. After a sufficiently high order of iterations a lobe may intersect its image and thus produce a Poincaré recurrence. Numerical results for a wide interval of energies are presented. The number of intersections changes through tangencies. In any finite interval of the energy between two tangencies of 1st order, an infinite number of higher order tangencies occur and thus, according to the Newhouse theorem, there exist nearby islands of stability.     

Restricted Three-Body Problem: An Approximation of its General Solution – Part One – the Manifold of Symmetric Periodic Solutions

This paper gives the results of a programme attempting to exploit ‘la seule bréche’ (Poincaré, 1892, p. 82) of non-integrable systems, namely to develop an approximate general solution for the three out of its four component-solutions of the planar restricted three-body problem. This is accomplished by computing a large number of families of ‘solutions précieuses’ (periodic solutions) covering densely the space of initial conditions of this problem. More specifically, we calculated numerically and only for μ = 0.4, all families of symmetric periodic solutions (1st component of the general solution) existing in the domain D:(x ₀ ∊ [−2,2],C ∊ [−2,5]) of the (x ₀, C) space and consisting of symmetric solutions re-entering after 1 up to 50 revolutions (see graph in Fig. 4). Then we tested the parts of the domain D that is void of such families and established that they belong to the category of escape motions (2nd component of the general solution). The approximation of the 3rd component (asymmetric solutions) we shall present in a future publication. The 4th component of the general solution of the problem, namely the one consisting of the bounded non-periodic solutions, is considered as approximated by those of the 1st or the 2nd component on account of the `Last Geometric Theorem of Poincaré' (Birkhoff, 1913). The results obtained provoked interest to repeat the same work inside the larger closed domain D:(x ₀ ∊ [−6,2], C ∊ [−5,5]) and the results are presented in Fig. 15. A test run of the programme developed led to reproduction of the results presented by Hénon (1965) with better accuracy and many additional families not included in the sited paper. Pointer directions construed from the main body of results led to the definition of useful concepts of the basic family of order n, n = 1, 2,… and the completeness criterion of the solution inside a compact sub-domain of the (x ₀, C) space. The same results inspired the ‘partition theorem’, which conjectures the possibility of partitioning an initial conditions domain D into a finite set of sub-domains D i that fulfill the completeness criterion and allow complete approximation of the general solution of this problem by computing a relatively small number of family curves. The numerical results of this project include a large number of families that were computed in detail covering their natural termination, the morphology, and stability of their member solutions. Zooming into sub-domains of D permitted clear presentation of the families of symmetric solutions contained in them. Such zooming was made for various values of the parameter N, which defines the re-entrance revolutions number, which was selected to be from 50 to 500. The areas generating escape solutions have being investigated. In Appendix A we present families of symmetric solutions terminating at asymptotic solutions, and in Appendix B the morphology of large period symmetric solutions though examples of orbits that re-enter after from 8 to 500 revolutions. The paper concludes that approximations of the general solution of the planar restricted problem is possible and presents such approximations, only for some sub-domains that fulfill the completeness criterion, on the basis of sufficiently large number of families.     

Hysteresis‐based analysis of overland metal transport

Introducing a concept of equivalent mass depth of flow, this study describes the phenomenon of non‐point source pollutant (metal) transport for pavement (or overland) flow in analogy with wave propagation in wide open channels. Hysteretic and normal mass rating curves are developed for runoff rate and mass of 12 dissolved and particulate‐bound metal elements (pollutants) using the rainfall‐runoff and water quality data of the 15 × 20 m² instrumented pavement in Cincinnati, USA. Normal mass rating curves developed for easy computation of pollutant load are found to be of a form similar to Manning's, which is valid for open channel flows. Based on the hysteresis analysis, wave types for dissolution and mixing of particulate‐bound metals are identified. The analysis finds that the second‐order partial‐differential equation normally used for metal transport does not have the efficacy to describe fully the strong non‐linear phenomena such as is described for various metal elements by dynamic waves. In addition, the proportionality concept of the popular SCS‐CN concept is extended for determining the potential maximum metal mass M_p of all the 12 elements transported by a rain storm and related to the antecedent dry period (ADP). For the primary metal zinc element, M_p is found to increase with the ADP and vice versa. Copyright © 2003 John Wiley & Sons, Ltd.     

A comparative study of four vulnerability mapping methods in a detritic aquifer under mediterranean climatic conditions

This paper presents the results of a comparative study relating to the application of four vulnerability mapping methods, GOD, AVI, DRASTIC and SINTACS, in a pilot detritic aquifer situated in NW Morocco, known as the Martil–Alila aquifer. The principal objective of this work is to determine the most suitable such methods for this aquifer type within a Mediterranean context, and to show the effect of the rainfall variations that are characteristic of the Mediterranean climate on the degree of vulnerability. The methods applied distinguish five classes of vulnerability, these being irregularly divided up in space, with the division varying according to the method in question. The vulnerability maps obtained by the different methods strongly suggest that the eastern half of the aquifer is more vulnerable to contamination than the western half, for all hydrological situations. The effect of climatic conditions on the degree of vulnerability is well represented by the DRASTIC, according to which the aquifer is moderately to strongly vulnerable during humid hydrological years and weakly to moderately vulnerable during dry ones. For the other methods, this climatic effect is limited to the area occupied by the two predominant classes (“High” and “Low” for GOD and “High” and “Moderate” for SINTACS) while it is null for AVI. In conclusion, DRASTIC appears the most suitable for mapping the vulnerability to contamination of Mediterranean coastal detritic aquifers such as the Martil–Alila aquifer.     

Reducing hazard vulnerability through local government engagement and action

The concept of a natural hazard is a human construct. It is the interaction with human communities and settlements that defines a natural phenomenon as a natural hazard. Thus the end point of hazard mitigation and hazard vulnerability assessment must involve an attempt to reduce, or mitigate, the impact of the natural hazard on human communities. The responsibility to mitigate hazard impact falls primarily upon governments and closely connected non-government and private institutional agencies. In particular, it is most often local government that takes the responsibility for safeguarding its own communities, infrastructure and people. Hazard vulnerability of specific local communities is best assessed by the local government or council, which then faces the responsibility to translate that assessment into community education and infrastructural safeguards for hazard mitigation. This paper illustrates the process of local government engagement in hazard mitigation in Australia, through the Natural Disaster Risk Management Studies, as a first step towards natural disaster reduction.     

河口Si输运过程研究

基于拉格朗日余流及其输运过程的一种三维空间弱非线性理论，进一步假定Ｓｉ在河口的无因次化学转移项量级为ｋ＾２，给出了Ｓｉ的长期输运方程。对于一种二维模型河口，数值求解了零阶天文潮，欧拉余流、斯托克斯漂移、拉格朗日余流和盐度的分布；分别计算了平均逼留时间为１３，３０，６０，１８０ｄ等所对应的Ｓｉ浓度分布；给出了Ｓｉ－ｓ相关图，讨论了河口中Ｓｉ的保守性问题。