Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Field Events”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

Big Data Synchrophasor Monitoring and Analytics for Resiliency Tracking (BDSMART)

This report contains key findings from a project titled Big Data Synchrophasor Monitoring and Analytics for Resiliency Tracking (BDSMART), which was carried out through a collaborative effort of a team of researchers from Texas A&M Engineering Experiment Station, Temple University, and Quanta Technology, LLC. The in-kind support came from OSIsoft (acquired by AVEVA), which provided their PI Historian software to demonstrate the use case of streaming PMU data. The first section of the report describes the project goals and objectives related to the development of Machine Learning (ML) models capable of detecting and classifying events by processing phasor measurements captured in the field by Phasor Measurement Units (PMUs). The data for this study was contributed by the utilities/ISOs from the Western and Eastern interconnects and ERCOT, further referred to as Interconnect B (IC B), Interconnect A (IC A), and Interconnect C (IC C), respectively. The approach that the BDSMART Research Team proposed and the key research tasks defined by the team are outlined in this section. The next section describes the technical approach. We first discuss the data constraints related to the PMU measurements and data interpretation constraints imposed by the data contributors. They provided neither the topological information of the grid nor PMU placement locations and captured recorded data at very few locations in the system with the reporting rate of either 30 or 60 fps. The recordings are mostly positive sequence voltage, frequency, and ROCOF, and in some limited cases, three-phase voltages and currents. We then reflect on the bad data issues that stem from poor recording practices and vague definitions of the PMU status bits to supposedly be used for bad data identification. Finally, the data discovery points to imprecise time stamps with incomplete event start/end time, as well as inconsistent and incomplete event labeling, which combined make the implementation of the data models using supervising learning quite challenging. Following the data discovery study, we hypothesize that because the IC B data has the most complete labels, we should focus our model development on that data and then test it on data from other interconnects. We also define the common metrics used to evaluate the results from the ML algorithm tests. We concluded this section by summarizing the common ML models we used and explaining how we implemented and tested them. The issues from this section are expanded in the Training Dataset Report from this project. The final section of this report deals with the accomplishments and conclusions. As the accomplishments, we formulate the problem we are solving and what is achieved by solving the problem. We then reflect on each of the analytics tools we developed and point out the performance of each tool when applied to solving the mentioned problems. We reference this work for further details to the papers we published on each tool. In the conclusions, we give recommendations on how to improve future PMU recording practices to facilitate the ML algorithm implementation and guidance for the future standardization work aimed at clarifying the ambiguities associated with the PMU status bits. We finally list future tasks that can bring about further improvements in the proposed algorithms. The issues from this section are expanded in the Training, and Test Dataset Report filed at the project completion date.

97 MATHEMATICS AND COMPUTING↗

Jet cone radius dependence of R AA and v 2 at PbPb 5.02 TeV from JEWEL+T R ENTo+v-USPhydro

We combine, for the first time, event-by-event T R ENTo initial conditions with the relativistic viscous hydrodynamic model v-USPhydro and the Monte Carlo event generator JEWEL to make predictions for the nuclear modification factor R AA and jet azimuthal anisotropies v n { 2 } in $\sqrt{s_{NN}}$ = 5.02 TeV PbPb collisions for multiple centralities and values of the jet cone radius R. The R-dependence of R AA and v 2 { 2 } strongly depends on the presence of recoiling scattering centers. We find a small jet v 3 { 2 } in mid-central collisions and consistent results in wide jet p T regions and centralities with ATLAS data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

PopSyCLE: A New Population Synthesis Code for Compact Object Microlensing Events

