Atmospheric Stability Research Articles

Observed Subdaily Variations in Air–Sea Turbulent Heat Fluxes under Different Marine Atmospheric Boundary Layer Stability Conditions in the Gulf Stream

Abstract Based on data collected from 14 buoys in the Gulf Stream, this study examines how hourly air–sea turbulent heat fluxes vary on subdaily time scales under different boundary layer stability conditions. The annual mean magnitudes of the subdaily variations in latent and sensible heat fluxes at all stations are 40 and 15 W m−2, respectively. Under near-neutral conditions, hourly fluctuations in air–sea humidity and temperature differences are the major drivers of subdaily variations in latent and sensible heat fluxes, respectively. When the boundary layer is stable, on the other hand, wind anomalies play a dominant role in shaping the subdaily variations in latent and sensible heat fluxes. In the context of a convectively unstable boundary layer, wind anomalies exert a strong controlling influence on subdaily variations in latent heat fluxes, whereas subdaily variations in sensible heat fluxes are equally determined by air–sea temperature difference and wind anomalies. The relative contributions by all physical quantities that affect subdaily variations in turbulent heat fluxes are further documented. For near-neutral and unstable boundary layers, the subdaily contributions are O(2) and O(1) W m−2 for latent and sensible heat fluxes, respectively, and they are less than O(1) W m−2 for turbulent heat fluxes under stable conditions. Significance Statement High-resolution buoy observations of air–sea variables in the Gulf Stream provide the opportunity to investigate the physical factors that determine subdaily variations in air–sea turbulent heat fluxes. This study addresses two key points. First, the observed subdaily amplitudes of heat fluxes are related to various processes, including wind fields and air–sea thermal effect differences. Second, the global sea surface heat budget is known to not be in near-zero balance and it ranges from several to tens of watts per square meter. Therefore, consideration of the relatively strong influence of subdaily variability in air–sea turbulent heat fluxes could provide a new strategy for solving the global heat budget balance problem.

Perbandingan Stabilitas Atmosfer di Permukaan Manifestasi Panas Bumi dengan Tanpa Manifestasi Panas Bumi menggunakan Bilangan Richardson

The stability of the atmosphere can affect the weather conditions in an area. This study aims to compare atmospheric stability carried out in Lake Linow which has two surfaces, namely the surface of geothermal manifestation and the surface without geothermal manifestation. In determining the stability of the atmosphere, researchers use the Richardson number method. To determine atmospheric stability using the Richardson number, data on wind speed and temperature at each altitude are needed. The measured height is 2 meters and 4 meters from the surface. The results of this study show that on the surface geothermal manifestations have unstable atmospheric stability with Richardson numbers -1.4729 to -0.1420, while on the surface without geothermal manifestations have 3 states of atmospheric stability, namely: stable, neutral, and unstable, with Richardson numbers -1.1316 to 0.9808. Therefore, it can be concluded that the state of the atmosphere on the surface of geothermal manifestations is more unstable, compared to surfaces without geothermal manifestations.

Investigating turbulence distribution in the lower atmosphere using time-lapse imagery from a camera bank

The atmosphere’s surface layer (first 50–100 m above the ground) is extremely dynamic and is influenced by surface radiative properties, roughness, and atmospheric stability. Understanding the distribution of turbulence in the surface layer is critical to many applications, such as directed energy and free space optical communications. Several measurement campaigns in the past have relied on weather balloons or sonic detection and ranging (SODAR) to measure turbulence up to the atmospheric boundary layer. However, these campaigns had limited measurements near the surface. We have developed a time-lapse imaging technique to profile atmospheric turbulence from turbulence-induced differential motion or tilts between features on a distant target, sensed between pairs of cameras in a camera bank. This is a low-cost and portable approach to remotely sense turbulence from a single site without the deployment of sensors at the target location. It is thus an excellent approach to study the distribution of turbulence in low altitudes with sufficiently high resolution. In the present work, the potential of this technique was demonstrated. We tested the method over a path with constant turbulence. We explored the turbulence distribution with height in the first 20 m above the ground by imaging a 30 m water tower over a flat terrain on three clear days in summer. In addition, we analyzed time-lapse data from a second water tower over a sloped terrain. In most of the turbulence profiles extracted from these images, the drop in turbulence with altitude in the first 15 m or so above the ground showed a h m dependence, where the exponent m varied from −0.3 to −1.0, quite contrary to the widely used value of −4/3.

A Systematic Investigation of the Applicability of Taylor’s Hypothesis in an Idealized Surface Layer

