Search

The Urban Heat Island (UHI) phenomenon exacerbates thermal discomfort in urban areas and significantly contributes to urban overheating when combined with climate change. This study investigates the spatiotemporal patterns of Surface Urban Heat Island Intensity (SUHII) in six major cities of Pakistan, focusing on the interplay between urban expansion, vegetation cover, and SUHII. To quantify SUHII dynamics, the impact of urban sprawl and vegetation changes was analyzed. The study offers critical insights into the implications for urban planning and policymaking in Pakistan. Using remote sensing data from Landsat satellites, analyzed with Geographic Information Systems (GIS) techniques, estimates of SUHII, urban expansion, and vegetation cover were derived. Specifically, imagery from Landsat-5 (2010−2013) and Landsat-8 (2014–2022), obtained from the US Geological Survey (USGS), was employed. Statistical analyses, including Pearson's correlation and linear regression, were conducted to assess relationships between these variables from 2010 to 2022. SUHII was found to increase annually by 0.18 °C in Islamabad and 0.19 °C in Peshawar, with corresponding urban expansion rates of 8.07 km2 (8967.75 pixels) and 1.67 km2 (1860.42 pixels) per year, respectively. Vegetation indices such as the Normalized Difference Vegetation Index (NDVI) and Fractional Vegetation Cover (FVC) were inversely correlated with SUHII, explaining up to 50 % of the variance in Peshawar. However, weaker correlations in Lahore suggest the presence of additional factors influencing SUHII. A distinct spatial relationship between increased vegetation and cooler areas was observed. For instance, Islamabad has greater vegetation cover and cool zones over 41.5 km2. In contrast, Lahore's hot spots spanned 127.1 km2, compared to Abbottabad's 10.4 km2, underscoring the thermal impact of reduced vegetation. The findings emphasize the effectiveness of urban greening, particularly in Islamabad's neutral thermal regions, in mitigating SUHII. This study offers important insights for urban planners in developing sustainable, climate-resilient cities within similar urban contexts. While the results are specific to Pakistani cities, the role of vegetation in mitigating SUHII may hold broader relevance for urban planning strategies in comparable settings.

This article presents an optimization-based approach to overcome redundancy arising from the multivariables enumeration process in multiple signal classification (MUSIC). By incorporating Broyden–Fletcher–Goldfarb–Shanno (BFGS) optimization, the computational speed of the MUSIC algorithm is significantly improved while maintaining mathematical accuracy. The optimization techniques require reasonable initial values to start the iteration, while for single target imaging purposes, the initial values can be acquired by the boundary between the near field and the far field. To generate suitable initial values for the optimization, we employ a modified convolutional neural network (CNN) to approximate the boundaries between the near and far fields, which vary with array system properties. Besides, the proposed method introduces a method for the Hessian matrix and gradient initialization for the BFGS method. Using simulation results as training samples, the modified CNN successfully establishes boundary approximations. Simulation and experimentation confirm the feasibility of our proposed method, showing its advantages in both accuracy and computation speed compared to existing techniques.

This work investigates the ionospheric response to the March 2023 geomagnetic storm over American and Asian sectors from total electron content (TEC), rate of TEC index, ionospheric heights, Swarm plasma density, radio occultation profiles of Formosat-7/Cosmic-2 (F7/C2), Fabry-Perot interferometer driven neutral winds, and E region electric field. During the storm’s main phase, post-sunset equatorial plasma bubbles (EPBs) extend to higher latitudes in the western American longitudes, showing significant longitudinal differences in the American sector. Over the Indian longitudes, suppression of post-sunset irregularities is observed, attributed to the westward prompt penetration electric field (PPEF). At the early recovery phase, the presence of post-midnight/near-sunrise EPBs till post-sunrise hours in the American sector is associated with the disturbance of dynamo-electric fields (DDEF). Additionally, a strong consistency between F7/C2 derived amplitude scintillation (S4) ≥ 0.5 and EPB occurrences is observed. Furthermore, a strong eastward electric field induced an increase in daytime TEC beyond the equatorial ionization anomaly crest in the American region, which occurred during the storm’s main phase. Both the Asian and American sectors exhibit negative ionospheric storms and inhibition of ionospheric irregularities at the recovery phase, which is dominated by the disturbance dynamo effect due to equatorward neutral winds. A slight increase in TEC in the Asian sector during the recovery phase could be explained by the combined effect of DDEF and thermospheric composition change. Overall, storm-time ionospheric variations are controlled by the combined effects of PPEF and DDEF. This study may further contribute to understanding the ionospheric responses under the influence of storm-phase and LT-dependent electric fields.

