The study, carried out by Liza Hovhannisyan of the A. I. Alikhanyan National Laboratory in Yerevan, Armenia, uses the CORSIKA Monte Carlo code, a well-established tool for simulating cosmic-ray air showers that includes the effects of atmospheric electric fields on particle transport. The simulations model a vertically oriented, uniform electric-field layer 2000 meters thick above the detector level, subdivided into twenty 100-meter slices so that the vertical growth of the electron population can be tracked in detail. Four high-altitude observational sites with distinct atmospheric densities were chosen, all of which have recorded thunderstorm-related particle enhancements: Aragats in Armenia at roughly 3000 meters above sea level, Nor Amberd at about 2000 meters, Lomnický Štít in the High Tatras at approximately 2600 meters, and the LHAASO site in China at around 4400 meters. Across the four configurations, the simulated field regions span altitudes of roughly 4000 to 6400 meters, where air densities range from 0.66 to 0.89 kilograms per cubic meter.

The electric fields applied in each configuration were tuned to the local runaway threshold, which scales with the normalized air density. Because thinner air at high altitude lowers the threshold field, the same absolute field strength means something quite different at LHAASO than it does at Nor Amberd. The simulations therefore used field ranges of 180 to 210 kilovolts per meter at LHAASO, 200 to 230 at Aragats, 210 to 240 at Lomnický Štít, and 230 to 260 at Nor Amberd, varied in steps of 10 kilovolts per meter. Each simulation injected a single seed electron with an initial energy spectrum drawn from the EXPACS model, following a power law between 1 and 300 megaelectronvolts, and between 500 and 1000 events were run for each field configuration to ensure statistical stability. Electrons and gamma rays were tracked until their energies fell below 0.05 megaelectronvolts, with particle numbers recorded at every layer boundary.

Within the supercritical region where the applied field exceeds the local runaway threshold, the number of runaway electrons grows approximately exponentially with distance, and the characteristic e-folding length of that growth is the central parameter of avalanche theory. A widely used empirical relation, introduced by J. R. Dwyer in 2003, links this avalanche length to the electric-field excess above the threshold through a single proportionality coefficient, with Dwyer reporting a value of about 7.2 times ten to the third kilovolts. When Hovhannisyan fitted this standard relation separately to the effective avalanche lengths extracted from the vertical electron profiles at each station, the results told a revealing story. The station-specific coefficients ranged from 8.6 to 11.2 times ten to the third kilovolts, all larger than the original value, and the spread between stations persisted even after fit-derived uncertainties were taken into account. A single global fit across all sixteen simulation points yielded a best-fit coefficient of about 9.58 times ten to the third kilovolts, but with a coefficient of determination of only about 0.90, leaving noticeable inter-station scatter.

The remedy was to let the coefficient itself depend on air density. By introducing a generalized relation with an additional density-dependent term, in which the avalanche length is multiplied by the normalized air density raised to a fitted exponent, the joint fit across all four stations produced a coefficient of about 6.39 times ten to the third kilovolts and an exponent of minus 0.923, with an uncertainty of plus or minus 0.083. That exponent, statistically indistinguishable from minus one, indicates an approximately inverse dependence on atmospheric density beyond what is already captured by the density-scaled runaway threshold. The improvement in descriptive power was striking: the coefficient of determination rose from roughly 0.90 to roughly 0.99, and plots comparing simulated and predicted avalanche lengths showed much better agreement, particularly at the larger avalanche lengths where the standard formulation deviated most from a one-to-one relation.

The study then turned to a second, complementary stage of avalanche evolution: what happens to the particles after they leave the accelerating field. Additional simulations for the Aragats configuration placed the detector level 25, 50, 100, and 200 meters below the lower boundary of the electric field, keeping all other parameters fixed. The results revealed sharply divergent transport behavior between the two particle species. The electron spectrum was progressively attenuated with increasing field-free distance, with a pronounced suppression of the high-energy component, and the electron population decreased rapidly and monotonically. Gamma rays, by contrast, persisted over larger distances, attenuated more gradually through stochastic interactions such as Compton scattering, photoelectric absorption, and pair production. At a separation of 100 meters, the electron component was already strongly depleted while a substantial gamma-ray component remained, making that distance a natural calibration point.

