The Advanced Regional Prediction System (ARPS) – A multi-scale nonhydrostatic atmospheric simulation and prediction model. Part I: Model dynamics and verification

Summary  A completely new nonhydrostatic model system known as the Advanced Regional Prediction System (ARPS) has been developed in recent years at the Center for Analysis and Prediction of Storms (CAPS) at the University of Oklahoma. The ARPS is designed from the beginning to serve as an effective tool for basic and applied research and as a system suitable for explicit prediction of convective storms as well as weather systems at other scales. The ARPS includes its own data ingest, quality control and objective analysis packages, a data assimilation system which includes single-Doppler velocity and thermodynamic retrieval algorithms, the forward prediction component, and a self-contained post-processing, diagnostic and verification package. The forward prediction component of the ARPS is a three-dimensional, nonhydrostatic compressible model formulated in generalized terrain-following coordinates. Minimum approximations are made to the original governing equations. The split-explicit scheme is used to integrate the sound-wave containing equations, which allows the horizontal domain-decomposition strategy to be efficiently implemented for distributed-memory massively parallel computers. The model performs equally well on conventional shared-memory scalar and vector processors. The model employs advanced numerical techniques, including monotonic advection schemes for scalar transport and variance-conserving fourth-order advection for other variables. The model also includes state-of-the-art physics parameterization schemes that are important for explicit prediction of convective storms as well as the prediction of flows at larger scales. Unique to this system are the consistent code styling maintained for the entire model system and thorough internal documentation. Modern software engineering practices are employed to ensure that the system is modular, extensible and easy to use. The system has been undergoing real-time prediction tests at the synoptic through storm scales in the past several years over the continental United States as well as in part of Asia, some of which included retrieved Doppler radar data and hydrometeor types in the initial condition. As the first of a two-part paper series, we describe herein the dynamic and numerical framework of the model, together with the subgrid-scale turbulence and the PBL parameterization. The model dynamic and numerical framework is then verified using idealized and realistic mountain flow cases and an idealized density current. Other physics parameterization schemes will be described in Part II, which is followed by verification against observational data of the coupled soil-vegetation model, surface layer fluxes and the PBL parameterization. Applications of the model to the simulation of an observed supercell storm and to the prediction of a real case are also found in Part II. In the latter case, a long-lasting squall line developed and propagated across the eastern part of the United States following a historical number of tornado outbreak in the state of Arkansas. Received April 14, 2000 Revised July 17, 2000     

The oxidation state diagram — a potential tool for studying redox chemistry in seawater

With minor modifications, the oxidation state diagram may be used to study the redox chemistry of elements in marine waters. In addition to indicating the most stable oxidation state of an element, the diagram also shows the possible existence of metastable species. Oxidation state diagrams of the nitrogen and the manganese systems are presented. The inferences drawn from these diagrams are consistent with our present understanding of the marine geochemistry of these elements.     

The stability of dissolved inorganic species of iodine in seawater

An oxidation state diagram was used to study the relative stability of inorganic iodine species in an aqueous system. It is shown that although iodate is the most stable form, iodide may exist as a metastable form in a basic solution. Molecular iodine may undergo disproportionation to form iodide and iodate. Results from laboratory studies suggest that molecular iodine is rapidly taken up by seawater, and hypoiodite is probably formed. Hypoiodite is also unstable in seawater, and may react with organic matter or undergo autodecomposition. Direct reactions between molecular iodine and organic matter were not observed.     

Modelling of vibratory pile driving in sand

A model for vibro-driving of rigid piles in sand is proposed incorporating the interaction of the vibrator–pile–soil system. The vibro-driver force and the non-linear soil resistance (dynamic f–w and q–w relationships) have been quantified in terms of in situ stress, relative density and particle size. The influence of relative density (0·65 and 0·90), particle size (effective grain size of 0·2 mm and 1·2 mm) and in situ stress (up to 20 psi) on the vibro-driving of a closed-ended pipe pile was investigated experimentally using a large-scale laboratory testing system. The vibro-driving model predicted the observed driving history reasonably well. A computer program (UHVIBRO) has been developed to analyse vibratory pile driving using the dynamic soil resistance relationships developed along with correlation factors from the systematic laboratory study.     

Total selenium and selenium (IV) in the James River estuary and southern Chesapeake Bay

The concentrations of total selenium (Se) and Se (IV) were determined in the surface waters of 30 stations located in the James River and southern Chesapeake Bay. The concentrations of total Se and Se (IV) ranged from 0·28 to 1·91 nM and from 0·07 to 1·36 nM, respectively, between salinities of 31·78 and 0·06‰. The concentration of Se (VI), calculated as the difference between the concentrations of total Se and Se (IV), ranged from 0·08 to 0·67 nM. While total Se seemed to be conservative in this study area at salinities above 0·36‰, Se (IV) might have been removed during estuarine mixing. The removal of Se (IV) occurred primarily at salinities below 4‰ possibly via the oxidation of Se (IV) to Se (VI).     