Earth observations from remotely sensed data have a substantial impact on natural hazard surveillance, specifically for earthquakes. The rapid emergence of diverse earthquake precursors has led to the exploration of different methodologies and datasets from various satellites to understand and address the complex nature of earthquake precursors. This study presents a novel technique to detect the ionospheric and atmospheric precursors using machine learning (ML). We examine the multiple precursors of different spatiotemporal nature from satellites in the ionosphere and atmosphere related to the Turkey earthquake on 6 February 2023 (Mw 7.8), in the form of total electron content (TEC), land surface temperature (LST), sea surface temperature (SST), air pressure (AP), relative humidity (RH), outgoing longwave radiation (OLR), and air temperature (AT). As a confutation analysis, we also statistically observe datasets of atmospheric parameters for the years 2021 and 2022 in the same epicentral region and time period as the 2023 Turkey earthquake. Moreover, the aim of this study is to find a synchronized and co-located window of possible earthquake anomalies by providing more evidence with standard deviation (STDEV) and nonlinear autoregressive network with exogenous inputs (NARX) models. It is noteworthy that both the statistical and ML methods demonstrate abnormal fluctuations as precursors within 6 to 7 days before the impending earthquake over the epicenter. Furthermore, the geomagnetic anomalies in the ionosphere are detected on the ninth day after the earthquake (Kp > 4; Dst < −70 nT; ap > 50 nT). This study indicates the relevance of using multiple earthquake precursors in a synchronized window from ML methods to support the lithosphere–atmosphere–ionosphere coupling (LAIC) phenomenon.

The ionospheric storm time responses during August 2018 are investigated over South American region using multiple observables, for example, Global Navigation Satellite System (GNSS) derived vertical total electron content (VTEC) from International GNSS Service, magnetic field data, geomagnetic indices, global ionospheric maps, thermospheric mass density (TMD), and [O/N2] ratio measurement. Strong-ionospheric and upper-atmospheric disturbances affected the ionospheric variables with long duration during the storm recovery phase and following after. First, daytime VTEC (9:00–20:00 UT) presented variations of >15 TECU during days 25 to 30 of August 2018 in low and middle latitudes of South America, this after sudden storm commencement (SSC). Furthermore, nighttime (21:00–24:00 and 00:00–05:00 UT) VTEC presented low values (58 TECU) in the recovery phase. Second, the ionospheric values during the storm main phase and following after, at low- and mid-latitudes, caused the equatorial ionization anomaly to expand due to prompt penetration electric field. Furthermore, VTEC enhancements are likely to occur few hours after the SSC of 25 August 2018, while enhancements of TMD and [O/N2] ratio started to appear later on 26 and 27 of August 2018.

Ionospheric anomalies through satellites can provide useful information about forthcoming earthquakes (EQs) over the epicentral regions. In this paper, we investigated seismo-ionospheric anomalies associated with the Mw 6.5 Sumatra earthquake that occurred in Indonesia on 06 December 2016 at 22:03 UT. We analyzed the total electron content (TEC) from Global Positioning System (GPS) signals received at the nearby stations around the epicenter. Furthermore, we also studied the TEC in local daytime and nighttime from Swarm satellites to confirm the EQ-induced ionospheric perturbations. The TEC showed significant perturbation within 5–10 days before the main shock in the form of positive anomalies beyond the upper bound. Similarly, Swarm satellites also validated the anomalies observed in the GPS TEC from nearby operating stations within 5–10 days over the epicentral region. The geomagnetic indices were quiet for the observed TEC anomalies within 5–10 days before the main shock, having Dst ≤ −40 nt and Kp ≤ 3. This study reveals the legitimate anomalies mainly associated with the EQ and suggests using the TEC from GPS and other satellites to look for possible future precursors with a more equipped satellite cluster.

