Integration of Landsat-8 and Sentinel-1 dataset to extract geological lineaments in complex formations of Tepal mountain area,Shahrood, north Iran

Recently, the detection and extraction of geological lineaments have become an essential analytical technique to find relationships between the characteristics and occurrence of hydrogeology, and tectonic studies. The use of remote sensing, with the progressive development of image enhancement techniques, provides an opportunity to produce more reliable and comprehensive lineament maps. In this paper, semi-automatic approach based on Landsat 8 and Sentinel 1 radar data is proposed for lineaments extraction and validation. The combined method of linear filtering and automatic line module ensures a high degree of accuracy resulting in a lineament map. Based on identified lineaments, Sentinel1 is more capable of detecting edges than Landsat8, but the primary orientation lineaments extracted from Landsat8 and Sentinel1 were different. So, by combining band6 of Landsat8, and VV and VH polarization of Sentinel1, the area lineaments were extracted with high accuracy. Rose diagram showed the extracted lineaments' orientation is in good compliance with the region's existing faults. Also, the formations' lineament length density has good consistent with the density of the faults in the geological map.     

Evaluation of MSS Landsat imagery over central Spain

Twenty four MSS Landsat I and II images have been selected to determine the most persistent lineaments within a given area in Central Spain. This area comprises four different geotectonic units: the Hercynian basement, an intermediate-type Alpine chain and two undeformed Tertiary basins. Lineaments appear to be well established in the Iberian Peninsula and their general pattern in good accordance with the faults system. Some of them correspond to simple or composite linears that extend over the entire peninsula, and must affect the whole crust. A generalized late-Hercynian fracture event is admitted to be the origin of these linears and their extension into the post-Hercynian areas is related to an Alpine reactivation. The main lineament directions are NS, N20, N70, N120 and N160. These trends seem to correlate well with arches and troughs inferred from gravimetric and aeromagnetic surveys carried out over both Tertiary basins.BPCG, Sn, W mineralizations within the Hercynian basement occur in relation to quartz and baryte dykes trending N20, N70, and N120.     

Automatic extraction of lineaments based on wavelet edge detection and aided tracking by hillshade

Lineaments refer to the linear or curvilinear textures on remote sensing image, whose general spatial distribution characteristics are often the response of deep geological sturcture at the surface. Firstly, we use wavelet modulus maximum transformation to detect the edges with 4 scales on Landsat – 8 OLI B5 image and analyze their multi-scale characteristics. As the result, it is determined that the optimal scale of edge detection is 4, and the outline that consist of the edge pixels is roughly corresponding to the geological structure of mine area. Thus the incomplete lineaments have been extracted by using the 2D otsu algorithm. Secondly, the hillshade map generated based on DEM is processed to generate binarized linear shadow. Finally, the linear shadow is superimposed on the lineaments preliminarily extracted to obtain the optimized lineaments. Experiment results show that, based on the method, there are some deformation and displacement between the lineaments extracted and the actual geological structure, and it fail to effectively extract Qilinchang Fault, but lineaments are in good correspondence with Kuangshanchang Fault, Dongtou Fault and Niulan River Fault, which are basically in accord with the geological structure framework of the mine area.     

Estimation of stress and its use in evaluation of landslide prone regions using remote sensing data

Recently, use of remote sensing data for determining the orientation of stress has been demonstrated. The present study deals with the estimation of stress pattern in the part of the Himalayan region which shows the ongoing neo-tectonic activities. The study area falls into a tectonically active zone of the Central-Himalaya, with a complex geotectonic set-up confined by a number of faults. Efforts have been made to evaluate the technique as a fast algorithm for quick and time limited analysis of linear feature from which the orientation of the lineaments are estimated by using remote sensing data. Further, the estimation of stress and the lineament analysis have been used in mapping of landslide prone areas. Terrain information such as land cover, geology, lineament, faults, mega faults, geomorphology and drainage has been derived from the satellite imageries, and the existing thematic information has been updated to enable the quantification of landslide causative parameters. Spatial and temporal multi-layered information have been used for landslides hazard susceptibility analysis. The qualitative hazard analysis has been carried out using the map overlying techniques in GIS environment along the central part of Himalayan region. It has been observed that the high potential zones have been found to have very high lineament density, moderate to low drainage density and high slope areas of the terrain. On the basis of the geological and morphological analysis, it is further suggested that the combined impacts of the crushed nature of bed rock (due to the neo-tectonic activities), heavy rainfall and lack of vegetation cover cause persistent recurrence of landslides along this region. The role of earthquake on induction of landslides will be presented.     

