Feed aggregator

Developing an Explainable Variational Autoencoder (VAE) Framework for Accurate Representation of Local Circulation in Taiwan

JGR–Atmospheres - Mon, 06/24/2024 - 19:34
Abstract

This study develops an explainable variational autoencoder (VAE) framework to efficiently generate high-fidelity local circulation patterns in Taiwan, ensuring an accurate representation of the physical relationship between generated local circulation and upstream synoptic flow regimes. Large ensemble semi-realistic simulations were conducted using a high-resolution (2 km) model, TaiwanVVM, where critical characteristics of various synoptic flow regimes were carefully selected to focus on the effects of local circulation variations. The VAE was constructed to capture essential representations of local circulation scenarios associated with the lee vortices by training on the ensemble data set. The VAE's latent space effectively captures the synoptic flow regimes as controlling factors, aligning with the physical understanding of Taiwan's local circulation dynamics. The critical transition of flow regimes under the influence of southeasterly synoptic flow regimes is also well represented in the VAE's latent space. This indicates that the VAE can learn the nonlinear characteristics of the multiscale interactions involving the lee vortex. The latent space within VAE can serve as a reduced-order model for predicting local circulation using synoptic wind speed and direction. This explainable VAE binds the physical reasoning to the predictions of the local circulation that ensures the physical examination of the uncertainty in accelerating the local weather assessments under various climate change scenarios.

Applying recession models for low-flow prediction: a comparison of regression and matching strip approaches

Natural Hazards and Earth System Sciences - Mon, 06/24/2024 - 18:14
Applying recession models for low-flow prediction: a comparison of regression and matching strip approaches
Michael Margreth, Florian Lustenberger, Dorothea Hug Peter, Fritz Schlunegger, and Massimiliano Zappa
Nat. Hazards Earth Syst. Sci. Discuss., https//doi.org/10.5194/nhess-2024-78,2024
Preprint under review for NHESS (discussion: open, 0 comments)
Recession models (RM) are crucial for observing the low flow behavior of a catchment. We developed two novel RM, which are designed to represent slowly draining catchment conditions. With a newly designed low flow prediction procedure we tested the prediction capability of these two models and three others from literature. One of our novel products delivered the best results, because it best represents the slowly draining catchment conditions.

Tsunami detection methods for Ocean-Bottom Pressure Gauges

Natural Hazards and Earth System Sciences - Mon, 06/24/2024 - 18:14
Tsunami detection methods for Ocean-Bottom Pressure Gauges
Cesare Angeli, Alberto Armigliato, Martina Zanetti, Filippo Zaniboni, Fabrizio Romano, Hafize Başak Bayraktar, and Stefano Lorito
Nat. Hazards Earth Syst. Sci. Discuss., https//doi.org/10.5194/nhess-2024-113,2024
Preprint under review for NHESS (discussion: open, 0 comments)
To issue precise and timely tsunami alerts, detecting the propagating tsunami is fundamental. The most used instruments are pressure sensors positioned at the ocean bottom, called Ocean-Bottom Pressure Gauges (OBPGs). In this work, we study four different techniques that allow to recognize a tsunami as soon as it is recorded by an OBPG and a methodology to calibrate them. The techniques are compared in terms of their ability to detect and characterize the tsunami wave in real time.

Evaluation of radiation schemes in the CMA-MESO model using high time-resolution radiation measurements in China: I. Long-wave radiation

Geoscientific Model Development - Mon, 06/24/2024 - 17:57
Evaluation of radiation schemes in the CMA-MESO model using high time-resolution radiation measurements in China: I. Long-wave radiation
Junli Yang, Weijun Quan, Li Zhang, Jianglin Hu, Qiying Chen, and Martin Wild
Geosci. Model Dev. Discuss., https//doi.org/10.5194/gmd-2024-74,2024
Preprint under review for GMD (discussion: open, 2 comments)
Due to the difficulties involved in the measurements of the Downward long-wave irradiance (DnLWI), the numerical weather prediction (NWP) models have been developed to obtain the DnLWI indirectly. In this study, a long-term high time-resolution (1 min) observational dataset of the DnLWI in China was used to evaluate the radiation scheme in the CMA-MESO model over various underlying surfaces and climate zones.

Periodicity in the Deccan Volcanism Modulated by Plume Perturbations at the Mid‐Mantle Transition Zone