This asymmetry matters for a phenomenological quantity called the free path distance, introduced in earlier work to characterize how far energetic particles can propagate beyond the electric-field region and still reach ground-based detectors. Because the maximum electron energy at the field boundary cannot be measured directly, the formulation estimates it from the maximum gamma-ray energy at the detector through an empirical electron-gamma energy coefficient. For the 100-meter calibration separation, the simulated gamma-ray spectra changed little over the propagation distance, with characteristic maximum energies of about 29 to 31 megaelectronvolts at both levels, validating the use of the detector value as a proxy. Electron spectra, however, shifted markedly toward lower energies: at Aragats, the characteristic maximum dropped from roughly 38 megaelectronvolts at the field exit to around 16 megaelectronvolts after 100 meters of propagation through air.

To quantify the electron-gamma coefficient, the study employed a bootstrap-based statistical procedure with 2000 Poisson resampling realizations of the energy histograms, applying a criterion that defines the characteristic maximum energy as the highest energy with at least three consecutive bins containing more than five counts, which suppresses sensitivity to isolated outliers. The station-specific median coefficients clustered tightly, from about 1.31 at Aragats to 1.46 at Lomnický Štít, and pooling all stations yielded a median of about 1.38, in close agreement with an independent global least-squares optimization that gave 1.375. The adopted value of 1.37 replaces the earlier empirical coefficient of 1.2, which the bootstrap analysis showed systematically underestimates the free path distance at all four stations. The second parameter, the altitude-dependent electron energy-loss rate, was obtained from the NIST ESTAR stopping-power database scaled by local air density, yielding values between roughly 0.17 and 0.23 megaelectronvolts per meter across the sites. Applying the recalibrated formulation reproduced free path distances of 85 to 105 meters, consistent with the 100-meter calibration scale, and an independent test at a 50-meter separation reconstructed a median of 69.9 meters, correctly identifying the short-distance regime.

The practical payoff of this work lies in the interpretation of observations. Mountain observatories such as Aragats, Lomnický Štít, and LHAASO routinely detect particle bursts during thunderstorms, and inferring the structure and strength of the electric fields inside thunderclouds from those detections depends on having reliable models of both avalanche growth and post-field transport. By showing that a single density-independent coefficient leaves systematic residuals across sites spanning a substantial range of atmospheric densities, and by demonstrating that a simple power-law correction removes nearly all of that scatter, the study provides a unified scaling that can be applied across different altitudes and storm conditions. The recalibrated free path distance formulation likewise gives experimentalists a firmer bridge between the gamma-ray spectra they measure and the electron populations that produced them at the edge of the accelerating field.

The author notes that the results are derived for idealized, uniform electric-field configurations and should not be generalized to arbitrary thunderstorm field structures without further validation, a caveat that reflects the persistent difficulty of measuring the spatial and temporal structure of fields inside storm clouds. Even so, the near-perfect fit achieved by the density-dependent formulation, together with the stability of the recalibrated transport coefficients across four very different atmospheric environments, marks a meaningful step toward a common physical framework for the growing family of high-energy atmospheric phenomena. All simulation inputs, analysis scripts, and processed outputs are publicly archived, allowing other groups to reproduce the workflow and extend the calibration to new observational sites and more realistic field geometries.

Subject of Research: Density-dependent scaling and transport of relativistic runaway electron avalanches in thunderstorm electric fields

Article Title: Relativistic runaway electron avalanches: unified density-dependent scaling and transport

Article References: Hovhannisyan, L. (2026). Relativistic runaway electron avalanches: unified density-dependent scaling and transport. Geoscientific Model Development, 19(19), 9289-9300. https://doi.org/10.5194/gmd-19-9289-2026

Image Credits: AI Generated

DOI: 10.5194/gmd-19-9289-2026

Keywords: relativistic runaway electron avalanches, thunderstorm electric fields, cosmic rays, gamma-ray glows, terrestrial gamma-ray flashes, CORSIKA simulations, air density, avalanche length, free path distance, thunderstorm ground enhancements, Monte Carlo modeling, high-altitude observatories