Use of the ASTER satellite images for evaluation of structural changes in the Popocatépetl volcano related to microseismicity

Mexico is one of the most volcanically active regions in North America. Volcanic activity in central Mexico is associated with the subduction of the Cocos and Rivera plates beneath the North American plate. Periods of enhanced microseismic activity, associated with the volcanic activity of the Popocatépetl volcano are compared with periods, during which the microseismic activity was low. We detected systematical changes in the number of lineaments, associated with the microseismic activity due to lineament analysis of a temporal sequence of high resolution satellite images of the Popocatépetl volcano, provided by the ASTER/VNIR instrument. The Lineament Extraction and Stripes Statistic Analysis (LESSA) software package was used for the lineament extraction. In the future it would allow develop a methodology for detection of possible elevation of pressure in volcano edifice.     

Automatic lineament extraction in a heavily vegetated region using Landsat Enhanced Thematic Mapper (ETM+) imagery

Lineament extraction from satellite remotely sensed data has been one of the widely used applications of remote sensing in geology. In fact, recent advances in digital image processing allow such lineament extraction to be accomplished in semi-automatic to fully automatic approaches. However, satellite remotely sensed data acquired in heavily vegetated regions such as tropical rainforest, are vulnerable to higher inherent noise levels attributed to the resultant effects of scattering by clouds and adjacency effects of highly inhomogeneous vegetation cover within the pixel dimension. In this study, we examined the effects of noise levels to lineament extraction using a fully automatic approach, consisting of a combination of edge-line detection algorithms. Ancillary information from a digitized topographic map and image classification was used to discriminate between cultural and natural lineaments from the extracted lineaments. Adapting the combination of edge detection and a line-linking algorithm, we have found the optimal parameters for automatic lineament extraction of such complex areas using Enhanced Thematic Mapper (ETM+) data. A noise level of 30% is the maximum threshold before artifacts are generated. It is therefore concluded that the combination of edge-based and line-linking digital image processing operations with the priori local optimal parameters is crucial in lineament feature extraction in heavily vegetated regions.     

RS-GIS based morphometrical and geological multi-criteria approach to the landslide susceptibility mapping in Gish River Basin,West Bengal,India

Darjeeling Himalaya is one of the several mountainous areas of India which is often suffered from landslide hazards. In this paper, a multi criteria evaluation is applied using 16 morphometric indicators, geology and lineaments to identify the areas vulnerable in respect to drainage and relief conditions. As both drainage and relief parameters exert strong influences on landslide intensity, both the diversity maps are integrated for final landslide susceptibility mapping. The obtained results show that 20.17?sq.?km (7.61%) area within the basin is highly susceptible for landslides, where average drainage density is 3.78?km/sq.?km, relative relief is greater than 408?m and slope is greater than 12°. The validation result shows that very high landslide susceptible zone is associated with very high frequency of landslide occurrence. Beside this, ROC curve also suggests good predicted rate (86.60%) for the model. So, the proposed method can be applied for predicting landslide susceptible zone.     

Lineament mapping and fractal analysis using SPOT-ASTER satellite imagery for evaluating the severity of slope weathering process

The present study describes a remote sensing approach for preparing lineament map that subsequently indicates the influence of lineament density in the severity of weathering development. In this study, SPOT-5 data, the integration of SPOT-ASTER and Digital Elevation Model (DEM) data were used and processed. The existence of an active fault system in the south of Mashhad city, NE Iran and presence of schistose rocks in this area result in the development of numerous lineament features. This region was selected for this research. Lineament features including fractures, bedding plane, cleavage, shear zones and schistosity were mapped in the study area. The results indicate that the highest concentration of lineaments occurred in the central-western and south-eastern parts of the study area, which coincide with metamorphic outcrops and NW-SE trending fault system. A comparison of lineament statistical analysis and field survey demonstrated that the structural discontinuities have a significant effect on forming and distribution of weathering profiles. It was observed that increasing the number, length and density of structural discontinuities led to strong severity in weathering, which can produce deep residual soils susceptible to landslide occurrence. The remote sensing approach developed in this study can be applicable for preparing lineament maps and evaluating the severity of weathering development in other active fault zones around the world.     

