Spatial Periodic Boundary Condition for MODFLOW     

Tariq Laattoe  Vincent E.A. Post  Adrian D. Werner 《Ground water》2014,52(4):606-612

Small‐scale hyporheic zone (HZ) models often use a spatial periodic boundary (SPB) pair to simulate an infinite repetition of bedforms. SPB's are common features of commercially available multiphysics modeling packages. MODFLOW's lack of this boundary type has precluded it from being effectively utilized in this area of HZ research. We present a method to implement the SPB in MODFLOW by development of the appropriate block‐centered finite‐difference expressions. The implementation is analogous to MODFLOW's general head boundary package. The difference is that the terms on the right hand side of the solution equations must be updated with each iteration. Consequently, models that implement the SPB converge best with solvers that perform both inner and outer iterations. The correct functioning of the SPB condition in MODFLOW is verified by two examples. This boundary condition allows users to build HZ‐bedform models in MODFLOW, facilitating further research using related codes such as MT3DMS and PHT3D.  相似文献   

Uncertainty in recurrence rates of large magnitude events due to short historic catalogs   总被引：2，自引：0，他引：2  

Tariq Mohammed  Gail M. Atkinson  Karen Assatourians 《Journal of Seismology》2014,18(3):565-573

Seismic hazard analysis requires knowledge of the recurrence rates of large magnitude earthquakes that drive the hazard at low probabilities of interest for seismic design. Earthquake recurrence is usually determined through studies of the historic earthquake catalogue for a given region. Reliable historic catalogues generally span time periods of 100–200 years in North America, while large magnitude events (M?≥?7) have recurrence rates on the order of hundreds or thousands of years in many areas, resulting in large uncertainty in recurrence rates for large events. Using Monte Carlo techniques and assuming typical recurrence parameters, we simulate earthquake catalogues that span long periods of time. We then split these catalogues into smaller catalogues spanning 100–200 years that mimic the length of historic catalogues. For each of these simulated “historic” catalogues, a recurrence rate for large magnitude events is determined. By comparing recurrence rates from one historic-length catalogue to another, we quantify the uncertainty associated with determining recurrence rates from short historic catalogues. The use of simulations to explore the uncertainty (rather than analytical solutions) allows us flexibility to consider issues such as the relative contributions of aleatory versus epistemic uncertainty, and the influence of fitting method, as well as lending insight into extreme-event statistics. The uncertainty in recurrence rates of large (M?>?7) events is about a factor of two in regions of high seismicity, due to the shortness of historic catalogues. This uncertainty increases greatly with decreasing seismic activity. Uncertainty is dependent on the length of the catalogue as well as the fitting method used (least squares vs. maximum likelihood). Examination of 90th percentile recurrence rates reveals that epistemic uncertainty in the true parameters may cause recurrence rates determined from historic catalogues to be uncertain by a factor greater than 50.  相似文献   

Two dimensional electromagnetic shock structures in dense electron-positron-ion magnetoplasmas     

W. Masood  H. Rizvi  S. Hussain 《Astrophysics and Space Science》2011,332(2):287-299

Linear and nonlinear analysis of low frequency magnetoacoustic waves propagating at an angle θ with the ambient magnetic field are investigated in dense electron-positron-ion (e-p-i) plasmas using the quantum magnetohydrodynamic (QMHD) model. In this regard, a quantum Kadomtsev-Petviashvili-Burgers (KPB) equation is derived in the small amplitude limit. The stability of KPB equation is also presented. The variation of the nonlinear fast and slow magnetoacoustic shock waves with the positron concentration, kinematic viscosity, obliqueness parameter θ, and the magnetic field, are also investigated. It is observed that the aforementioned plasma parameters significantly modify the propagation characteristics of two dimensional nonlinear magnetoacoustic shock waves in dissipative quantum magnetoplasmas. The relevance of the present investigation with regard to dense astrophysical environments is also pointed out.  相似文献   

Sea Level Variations and Geomorphological Changes in the Coastal Belt of Pakistan     

Tariq Masood Ali Khan  D. A. Razzaq  Qamar-Uz-Zaman Chaudhry  Dewan Abdul Quadir  Anwarul Kabir  Majajul Alam Sarker 《Marine Geodesy》2013,36(1-2):159-174

The UNEP in its regional seas program in 1989 has included Pakistan in a group of countries which are vulnerable to the impact of rising sea level. If the present trend of sea level rise (SLR) at Karachi continues, in the next 50 years the sea level rise along the Pakistan Coast will be 50 mm (5 cm). Since the rising rates of sea level at Karachi are within the global range of 1-2 mm/year, the trends may be treated as eustatic SLR. Historical air temperature and sea surface temperature (SST) data of Karachi also show an increasing pattern and an increasing trend of about 0.67°C has been registered in the air temperature over the last 35 years, whereas the mean SST in the coastal waters of Karachi has also registered an increasing trend of about 0.3°C in a decade. Sindh coastal zone is more vulnerable to sea level rise than Baluchistan coast, as uplifting of the coast by about 1-2 mm/year due to subduction of Indian Ocean plate is a characteristic of Baluchistan coast. Within the Indus deltaic creek system, the area nearby Karachi is more vulnerable to coastal erosion and accretion than the other deltaic region, mainly due to human activities together with natural phenomena such as wave action, strong tidal currents, and rise in sea level. Therefore, The present article deals mainly with the study of dynamical processes such as erosion and accretion associated with sea level variations along the Karachi coast and surrounding Indus deltaic coastline. The probable beach erosion in a decade along the sandy beaches of Karachi has been estimated. The estimates show that 1.1 mm/year rise in sea level causes a horizontal beach loss of 110 mm per year. Therefore, coast eroded with rise in sea level at Karachi and surrounding sandy beaches would be 1.1 m during a period of next 10 years. The northwestern part of Indus delta, especially the Gizri and Phitti creeks and surrounding islands, are most unstable. Historical satellite images are used to analyze the complex pattern of sediment movements, the change in shape of coastline, and associated erosion and accretion patterns in Bundal and Buddo Islands. The significant changes in land erosion and accretion areas at Bundal and Buddo Islands are evident and appear prominently in the images. A very high rate of accretion of sediments in the northwestern part of Buddo Island has been noticed. In the southwest monsoon season the wave breaking direction in both these islands is such that the movement of littoral drift is towards west. Erosion is also taking place in the northeastern and southern part of Bundal Island. The erosion in the south is probably due to strong wave activities and in the northeast is due to strong tidal currents and seawater intrusion. Accretion takes place at the northwest and western parts of Bundal Island. By using the slope of Indus delta, sea encroachment and the land area inundation with rising sea level of 1 m and 2 m have also been estimated.  相似文献   

Investigation of the groundwater modelling component of the Integrated Water Flow Model (IWFM)     

Ali Ercan  Emin C. Dogrul  Tariq N. Kadir 《水文科学杂志》2013,58(16):2834-2848

ABSTRACT

The Integrated Water Flow Model (IWFM), developed by the California Department of Water Resources, is an integrated hydrological model that simulates key flow processes including groundwater flows, streamflow, stream–aquifer interactions, rainfall–runoff and infiltration. It also simulates the agricultural water demand as a function of soil, crop and climatic characteristics, as well as irrigation practices, and allows the user to meet these demands through pumping and stream diversions. This study investigates the modelling performance of the groundwater module of IWFM using several hypothetical test problems that cover a wide range of settings and boundary conditions, by comparing the simulation results with analytical solutions, field and laboratory observations, or with results from MODFLOW outputs. The comparisons demonstrate that IWFM is capable of simulating various hydrological processes reliably.