Parallel Simulation in Subsurface Hydrology: Evaluating the Performance of Modeling Computers

Monte Carlo uncertainty analysis, model calibration and optimization applications in hydrology, usually involve a very large number of forward transient model solutions, often resulting in computational bottlenecks. Parallel processing can significantly reduce overall simulation time, benefiting from the architecture of modern computers. This work investigates system performance using two realistic flow and transport modeling scenarios, applied to various modeling hardware, to provide information on the expected performance of parallel simulations and inform investment decisions. We investigate how performance, measured in terms of speedup and efficiency, changes with increasing number of parallel processes. We conclude that the maximum performance achieved by parallelization can range from 40% to 100% of the theoretical limit, with the lower increases associated with multi-CPU servers. The number of parallel processes required to maximize performance is application dependent, and in contrast to common practice, often needs to be significantly larger than the total number of system CPU cores. Further testing is required to better understand how the physical problem being simulated affects the optimal number of parallel processes needed. Finally, when laptops are considered for modeling applications, careful consideration should be given not only to the specifications but also to the intended use designated by the manufacturer.     

Multistop-based measurements of accessibility in a GTS environment

Abstract

The opportunities available at a demand location are usually measured as the costs of reaching a specified critical number of facilities from that location. This method does not however, account for multistop trips nor for differences in the diversity of supply at the level of individual facilities. In this paper we introduce an alternative measurement method that overcomes these shortcomings. In this method the probability of successfully visiting a specific facility is assumed to be a function of the diversity of supply provided. Trip routes are constructed that have an acceptable probability of success. Then, the expected costs of travelling the optimum route are determined as an indicator of spatial opportunities. The proposed method has been implemented in a GIS environment, using typical GIS data and GIS tools for spatial analysis and display. The results of a case study indicate that the new method, compared to current methods, may lead to different evaluations of the level of opportunities at demand locations.     

Active and passive ocean regimes in a low‐order climate model

A low‐order climate model is studied which combines the Lorenz‐84 model for the atmosphere on a fast time scale and a box model for the ocean on a slow time scale. In this climate model, the ocean is forced strongly by the atmosphere. The feedback to the atmosphere is weak. The behaviour of the model is studied as a function of the feedback parameters. We find regions in parameter space with dominant atmospheric dynamics, i.e., a passive ocean, as well as regions with an active ocean, where the oceanic feedback is essential for the qualitative dynamics. The ocean is passive if the coupled system is fully chaotic. This is illustrated by comparing the Kaplan–Yorke dimension and the correlation dimension of the chaotic attractor to the values found in the uncoupled Lorenz‐84 model. The active ocean behaviour occurs at parameter values between fully chaotic and stable periodic motion. Here, intermittency is observed. By means of bifurcation analysis of periodic orbits, the intermittent behaviour, and the rôle played by the ocean model, is clarified. A comparison of power spectra in the active ocean regime and the passive ocean regime clearly shows an increase of energy in the low frequency modes of the atmospheric variables. The results are discussed in terms of itinerancy and quasi‐stationary states observed in realistic atmosphere and climate models.     

Toxic metal adsorption in the triassic sandstone aquifer of the English Midlands

Laboratory experimental results show that adsorption can be an important mechanism in taking up heavy metals contained in groundwaters in the sandstones. The adsorption can be described in terms of Freundlich isotherms and can be simply modeled mathematically. The adsorption function has been integrated into a mass balance equation, and tested using a laboratory sand tank model and in a small field study. The results show that the transport-adsorption models are not easily applied in the field because of hydrogeological complications.     

A MODIFIED WEATHER GENERATOR(WGEN)FOR THE GENERATION OF DAILY WEATHER VARIABLES

Long series of daily weather data are frequently needed to evaluate diachronic climatic impacton water resources,the effects of watershed changes on hydrology and to use in a variety ofweather and general circulation models.A computer generation model called WGEN that wasdeveloped in the United States of America,was modified in this study and applied to Nigerianweather data spanning the period between 1969 to 1988 and covering 17 sites located in thedifferent climatic zones in Nigeria.The model generated the monthly mean maximum and minimumtemperatures,solar radiation,total rainfall,and number of wet and dry days with high accuracy,95% of the times.The Mann-Whitney U-test revealed that the number of months per year forwhich observed and generated weather variables were significant,was less than 4 in majority of thesites.     