JGR–Solid Earth - Mon, 06/24/2024 - 16:38
Abstract

In peninsular India, the Deccan Traps record massive, continental-scale volcanism in a sequence of magmatic events that corresponds with the timing of mass extinction at the Cretaceous-Paleogene boundary. Although the Deccan volcanism is linked with the Réunion hotspot, the origin of its periodic magmatic pulses is still debated. We developed a numerical model replicating the geodynamic scenario of the African superplume underneath a moving Indian plate to explore the mechanism of magmatic pulse generation during the Deccan volcanism. Our model results revealed a connection between the Réunion hotspot and the African large low shear-wave velocity province (LLSVP), suggesting that the pulses were produced from a thermochemical plume originated in the lower mantle. The ascending plume had stagnation at 660 km due to phase changes in the transition zone, and its head eventually underwent detachment from the tail under the influence of Indian plate movement to produce sequentially four major pulses (periodicity: 5–8 Ma), each giving rise to multiple secondary magmatic pulses at a time interval of ∼0.15–0.4 Ma. This study sheds a new light on the mechanism of periodic hotspot activities from a global perspective.

Revealing halos concealed by cirrus clouds

Atmos. Meas. techniques - Mon, 06/24/2024 - 15:57
Revealing halos concealed by cirrus clouds
Yuji Ayatsuka
Atmos. Meas. Tech., 17, 3739–3750, https://doi.org/10.5194/amt-17-3739-2024, 2024
Many types of halos appear in the sky. Each type of halo reflects the state of the atmosphere; therefore observing them from the ground greatly helps in understanding the state of the atmosphere. However, halos are easily obscured by the contrast of the cloud itself, making it difficult to observe them. This study describes the construction of a sky-color model for halos and a new effective algorithm to reveal halos in images.

Algorithm to retrieve aerosol optical properties using lidar measurements on board the EarthCARE satellite

Atmos. Meas. techniques - Mon, 06/24/2024 - 15:57
Algorithm to retrieve aerosol optical properties using lidar measurements on board the EarthCARE satellite
Tomoaki Nishizawa, Rei Kudo, Eiji Oikawa, Akiko Higurashi, Yoshitaka Jin, Nobuo Sugimoto, Kaori Sato, and Hajime Okamoto
Atmos. Meas. Tech. Discuss., https://doi.org/10.5194/amt-2024-100,2024
Preprint under review for AMT (discussion: open, 0 comments)
We developed algorithms to produce JAXA ATLID L2 aerosol products using ATLID L1 data. The algorithms estimate layer identifiers such as aerosol or cloud layers, (2) particle optical properties at 355 nm, (3) particle type identifiers, and (4) planetary boundary layer height. We demonstrated the algorithm performance using the simulated ATLID L1 data and found the algorithm’s capability to provide valuable insights into the global distribution of aerosols and clouds.

Hyper-resolution flood hazard mapping at the national scale

Natural Hazards and Earth System Sciences - Mon, 06/24/2024 - 12:41
Hyper-resolution flood hazard mapping at the national scale
Günter Blöschl, Andreas Buttinger-Kreuzhuber, Daniel Cornel, Julia Eisl, Michael Hofer, Markus Hollaus, Zsolt Horváth, Jürgen Komma, Artem Konev, Juraj Parajka, Norbert Pfeifer, Andreas Reithofer, José Salinas, Peter Valent, Roman Výleta, Jürgen Waser, Michael H. Wimmer, and Heinz Stiefelmeyer
Nat. Hazards Earth Syst. Sci., 24, 2071–2091, https://doi.org/10.5194/nhess-24-2071-2024, 2024
A methodology of regional flood hazard mapping is proposed, based on data in Austria, which combines automatic methods with manual interventions to maximise efficiency and to obtain estimation accuracy similar to that of local studies. Flood discharge records from 781 stations are used to estimate flood hazard patterns of a given return period at a resolution of 2 m over a total stream length of 38 000 km. The hazard maps are used for civil protection, risk awareness and insurance purposes.

An Empirical Model of the Occurrence Rate of Low Latitude Post‐Sunset Plasma Irregularities Derived From CHAMP and Swarm Magnetic Observations

Space Weather - Mon, 06/24/2024 - 07:00
Abstract