Landslides are a recurrent environmental hazard in hilly regions and affect the socioeconomic development in Pakistan. The current study area is the tourism and hydro energy hub of Pakistan and is affected by environmental hazard. A landslide susceptibility mapping (LSM) of the Hindu Kush Himalayan, Swat District, Pakistan, can be created to reduce demographic losses due to landslides. This current study is conducted to apply three bivariate models, including weights of evidence (WOE), frequency ratio (FR), and information value (IV) for an LSM that has not been explored or applied in the current study area. For this purpose, first, an inventory map of 495 landslides was constructed from both ground and satellite data and randomly divided into training (70%) and testing (30%) datasets. Furthermore, 10 conditioning factors (elevation, slope, aspect, curvature, fault, rainfall, land use land cover (LULC), lithology, road, and drainage) used for the mapping of landslides were prepared in ArcGIS 10.8. Finally, LSM is generated based on WOE, FR, and IV models and validated the performance of LSM models using the area under receiver operating characteristic curve (AUROC). The findings of success rate curve (SRC) of the WOE, FR, and IV models were 67%, 93%, and 64%, respectively, while the prediction rate curves (PRCs) of the three models were 87%, 95%, and 73%, respectively. The validation results for WOE, FR, and IV justified that the FR model is the most reliable technique of all three of these models to produce the highest accuracy LSM for the present study area. Policymakers can use the findings of current research work to mitigate the loss due to landslide hazard.

The identification of atmospheric and ionospheric variations through multiple remote sensing and global navigation satellite systems (GNSSs) has contributed substantially to the development of the lithosphere-atmosphere-ionosphere coupling (LAIC) phenomenon over earthquake (EQ) epicenters. This study presents an approach for investigating the Petrolia EQ (Mw 6.2; dated 20 December 2021) and the Monte Cristo Range EQ (Mw 6.5; dated 15 May 2020) through several parameters to observe the precursory signals of various natures. These parameters include Land Surface Temperature (LST), Air Temperature (AT), Relative Humidity (RH), Air Pressure (AP), Outgoing Longwave Radiations (OLRs), and vertical Total Electron Content (TEC), and these are used to contribute to the development of LAIC in the temporal window of 30 days before and 15 days after the main shock. We observed a sharp increase in the LST in both the daytime and nighttime of the Petrolia EQ, but only an enhancement in the daytime LST for the Monte Cristo Range EQ within 3–7 days before the main shock. Similarly, a negative peak was observed in RH along with an increment in the OLR 5–7 days prior to both impending EQs. Furthermore, the Monte Cristo Range EQ also exhibited synchronized ionospheric variation with other atmospheric parameters, but no such co-located and synchronized anomalies were observed for the Petrolia EQ. We also applied machine learning (ML) methods to confirm these abrupt variations as anomalies to further aid certain efforts in the development of the LAIC in order to forecast EQs in the future. The ML methods also make prominent the variation in the different data.

We analyze vertical total electron content (vTEC) variations from the Global Navigation Satellite System (GNSS) at different latitudes in different continents of the world during the geomagnetic storms of June 2015, August 2018, and November 2021. The resulting ionospheric perturbations at the low and mid-latitudes are investigated in terms of the prompt penetration electric field (PPEF), the equatorial electrojet (EEJ), and the magnetic H component from INTERMAGNET stations near the equator. East and Southeast Asia, Russia, and Oceania exhibited positive vTEC disturbances, while South American stations showed negative vTEC disturbances during all the storms. We also analyzed the vTEC from the Swarm satellites and found similar results to the retrieved vTEC data during the June 2015 and August 2018 storms. Moreover, we observed that ionospheric plasma tended to increase rapidly during the local afternoon in the main phase of the storms and has the opposite behavior at nighttime. The equatorial ionization anomaly (EIA) crest expansion to higher latitudes is driven by PPEF during daytime at the main and recovery phases of the storms. The magnetic H component exhibits longitudinal behavior along with the EEJ enhancement near the magnetic equator.

