Geochronological Constraints on the Haftcheshmeh Porphyry Cu-Mo-Au Ore Deposit, Central Qaradagh Batholith, Arasbaran Metallogenic Belt, Northwest Iran

The Haftcheshmeh porphyry Cu-Mo-Au deposit in the Arasbaran metallogenic belt(AMB) of NW Iran contains more than 185 Mt of ore, with a grade ranging from 0.3% to 0.4%. It is hosted within a porphyritic diorite to granodiorite intruded into an older gabbro-diorite intrusion. 40 Ar/39 Ar analyses of primary magmatic hornblende from the granodiorite porphyry and gabbro-diorite show plateau ages of 26.41 ± 0.59 Ma, with an inverse isochron age of 25.9 ± 1.0 Ma and a plateau age of 27.47 ± 0.17 Ma, with an inverse isochron age of 27.48 ± 0.35 Ma for these two rock types, respectively. Comparing these new age data with those from the nearby Sungun(20.69 ± 0.35 Ma) and Kighal porphyry deposits defines a northwest-southeast Cu-Mo-Au mineralization zone extending for 20 km over the time span of ~27 to 20 Ma. Geochemically, Haftcheshmeh rocks are calc-alkaline with high potassium affinities with tectonic setting in relation to volcanic arc setting. Large ion lithophile elements(LILE) such as Th, U and K show enrichment on a primitive mantle normalized diagram(specially Pb), and are depleted in high field strength elements(HFSE) such as Ti and Nb, pointing to a mantle magma source contamination with crustal materials by subducted oceanic crust.     

On inversion of the second- and third-order gravitational tensors by Stokes’ integral formula for a regional gravity recovery

A regional recovery of the Earth’s gravity field from satellite observables has become particularly important in various geoscience studies in order to better localize stochastic properties of observed data, while allowing the inversion of a large amount of data, collected with a high spatial resolution only over the area of interest. One way of doing this is to use observables, which have a more localized support. As acquired in recent studies related to a regional inversion of the Gravity field and steady-state Ocean Circulation Explorer (GOCE) data, the satellite gravity-gradient observables have a more localized support than the gravity observations. Following this principle, we compare here the performance of the second- and third-order derivatives of the gravitational potential in context of a regional gravity modeling, namely estimating the gravity anomalies. A functional relation between these two types of observables and the gravity anomalies is formulated by means of the extended Stokes’ integral formula (or more explicitly its second- and third-order derivatives) while the inverse solution is carried out by applying a least-squares technique and the ill-posed inverse problem is stabilized by applying Tikhonov’s regularization. Our results reveal that the third-order radial derivatives of the gravitational potential are the most suitable among investigated input data types for a regional gravity recovery, because these observables preserve more information on a higher-frequency part of the gravitational spectrum compared to the vertical gravitational gradients. We also demonstrate that the higher-order horizontal derivatives of the gravitational potential do not necessary improve the results. We explain this by the fact that most of the gravity signal is comprised in its radial component, while the horizontal components are considerably less sensitive to spatial variations of the gravity field.     

Water use at Luxor,Egypt: consumption analysis and future demand forecasting

Water resources in Egypt are becoming scarce and the demand for clean drinking water supply is one of the most important priorities of the Egyptian government in recent years. Analyzing water use and future demand forecasting is a primitive clue for water demand management. Water in Luxor is used for agricultural, residential, institutional, commercial, and touristic purposes. The results of water use analysis indicated that for the time period from 1983 to 2012, agriculture is the highest consumer of water which reached about 94.76–97.38 % followed by residential water consumption (1.90–3.05 %), institutional water consumption (0.71–1.75 %), and touristic water consumption (0.02–0.43 %), respectively. The future demand forecasting results revealed that the present situation may continue to rise in the next 50 years which will increase the water demand with a water deficit ranging between 15 and 114 MCM/year. To fill the gap between the present water consumption and future residential, institutional, commercial, and touristic water demand, additional municipal facilities, and improvement and management of water supply/demand are needed. To cope with the predicted future water demand, it is recommended to improve the on-farm irrigation, reduce the demand for irrigation water, rationalize the irrigation water use, and enhance the integrated role of water users in integrated water resources management.     

Introducing an effective coherence function to generate non-uniform ground motion on topographic site using time-domain boundary element method

In this study,a comprehensive parametric analysis was performed on non-uniform excitation of V-shaped topography using the boundary element method in time domain.For this purpose,wave scattering analysis was carried out on a topography subjected to the SV-wave for different predominant frequencies and shape ratios.Based on the numerical results,new coherence and time delay functions are proposed to generate non-uniform ground motion for topographic irregularities.The efficiency and accuracy of the proposed functions for real engineering problems are indicated by comparison with observations reported in previous literature.     

Multicomponent theory of buoyancy instabilities in magnetized plasmas: the case of magnetic field parallel to gravity

We investigate electromagnetic buoyancy instabilities of the electron-ion plasma with the heat flux based on not the magnetohydrodynamic (MHD) equations, but using the multicomponent plasma approach when the momentum equations are solved for each species. We consider a geometry in which the background magnetic field, gravity, and stratification are directed along one axis. The nonzero background electron thermal flux is taken into account. Collisions between electrons and ions are included in the momentum equations. No simplifications usual for the one-fluid MHD-approach in studying these instabilities are used. We derive a simple dispersion relation, which shows that the thermal flux perturbation generally stabilizes an instability for the geometry under consideration. This result contradicts to conclusion obtained in the MHD-approach. We show that the reason of this contradiction is the simplified assumptions used in the MHD analysis of buoyancy instabilities and the role of the longitudinal electric field perturbation which is not captured by the ideal MHD equations. Our dispersion relation also shows that the medium with the electron thermal flux can be unstable, if the temperature gradients of ions and electrons have the opposite signs. The results obtained can be applied to the weakly collisional magnetized plasma objects in laboratory and astrophysics.     

