A spatio-temporal three-dimensional conceptualization and simulation of Dera Ismail Khan alluvial aquifer in visual MODFLOW: a case study from Pakistan

Dera Ismail Khan (DIK) is situated in the Lower Indus Basin of Pakistan. The land use has been changed in the canal command area due to irrigation activities near the Indus River. To check the current status and predict the groundwater levels in the area, the unconfined aquifer has been simulated in Visual MODFLOW for a period of 35 years, i.e., from 1985 to 2020. The 2900-km² area has been modeled with a grid of 500 by 500 m and the depth set to 100 m. The aquifer in the study area has been divided vertically and laterally into three and ten zones, respectively, for the characterization. Water wells and streams were used as the sinks and hydrologic boundaries, respectively. The model was successfully calibrated in steady and the non-steady state. The simulation revealed that the whole simulation can be divided into two phases, i.e., before and after the construction of the Chashma Right Bank Canal (CRBC), whereas the results were summarized in the form of water table depth maps and groundwater budget calculations. To determine the groundwater sustainability, a conjunctive use scenario has been employed to simulate the aquifer dynamics till 2020. The simulation revealed incremental drawdowns till the end.     

Subduction-related mafic to felsic magmatism in the Malayer–Boroujerd plutonic complex,western Iran

The Malayer–Boroujerd plutonic complex (MBPC) in western Iran, consists of a portion of a magmatic arc built by the northeast verging subduction of the Neo-Tethys plate beneath the Central Iranian Microcontinent (CIMC). Middle Jurassic-aged felsic magmatic activity in MBPC is manifested by I-type and S-type granites. The mafic rocks include gabbroic intrusions and dykes and intermediate rocks are dioritic dykes and minor intrusions, as well as mafic microgranular enclaves (MMEs). MBPC Jurassic-aged rocks exhibit arc-like geochemical signatures, as they are LILE- and LREE-enriched and HFSE- and HREE-depleted and display negative Nb–Ta anomalies. The gabbro dykes and intrusions originated from metasomatically enriched garnet-spinel lherzolite [Degree of melting (f_mel) ~ 15%] and exhibit negative Nd and positive to slightly negative ε_Hf(T) (+ 3.0 to ? 1.6). The data reveal that evolution of Middle Jurassic magmatism occurred in two stages: (1) deep mantle-crust interplay zone and (2) the shallow level upper crustal magma chamber. The geochemical and isotopic data, as well as trace element modeling, indicate the parent magma for the MBPC S-type granites are products of upper crustal greywacke (f_mel: 0.2), while I-type granites formed by partial melting of amphibolitic lower crust (f_mel: 0.25) and mixing with upper crustal greywacke melt in a shallow level magma chamber [Degree of mixing (f_mix): 0.3]. Mixing between andesitic melt leaving behind a refractory dense cumulates during partial crystallization of mantle-derived magma and lower crustal partial melt most likely produced MMEs (f_mix: 0.2). However, enriched and moderately variable ε_Nd(T) (? 3.21 to ? 4.33) and high (⁸⁷Sr/⁸⁶Sr)_i (0.7085–0.7092) in dioritic intrusions indicate that these magmas are likely experienced assimilation of upper crustal materials. The interpretations of magmatic activity in the MBPC is consistent with the role considered for mantle-derived magma as heat and mass supplier for initiation and evolution of magmatism in continental arc setting, elsewhere.     

Wetting Induced Deformation of Soils Triggering Landslides in Pakistan

Geotechnical and Geological Engineering - Shallow and deep seated landslides in natural slopes are often induced by rainfall. The cause of the failure is usually considered to be due to the...     

Cathodoluminescence of high-pressure feldspar minerals as a shock barometer

Cathodoluminescence (CL) analyses were carried out on maskelynite and lingunite in L6 chondrites of Tenham and Yamato-790729. Under CL microscopy, bright blue emission was observed in Na-lingunite in the shock veins. Dull blue-emitting maskelynite is adjacent to the shock veins, and aqua blue luminescent plagioclase lies farther away. CL spectroscopy of the Na-lingunite showed emission bands centered at ~330, 360–380, and ~590 nm. CL spectra of maskelynite consisted of emission bands at ~330 and ~380 nm. Only an emission band at 420 nm was recognized in crystalline plagioclase. Deconvolution of CL spectra from maskelynite successfully separated the UV–blue emission bands into Gaussian components at 3.88, 3.26, and 2.95 eV. For comparison, we prepared K-lingunite and experimentally shock-recovered feldspars at the known shock pressures of 11.1–41.2 GPa to measure CL spectra. Synthetic K-lingunite has similar UV–blue and characteristic yellow bands at ~550, ~660, ~720, ~750, and ~770 nm. The UV–blue emissions of shock-recovered feldspars and the diaplectic feldspar glasses show a good correlation between intensity and shock pressure after deconvolution. They may be assigned to pressure-induced defects in Si and Al octahedra and tetrahedra. The components at 3.88 and 3.26 eV were detectable in the lingunite, both of which may be caused by the defects in Si and Al octahedra, the same as maskelynite. CL of maskelynite and lingunite may be applicable to estimate shock pressure for feldspar-bearing meteorites, impactites, and samples returned by spacecraft mission, although we need to develop more as a reliable shock barometer.     