Identifying atmospheric and ionospheric anomalies based on remote sensing satellites has contributed highly to develop the hypothesis of lithosphere-atmosphere-ionosphere coupling over the earthquake (EQ) epicenter during the seismic preparation period. This article has investigated the variations of potential EQ precursor in daytime and nighttime land surface temperature (LST) before and after the 2019 Pakistan EQ from Moderate Resolution Imaging Spectroradiometer (MODIS) satellite. The nighttime LST values of MODIS exhibit temporal anomalies during nighttime period within a time window of five days before and after the main shock day. Furthermore, the LST values predicted by artificial neural network (ANN) validate the significant enhancement in nighttime time series of MODIS. The nighttime LST anomalies obtained from the observation and ANN prediction are more than 20% and 7% of normal distribution beyond the confidence bounds, respectively, within five days after the main shock. Likewise, the ionospheric anomaly from daily total electron content (TEC) values at Sukkur Global Positioning System (GPS) station confirms the EQ associated ionospheric perturbations on the day after the main shock. The Global Ionospheric Maps (GIMs) also show the TEC anomalies during 1000–1400 LT on September 25, 2019.

Global Navigation Satellite System (GNSS)- and Remote Sensing (RS)-based Earth observations have a significant approach on the monitoring of natural disasters. Since the evolution and appearance of earthquake precursors exhibit complex behavior, the need for different methods on multiple satellite data for earthquake precursors is vital for prior and after the impending main shock. This study provided a new approach of deep machine learning (ML)-based detection of ionosphere and atmosphere precursors. In this study, we investigate multi-parameter precursors of different physical nature defining the states of ionosphere and atmosphere associated with the event in Japan on 13 February 2021 (Mw 7.1). We analyzed possible precursors from surface to ionosphere, including Sea Surface Temperature (SST), Air Temperature (AT), Relative Humidity (RH), Outgoing Longwave Radiation (OLR), and Total Electron Content (TEC). Furthermore, the aim is to find a possible pre-and post-seismic anomaly by implementing standard deviation (STDEV), wavelet transformation, the Nonlinear Autoregressive Network with Exogenous Inputs (NARX) model, and the Long Short-Term Memory Inputs (LSTM) network. Interestingly, every method shows anomalous variations in both atmospheric and ionospheric precursors before and after the earthquake. Moreover, the geomagnetic irregularities are also observed seven days after the main shock during active storm days (Kp > 3.7; Dst < −30 nT). This study demonstrates the significance of ML techniques for detecting earthquake anomalies to support the Lithosphere-Atmosphere-Ionosphere Coupling (LAIC) mechanism for future studies.

Global Navigation Satellite System (GNSS)-based ionospheric anomalies are nowadays used to identify a possible earthquake (EQ) precursor and hence a new research topic in seismic studies. The current study also aims to provide an investigation of ionospheric anomalies associated to EQs. In order to study possible pre-and post-seismic perturbations during the preparation phase of large-magnitude EQs, statistical and machine learning algorithms are applied to Total Electron Content (TEC) from the Global Positioning System (GPS) and Global Ionosphere Maps (GIMs). We observed TEC perturbation from the Sukkur (27.8° N, 68.9° E) GNSS station near the epicenter of Mw 5.4 Mirpur EQ within 5–10 days before the main shock day by implementing machine learning and statistical analysis. However, no TEC anomaly occurred in GIM-TEC over the Mirpur EQ epicenter. Furthermore, machine learning and statistical techniques are also implemented on GIM TEC data before and after the Mw 7.7 Awaran, where TEC anomalies can be clearly seen within 5–10 days before the seismic day and the subsequent rise in TEC during the 2 days after the main shock. These variations are also evident in GIM maps over the Awaran EQ epicenter. The findings point towards a large emission of EQ energy before and after the main shock during quiet storm days, which aid in the development of lithosphere ionosphere coupling. However, the entire analysis can be expanded to more satellite and ground-based measurements in Pakistan and other countries to reveal the pattern of air ionization from the epicenter through the atmosphere to the ionosphere.