We present a new Milky Way microlensing simulation code, dubbed PopSyCLE (Population Synthesis for Compact object Lensing Events). PopSyCLE is the first resolved microlensing simulation to include a compact object distribution derived from numerical supernovae explosion models and both astrometric and photometric microlensing effects. We demonstrate the capabilities of PopSyCLE by investigating the optimal way to find black holes (BHs) with microlensing. Candidate BHs have typically been selected from wide-field photometric microlensing surveys, such as OGLE, by selecting events with long Einstein crossing times (tE > 120 days). These events can be selected at closest approach and monitored astrometrically in order to constrain the mass of each lens; PopSyCLE predicts a BH detection rate of ~40% for such a program. We find that the detection rate can be enhanced to ~ 85% by selecting events with both tE > 120 days and a microlensing parallax of πE < 0.08. Unfortunately, such a selection criterion cannot be applied during the event as πE requires both pre- and post-peak photometry. However, historical microlensing events from photometric surveys can be revisited using this new selection criteria in order to statistically constrain the abundance of BHs in the Milky Way. The future WFIRST microlensing survey provides both precise photometry and astrometry and will yield individual masses of O(100 - 1000) black holes, which is at least an order of magnitude more than is possible with individual candidate follow-up with current facilities. The resulting sample of BH masses from WFIRST will begin to constrain the shape of the black hole present-day mass function, BH multiplicity, and BH kick velocity distributions

79 ASTRONOMY AND ASTROPHYSICS↗

Evolution of the ATLAS event data model for the HL-LHC

The upcoming high-luminosity run of the CERN Large Hadron Collider (HL-LHC) will yield an unprecedented volume of data. In order to process this data, the ATLAS collaboration is evolving its offline software to be able to use heterogeneous resources such as graphical processing units (GPUs) and field-programmable gate arrays (FPGAs). To reduce conversion overheads, the event data model (EDM) should be compatible with the requirements of these resources. While the ATLAS EDM has long allowed representing data as a structure of arrays, further evolution of the EDM can enable more efficient sharing of data between CPU and GPU resources. Some of this work will be summarized here, including extensions to allow controlling how memory for event data is allocated and the implementation of jagged vectors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Predicting cutoff L-shells of solar protons using the GPPSn particle dataset

Solar energetic protons (SEPs) arriving at the Earth trigger severe radiation storms in the near-Earth space, directly impacting space missions operating at various altitudes. Therefore, monitoring SEP events and predicting the penetration depths of solar protons are critical for aerospace sectors. Building on previous efforts, here we demonstrate the feasibility of using proton measurements from the Global Prompt Proton Sensor network (GPPSn), enabled by Los Alamos National Laboratory developed combined X-ray dosimeters aboard GPS satellites, to characterize and predict the penetration of solar protons into the geomagnetic field. The inclined medium-Earth-orbits (MEOs) of the global GPS constellation offer a unique advantage of allowing simultaneous measurements of penetrating solar protons inside both open- and closed-field line regions. Therefore, the L-profiles of ∼10s–100 MeV solar protons and their associated cutoff L-shells can be determined from the GPPSn dataset, using predefined threshold proton flux values rather than traditional flux ratios. After examining a list of SEP event intervals across solar cycles 23, 24 and 25—including the 2024 Mother’s Day superstorm, we showcase how the latest GPPSn proton dataset (release v1.10), reprocessed and calibrated, can not only be used to monitor solar proton distributions inside the dynamic geomagnetic field for individual events, but also to derive a new empirical model linking cutoff L-shells with several key space weather parameters. This newly developed SEPCL-MEO model demonstrates high predictive performance; for example, predictions for > 30 MeV solar protons yield a correlation coefficient of 0.85 and performance efficiency of 0.67 when validated against GPPSn observations. Results from this pilot study underscores the scientific and operational value of the GPPSn dataset, and this dataset—when paired with machine-learning techniques—can play a critical role in observing and predicting the effects of future incoming SEP events, including extreme ones.

58 GEOSCIENCES↗

Performance of the ATLAS Level-1 topological trigger in Run 2

During LHC Run 2 (2015–2018) the ATLAS Level-1 topological trigger allowed efficient data-taking by the ATLAS experiment at luminosities up to 2.1x10 34 cm -2 s -1 , which exceeds the design value by a factor of two. The system was installed in 2016 and operated in 2017 and 2018. It uses Field Programmable Gate Array processors to select interesting events by placing kinematic and angular requirements on electromagnetic clusters, jets, τ-leptons, muons and the missing transverse energy. It allowed to significantly improve the background event rejection and signal event acceptance, in particular for Higgs and B-physics processes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

“Godzilla,” the Extreme African Dust Event of June 2020: Origins, Transport, and Impact on Air Quality in the Greater Caribbean Basin