Implications for behavior of volatile elements during impacts—Zinc and copper systematics in sediments from the Ries impact structure and central European tektites

Moldavites are tektites genetically related to the Ries impact structure, located in Central Europe, but the source materials and the processes related to the chemical fractionation of moldavites are not fully constrained. To further understand moldavite genesis, the Cu and Zn abundances and isotope compositions were measured in a suite of tektites from four different substrewn fields (South Bohemia, Moravia, Cheb Basin, Lusatia) and chemically diverse sediments from the surroundings of the Ries impact structure. Moldavites are slightly depleted in Zn (~10–20%) and distinctly depleted in Cu (>90%) relative to supposed sedimentary precursors. Moreover, the moldavites show a wide range in δ⁶⁶Zn values between 1.7 and 3.7‰ (relative to JMC 3‐0749 Lyon) and δ⁶⁵Cu values between 1.6 and 12.5‰ (relative to NIST SRM 976) and are thus enriched in heavy isotopes relative to their possible parent sedimentary sources (δ⁶⁶Zn = ?0.07 to +0.64‰; δ⁶⁵Cu = ?0.4 to +0.7‰). In particular, the Cheb Basin moldavites show some of the highest δ⁶⁵Cu values (up to 12.5‰) ever observed in natural samples. The relative magnitude of isotope fractionation for Cu and Zn seen here is opposite to oxygen‐poor environments such as the Moon where Zn is significantly more isotopically fractionated than Cu. One possibility is that monovalent Cu diffuses faster than divalent Zn in the reduced melt and diffusion will not affect the extent of Zn isotope fractionation. These observations imply that the capability of forming a redox environment may aid in volatilizing some elements, accompanied by isotope fractionation, during the impact process. The greater extent of elemental depletion, coupled with isotope fractionation of more refractory Cu relative to Zn, may also hinge on the presence of carbonyl species of transition metals and electromagnetic charge, which could exist in the impact‐induced high‐velocity jet of vapor and melts.     

Pressure and temperature dependence of cation distribution in Mg-Mn olivine

Synthetic (Mg_0.51, Mn_0.49)₂SiO₄ olivine samples are heat-treated at three different pressures; 0, 8 and 12 GPa, all at the same temperature (～500° C). X-ray structure analyses on these single crystals are made in order to see the pressure effect on cation distribution. The intersite distribution coefficient of Mg and Mn in M1 and M2 sites, K _D = (Mn/Mg) _M1/(Mn/Mg) _M2, of these samples are 0.192 (0 GPa), 0.246 (8 GPa) and 0.281 (12 GPa), indicating cationic disordering with pressure. The small differences of cell dimensions between these samples are determined by powder X-ray diffraction. Cell dimensions b and c decrease, whereas a increases with pressure of equilibration. Cell volume decreases with pressure as a result of a large contraction of the b cell dimension. The effect of pressure on the free energy of the cation exchange reaction is evaluated by the observed relation between the cell volume and the site occupancy numbers. The magnitude of the pressure effect on cation distribution is only a fifth of that predicted from the observed change in volume combined with thermodynamic theory. This phenomenon is attributed to nonideality in this solid solution, and nonideal parameters are required to describe cation distribution determined in the present and previous experiments. We use a five-parameter equation to specify the cationic equilibrium on the basic of thermodynamic theory. It includes one energy parameter of ideal mixing, two parameters for nonideal effects, one volume parameter, and one thermal parameter originated from the lattice vibrational energy. The present data combined with some of the existing data are used to determine the five parameters, and the cation distribution in Mg-Mn olivine is described as a function of temperature, pressure, and composition. The basic framework of describing the cationic behavior in olivine-type mineral is worked out, although the result is preliminary: each of the determined parameters is not accurate enough to enable us to make a reliable prediction.     

High time resolution monitoring of tropospheric temperature with a Radio Acoustic Sounding System (RASS)

We have observed the time-height variation of the temperature field in the upper troposphere using a Radio Acoustic Sounding System (RASS) which consists of the MU radar and a high-power acoustic transmitter. The fast beam steerability of the MU radar has made it possible to measure temperature profiles in a fairly wide height range in the upper troposphere (5–11 km), even under intense wind conditions. Observations were continued for about 32 hr on 24–26 December, 1986 with a time-height resolution of 30 min and 150 m. During the observation period, the tropospheric jet was so intense that the acoustic wavefronts were severely distorted. Using wind velocity profiles observed by the MU radar we have numerically estimated the propagation of acoustic wavefronts, and further determined favorable pointing directions for the MU radar to receive significant backscattering from refractive index fluctuations produced by the acoustic waves. Conventional radiosonde soundings were carried out every 6 hr, which showed a temperature decrease of 4 K/day in the upper troposphere during the observation period. Temperature profiles taken by RASS agree well with the radiosonde results.     