Taylor’s Frozen Turbulence Hypothesis (TH) is a critical assumption in turbulent theory and practice which allows time series of point measurements of turbulent variables to be translated to the spatial domain via the mean wind. Using a 3D array of fibre-optic distributed temperature sensing in the atmospheric surface layer over an idealized desert site we present a systematic investigation of the applicability of Taylor’s Hypothesis to atmospheric surface layer flows over a variety of conditions: unstable, near-neutral, and stable atmospheric stabilities; and multiple measurement heights between the surface and 3 m above ground level. Both spatially integrated and spatially scale-dependent eddy velocities are investigated by means of time-lagged streamwise two-point correlations and compared to the mean Eulerian wind. We find that eddies travel slower than predicted by TH at small spatial separations, as predicted by TH at separations typically between 5 and 16 m, and faster than predicted by TH at larger spatial separations. In unstable atmospheric conditions the spatial separation at which eddy velocity is larger than Eulerian velocity decreases with height.

Surface heat fluxes at coarse blocky Murtèl rock glacier (Engadine, eastern Swiss Alps)

Abstract. We estimate the surface energy balance (SEB) of the Murtèl rock glacier, a seasonally snow-covered permafrost landform with a ventilated coarse blocky active layer (AL) located in the eastern Swiss Alps. We focus on the parameterisation of the turbulent heat fluxes. Seasonally contrasting atmospheric conditions occur in the Murtèl cirque, with downslope katabatic jets in winter and a strongly unstable atmosphere over the heated blocky surface in summer. We use a novel comprehensive sensor array both above the ground surface and in the coarse blocky AL to track the rapid coupling by convective heat and moisture fluxes between the atmosphere, the snow cover, and the AL for the time period September 2020–September 2022. The in situ sensor array includes a sonic anemometer for eddy-covariance flux above-ground and sub-surface long-wave radiation measurements in a natural cavity between the AL blocks. During the thaw seasons, the measurements suggest an efficient (∼ 90 %) export of the available net radiation by sensible and latent turbulent fluxes, thereby strongly limiting the heat available for melting ground ice. Turbulent export of heat and moisture drawn from the permeable AL contributes to the well-known insulating effect of the coarse blocky AL and partly explains the climate resiliency of rock glaciers. This self-cooling capacity is counteracted by an early snow melt-out date, exposing the low-albedo blocky surface to the intense June–July insolation and causing reduced evaporative cooling due to exacerbated moisture scarcity in the near-surface AL during dry spells. With climate change, earlier snowmelt and increased frequency, duration, and intensity of heat waves and droughts are projected. Regarding the parameterisation of the turbulent fluxes, we estimated the year-round turbulent fluxes using a modified Louis (1979) scheme. The monthly SEB is closed within 20 W m−2 except during the snowmelt months and under katabatic drainage winds in winter. Detected sensible turbulent fluxes from nocturnal ventilation processes, although a potentially important ground cooling mechanism, are within our 20 W m−2 uncertainty because nighttime wind speeds are low. Wintertime katabatic wind speeds needed to be scaled to close the SEB, which hints at the limits of parameterisations based on the Monin–Obukhov similarity theory in complex mountain terrain and katabatic drainage winds. The present work contributes to the process understanding of the SEB and climate sensitivity of coarse blocky landforms.

Chemical Reaction of FA Cations Enables Efficient and Stable Perovskite Solar Cells.

Organometal halide perovskite solar cells (PSCs) have received great attention owing to a rapid increase in power conversion efficiency (PCE) over the last decade. However, the deficit of long-term stability is a major obstacle to the implementation of PSCs in commercialization. The defects in perovskite films are considered as one of the primary causes. To address this issue, isocyanic acid (HNCO) is introduced as an additive into the perovskite film, in which the added molecules form covalent bonds with FA cations via a chemical reaction. This chemical reaction gives rise to an efficient passivation on the perovskite film, resulting in an improved film quality, a suppressed non-radiation recombination, a facilitated carrier transport, and optimization of energy band levels. As a result, the HNCO-based PSCs achieve a high PCE of 24.41% with excellent storage stability both in an inert atmosphere and in air. Different from conventional passivation methods based on coordination effects, this work presents an alternative chemical reaction for defect passivation, which opens an avenue toward defect-mitigated PSCs showing enhanced performance and stability.

Fabrication of stable monolayer liquid marbles with reduced particle coverage and locomotion on hydrophilic surface