Venus atmospheric circulation: Observations and implications of the thermal structure

Continued analysis of Pioneer Venus imaging and polarimetry data indicates that the average cloud-top level circulation is mainly zonal (east to west) with a small meridional component. Presence of planetary scale waves and a possible sun-related component are evident in the data. If the tracked features refer to the same vertical level, then some variability of the circulation would have to be present to account for the Pioneer and Mariner 10 cloud-tracking results. However, the implied balanced flow from the observed thermal structure analysis strongly suggests that at least some of the variations in these observations is due to apparent cloud-top variations and that the circulation itself is relatively stable.Direct cyclostrophic calculations based on the observed thermal structure of the atmosphere yield a balanced zonal circulation with distinct mid-latitude jets (peak velocities about 110–120 ms^?1) located between 50 and 40 mb in each hemisphere of the planet near 45° latitude. The calculations which extend to about 40 km altitude from 80 km above the surface agree well with the observed entry probe zonal components and indicate breakdown of the balance condition near the upper and lower boundaries at low latitudes.The balanced flow results are consistent with the Mariner 10 and Pioneer cloud tracked estimates of the zonal circulation provided the effective altitude of the tracked features is slightly different at different observation periods. The features in the Pioneer Venus data would then lie on a sloping surface that extends from about 68 km (40 mb) at low latitudes to about 75 km (10 mb) in mid-latitudes. The polarization features would occur on a roughly parallel surface that is 1–2 km above the effective cloud-height surface, and Mariner 10 features would have effective altitudes somewhat lower than the Pioneer ultraviolet features. A slight asymmetry is evident in the balanced zonal circulation arising out of an asymmetry in the thermal field.Finally, the solenoids formed by intersecting isobaric and isosteric (constant specific volume) surfaces deduced from the Pioneer Venus radio occultation data show distinct evidence of a direct meridional circulation that may be important in sustaining the Venus atmospheric circulation.     

太阳耀斑和地磁扰动的统计分析(1966—1978)

本文给出了不同等级、不同持续时间、发生在日面不同日心经距上的太阳耀斑所引起的各种地磁指数(Ap、AE、Dst和Kp指数)的变化,讨论了各类耀斑所引起的地磁扰动的特点以及耀斑的日心经距及其喷发物在日地空间的传播速度对磁扰强度的影响。      

Neotectonic investigations in the Pannonian basin based on satellite images

Analysis of Landsat images of the Mura depression in the Pannonian basin enabled identification of a number of lineaments that, in most cases, were not known from previous geological investigations. Many of these lineaments act as real faults and they have neotectonic significance. Combination of Landsat and field data made possible construction of a structural model that can explain the neotectonic evolution of this area. The model is based on wrench-tectonics, where Sava, ?o?tanje and Labod faults represent the first order wrench faults, and Donat, Ormo? and Ljutomer the second order wrench faults. During this wrenching only a few accompanying structural forms were newly generated. In most cases they follow a pre-existing regional fracture pattern.     

Modeling the ionospheric gradient by using VSS base functions and GPS data over Iran

The spatial and temporal variations of ionosphere play an important role in positioning and navigation by the space geodetic techniques. Therefore, the ionospheric gradient should be calculated, analyzed, and applied in both space and time. Spatial gradients of the ionosphere have remarkable delay on the propagation of electromagnetic waves. This study intends to propose a new method for simultaneous modeling of the spatial gradients of ionosphere and VTECs in the local scale for Iran. Vector Spherical Slepian (VSS) base functions are used for the development of this method.Five VSS models with the maximal degrees of L = 30, 35, 40, 45 and 50 are taken into account. For implementing the VSS method, 24 permanent GPS stations from the Iranian Permanent GPS Network (IPGN) have been used. The unknown coefficients are estimated with the observations of these stations with least squares technique. Four other stations are used for evaluating the accuracy of the models. Repeatability of baselines is the measure that is used for this purpose. Based on the results obtained, L = 40 is the optimum degree for the VSS model with this input data over Iran.The baselines’ repeatability showed that ionospheric gradients have more influence on the north–south component. Moreover, the spatial gradient is negligible in the east–west and up-down component when a short baseline is processed. In other words, this kind of ionospheric modeling has significant application for baseline, which is longer than 1000 km. In the study, proposed method has improved the long baselines' solution by more than 12%, 18% and 10% in east–west, north–south and up-down components, respectively.     