In June 2020, the tropical Atlantic and the Caribbean Basin were affected by a series of African dust outbreaks unprecedented in size and intensity. These events, informally named “Godzilla,” coincided with CALIMA, a large field campaign, offering a rare opportunity to assess the impact of African dust on air quality in the Greater Caribbean Basin. Network measurements of respirable particles (i.e., PM 10 and PM 2.5 ) showed that dust significantly degraded regional air quality and increased the risk to public health in the Caribbean, the southern United States, northern South America, and Central America. CALIMA examined the meteorological context of Godzilla dust events over North Africa and how these conditions might relate to the greatly increased dust emissions and enhanced transport to the Americas. Godzilla was linked to strong pressure anomalies over West Africa, resulting in a large-scale geostrophic wind anomaly at 700 hPa over North Africa. We used surface-based and columnar measurements to test the performance of two frequently used aerosol forecast models: the NASA Goddard Earth Observing System (GEOS) and Weather Research and Forecasting Model coupled with Chemistry (WRF-Chem) models. The models showed some skills but differed substantially between their forecasts, suggesting large uncertainties in these forecasts that are critical for issuing early warnings of health-threatening dust events. Our results demonstrate the value of an integrated approach in characterizing the spatial and temporal variability of African dust transport and assessing its impact on regional air quality. Future studies are needed to improve models and to track the long-term changes in dust transport from Africa under a changing climate.

Aerosols/particulates↗

A Closed Bay-Breeze Circulation and Its Lifecycle From TRACER With a New Orienteering Tape Recorder Diagram

A closed bay-breeze circulation (BBC) followed by a gulf-breeze front (GBF) was observed on 10 September 2022 during the Department of Energy's (DOE) Atmospheric Radiation Measurement (ARM) TRACER campaign. Using high-resolution X-band Scanning Cloud Radar (XSACR) and a newly developed orienteering tape recorder diagram, the study analyzed radar reflectivity and Doppler velocity to identify anomalies and track the evolution of these circulations. The BBC, a mesoscale system approximately 30-km long, 30-km wide, and 2-km deep, formed from enhanced horizontal convective rolls along the Galveston Bay coast, progressing northwestward under a 6 m s−1 onshore flow anomaly and reaching 1.5 km depth with return flow aloft. The inland penetration speed was 2 m s−1, driven by an observed 6-9 K land-water temperature contrast. The GBF, coupled to the BBC, intensified with additional southerly flow, penetrating further inland after the BBC collapsed. Passing over the TRACER field site, both fronts significantly impacted boundary layer thermodynamics, dynamics, and aerosol concentration. The BBC event exhibited four lifecycle stages-formation, development, maturation, and dissipation-driven by solar heating, wind field rotation, and interactions with convective eddies and the GBF. This study provides insights into the inland evolution of coastal breeze circulations and their interactions with environmental processes.

54 ENVIRONMENTAL SCIENCES↗

Self-Propulsion by Directed Explosive Emulsification

An active droplet system, programmed to repeatedly move autonomously at a specific velocity in a well-defined direction, is demonstrated. Coulombic energy is stored in oversaturated interfacial assemblies of charged nanoparticle-surfactants by an applied DC electric field and can be released on demand. Spontaneous emulsification is suppressed by an increase in the stiffness of the oversaturated assemblies. Rapidly removing the field releases the stored energy in an explosive event that propels the droplet, where thousands of charged microdroplets are ballistically ejected from the surface of the parent droplet. The ejection is made directional by a symmetry breaking of the interfacial assembly, and the combined interaction force of the microdroplet plume on one side of the droplet propels the droplet distances tens of times its size, making the droplet active. The propulsion is autonomous, repeatable, and agnostic to the chemical composition of the nanoparticles. The symmetry-breaking in the nanoparticle assembly controls the microdroplet velocity and direction of propulsion. Finally, this mechanism of droplet propulsion will advance soft micro-robotics, establishes a new type of active matter, and introduces new vehicles for compartmentalized delivery.

36 MATERIALS SCIENCE↗

Data-driven picosecond X-ray imaging for quantitative plasma-induced shock characterization