The prediction of post-sunset equatorial plasma depletions (EPDs), often called ionospheric plasma bubbles, has remained a challenge for decades. In this study, we introduce the Ionospheric Bubble Probability (IBP) model, an empirical model to predict the occurrence probability of EPDs derived from 9 years of CHAMP and 9 years of Swarm magnetic field measurements. The model predicts the occurrence probability of EPDs for a given longitude, day of year, local time and solar activity, for the altitude range of about 350–510 km, and low geographic latitudes of ±45°. IBP has been found to successfully reconstruct the distribution of EPDs as reported in previous studies from independent data. IBP has been further evaluated using 1-year of untrained data of the Ionospheric Bubble Index (IBI). IBI is a Level 2 product of the Swarm satellite mission used for EPD identification. The relative operating characteristics (ROC) curve shows positive excursion above the no-skill line with Hanssen and Kuiper's Discriminant (H&KSS) score of 0.52, 0.51, and 0.55 at threshold model output of 0.16 for Swarm A, B, and C satellites. Additionally, the reliability plots show proximity to the diagonal line with a decent Brier Skill Score (BSS) of 0.249, 0.210, and 0.267 for Swarm A, B, and C respectively at 15% climatological occurrence rate. These tests indicate that the model performs significantly better than a no-skill forecast. The IBP model offers compelling glimpses into the future of EPD forecasting, thus demonstrating its potential to reliably predict EPD occurrences. The IBP model is publicly available.

Evolution of Energetic Proton Parallel Pressure Anisotropy at Geosynchronous Altitudes: Potential Role in Triggering Substorm Expansion Phase Onset

GRL - Mon, 06/24/2024 - 06:05
Abstract