The influence of winds on the time-dependent behaviour of self-gravitating accretion discs

We study the effects of winds on the time evolution of isothermal, self-gravitating accretion discs by adopting a radius-dependent mass-loss rate because of the existence of the wind. Our similarity and semi-analytical solution describes time evolution of the system in the slow accretion limit. The disc structure is distinct in the inner and outer parts, irrespective of the existence of the wind. We show that the existence of wind will lead to a reduction of the surface density in the inner and outer parts of the disc in comparison to a no-wind solution. Also, the radial velocity significantly increases in the outer part of the disc, however, the accretion rate decreases due to the reduced surface density in comparison to the no-wind solution. In the inner part of the disc, mass loss due to the wind is negligible according to our solution. But the radial size of this no-wind inner region becomes smaller for stronger winds.     

Gravitational instability of discs with dissipative coronae around supermassive black holes

We study the dynamical structure of a self-gravitating disc with coronae around a supermassive black hole. Assuming that the magnetorotational instability responsible for generating the turbulent stresses inside the disc is also the source for a magnetically dominated corona, a fraction of the power released when the disc matter accretes is transported to and dissipated in the corona. This has a major effect on the structure of the disc and its gravitational (in)stability according to our analytical and self-consistent solutions. We determine the radius where the disc crosses the inner radius of gravitational instability and forms the first stars. Not only the location of this radius which may extend to very large distances from the central black hole, but also the mass of the first stars highly depends on the input parameters, notably the viscosity coefficient, the mass of the central object and the accretion rate. For accretion discs around quasi-stellar objects (QSOs) and the Galactic Centre, we determine the self-gravitating radius and the mass of the first clumps. Comparing the cases with a corona and without a corona for typical discs around QSOs or the Galactic Centre, when the viscosity coefficient is around 0.3, we show that the self-gravitating radius decreases by a factor of approximately 2, but the mass of the fragments increases with more or less the same factor. The existence of a corona implies a more gravitationally unstable disc according to our results. The effect of a corona on the instability of the disc is more effective when the viscosity coefficient increases.     

Formation of fluctuations in a molecular slab via isobaric thermal instability

Frictional heating by the ion-neutral drift is calculated and its effect on the isobaric thermal instability is studied. Ambipolar drift heating of a one-dimensional self-gravitating magnetized molecular slab is used under the assumptions of quasi-magnetohydrostatic and local ionization equilibrium. We see that ambipolar drift heating is inversely proportional to density and its value in some regions of the slab can be significantly larger than the average heating rates of cosmic rays and turbulent motions. The results show that isobaric thermal instability can occur in some regions of the slab, and thus it may produce slab fragmentation and formation of astronomical unit scale condensations.     

Titanium‐in‐biotite thermometry in porphyry copper systems: Challenges to application of the thermometer

Empirical geothermometer dealing with Ti solubility in the Fe‐Mg biotites was originally proposed for biotites in graphitic, peraluminous metapelites containing ilmenite or rutile that equilibrated roughly at 4–6 kbar. Given that biotites are abundant in the porphyry copper systems, this geothermometer has frequently been used for the determination of magmatic–hydrothermal temperatures in the porphyry copper systems. Common associations of porphyry copper deposits (PCDs), that is, low Al content of biotite, biotite chloritization (causes the biotite to become more magnesian and to lose Ti), and biotite formation by amphibole replacement, as well as disequilibrium, local equilibrium, or re‐equilibration of biotites, especially through potassic alteration, may provide significant uncertainty in the temperatures estimated a by Ti‐in‐biotite geothermometer. In addition, besides the calibration range of thermometer for pressure (400–600 MPa), the temperatures of major sulfide precipitation in PCDs (>~400°C) does not fit with the temperature range of thermometer calibration (480–800°C). Worth noting, as confirmed by fluid inclusion data in the Sarkuh PCD, regardless of presence of mineralogical requirements, obtained temperatures of sulfide mineralization using Ti in biotite thermometer could be overestimated. This may be due to the difference between general conditions of sulfide mineralization and calibration range of Ti in the biotite thermometer for pressure and temperature, as well as the metaluminous nature of biotites in PCDs.     

Seismic analysis of dam-foundation-reservoir system including the effects of foundation mass and radiation damping

One of the main concerns in using commercial software for finite element analyses of dam-foundation-reservoir systems is that the simplifying assumptions of the massless foundation are unreliable. In this study, an appropriate direct finite element method is introduced for simulating the mass, radiation damping and wave propagation effect in foundations of damfoundation-reservoir systems using commercial software ABAQUS. The free-field boundary condition is used for modeling the semi-infinite foundation and radiation damping, which is not a built-in boundary condition in most of the available commercial software for finite element analysis of structures such as ANSYS or ABAQUS and thus needs to be implemented differently. The different mechanism for modeling of the foundation, earthquake input and far-field boundary condition is described. Implementation of the free-field boundary condition in finite element software is verified by comparing it with analytical results. To investigation the feasibility of the proposed method in dam-foundation-reservoir system analysis, a series of analyses is accomplished in a variety of cases and the obtained results are compared with the substructure method by using the EAGD-84 program. Finally, the massed and massless foundation results are compared and it is concluded that the massless foundation approach leads to the overestimation of the displacements and stresses within the dam body.