The core challenge the researchers confronted is one of boundary characterization. Every battery thermal model depends on assumptions about how heat leaves the cell’s surface, and at sea level those assumptions are well rehearsed. Air density, dynamic viscosity, and thermal conductivity are treated as constants or simple functions of temperature, and convective heat-transfer coefficients follow familiar correlations. But when atmospheric pressure drops from 100 kilopascals to 50 kilopascals, the air density falls by roughly half, and the viscosity and thermal conductivity of the gas shift in ways described by Sutherland’s law and related temperature-correction formulations. A simulation that ignores these corrections will systematically misjudge how effectively the surrounding air can carry heat away from the cell, and the error compounds at higher discharge rates when more heat is being generated in the first place.

To build a model adequate to this environment, the team anchored their electrochemical description in the Newman–Tiedemann–Gu–Kim, or NTGK, framework, a semi-empirical approach that represents the battery through parameters related to the equilibrium potential and the equivalent conductance of the cell. The NTGK model is prized in battery engineering because it captures the essential electrochemistry without the computational burden of full porous-electrode theory, making it practical for the multi-scale multi-domain simulations needed to resolve temperature fields across an entire cell. The researchers embedded this electrochemical core within a coupled thermal-fluid simulation, and then made the crucial modification: at the convective heat-transfer boundary, they corrected the air density and viscosity for the actual ambient pressure and temperature, so that the fluid exchanging heat with the battery behaved like the thin air it truly was.

Validation came from experiment. The team compared their simulation predictions against measured battery surface temperatures and found an absolute error of less than 1.5 kelvin and a relative error of less than 5 percent in predicting the maximum surface temperature. The simulated temperature-rise trends tracked the experimental curves closely across the tested conditions. In a field where thermal models of batteries under nonstandard atmospheres often carry uncertainties several times larger, this level of agreement matters. It means engineers can now trust the model’s output when designing thermal management systems for aircraft battery packs or high-altitude installations, rather than treating the simulations as rough qualitative guides requiring extensive physical retesting.

The quantitative findings carry a message that is both reassuring and cautionary. Under identical atmospheric pressure, when the discharge rate increased from 0.5C to 2C, the battery’s temperature rise grew by 7.42 to 8.21 kelvin. That is a substantial thermal penalty for pushing a cell harder, and it confirms that discharge rate is the dominant lever on battery thermal response. By contrast, under the same discharge rate, when atmospheric pressure fell from 100 to 50 kilopascals, the overall battery temperature increased by only 0.49 to 1.37 kelvin. The effect of halving the ambient pressure is real but modest, roughly an order of magnitude smaller than the effect of quadrupling the discharge current.

Perhaps more intriguing is what the pressure change did not do. The temperature difference between the interior and exterior of the battery remained essentially constant at approximately 1.5 kelvin across the tested pressures, showing minimal impact from the reduced atmosphere. This is a technically significant observation. It suggests that low pressure does not dramatically alter the internal thermal gradients within the cell, which are governed mainly by the cell’s own thermal conductivity and the distribution of heat generation. Instead, the primary mechanism by which reduced pressure harms battery thermal performance is external: thinner air convects heat away less effectively, so heat accumulates at the surface and the whole cell runs slightly hotter, but the interior-to-surface temperature difference stays largely intact.

The physics behind this asymmetry is worth unpacking. Convective heat transfer at a battery surface scales with the properties of the boundary-layer air, and as density decreases, the mass flow of air available to absorb and transport heat diminishes even if the flow velocity is unchanged. The specific heat capacity per unit volume of the air drops in proportion to its density, so each cubic meter of thin air carries away less thermal energy than the same volume at sea level. Meanwhile, the heat generated inside the cell, which the NTGK model partitions into irreversible contributions from electrochemical resistance and reversible contributions from entropy change, depends on the discharge current and the cell’s state of charge, not on the ambient pressure. Discharge rate therefore controls the heat input directly, while pressure controls only the efficiency of heat rejection, which explains why the former dominates the thermal response.

For the emerging field of electric aviation, these results arrive at a critical moment. Battery thermal runaway hazards at cruise altitudes have become a focus of intense research, with prior studies documenting how cylindrical and pouch cells behave under the low pressures typical of civil aircraft flight, and how thermal runaway propagates differently in modules at reduced ambient pressure. The new work complements that hazard-focused literature by providing a validated predictive tool for normal operating conditions, the regime in which thermal management systems must function every flight. A designer of an air-cooled battery pack for an electric aircraft can now use the corrected convective boundary conditions to estimate how much extra temperature margin is needed at altitude, and the answer, according to this study, is on the order of one kelvin per halving of pressure at moderate discharge rates, a manageable but non-negligible correction.