Liquid marbles are non-wetting, particle-covered microdroplets with a core-shell structure that are used in sample transport, material synthesis, and real-time sensing. Optimizing the distribution of shell particles remains a challenge, due to a tendency for aggregation via spontaneous assembly, which often leads to multilayered structures. Here, we outline a simple method for fabricating water-filled, monolayer liquid marbles with adjustable particle coverage rates, greatly reducing particle consumption. The soft liquid marbles are enclosed by a small quantity of modified polystyrene microspheres and display good atmospheric stability. The rolling behavior of flexible liquid marbles with wide coverage rates is then characterized. Contrary to common perception, the marbles with transparent openings exhibit high maneuverability on hydrophilic surfaces, and also excel in fusion, reaction and surface cleaning, with an elongated operational duration and a wide visualization range. The study provides new insights into the implementation of liquid marble-based miniaturized platforms.

Intrinsically Stretchable Polymer Semiconductor with Regional Conjugation for Stretchable Electronics

AbstractThe development of intrinsically stretchable polymer semiconductor holds substantial promise in the field of wearable electronics. However, charge transport mobility is typically compromised in existing stretchable semiconductors to achieve the desired stretchability. Herein, a novel “regional conjugation” strategy is proposed to design an intrinsically stretchable polymer semiconductor oligo‐diketopyrrolopyrrole‐thieno[3,2‐b]thiophene (DPPTT)–urethane, in which oligo‐DPPTT conjugated units and alkyl urethane nonconjugated units are introduced. The regional conjugation of oligo‐DPPTT in the polymer backbone endows DPPTT–urethane with good molecular packing, leading to a high mobility of up to 1.7 cm2 V−1 s−1. Additionally, incorporating alkyl urethane nonconjugated units in the backbone can reduce film crystallinity and chain aggregation, which contribute to the stretchability of the polymer thin film. Consequently, fully stretchable transistors retain carrier mobility even at 100% biaxial tensile strain. Furthermore, the fully stretchable organic field‐effect transistor arrays show remarkable charge transport reversibility and durability after 1000 stretch–release cycles at 25% strain. Additionally, the device exhibits extraordinary electrical stability in air atmosphere. Overall, these results indicate that the “regional conjugation” strategy provides an effective and promising methodology to design intrinsically stretchable and high‐performance polymer semiconductor that can advance the development of soft and wearable electronics.

Boosting the Carrier Lifetime and Optical Activity of CsPbX3 Nanocrystals through Aromatic Ligand Passivation.

Ligand engineering is crucial for tuning the structural and optoelectronic properties of perovskite nanocrystals (NCs), which also improves their stability. In contrast to the typically used long-chain alkylamine ligands, we successfully introduced an aromatic 1-(p-tolyl)ethylamine (PTEA) ligand to synthesize the CsPbX3 (X = Br or I) NCs. The CsPbI3 and CsPbBr3 NCs demonstrated long carrier lifetimes of ∼877 and 49 ns, respectively, as well as high photoluminescence quantum yields (PLQYs) of ∼99% and 95%, respectively. Furthermore, such NCs realized excellent long-term stability in an ambient atmosphere without obvious degradation over one month. All of these properties were better than the properties of NCs coated with the conventional alkylamine ligands. The high performance of these NCs was discussed with the effective surface passivation by PTEA. Our finding suggests a facile and effective ligand (PTEA) for modulating perovskites, achieving enhancement of both the carrier lifetime and the PLQY.

Air-Stable and Eco-Friendly Symmetrical Imine with Thiadiazole Moieties in Neutral and Protonated form for Perovskite Photovoltaics.

This paper proposes molecular and supramolecular concepts for potential application in perovskite solar cells. New air-stable symmetrical imine, with thiadiazole moieties PPL2: (5E,6E)-N2,N5-bis(4-(diphenylamino)benzylidene)-1,3,4-thiadiazole-2,5-diamine), as a hole-transporting material was synthesised in a single-step reaction, starting with commercially available and relatively inexpensive reagents, resulting in a reduction in the cost of the final product compared to Spiro-OMeTAD. Moreover, camphorsulfonic acid (CSA) in both enantiomeric forms was used to change the HOMO-LUMO levels and electric properties of the investigated imine-forming complexes. Electric, optical, thermal, and structural studies of the imine and its complexes with CSA were carried out to characterise the new material. Imine and imine/CSA complexes were also characterised in depth by the proton Nuclear Magnetic Resonance 1H NMR method. The position of nitrogen in the thidiazole ring influences the basicity of donor centres, which results in protonation in the imine bond. Simple devices of ITO/imine (with or without CSA(-) or CSA(+))/Ag/ITO architecture were constructed, and a thermographic camera was used to find the defects in the created devices. Electric behaviour was also studied to demonstrate conductivity properties under the forward current. Finally, the electrical properties of imine and its protonated form with CSA were compared with Spiro-OMeTAD. In general, the analysis of thermal images showed a very similar response of the samples to the applied potential in terms of the homogeneity of the formed organic layer. The TGA analysis showed that the investigated imine exhibits good thermal stability in air and argon atmospheres.