The search for Earthquake (EQ) precursors in the ionosphere and atmosphere from satellite data has provided significant information about the upcoming main shock. This study presents the abnormal atmospheric and ionospheric perturbations associated with the Mw 7.2 Haiti EQ on 14 August 2021 at geographical coordinates (18° N, 73° W) and shallow hypocentral depth of 10 km from the data of permanent Global Navigation Satellite System (GNSS) stations near the epicenter, followed by Swarm satellites data. The total vertical electron (VTEC) anomalies occur within a 5-day window before the main shock in the analysis of nearby operation stations, followed by Swarm (A and C satellites) ionospheric anomalies in the same 5-day window before the main shock. Moreover, the geomagnetic activities are completely quiet within 10 days before and 10 days after the main shock. Similarly, the atmospheric parameters endorse the EQ anomalies within 5 days before the main shock day. The evolution of gases from the lithosphere at the epicentral region possessed significant atmospheric and ionospheric perturbations within the EQ preparation period of 5-day before the main shock under the hypothesis of Lithosphere-Atmosphere-Ionosphere Coupling (LAIC).

Due to the scarcity of studies linking the variability of rainfall and population growth in the capital cities of Northeastern Brazil (NEB), the purpose of this study is to evaluate the variability and multiscale interaction (annual and seasonal), and in addition, to detect their trends and the impact of urban growth. For this, monthly rainfall data between 1960 and 2020 were used. In addition, the detection of rainfall trends on annual and seasonal scales was performed using the Mann–Kendall (MK) test and compared with the phases of El Niño-Southern Oscillation (ENSO) and Pacific Decadal Oscillation (PDO). The relationship between population growth data and rainfall data for different decades was established. Results indicate that the variability of multiscale urban rainfall is directly associated with the ENSO and PDO phases, followed by the performance of rain-producing meteorological systems in the NEB. In addition, the anthropic influence is shown in the relational pattern between population growth and the variability of decennial rainfall in the capitals of the NEB. However, no capital showed a significant trend of increasing annual rainfall (as in the case of Aracaju, Maceió, and Salvador). The observed population increase in the last decades in the capitals of the NEB and the notable decreasing trend of rainfall could compromise the region’s water security. Moreover, if there is no strategic planning about water bodies, these changes in the rainfall pattern could be compromising.

The ionospheric response and the associated mechanisms to geomagnetic storms are very complex, particularly during the February 2014 multiphase geomagnetic storm. In this paper, the low-latitude ionosphere responses and their coupling mechanisms, during the February 2014 multiphase geomagnetic storm, are investigated from ground-based magnetometers and global navigation satellite system (GNSS), and space weather data. The residual disturbances between the total electron content (TEC) of the International GNSS Service (IGS) global ionospheric maps (GIMs) and empirical models are used to investigate the storm-time ionospheric responses. Three clear sudden storm commencements (SSCs) on 15, 20, and 23 February are detected, and one high speed solar wind (HSSW) event on 19 February is found with the absence of classical SSC features due to a prevalent magnetospheric convection. The IRI-2012 shows insufficient performance, with no distinction between the events and overestimating approximately 20 TEC units (TECU) with respect to the actual quiet-time TEC. Furthermore, the median average of the IGS GIMs TEC during February 2014 shows enhanced values in the southern hemisphere, whereas the IRI-2012 lacks this asymmetry. Three low-latitude profiles extracted from the IGS GIM data revealed up to 20 TECU enhancements in the differential TEC. From these profiles, longer-lasting TEC enhancements are observed at the dip equator profiles than in the profiles of the equatorial ionospheric anomaly (EIA) crests. Moreover, a gradual increase in the global electron content (GEC) shows approximately 1 GEC unit of differential intensification starting from the HSSW event, while the IGS GIM profiles lack this increasing gradient, probably located at higher latitudes. The prompt penetration electric field (PPEF) and equatorial electrojet (EEJ) indices estimated from magnetometer data show strong variability after all four events, except the EEJ’s Asian sector. The low-latitude ionosphere coupling is mainly driven by the variable PPEF, DDEF (disturbance dynamo electric fields), and Joule heating. The auroral electrojet causing eastward PPEF may control the EIA expansion in the Asian sector through the dynamo mechanism, which is also reflected in the solar-quiet current intensity variability.