Melting relationships in the system Fe-Feo at high pressures: Implications for the composition and formation of the earth's core

A reconnaissance investigation has been carried out on melting relationships in the system Fe-FeO at pressures up to 25 GPa and temperatures up to 2200° C using an MA-8 apparatus. Limited studies were also made of the Co-CoO and Ni-NiO systems. In the system FeFeO, the rapid exsolution of FeO from liquids during quenching causes some difficulties in interpretation of textures and phase relationships. The Co-CoO and Ni-NiO systems are more tractable experimentally and provide useful analogues to the Fe-FeO system. It was found that the broad field of liquid immiscibility present at ambient pressure in the Co-CoO system had disappeared at 18 GPa, 2200° C and that the system displayed complete miscibility between molten Co and CoO, analogous to the behaviour of the Ni-NiO system at ambient pressure. The phase diagram of the system Fe-FeO at 16 GPa and from 1600–2200° C was constructed from interpretations based on the textures of quenched run products. The solubility of FeO in molten iron is considerably enhanced by high pressures. At 16 GPa, the Fe-FeO eutectic contains about 10–15 mol percent FeO and the eutectic temperature in this iron-rich region of the system occurs at 1700±25° C, some 350° C below the melting point of pure iron at the same pressure. The solubility of FeO in molten Fe increases rapidly as temperature increases from 1700 to 2200° C. A relatively small liquid immiscibility field is present above 1900° C but is believed to be eliminated above 2200° C. This inference is supported by thermodynamic calculations on the positions of key phase boundaries. A single run carried out on an Fe₅₀ FeO₅₀ composition at 25 GPa and 2200° C demonstrated extensive and probably complete miscibility between Fe and FeO liquids under these conditions. The melting point of iron is decreased considerably by solution of FeO at high pressures; moreover, the melting point gradient (dP/dT) of the Fe-FeO eutectic is much smaller than that of pure iron and is also smaller than that of mantle pyrolite under the P, T conditions studied. These characteristics make it possible for melting of metal phase and segregation of the core to proceed within the Earth under conditions where most of the mantle remains below solidus temperatures. Under these conditions, the core would inevitably contain a large proportion of dissolved FeO. It is concluded therefore, that oxygen is likely to be the principal light element in the core. The inner core may not be composed of pure iron, as often proposed. Instead, it may consist of a crystalline oxide solid solution (Ni, Fe)₂O.     

Geomorphological study on a clastic accretionary prism: The Nankai Trough

Abstract The Nankai Trough, off southwest Japan, is one of the best sites for the study of geomorphic characteristics of a clastic accretionary prism. A recent multibeam survey over the central and eastern parts of the Nankai accretionary prism has revealed a large variation of the topography along the trough axis. Analysis of the bathymetric data suggests the existence of prism deformational features of different scales, such as depressions, embayment structures and cusps. These structures are the results of slope instability caused by basement relief of subducted oceanic plate. Unstable slopes recover by new accretion and development of a low angle thrust. Small-scale deformation due to the subduction of a small isolated seamount is then adjusted to the regional trend. By contrast, a 30 km indentation of the wedge observed in the eastern part of the Nankai Trough, the Tenryu Cusp, has seemed to retain its geometry. The subducted Philippine Sea plate has deformed greatly near the eastern end of the Nankai Trough, because of the collision between the Izu-Ogasawara (Bonin) arc and central Japan. Therefore, the indentation may be the result of the continuous subduction of a basement high, such as the Zenisu Ridge, which has been formed under north-south compression due to the arc-arc collision.     

High-angle reverse faulting associated with the 1946 Nankaido earthquake

Analysis of vertical crustal deformation data in the southwestern part of Shikoku, southwest Japan, suggests that the Nankaido earthquake of 1946 (M_w = 8.1), which is a principal interplate thrust earthquake, was accompanied by subsidiary faulting on a splay fault adjacent to the coast of Shikoku. Discarding crustal movement resulting from the main thrusting of the Nankaido earthquake, local leveling data are explained by slip on a simple rectangular thrust fault located just offshore of Shikoku. Although it is difficult to constrain the fault location, a possible result is a high-angle thrust dipping landward at an angle of about 70°, with a dislocation of about 1.5 m, and source dimensions of 30 × 13 km along strike and dip. respectively. This result indicates that the fault may be one of the steeply dipping subsidiary faults branching from the main low-angle thrust, as was the case in the Alaska earthquake of 1964. Although several lines of evidence suggest that this faulting occurred as slow aseismic slip, its discrimination from the main seismic event is extremely difficult. This kind of high-angle thrusting just offshore of the coast would play an important role for the formation of the marine terraces during the late Quaternary period.