Unraveling the Role of 2D Ti3C2Tx MXene Nanosheets in Cu-Based Double Perovskite Active Layer for Enhanced Photovoltaic Performance.

Although the atmospheric stability of lead-free inorganic double perovskite (DP) solar cells (PSCs) looks promising, their further development is hampered by inadequate film quality and non-radiative carrier recombination at the interfaces. Herein, the incorporation of a newly developed intriguing class of 2D material Ti3C2Tx MXene nanosheets with the photo-absorbing Cu2AgBiI6 (CABI) active layer of a fully inorganic solar cell is reported. The highly conductive Ti3C2Tx nanosheets work as a multi-functional additive by tuning the band gap, reducing the non-radiative carrier recombination, and inhibiting carrier accumulation. In addition, the presence of Ti3C2Tx MXene increases the surface free energy of the perovskite film, which elevates the energy barrier for nucleation and realizes a highly crystalline CABI perovskite film. Primarily, the MXene modification accelerates the charge extraction and transport at the interfaces of the active layer, utilizing energy level alignment with the charge transport layers. Consequently, the photo-conversion efficiency (PCE) of the device with MXene is substantially enhanced to 1.50%. Moreover, the 2D Ti3C2Tx nanosheets increased the long-term stability of the devices by retaining 70% of the initial PCE after 1680 h. With regard to relieving the severe carrier recombination at the interfaces, this work sets a new paradigm toward imminent solar energy conversion.

Thermal Energy Storage Using Phase Change Materials in High-Temperature Industrial Applications: Multi-Criteria Selection of the Adequate Material.

Thermal energy storage (TES) plays an important role in industrial applications with intermittent generation of thermal energy. In particular, the implementation of latent heat thermal energy storage (LHTES) technology in industrial thermal processes has shown promising results, significantly reducing sensible heat losses. However, in order to implement this technology, a proper selection of materials is important. In this study, a new multi-criteria phase change material (PCM) selection methodology is presented, which considers relevant factors from an application and material handling point of view, such as hygroscopicity, metal compatibility (corrosion), level hazard, cost, and thermal and atmospheric stability. The methodology starts after setting up the system requirements where the PCM will be used, then a material screening is able to find all possible candidates that are listed with all available properties as listed before. Then, a color map is produced, with a qualitative assessment of material properties drawbacks, hazard level, melting enthalpy, and price. The experimentation starts with a preliminary set of tests on hygroscopicity and one-week corrosion test, which allows disregarding PCMs and selecting a short list of potential PCMs that would need further characterization before the final selection.

Evaluation of wind farm parameterizations in the WRF model under different atmospheric stability conditions with high-resolution wake simulations

Abstract. Wind farm parameterizations (WFPs) are used in mesoscale models for predicting wind farm power production and its impact on wind resources while considering the variability of the regional wind climate. However, the performance of WFPs is influenced by various factors including atmospheric stability. In this study, we compared two widely used WFPs in the Weather Research and Forecasting (WRF) model to large-eddy simulations (LES) of turbine wakes performed with the same model. The Fitch WFP and the explicit wake parameterization were evaluated for their ability to represent wind speed and turbulent kinetic energy (TKE) in a two-turbine wind farm layout under neutral, unstable, and stable atmospheric stability conditions. To ensure a fair comparison, the inflow conditions were kept as close as possible between the LES and mesoscale simulations for each type of stability condition, and the LES results were spatially aggregated to align with the mesoscale grid spacing. Our findings indicate that the performance of WFPs varies depending on the specific variable (wind speed or TKE) and the area of interest downwind of the turbine when compared to the LES reference. The WFPs can accurately depict the vertical profiles of the wind speed deficit for either the grid cell containing the wind turbines or the grid cells in the far wake, but not both simultaneously. The WFPs with an explicit source of TKE overestimate TKE values at the first grid cell containing the wind turbine; however, for downwind grid cells, agreement improves. On the other hand, WFPs without a TKE source underestimate TKE in all downwind grid cells. These agreement patterns between the WFPs and the LES reference are consistent under the three atmospheric stability conditions. However, the WFPs resemble less the wind speed and TKE from the LES reference under stable conditions than that under neutral or unstable conditions.