Seasonal variability of atmospheric aerosol over the North Indian region during 2005–2009

The Indo-Gangetic basin (IGB) extends 2000 km in length along NW–SE and has 400 km width, in the north the basin is bounded by towering Himalaya. High aerosol optical depth (AOD) is observed over the IGB throughout the year. The Himalaya restricts the transport of aerosols across Tibet and China. We have used ground based Kanpur and Gandhi College Aerosol Robotic Network (AERONET) stations and Multiangle Imaging SpectroRadiometer (MISR) and Moderate Resolution Imaging Spectroradiometer (MODIS) Terra level-3 AOD products for the years 2005–2009 to study the variability of aerosol over the Indo-Gangetic (IG) plains. An increase in both satellite-derived as well as ground observed aerosol loading during 2005–2009 has been found over major cities located in the IG plains. The correlation coefficients between AERONET and MISR data are found to be 0.70, 0.36 0.82, in contrast the correlation coefficients between AERONET and MODIS 0.49, 0.68, and 0.43, respectively during summer, winter and monsoon seasons. The AOD estimation using MISR is found to be close to AERONET data during summer and monsoon seasons, in contrast MODIS estimation is better during winter season.     

Circulation of the atmosphere from the surface to 100 km

A reference model of the atmospheric circulation on Venus based on available observations is presented. The reference atmosphere has the following main features: (i) the entire atmosphere below 85 km moves predominantly from east to west in the planet's reference system (in the same direction as the rotation of the solid planet itself) with the possible exception of the lowest 10 km where velocities are low, (ii) a jet is present near the cloud-top level at 45° latitude in both hemispheres, with a magnitude of approximately 100 ms^?1, (iii) a weak meridional (north-south component) flow directed towards either pole is superimposed on the zonal (east-west) motion at cloud-top level (about 68 km) altitudes, and, (iv) eddies or wave motions are present in the atmosphere, with amplitudes of less than 15 ms^?1 in the upper atmosphere.     

Study of seasonal and long-term vertical deformation in Nepal based on GPS and GRACE observations

Lithospheric deformation signal can be detected by combining data from continuous global positioning system (CGPS) and satellite observations from the Gravity Recovery and Climate Experiment (GRACE). In this paper, we use 2.5- to 19-year-long time series from 35 CGPS stations to estimate vertical deformation rates in Nepal, which is located in the southern side of the Himalaya. GPS results were compared with GRACE observations. Principal component analysis was conducted to decompose the time series into three-dimensional principal components (PCs) and spatial eigenvectors. The top three high-order PCs were calculated to correct common mode errors. Both GPS and GRACE observations showed significant seasonal variations. The observed seasonal GPS vertical variations are in good agreement with those from the GRACE-derived results, particularly for changes in surface pressure, non-tidal oceanic mass loading, and hydrologic loading. The GPS-observed rates of vertical deformation obtained for the region suggest both tectonic impact and mass decrease. The rates of vertical crustal deformation were estimated by removing the GRACE-derived hydrological vertical rates from the GPS measurements. Most of the sites located in the southern part of the Main Himalayan Thrust subsided, whereas the northern part mostly showed an uplift. These results may contribute to the understanding of secular vertical crustal deformation in Nepal.     

Induced geomagnetic variation and crustal evolution of the southern peninsula of India

The current thinking on crustal evolution and the influence of Precambrian lineaments on tectonism and mineralization provide support for the conjecture that zones of high electrical conductivity exist in the crust which are related to the tectonic and mineralized zones. Studies by the present authors using geological and tectonic data from remote sensing techniques, in particular, together with relevant data on anomalous geomagnetic variations reported in this paper would serve as additional contribution towards the global model.     

Statistical characteristics of low latitude Pc4s: Influence of solar wind and IMF parameters