Numerical simulation of unsteady fluid flow and propagation of a circular hydraulic fracture

A computational procedure is developed for solving the problem of a circular hydraulic fracture propagating under the action of frac-0fluid being pumped in at a central wellbore. The crack is modelled as continuous distributions of ring dislocations and the resulting elasticity singular integral equation is solved numerically. The fluid flow equations are approximated by local and global interpolation finite difference schemes. The coupling between elasticity and fluid flow is handled numerically, by, two different algorithms: one iterates on crack tip velocity whereas the other varies the time step size until it agrees with the chosen increment in crack length. Sample results are given; it is found that the velocity algorithm is computationally more, efficient and more stable. The model allows detailed tracing of pressure distribution and fluid flow in the fracture, even under complex conditions of cyclic injection and fluid rheology. It may serve as a stand-alone model of (horizontal) hydrafracs–especially at shallow depths–or it may be used as a reference frame to test the various numerical formulation/algorithms required for the ongoing development of a fully 3-D hydrafrac simulator.     

Cesium-137, metals and organic carbon in the sediments of the James River estuary,Virginia

A wide variety of sedimentary subenvironments are found within a 10-km stretch of James River including a flood dominated channel (Rocklanding Channel) and its bank (Rocklanding Shoal), a shoal with a water depth of 1 m separating two channels (Point of Shoals), an ebb-dominated channel (Burwell Bay Channel) and its bank (Burwell Bay Bank) and a tributary (Warwick River). The concentrations of Cs-137, Cu, Pb, Zn and organic carbon in the fine-grained sediments (i.e. < 63 μm) and the amount of fine-grained sediments in eight cores covering these subenvironments were determined. The sedimentation rates, estimated by Cs-137 geochronology, range from 0·4 to > 3 cm year^?1. The sedimentation rates in the Burwell Bay region are two or more, times those in the Point of Shoals and in the Rocklanding Channel and Shoal, reflecting the weaker currents in the Burwell Bay region. These sedimentation rates agree well with those obtained independently by measuring changes in the bathymetry of this area between 1873 and 1943. The concentrations of Cs-137, Cu, Pb, Zn and organic carbon in surface sediments vary by a factor of two to three. The concentrations are higher in the Burwell Bay region, probably as a result of the higher rates of accumulation of recently formed sediments in these subenvironments. The inventories of fine-grained sediments and of Cs-137, Cu, Pb, Zn and organic carbon accumulated since 1954 are also up to an order of magnitude higher in the Burwell Bay region. Although the concentrations of fine-grained sediments in three cores obtained in this region are similar, the inventories still vary by a factor of two to three. The inventories of Cs-137, Cu, Pb, Zn, organic carbon and fine-grained sediments correlate well with each other indicating that Cs-137 can be a useful tracer for studying the fate of these metals and organic carbon in estuarine environments. The inhomogeneity of the concentrations and inventories of the different elements along a 10-km segment of a river suggests that a closely-spaced sampling programme is essential for characterizing the sedimentary provinces within an estuary. The concentrations of Cs-137, metals and organic carbon in the coarse-grained sediments (i.e. > 63 μm) are considerably lower than those in fine-grained sediments. Thus, the contribution of coarse-grained sediments to the total inventory of these elements is small.     

The old alluvium of Singapore and the extinct drainage system to the South China Sea

In both sand quarry faces and boreholes the Old Alluvium of Singapore is a matrix supported pebbly sand with minor beds of better sorted sands and clays. The beds can be grouped into four textural classes (pebbles, coarse sand with fine pebbles, medium to coarse sand, clay, and silt), each with characteristic sedimentary structures. The deposit appears to be the proximal facies of an ancient braided river alluvium of possible Pleistocene age. This widespread but scattered alluvium is found both on land and offshore on the Sunda Shelf in Southeast Asia. It is believed to have been deposited during low sea levels. In Singapore, the mineralogy of the grains forming the Old Alluvium suggests a mixed provenance of granitic and low-grade metamorphic origin. The volume and freshness of the deposited material indicate an environment of considerable relief, seasonal rainfall, and extremely active erosional processes, conditions considerably different from the current ones. The Old Alluvium was probably deposited by seasonal rivers which experienced periodic large floods. Such regional conditions could have prevailed over much of Southeast Asia at the time of deposition of the Old Alluvium.