Imaging dynamic events, especially shockwave behavior, is key to advancing high-energy-density (HED) research. Recent advances in fourth- and fifth-generation X-ray light sources allow for high-resolution imaging of fast phenomena, but limited beam time necessitates maximizing data acquisition. We present a benchtop-scale pulsed plasma device submerged in liquid heptane, capable of generating dynamic events at rates exceeding 10 Hz, supporting the field’s data-driven goals by producing large, high-quality imaging datasets. Using X-ray phase contrast imaging (XPCI) at the Advanced Photon Source, we imaged weak shockwaves (Mach ~ 1.2) in heptane interacting with plasma-induced cavitation bubbles, causing deviation from Rankine-Hugoniot behavior; to our knowledge, this represents the first direct imaging of such interaction. Our quantitative analysis offers insight into weak shock phenomena and energy-focusing applications in pulsed plasmas. These results highlight the potential for large datasets to advance dynamic HED research at current light source facilities, and have implications for fields such as inertial confinement fusion, plasma-enhanced chemical processing, and biomedical applications.

36 MATERIALS SCIENCE↗

Fail-Safe Logic Design Strategies Within Modern FPGA Architectures

Fail-safe computing refers to computing systems that revert to a non-operational safe state when a fault occurs. In this paper, we investigate a circuit level technique as mitigation for single event upsets (SEUs) and fault injection attacks on field programmable gate arrays (FPGAs), and analyze the effectiveness of the technique as a fail-safe monitor for an encryption algorithm. The propagation of fault effects through FPGA primitives including lookup tables (LUTs) and programmable interconnect points (PIPs) is assessed within an FPGA architecture created using an open source tool, and validated using fault injection experiments on an FPGA. The analysis reveals additional vulnerabilities exist within reconfigurable architectures over those in equivalent fail-safe application specific integrated circuit (ASIC), thus requiring a more elaborate network of redundant circuits and checking logic. The configuration memory bits (CMBs), which configure routing and designate logic functions within the LUTs of the FPGA, add complexity to fail-safe design strategies by introducing additional fault conditions and fault propagation paths. A resource-efficient fail-safe circuit design technique called DEsign for Fail-safe in reCONfigurable systems (DEFCON) is proposed. The benefits and limitations associated with DEFCON are described in the context of fault injection experiments carried out as simulations and in FPGA hardware.