Pc4 signatures for the year 2013, extracted from geomagnetic north–south and east–west components of induction coil magnetometer (LEMI 30) from low latitude station Desalpar (DSP), operated by Institute of Seismological Research (ISR), India have been investigated vis-à-vis the prevalent interplanetary parameters (IMF) as well as the geomagnetic activity indices. A clear dominance of Pc4-5 (467 events) over Pc3 (17 events) is observed. Local time variation of Pc4 shows a peak in the noon sector in both X and Y components. Our investigations show that the dominant peak frequency is 10 mHz at low latitude region. Correlations with solar wind and IMF parameters illustrate highest occurrence of Pc4 for a solar wind speed of 300–400 km/s and average IMF B field of 3–6 nT. The amplitude of Pc4s at DSP shows an increase with increasing solar wind speed, plasma density, solar wind dynamic pressure and average B field which is also reflected in the trend of frequency variation of these pulsations. We report that IMF clock angle at low latitude does not have influence on Pc4 occurrence. Based on the characteristics of these events, detected in latitudinally distributed stations from low and mid-latitudes from northern and southern hemisphere, we infer that modes were compressional, which could be driven by K-H instability or solar wind dynamic pressure, as compressional modes can propagate to low latitude with little attenuation.     

Effects of variations of geomagnetic field on VLF waves induced by beating of two HF waves

Beat wave (BW) high frequency (HF) ionospheric heating experiments were conducted to generate very low frequency (VLF) waves. The VLF waves were registered with a VLF receiver located ～15?km east of the European Incoherent Scatter (EISCAT) heating facility in Tromsø, Norway. A fluxgate magnetometer was used to monitor auroral electrojet current, and ionospheric conditions were measured using a Dynasonde. Correlation coefficients between VLF amplitudes and the deviation of geomagnetic north–south components were calculated. Experimental results show that strong and positive correlation exists the majority of the time, but sometimes no correlation or even a negative correlation occurred. This is consistent with similar past experiments that took place with exclusively AM generation. These results therefore support the conclusion that BW generation of VLF waves is no different than with AM, likely occurring in the D or lower E ionospheric region.     

月球的构造格架及其演化差异

根据以GRAIL重力场数据和LOLA地形数据计算的月壳厚度,将月球构造格架初步划分为三个构造单元:主要覆盖月球正面风暴洋区域的月海构造单元、主要覆盖月球背面高地的月陆构造单元以及主要位于南半球背面的南极艾肯盆地构造单元。结合最新的研究成果,对各构造单元上的重大地质事件包括岩浆事件、火山事件和撞击事件分别进行了简单的梳理,结果表明月球不同构造单元的演化事件具有明显的差异。     

Pre-earthquake ionospheric anomalies before three major earthquakes by GPS-TEC and GIM-TEC data during 2015–2017

The dual-frequency satellite-based measurements from Global Positioning System (GPS) may provide feasible ways of studying and potentially detecting of earthquake (EQ) related anomalies in the ionosphere. In this paper, GPS based Total Electron Content (TEC) data are studied for three major M?>?7.0 EQs in Nepal and Iran-Iraq border during 2015–2017 by implementing statistical procedures on temporal and spatial scale. Previous studies presented different time interval of pre-seismic ionospheric anomalies, however, this study showed that EQs ionospheric precursors may occur within 10?days. Furthermore, the ionospheric anomalies on the suspected day occurred during UT?=?08:00–12:00?h before the main shock. The Global Ionospheric Map TEC (GIM-TEC) data retrieved over the epicenter of M7.8 (Nepal EQ) showed a significant increase of 6 TECU on April 24, 2015 (one day before the main shock), which is recorded by the ground GPS station data of Islamabad (station lies within the EQ preparation zone). Furthermore, the spatial GIM-TEC result imply significant anomalies over the epicenter during the time interval between UT?=?08:00–12:00?h (LT?=?13:00–17:00). For M7.3 (Nepal EQ), the TEC anomalies were detected on May 10, 2015 (2?days before the event) in the temporal data. The spatial TEC data imply the huge clouds over the epicenter at about UT?=?08:00–12:00?h on May 10, 2015, that may be associated with this EQ in the quiet geomagnetic storm conditions. Similarly, temporal and spatial TEC showed anomaly on November 3, 2017, during UT?=?08:00–12:00 (9?days before the Iran-Iraq border EQ) after implementing the statistical method on it. Conversely, there exists a short-term but low magnitude TEC anomaly synchronized with a geomagnetic storm on November 7–8, 2017 (4 to 5?days prior to M7.3 Iran-Iraq border EQ). The diurnal and hourly GIM-TEC and VTEC data also imply the execution of ionospheric anomalies within 10?days prior to all events. All these positive anomalies in TEC may be due to the existence of a huge energy from the epicenter during the EQ preparation period.