The study also speaks to a broader methodological debate in battery modeling. Computational fluid dynamics approaches to battery thermal management have proliferated, covering pouch cells, prismatic cells, and full modules with liquid cooling and phase-change materials. But many of these studies implicitly assume sea-level atmospheres, and the authors’ demonstration that pressure-corrected air properties must enter the convective boundary to achieve sub-1.5-kelvin accuracy is a warning that extends well beyond low-pressure applications. Any simulation of an air-cooled battery system operating in a hot, dry climate, at elevation, or inside a pressurized cabin with altered air composition should apply the same rigor to its boundary conditions. The framework’s modular structure, combining the NTGK electrochemical core with corrected fluid properties, makes such extensions straightforward.

Funded by the National Natural Science Foundation of China and the Natural Science Foundation of Hunan Province, the work from the Hunan University of Technology and its collaborators provides a quantitative foundation for a question that has lingered as batteries have taken to the skies. The takeaway is crisp: if you want to keep a lithium-ion battery cool, the discharge rate you demand from it matters far more than the altitude at which you fly it, but thin air will quietly erode your cooling margin, and now you can calculate exactly by how much. As electric aircraft move from prototypes to passenger service, models like this one, validated against experiment and honest about their boundary conditions, will be part of the safety case that makes that transition possible.

Subject of Research: Thermal behavior and convective heat-transfer modeling of lithium-ion batteries under low-pressure conditions

Article Title: Convective heat-transfer boundary characterization and numerical simulation of lithium-ion batteries under low-pressure conditions

Article References: Zhang, Y., Huang, Z., Liu, S., Jiang, Y., & Shen, S. (2026). Convective heat-transfer boundary characterization and numerical simulation of lithium-ion batteries under low-pressure conditions. Ionics. https://doi.org/10.1007/s11581-026-07553-5

Image Credits: AI Generated

DOI: 10.1007/s11581-026-07553-5

Keywords: lithium-ion batteries, low pressure, convective heat transfer, NTGK model, multiphysics simulation, thermal management, electric aircraft, discharge rate, air density, battery temperature, numerical simulation, heat dissipation

The motivation for the study stems from a persistent blind spot in climate scenario databases. Future projections produced by coupled climate models under the Shared Socioeconomic Pathways are readily available for some variables, such as air temperature and precipitation, but are not directly published for several key atmospheric quantities. Dew point temperature, wet-bulb temperature, humidity ratio, water vapor partial pressure, and total column precipitable water all fall into this category, even though each plays a decisive role in human heat stress, agriculture, water resources, and the energy balance of the climate system. Without forward projections of these variables, comprehensive assessments of future climate risk remain incomplete, and adaptation planners are left working with half the picture.

To overcome this limitation, the researcher developed an integrated forecasting framework that combines three complementary ingredients: externally available climate scenario data, a statistical time-series forecasting model known as ARIMAX, and physically based thermodynamic equations. The logic is elegant in its simplicity. Where future scenario outputs exist, they are used directly. Where they do not, the ARIMAX model—a variant of the classical autoregressive integrated moving average approach that incorporates exogenous predictor variables—generates statistically robust projections using air temperature as the driving input. Finally, established psychrometric relationships, of the kind used for decades in engineering handbooks, translate the projected temperatures and humidity measures into the remaining indicators, ensuring that every number in the resulting dataset is physically consistent with every other.

Specifically, future projections of dew point temperature and total column precipitable water were obtained using ARIMAX with air temperature as the exogenous variable, while the remaining indicators—wet-bulb temperature, dew point depression, relative humidity, humidity ratio, water vapor partial pressure, air density, and total precipitation—were derived through established thermodynamic relationships. This hybrid design means the framework does not attempt to conjure atmospheric physics from statistics alone. Instead, it anchors the statistical forecast to a well-constrained predictor and lets the known equations of moist air physics do the rest, a strategy that dramatically reduces the risk of producing projections that violate conservation laws or known vapor-pressure behavior.