Bhakta, Priya A. [Univ. of New Mexico, Albuquerque↗

Optimization of Carbon-Nanotube-Reinforced Composite Gridlines Towards Commercialization

Solar cells in PV modules crack during field operation from environmental stressors, including extreme weather events, such as hailstorms and hurricanes. These cracks can lead to gradual or immediate acute power degradation. To directly address cell-crack-induced degradation, we have formulated a carbon nanotube additive for commercial screen printed silver pastes. We have shown in previous work that these metal matrix composites have little to no effect on the cell’s efficiency while enhancing the metallization’s fracture toughness and electrical gap-bridging capability. In this work, we focus on optimizing the composite metallization to achieve the best possible performance. We discover that reducing carbon nanotube agglomerations in the paste not only improves the printability for narrow gridlines, but also increases the modulus of toughness of the metallization by over 500%.

carbon nanotubes↗

Beryllium Isotopes in Rainfall for Stratospheric-Tropospheric Exchange Dynamics Experiment (BIRDIE) (Final Report)

The Beryllium Isotopes in Rainfall for Stratospheric-Tropospheric Exchange Dynamics Experiment or BIRDIE was a precipitation collection campaign designed to measure cosmogenic isotopes of beryllium in rainfall and other forms of precipitation across the DOE Southern Great Plains (SGP) Atmospheric Radiation Measurement (ARM) site. Precipitation samples were collected over three weeks (May 21-June 10, 2022) at 11 SGP sites across a 120 km 2 domain in northern and central Oklahoma. The field campaign targeted summer-time convection rain events. We aimed to capture any spatial and/or temporal variability in the concentrations of beryllium-10 (10Be) and beryllium-7 (7Be) across the SGP ARM site during these rainfall events. The ratio of 10Be to 7Be (10/7Be) has been shown to be a powerful tracer for transport processes in the stratosphere-troposphere exchange (STE) system (e.g., Koch and Rind 1998). Any observed variability is analyzed to see if it could be linked to STE processes, including overshooting of deep convection clouds which penetrate the lower stratosphere. Six wet weather system events were collected within our three-week span. Five out of six of these rainfall events are classified as moderately to strongly convective. Those events are studied in more detail. The five systems also follow each other over a short time period of eleven days. In our rainfall samples we observed spatial and event variability in the concentrations of beryllium-10 and -7 which resulted in variability of 10/7 Be ratios. Each weather event appeared to exhibit a unique 10/7Be pattern across the domain that was not directly linked to storm total precipitation amount. To better understand the causes for this spatial variability, cloud characteristics in the domain were looked at to assess cloud type, cloud height and depth information. We also ran Hysplit back trajectories to identify possible air mass and moisture origins across our events. Detailed analysis is still on-going. The purpose of this technical report is to provide background on using cosmogenic isotopes as atmospheric tracers, describe the experimental design, and present the beryllium concentration and ratio data.

58 GEOSCIENCES↗

Regional inertia dynamics of U.S. interconnections: An event-based measurement approach

Power grid inertia plays a vital role in frequency stability following large disturbances, yet its distribution across the U.S. grid is highly uneven. While interconnection-wide inertia benchmarks are useful, they can mask regional variability driven by resource mix, network coupling, and geographic separation. This paper extends event-driven inertia estimation to the regional scale using field measurements from the Frequency Monitoring Network (FNET/GridEye). Starting from balancing authority and independent system operator footprints, candidate regions are refined using a composite coherency score that combines frequency-trajectory shape similarity, timing spread, and lead/lag behavior to ensure dynamic consistency. A filtered sliding difference method (FSDM) is then used to construct regional frequency trajectories, detect disturbance onset, and compute robust regional rate-of-change of frequency (RoCoF). Regional, local, and interconnection inertia are estimated by combining RoCoF with event power imbalance, and additional indicators (regional-to-system inertia ratio and inertial-support arrival time) quantify regional-to-interconnection coupling and relative regional contributions. The method is demonstrated on eleven regions across the Eastern Interconnection (EI) and the Western Electricity Coordinating Council (WECC), with the Electric Reliability Council of Texas (ERCOT) used for validation. In ERCOT, estimates compared against energy management system (EMS) values achieve a mean absolute percentage error of 17.94%. WECC exhibits consistently shorter inertial-support arrival times (0.15–0.3 s) than EI (0.7–1.1 s), highlighting contrasting coupling and disturbance-propagation behavior. Overall, the results reveal pronounced spatial heterogeneity in inertia and coupling, underscoring the value of regional monitoring for both operational decision-making and long-term system planning.

Disturbance events↗

Understanding and improving model representation of aerosol optical properties for a Chinese haze event measured during KORUS-AQ

KORUS-AQ was an international cooperative air quality field study in South Korea that measured local and remote sources of air pollution affecting the Korean Peninsula during May–June 2016. Some of the largest aerosol mass concentrations were measured during a Chinese haze transport event (24 May). Air quality forecasts using the WRF-Chem model with aerosol optical depth (AOD) data assimilation captured AOD during this pollution episode but overpredicted surface particulate matter concentrations in South Korea, especially PM 2.5 , often by a factor of 2 or larger. Analysis revealed multiple sources of model deficiency related to the calculation of optical properties from aerosol mass that explain these discrepancies. Using in situ observations of aerosol size and composition as inputs to the optical properties calculations showed that using a low-resolution size bin representation (four bins) underestimates the efficiency with which aerosols scatter and absorb light (mass extinction efficiency). Besides using finer-resolution size bins (8–16 bins), it was also necessary to increase the refractive indices and hygroscopicity of select aerosol species within the range of values reported in the literature to achieve better consistency with measured values of the mass extinction efficiency (6.7 m 2 g –1 observed average) and light-scattering enhancement factor (ƒ(RH)) due to aerosol hygroscopic growth (2.2 observed average). Furthermore, an evaluation of the optical properties obtained using modeled aerosol properties revealed the inability of sectional and modal aerosol representations in WRF-Chem to properly reproduce the observed size distribution, with the models displaying a much wider accumulation mode. Other model deficiencies included an underestimate of organic aerosol density (1.0 g cm –3 in the model vs. observed average of 1.5 g cm –3 ) and an overprediction of the fractional contribution of submicron inorganic aerosols other than sulfate, ammonium, nitrate, chloride, and sodium corresponding to mostly dust (17 %–28 % modeled vs. 12 % estimated from observations). These results illustrate the complexity of achieving an accurate model representation of optical properties and provide potential solutions that are relevant to multiple disciplines and applications such as air quality forecasts, health impact assessments, climate projections, solar power forecasts, and aerosol data assimilation.

54 ENVIRONMENTAL SCIENCES↗

A Frequency-Domain-Based Algorithm for Detecting Microseismicity Using Dense Surface Seismic Arrays

We propose a new frequency-domain-based algorithm for detecting small-magnitude seismic events using dense surface seismic arrays. Our proposed method takes advantage of the high energy carried by S waves, and approximate known source locations, which are used to rotate the horizontal components to obtain the maximum amplitude. By surrounding the known source area with surface geophones, we achieve a favorable geometry for locating the detected seismic events with the backprojection method. To test our new detection method, we used a dense circular array, consisting of 151 5 Hz three-component geophones, over a 5 km aperture that was in operation at the Utah Frontier Observatory for Research in Geothermal Energy (FORGE) in southcentral Utah. We apply the new detection method during a small-scale test injection phase at FORGE, and during an aftershock sequence of an Mw 4.1 earthquake located ~30 km north of the geophone array, within the Black Rock volcanic field. We are able to detect and locate microseismic events (Mw<0) during injections, despite the high level of anthropogenic activity, and several aftershocks that are missing from the regional catalog. By comparing our method with known algorithms that operate both in the time and frequency domain, we show that our proposed method performs better in the case of the FORGE injection monitoring, and equally well for the off-array aftershock sequence. Our new method has the potential to improve microseismic event detections even in extremely noisy environments, and the proposed location scheme serves as a direct discriminant between true and false detections.

Geochemistry & Geophysics↗

Waveform retrieval for ultrafast applications based on convolutional neural networks

Electric field waveforms of light carry rich information about dynamical events on a broad range of timescales. The insight that can be reached from their analysis, however, depends on the accuracy of retrieval from noisy data. In this article, we present a novel approach for waveform retrieval based on supervised deep learning. We demonstrate the performance of our model by comparison with conventional denoising approaches, including wavelet transform and Wiener filtering. The model leverages the enhanced precision obtained from the nonlinearity of deep learning. The results open a path toward an improved understanding of physical and chemical phenomena in field-resolved spectroscopy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling the role of lateral surface flow in low‐relief polygonal tundra

Abstract Ice‐wedge polygon troughs play an important role in controlling the hydrology of low‐relief polygonal tundra regions. Lateral surface flow is confined to troughs only, but it is often neglected in model projections of permafrost thermal hydrology. Recent field and modeling studies have shown that, after rain events, increases in trough water levels are significantly more than the observed precipitation, highlighting the role of lateral surface flow in the polygonal tundra. Therefore, understanding how trough lateral surface flow can influence polygonal tundra thermal hydrology is important, especially under projected changes in temperatures and rainfall in the Arctic regions. Using an integrated cryohydrology model, this study presents plot‐scale end‐of‐century projections of ice‐wedge polygon water budget components and active layer thickness with and without trough lateral surface flow under the Representative Concentration Pathway 8.5 scenario. Trough lateral surface flow is incorporated through a newly developed empirical model, evaluated against field measurements. The numerical scenario that includes trough lateral surface flow simulates discharge (outflow from a polygon) and recharge (rain‐induced inflow to a polygon trough from upslope areas), while the scenario that does not include trough lateral surface flow ignores recharge. The results show considerable reduction (about 100–150%) in evapotranspiration and discharge in rainy years in the scenarios ignoring trough lateral surface flow, but less effect on soil water storage, in comparison with the scenario with trough lateral surface flow. In addition, the results demonstrate long‐term changes (~10–15 cm increase) in active layer thickness when trough lateral surface flow is modeled. This study highlights the importance of including lateral surface flow processes to better understand the long‐term thermal and hydrological changes in low‐relief polygonal tundra regions under a changing climate.

Jan, Ahmad↗