The African Climate as Predicted by the IAP Grid-Point Nine-Layer Atmospheric General Circulation Model (IAP-9L-AGCM)

ＴｈｅＡｆｒｉｃａｎＣｌｉｍａｔｅａｓＰｒｅｄｉｃｔｅｄｂｙｔｈｅＩＡＰＧｒｉｄ－ＰｏｉｎｔＮｉｎｅ－ＬａｙｅｒＡｔｍｏｓｐｈｅｒｉｃＧｅｎｅｒａｌＣｉｒｃｕｌａｔｉｏｎＭｏｄｅｌ（ＩＡＰ－９Ｌ－ＡＧＣＭ）ＣｈｉｎｅｋｅＴｈｅｏＣｈｉｄｉｅｚｉｅ①，...     

Intercomparison of Two Hydroxyl Radical Measurement Techniques at the Atmosphere Simulation Chamber SAPHIR

At the atmosphere simulation chamber SAPHIR in Jülich both Laser-Induced Fluorescence Spectroscopy (LIF) and Long-Path Differential Optical Laser Absorption Spectroscopy (DOAS) are operational for the detection of OH radicals at tropospheric levels. The two different spectroscopic techniques were compared within the controlled environment of SAPHIR based on all simultaneous measurements acquired in 2003 (13 days). Hydroxyl radicals were scavenged by added CO during four of these days in order to experimentally check the calculated precisions at the detection limit. LIF measurements have a higher precision (σ= 0.88×10⁶ cm^–3) and better time resolution (Δt = 60 s), but the DOAS method (σ= 1.24×10⁶ cm^–3, Δt = 135 s) is regarded as primary standard for comparisons because of its good accuracy. A high correlation coefficient of r = 0.95 was found for the whole data set highlighting the advantage of using a simulation chamber. The data set consists of two groups. The first one includes 3 days, where the LIF measurements yield (1 – 2) ×10⁶ cm^–3 higher OH concentrations than observed by the DOAS instrument. The experimental conditions during these days are characterized by increased NO_x concentration and a small dynamic range in OH. Excellent agreement is found within the other group of 6 days. The regression to the combined data of this large group yields unity slope without a significant offset.     

Application of SAXS and SANS in evaluation of porosity, pore size distribution and surface area of coal

This paper discusses the applicability of small angle X-ray scattering (SAXS) and small angle neutron scattering (SANS) techniques for determining the porosity, pore size distribution and internal specific surface area in coals. The method is noninvasive, fast, inexpensive and does not require complex sample preparation. It uses coal grains of about 0.8 mm size mounted in standard pellets as used for petrographic studies.Assuming spherical pore geometry, the scattering data are converted into the pore size distribution in the size range 1 nm (10 Å) to 20 μm (200,000 Å) in diameter, accounting for both open and closed pores. FTIR as well as SAXS and SANS data for seven samples of oriented whole coals and corresponding pellets with vitrinite reflectance (R_o) values in the range 0.55% to 5.15% are presented and analyzed. Our results demonstrate that pellets adequately represent the average microstructure of coal samples.The scattering data have been used to calculate the maximum surface area available for methane adsorption. Total porosity as percentage of sample volume is calculated and compared with worldwide trends. By demonstrating the applicability of SAXS and SANS techniques to determine the porosity, pore size distribution and surface area in coals, we provide a new and efficient tool, which can be used for any type of coal sample, from a thin slice to a representative sample of a thick seam.     

The role of antelope trampling on scarp erosion in a hyper-arid environment, Skeleton Coast, Namibia

The role of animal trampling in a hyper-arid environment is evaluated. Field observations on fluvial terraces of the Uniab fan in Namibia suggest that animal impact on level surfaces with gravel armouring is minimal. However, animal trampling along tracks across scarp slopes results in scarp recess and initiate cutback development. With continued growth ephemeral fluvial processes become more dominant and obscure the animal impact. Animal trampling effects are pronounced in hyper-arid environments due to the low rates of other slope processes. However, it is argued that large herd migration may have contributed significantly to cutback development in other environments, especially in Africa, but may be very difficult to substantiate by landform analysis.