The sequence of events associated with the triggering of energy release during substorm expansion phase onset is still not well-understood. Oberhagemann and Mann (2020b, https://doi.org/10.1029/2019gl085271) proposed a new substorm onset mechanism, where the transition toward parallel proton pressure anisotropy during tail stretching in the late growth phase could trigger a pressure anisotropic ballooning instability. Here we examine the evolution of energetic proton parallel pressure anisotropy at geosynchronous altitudes, seeking evidence in support of the proposed substorm onset mechanism. We use the Geostationary Operational Environment Satellite (GOES) proton flux and magnetometer data combined with substorm onset indicators derived from ground-based magnetometers. Superposed epoch analysis of substorm onset times for 2014 using the isolated substorm list (Ohtani & Gjerloev, 2020, https://doi.org/10.1029/2020ja027902) clearly shows signatures of energetic proton parallel pressure anisotropy immediately before substorm onset, potentially supportive of the Oberhagemann and Mann theory.

First Observations of Severe Scintillation Over Low‐to‐Mid Latitudes Driven by Quiet‐Time Extreme Equatorial Plasma Bubbles: Conjugate Measurements Enabled by Citizen Science Initiatives

JGR:Space physics - Mon, 06/24/2024 - 04:59
Abstract

Low-cost instrumentation combined with volunteering and citizen science educational initiatives allowed the deployment of L-band scintillation monitors to remote sense areas that are geomagnetically conjugated and located at low-to-mid latitudes in the American sector (Quebradillas in Puerto Rico and Santa Maria in Brazil). On 10 and 11 October, 2023, both monitors detected severe scintillations, some reaching dip latitudes beyond 26°N. The observations show conjugacy in the spatio-temporal evolution of the scintillation-causing irregularities. With the aid of collocated all-sky airglow imager observations, it was shown that the observed scintillation event was caused by extreme equatorial plasma bubbles (EPBs) reaching geomagnetic apex altitudes exceeding 2,200 km. The observations suggest that geomagnetic conjugate large-scale structures produced conditions for the development of intermediate scale (few 100 s of meters) in both hemispheres, leading to scintillation at conjugate locations. Finally, unlike previous reports, it is shown that the extreme EPBs-driven scintillation reported here developed under geomagnetically quiet conditions.

The Dynamics of the India‐Eurasia Collision: Faulted Viscous Continuum Models Constrained by High‐Resolution Sentinel‐1 InSAR and GNSS Velocities

JGR–Solid Earth - Mon, 06/24/2024 - 03:41
Abstract

The distribution and magnitude of forces driving lithospheric deformation in the India-Eurasia collision zone have been debated over many decades. Here we test a two-dimensional (2-D) Thin Viscous Shell approach that has been adapted to explicitly account for displacement on major faults and investigate the impact of lateral variations in depth-averaged lithospheric strength. We present a suite of dynamic models to explain the key features from new high-resolution Sentinel-1 Interferometric Synthetic Aperture Radar as well as Global Navigation Satellite System velocities. Comparisons between calculated and geodetically observed velocity and strain rate fields indicate: (a) internal buoyancy forces from Gravitational Potential Energy acting on a relatively weak region of highest topography (>2,000 m) contribute to dilatation of the high plateau and contraction on the margins; (b) a weak central Tibetan Plateau (∼1021 Pa s compared to far-field depth-averaged effective viscosity of at least 1022–1023 Pa s) is required to explain the observed long-wavelength eastward velocity variation; (c) localized displacement on fault systems enables strain concentration and clockwise rotation around the Eastern Himalayan Syntaxis. We discuss the tectonic implications for rheology of the lithosphere, distribution of geodetic strain, and partitioning of active faulting and seismicity.

Toward Enhanced Prediction of High‐Impact Solar Energetic Particle Events Using Multimodal Time Series Data Fusion Models

Space Weather - Sun, 06/23/2024 - 07:00
Abstract

Solar energetic particle (SEP) events, originating from solar flares and Coronal Mass Ejections, present significant hazards to space exploration and technology on Earth. Accurate prediction of these high-energy events is essential for safeguarding astronauts, spacecraft, and electronic systems. In this study, we conduct an in-depth investigation into the application of multimodal data fusion techniques for the prediction of high-energy SEP events, particularly ∼100 MeV events. Our research utilizes six machine learning (ML) models, each finely tuned for time series analysis, including Univariate Time Series (UTS), Image-based model (Image), Univariate Feature Concatenation (UFC), Univariate Deep Concatenation (UDC), Univariate Deep Merge (UDM), and Univariate Score Concatenation (USC). By combining time series proton flux data with solar X-ray images, we exploit complementary insights into the underlying solar phenomena responsible for SEP events. Rigorous evaluation metrics, including accuracy, F1-score, and other established measures, are applied, along with K-fold cross-validation, to ensure the robustness and generalization of our models. Additionally, we explore the influence of observation window sizes on classification accuracy.

Dependence of Cloud Macrophysical Properties and Phase Distributions on Environmental Conditions Over the North Atlantic and Southern Ocean: Results From COMBLE and MARCUS

JGR–Atmospheres - Sat, 06/22/2024 - 20:34
Abstract

The accurate representation of Cold Air Outbreaks (CAOs) and affiliated mixed-phase boundary layer (BL) clouds in models is challenging. How BL cloud properties evolve during CAOs and their dependence on meteorological conditions is not well understood but is important for the simulation of Earth's energy budgets. Here the properties of polar BL clouds over the North Atlantic (NA) and Southern Ocean (SO) are compared using observations from the Measurements of Aerosol Radiation and CloUds over the SO (MARCUS) and CAOs in the Marine BL Experiment (COMBLE) conducted over the NA. MARCUS observations show a stronger BL inversion than COMBLE, with a higher mean EIS (estimated inversion strength)/LTS (lower tropospheric stability) of −0.03 K/13 K compared to COMBLE’s −3.2 K/9.3 K. 39% of CAOs observed during COMBLE were intense with M > 5 K, while MARCUS only had 1.3%. 78%/72% of clouds sampled in CAOs during COMBLE/MARCUS had cloud top heights <4 km. The mean BL cloud top height was over 400 m higher, and the BL was over 500 m deeper for M of 10 K compared to 0 K for both regions. MARCUS observed a 27% moister BL structure than COMBLE when M > 5 K due to stronger BL inversion trapping more moisture within the BL. Under the same LTS, EIS, and M conditions, MARCUS observed a 12% drier BL structure, and clouds were 46% more turbulent than COMBLE. During CAOs, 54% of single-layer BL clouds sampled during MARCUS had liquid-dominated bases compared to 39% during COMBLE.

Aitken Mode Aerosols Buffer Decoupled Mid‐Latitude Boundary Layer Clouds Against Precipitation Depletion

JGR–Atmospheres - Sat, 06/22/2024 - 20:18
Abstract

Aerosol-cloud-precipitation interactions are a leading source of uncertainty in estimating climate sensitivity. Remote marine boundary layers where accumulation mode (∼100–400 nm diameter) aerosol concentrations are relatively low are very susceptible to aerosol changes. These regions also experience heightened Aitken mode aerosol (∼10–100 nm) concentrations associated with ocean biology. Aitken aerosols may significantly influence cloud properties and evolution by replenishing cloud condensation nuclei and droplet number lost through precipitation (i.e., Aitken buffering). We use a large-eddy simulation with an Aitken-mode enabled microphysics scheme to examine the role of Aitken buffering in a mid-latitude decoupled boundary layer cloud regime observed on 15 July 2017 during the Aerosol and Cloud Experiments in the Eastern North Atlantic flight campaign: cumulus rising into stratocumulus under elevated Aitken concentrations (∼100–200 mg−1). In situ measurements are used to constrain and evaluate this case study. Our simulation accurately captures observed aerosol-cloud-precipitation interactions and reveals time-evolving processes driving regime development and evolution. Aitken activation into the accumulation mode in the cumulus layer provides a reservoir for turbulence and convection to carry accumulation aerosols into the drizzling stratocumulus layer above. Further Aitken activation occurs aloft in the stratocumulus layer. Together, these activation events buffer this cloud regime against precipitation removal, reducing cloud break-up and associated increases in heterogeneity. We examine cloud evolution sensitivity to initial aerosol conditions. With halved accumulation number, Aitken aerosols restore accumulation concentrations, maintain droplet number similar to original values, and prevent cloud break-up. Without Aitken aerosols, precipitation-driven cloud break-up occurs rapidly. In this regime, Aitken buffering sustains brighter, more homogeneous clouds for longer.

Contributions of Anthropogenic Greenhouse Gases and Aerosol Emissions to Changes in Summer Precipitation Over Southern China

JGR–Atmospheres - Sat, 06/22/2024 - 20:04
Abstract

This study investigates the anthropogenic contribution to the summer precipitation changes over southern China and the underlying physical mechanisms. Observations show a wetting trend over southeastern China (SEC) but a drying trend over southwestern China (SWC) in summer of 1961–2014. The dipole pattern can be reasonably reproduced by the anthropogenic forcing simulations of CMIP6 models but with weak trends under the external natural forcing simulations, suggesting the vital human contribution to the observed changes. Particularly, anthropogenic greenhouse gases (GHG) dominate the wetting trend over SEC, while the drying trend over SWC is primarily attributed to anthropogenic aerosol (AA) emissions. Further analysis shows that the GHG concentrations enhance the subtropical high over the western North Pacific (WNP) via the heterogeneous warming of the sea surface temperature, decrease the sea level pressure over eastern China, and increase the atmospheric moisture, facilitating the moisture flux convergence (MFC) and the precipitation over SEC. The GHG-induced wetting trend is somewhat offset by the inhibited evaporation due to the AA forcing. For SWC, the decreased precipitation is influenced by the anomalous high pressure from India to WNP, which is closely associated with the enhanced Asian AA emission and the interhemispheric asymmetrical distribution of AA emissions. In the upper troposphere, the uneven AA emissions between South and East Asia and Europe weaken the East Asian summer subtropical jet, resisting the western moisture to SWC. Both factors in the low-and-high levels suppress the MFC and precipitation over SWC, counteracted by the thermodynamical effects of GHG forcing.

Dynamic Variations of Agricultural Drought and Its Response to Meteorological Drought: A Drought Event‐Based Perspective

JGR–Atmospheres - Sat, 06/22/2024 - 19:48
Abstract

The dynamic variations of agricultural drought can reflect the water shortage situation of the agricultural system, and there is a progressive relationship in the response of agricultural drought to meteorological drought on a spatiotemporal scale. In this study, the vegetation health index and the standardized precipitation evapotranspiration index were adopted as agricultural and meteorological drought indicators, respectively. Additionally, using the three-dimensional spatiotemporal clustering technology, the dynamic evolutions of typical drought events were clarified, and the spatiotemporal response characteristics of agricultural drought to meteorological drought were revealed. The results indicated that: (a) there were 81 agricultural drought events in the North China Plain (NCP) during 1982–2020, with a largest drought severity (12.82 × 104 month km2), a 6-month duration, and a 23.24 × 104 km2 affected area occurring in the No. 4 event; (b) from the 1980s to the 2010s, the agricultural drought gradually decreased and large-scale droughts mainly concentrated in the border areas of Hebei, Shandong, and Henan; and (c) a total of 13 drought event pairs were successfully matched in the NCP, including 7 pairs of “one-to-one,” 4 pairs of “one-to-many,” 1 pair of “many-to-one,” and 1 pair of “many-to-many.” The spatiotemporal responses of agricultural drought were elucidated in a three-dimensional perspective, which can propose a new approach for establishing drought propagation model, predicting future agricultural drought conditions, improving ecological environment quality, and can also be applied for the investigation of other drought types.

An Improved Non‐Local Planetary Boundary Layer Parameterization Scheme in Weather Forecasting and Research Model Based on a 1.5‐Order Turbulence Closure Model

JGR–Atmospheres - Sat, 06/22/2024 - 18:54
Abstract

Planetary boundary layer (PBL) modeling is a primary contributor to uncertainties in a numerical weather prediction (NWP) model due to difficulties in modeling the turbulent transport of surface fluxes. The Weather Research and Forecasting model (WRF) has provided many PBL schemes that may feature a non-local transport component driven by large eddies or a one-and-half order turbulence closure model, but few of them possess the two features at once. In the present study, a turbulent kinetic energy (TKE)-based eddy diffusivity/viscosity method is integrated into the non-local Asymmetric Convective Model version 2 (ACM2) PBL scheme and implemented in WRF. The original first-order eddy-diffusivity term in ACM2 is discarded and an extra prognostic equation for TKE, which considers the tendency of TKE by buoyancy, wind shear, vertical transport, and dissipative processes, is supplied to calculate the diffusivity/viscosity. Non-local transport is modeled the same as ACM2 using the transilient matrix method. Idealized tests using prescribed surface heat flux and roughness length are performed. TKE-ACM2 displays advantages over the PBL scheme developed by Bougeault and Lacarrère (hereinafter referred to as Boulac) and ACM2 in the wind speeds (WS) profile because it better matches large-eddy simulations results in the surface momentum flux. Real case simulations show that TKE-ACM2 generally outperforms in the diurnal vertical profiles of WS under stable conditions. TKE-ACM2 also produces a better alignment under moderately unstable conditions in the early nighttime at the urban LiDAR station. However, the model exhibits discrepancies more apparently under a more unstable condition during the winter daytime.

Partitioning of Iron Between (Mg,Fe)SiO3 Liquid and Bridgmanite

GRL - Sat, 06/22/2024 - 17:54
Abstract

The evolution of the magma ocean that occupied the early Earth is influenced by the buoyancy of crystals in silicate liquid. At lower mantle pressures, silicate crystals are denser than the iso-chemical liquid, but heavy elements like iron can cause crystals to float if they partition into the liquid phase. Crystal flotation allows for a basal magma ocean, which might explain geochemical anomalies in mantle-derived magmas, seismic anomalies in the lower mantle, and the source of the Earth's early magnetic field. To examine whether a basal magma ocean is gravitationally stable, we investigate the degree of iron partitioning between (Mg,Fe)SiO3 liquid and bridgmanite. By utilizing ab initio molecular dynamics simulations coupled with thermodynamic integration, we find that iron partitions into the liquid, and increasingly so with increasing pressure. Bridgmanite crystals are found to be buoyant at lower mantle conditions, stabilizing the basal magma ocean.

Elevation Regulates the Response of Climate Heterogeneity to Climate Change

GRL - Sat, 06/22/2024 - 17:44
Abstract

Climate change represents a profound threat to the diversity and stability of global climate zones. However, the complex interplay between climate change and elevation in shaping climate heterogeneity is not yet fully understood. Here, we combine Shannon's diversity index (SHDI) with the Köppen-Geiger climate classification to explore the altitudinal distributions of global climate heterogeneity; and their responses to climate change. The study reveals a distinctive pattern: SHDI, a proxy for climate heterogeneity tends to slow down or decline at lower elevations with increasing temperatures, while at higher elevations, it continues to rise due to continuing cold conditions. Examination of climate simulations, both with and without anthropogenic forcing, confirms that observed changes in climate heterogeneity are primarily attributable to anthropogenic climate change within these high-elevation regions. This study underscores the importance of high-elevation regions as not only custodians of diverse climate types but also potential refuges for species fleeing warmer climates.

Theme by Danetsoft and Danang Probo Sayekti inspired by Maksimer