The seismo-ionospheric anomalies may provide some insights about the earthquake. However, pre-seismic ionospheric anomalies are still challenging. In this paper, seismo-ionospheric anomalies are investigated before the September 24, 2013 (Mw = 7.7) Awaran (Pakistan) earthquake from GPS TEC (Total Electron Content) and COSMIC (Constellation Observing System for Meteorology, Ionosphere and Climate) data. The TEC data are showing anomalies on September 21, 2013. The abnormality detected in the temporal data is about 10 TECU beyond the 30-day running median. The percentage deviation of the TEC anomaly on September 21, is 30% above the upper confidence interval. The anomalies prevail 5° in Latitude and 10° in Longitude over the epicenter. The spatial and temporal data of TEC showed anomalies in TEC from UT = 08 to UT = 12. In addition, the enhancement on September 21, 2013 is also very significant in COSMIC data. The results of COSMIC completely agreed with GPS TEC anomalies on September 21, 2013. The percentage deviation of the peak plasma frequency on September 21 is 5% of the normal distribution. The storm indices are quiet before and after the earthquake. The pre-sesimic ionospheric anomalies are most probably associated with the 2013 Mw = 7.7 Awaran (Pakistan) earthquake. Keywords: Seismo-ionospheric disturbance, GPS TEC, COSMIC, Earthquake

Earthquakes can generate thermal anomalies in the atmosphere at low altitudes. Pending well-focused detailed studies, such phenomenon may be referred to as a precursor for earthquake prediction. However, today the pre-earthquake thermal anomalies are not clear enough. In this paper, the thermal anomalies prior to the April 25, 2015 Mw 7.8 Gorkha (Nepal) earthquake are investigated from the Moderate Resolution Imaging Spectroradiometer (MODIS) land surface temperature (LST), air temperature and Outgoing Longwave Radiations (OLR) data. The 2D and 3D wavelet transformation techniques are used to interpret the real time enhancement of the daily MODIS and OLR data before the impending earthquake. Using the wavelet density spectrum, pre-earthquake anomalies in MODIS and OLR are found in connection to the impending earthquake. The spatial images of MODIS and OLR show the evolutionary pattern of the emanation of ions from the epicenter and the surrounding area. The most important feature revealed by the spatial analysis is the eastward migration of temperature clouds due to a strong electric field. The satellite based LST data showed deviation, which crosses the upper bound by 5 °C. All the observations in our case study strongly support the notion of pre-earthquake thermal anomalies. Based on the analysis of the results, it can be concluded that the overabundance of ions from the seismogenic zone is responsible for prompting large temperature perturbations in atmospheric layers.

Rice is a primary food for more than three billion people worldwide and cultivated on about 12% of the world’s arable land. However, more than 88% production is observed in Asian countries, including Pakistan. Due to higher population growth and recent climate change scenarios, it is crucial to get timely and accurate rice yield estimates and production forecast of the growing season for governments, planners, and decision makers in formulating policies regarding import/export in the event of shortfall and/or surplus. This study aims to quantify the rice yield at various phenological stages from hyper-temporal satellite-derived-vegetation indices computed from time series Sentinel-II images. Different vegetation indices (viz. NDVI, EVI, SAVI, and REP) were used to predict paddy yield. The predicted yield was validated through RMSE and ME statistical techniques. The integration of PLSR and sequential time-stamped vegetation indices accurately predicted rice yield (i.e., maximum R2 = 0.84 and minimum RMSE = 0.12 ton ha−1 equal to 3% of the mean rice yield). Moreover, our results also established that optimal time spans for predicting rice yield are late vegetative and reproductive (flowering) stages. The output would be useful for the farmer and decision makers in addressing food security.