Green buffers near industrial plants, examples in Jefferson County TX

Green buffers between industrial plants and residential areas have multiple purposes in mitigating air pollution and protecting the environment. This study focuses on providing information on examples of green space development by industries near residential areas in Jefferson County, Texas. One represents a full buffer where residences were totally removed, another is a partial buffer where some residences were removed in a neighborhood and the third represents vegetation of a formerly industrial facility. Additional estimates on benefits of buffer spaces with respect to a pollutant plume are also considered. The study used geospatial and property data from 1966 to 2020. The transition from a residential area to a full or partial buffers highlights possible mitigation of air pollution impacts while serving additional functions, such as boosting stormwater infiltration and improving the aesthetic attractiveness of the area. The transformation from oil tanks to forested space may display the potential for environmental remediation, but not all spaces are suitable for fully forested buffers, as safety and other limitations should be considered. With the analysis of the classic Plume Equation insights can be provided regarding the possible impact of green buffers on ground-level pollutant concentrations. The calculations for atmospheric stability classes A-F and for various stack heights provide an understanding of how different conditions may affect the pollutant concentrations. The presence of a green buffer can play a crucial role in controlling ground-level pollution concentrations and reduce a community’s perception of risk.

Variability and trends of near-surface wind speed over the Tibetan Plateau: The role played by the westerly and Asian monsoon

Near-surface wind speed exerts profound impacts on many environmental issues, while the long-term (≥60 years) trend and multidecadal variability in the wind speed and its underlying causes in global high-elevation and mountainous areas (e.g., Tibetan Plateau) remain largely unknown. Here, by examining homogenized wind speed data from 104 meteorological stations over the Tibetan Plateau for 1961–2020 and ERA5 reanalysis datasets, we investigated the variability and long-term trend in the near-surface wind speed and revealed the role played by the westerly and Asian monsoon. The results show that the homogenized annual wind speed displays a decreasing trend (−0.091 m s−1 per decade, p < 0.05), with the strongest in spring (−0.131 m s−1 per decade, p < 0.05), and the weakest in autumn (−0.071 m s−1 per decade, p < 0.05). There is a distinct multidecadal variability of wind speed, which manifested in an prominent increase in 1961–1970, a sustained decrease in 1970–2002, and a consistent increase in 2002–2020. The observed decadal variations are likely linked to large-scale atmospheric circulation, and the correlation analysis unveiled a more important role of westerly and East Asian winter monsoon in modulating near-surface wind changes over the Tibetan Plateau. The potential physical processes associated with westerly and Asian monsoon changes are in concordance with wind speed change, in terms of overall weakened horizontal air flow (i.e., geostrophic wind speed), declined vertical thermal and dynamic momentum transfer (i.e., atmospheric stratification thermal instability and vertical wind shear), and varied Tibetan Plateau vortices. This indicates that to varying degrees these processes may have contributed to the changes in near-surface wind speed over the Tibetan Plateau. This study has implications for wind power production and soil wind erosion prevention in the Tibetan Plateau.

Impact of Atmospheric Compressibility and Stokes Drift on the Vertical Transport of Heat and Constituents by Gravity Waves

AbstractVertical transport of heat and atmospheric constituents by gravity waves plays a crucial role in shaping the thermal and constituent structure of the middle atmosphere. We show that atmospheric mixing by non‐breaking waves can be described as a diffusion process where the potential temperature (KH) and constituent (KWave) diffusivities depend on the compressibility of the wave fluctuations and the vertical Stokes drift imparted to the atmosphere by the wave spectrum. KH and KWave are typically much larger than the eddy diffusivity (Kzz), arising from the turbulence generated by breaking waves, and can exceed several hundred m2s−1 in regions of strong wave dissipation. We also show that the total diffusion of heat and constituents caused by waves, turbulence, and the thermal motion of molecules, is enhanced in the presence of non‐breaking waves by a factor that is proportional to the variance of the wave‐driven lapse rate fluctuations. Diffusion enhancements of both heat and constituents of 50% or more can be experienced in regions of low atmospheric stability, where the lapse rate fluctuations are large. These important transport effects are not currently included in most global chemistry‐climate models, which typically only consider the eddy diffusion that is induced when the unresolved, but parameterized waves, experience dissipation. We show that the theoretical results compare favorably with observations of the mesopause region at midlatitudes and describe how the theory may be used to more fully account for the unresolved wave transport in global models.