The historical reconstruction, spanning 1940 to 2024, establishes the baseline against which the future is measured. By 2024, the study records global mean dry-bulb temperature at 15.1 degrees Celsius, dew point temperature at 10.0 degrees Celsius, and wet-bulb temperature at 12.2 degrees Celsius. Dew point depression—the gap between air temperature and dew point, a measure of how much additional moisture the air can hold before saturation—stood at 5.1 degrees Celsius. Total column precipitable water, the depth of liquid water that would result if all vapor in a vertical column of atmosphere were condensed, reached 25.6 kilograms per square meter. Water vapor partial pressure registered 12.3 hectopascals, the humidity ratio 7.8 grams of water vapor per kilogram of dry air, relative humidity 71.4 percent, air density 1.191 kilograms per cubic meter, and total annual precipitation 1.073 meters.

Projected forward under SSP5-8.5, the high-emissions pathway in which fossil fuel use continues largely unabated through the century, these indicators shift substantially by 2100. Dry-bulb temperature climbs to 18.9 degrees Celsius, a rise of 3.8 degrees over the 2024 value. Dew point temperature reaches 13.4 degrees Celsius and wet-bulb temperature 15.6 degrees Celsius, while dew point depression widens slightly to 5.5 degrees Celsius. Total column precipitable water rises to 30.0 kilograms per square meter, water vapor partial pressure increases to 15.5 hectopascals, and the humidity ratio grows to 9.9 grams per kilogram. Relative humidity edges down marginally to 70.5 percent, air density thins to 1.174 kilograms per cubic meter, and total precipitation increases to 1.128 meters per year.

Several of these projected changes carry implications that extend well beyond academic climatology. The rise in wet-bulb temperature is perhaps the most consequential for human societies, because this quantity effectively measures the body’s ability to cool itself through sweating; sustained wet-bulb values approaching the mid-thirties Celsius are incompatible with human survival outdoors, and even far lower values impose severe health and labor-productivity penalties. A global mean wet-bulb temperature of 15.6 degrees Celsius by 2100 implies that regional extremes in already hot, humid parts of the world will push far closer to dangerous thresholds than they do today. Meanwhile, the roughly 17 percent increase in precipitable water is a direct manifestation of the Clausius–Clapeyron relationship, which dictates that the atmosphere holds about seven percent more water vapor for every degree of warming—a feedback that intensifies heavy rainfall events even as some regions grow drier overall.

The projected decline in air density, from 1.191 to 1.174 kilograms per cubic meter, is a subtler but technically significant outcome. Warmer, more humid air is less dense, and this thinning has measurable consequences for wind turbine power output, aircraft lift and engine performance, and the combustion efficiency of engines and turbines. As the atmosphere warms, engineers and energy planners will increasingly need to account for these density effects in siting and design decisions. The slight decrease in relative humidity alongside rising absolute moisture content reflects the fact that warmer air’s capacity to hold vapor grows faster than actual vapor supply in many regions, a distinction that matters for evaporation rates, wildfire risk, and agricultural water demand.

The study’s methodological contribution may prove as influential as its numerical results. By demonstrating that ARIMAX modeling anchored to temperature projections can generate physically plausible, long-horizon estimates of variables missing from standard scenario archives, the framework provides a practical, reproducible template that other researchers can apply at regional scales or with different scenario pathways. The author acknowledges the Climate Change Institute at the University of Maine for supplying the reanalysis data underpinning the historical record. As climate impact assessment, adaptation planning, and long-term environmental decision-making increasingly demand multi-variable, internally consistent futures, tools of this kind bridge the gap between what climate models currently publish and what practitioners actually need—turning the full thermodynamic state of the future atmosphere from an afterthought into a forecastable, actionable dataset.

Subject of Research: Historical analysis and future projections of ten fundamental global climate indicators using statistical forecasting and thermodynamic modeling

Article Title: Historical analysis and future projections of fundamental global climate indicators

Article References: Historical analysis and future projections of fundamental global climate indicators. (n.d.). https://doi.org/10.1007/s00704-026-06543-x

Image Credits: AI Generated

DOI: 10.1007/s00704-026-06543-x

Keywords: climate change, climate indicators, ARIMAX forecasting, dew point temperature, wet-bulb temperature, relative humidity, precipitable water, air density, SSP5-8.5 scenario, thermodynamic modeling, climate projections, global warming