Modulations of dust aerosols on precipitation: Evidence from a typical heavy sandstorm event

Dust aerosols pose significant modulations on precipitation, varying with precipitation intensity. In particular, the microphysical influences of dust aerosols on the clouds and diurnal variations of heavy precipitation exhibit substantial uncertainties. Taking the typical heavy sandstorm event occurring on March 15, 2021 in China as an example, we investigate the distinct modulations of dust aerosols on precipitation by the Weather Research and Forecasting model coupled to Chemistry (WRF-Chem) model and multiple observational and reanalysis data. Note that dust aerosols exert opposite influences on moderate and heavy precipitation. The inhibiting effect of dust aerosols on the intensity of light precipitation was observed as dust aerosols increase atmospheric stability by cooling in the lower layers, warming in the upper layers, reducing CAPE (Convective Available Potential Energy), increasing CIN (Convective Inhibition). Additionally, dust aerosols hinder the transformation of cloud droplets to raindrops, leading to decreased mass of rainwater (Qr) and snow (Qs) in the lower and middle atmosphere. In contrast, dust aerosols delay but intensify heavy rainfalls. During the night of heavy rain, dust aerosols can can warm the lower layers and cool the upper layers, increase CAPE and decrease CIN, thereby enhancing convection, making it easier for small cloud droplets to grow into larger raindrops through collision and coalescence. Simultaneously, small cloud droplets in the mid-levels ascend to higher altitudes where they freeze, transforming into ice particles. This process, coupled with the Bergeron-Findeisen mechanism, enhances the intensity of nighttime. These findings provide novel evidence to improve our understanding of dust aerosols-precipitation interaction.

Aviation accidents related to atmospheric instability in the United States (2000–2020)

Aviation accidents related to atmospheric instability in the United States (2000–2020)

A novel halide hybrid white light emitter with remarkable structural and emission stability

Addressing the issues of self-absorption effect among different components in traditional white light-emitting diodes (WLEDs) and the high cost associated with rare earth elements, we propose the potential utility of zero-dimensional organic–inorganic metal halide hybrids (0D OIMHs) in white WLEDs fabrication. In this study, we synthesized a novel 0D OIMH, (C5N2H14)ZnCl4, characterized by broadband white light emission. The (C5N2H14)ZnCl4 crystal features a short ultraviolet cutoff situated at 190 nm and emits cold white light with a Color Rendering Index (CRI) of 73 under an excitation wavelength of 260 nm. Furthermore, (C5N2H14)ZnCl4 demonstrates exceptional structural and emission stability in an air atmosphere over a span of three months.

Chemical release risk assessment in earthquake: Natech event scenario

IntroductionIn recent decades, an increase in natural disasters has led to the Natech events occurrence, which are technological accidents triggered by natural hazards. This study aims to assess the risk of chemical release of the H2S toxic gas following an earthquake from the refinery. MethodThe Natech risk assessment was carried out in potential earthquake risk scenarios via a semi-quantitative method using Rapid-N software. The reference scenario was obtained using the Australian SMUG Model, which has presented precise prioritization. FindingsThe Natech risk assessment has shown 40 possible earthquake risk scenarios from 5 enclosing Faults in Tehran (North Tehran, Mosha, North Ray, South Ray, and Kahrizak) with an extension of more than 15 km. The earthquake moment magnitudes of 7.5, 7.1, 6.5, and 5.9 were obtained on the Richter scale and at two Focal depths of 5 and 10 km. The South-Ray Fault (HSR-1) was selected as the reference scenario with a moment magnitude of 7.5 on the Richter scale at a Focal depth of 5 km and a distance of 5.5 km from the earthquake epicenter. The highest probability of Natech risk has shown the release of H2S toxic gas in an area with a radius of 6.59 km from the studied vessel in atmospheric stability conditions. ConclusionLegislation in line with risk reduction, planning for the transfer of potentially dangerous industries to outside urban areas, monitoring land use laws, and promoting coordination between the government, industries, and people in the management of risks caused by industries containing hazardous substances and processes, assessing the risk of Natech events and their consequences, strengthening risk communication and the installation of early warning systems, and public education to adopt personal and socially safe behaviors when facing natural disasters are recommended.