Every time a power grid operator decides how much electricity each generator should produce, a quiet balancing act takes place. Burning less fuel saves money, but it often means allowing more pollution, and cutting emissions usually drives costs up. Engineers have long searched for a way to satisfy both goals at once, a challenge known as combined economic emission dispatch, or CEED. Now, a team of researchers has turned to one of the ocean’s most graceful feeders, the manta ray, for inspiration, and the result is an optimization algorithm that outperforms a crowded field of competitors on both cost and cleanliness.

The study, published in the journal Results in Engineering by Yu-Feng Sun, Jie-Sheng Wang and colleagues, describes an enhanced version of the Manta Ray Foraging Optimization algorithm, or MRFO. The original algorithm, first introduced a few years ago, mimics the cooperative hunting strategies that manta rays use to harvest plankton. Because plankton drifts unpredictably with seasons and weather, manta rays have evolved three distinct foraging behaviors: chain foraging, in which animals line up in arcs so that each one catches food missed by the one ahead; spiral foraging, where the group coils around a dense patch of food; and somersaulting, in which individuals flip back and forth around a promising spot. Translated into mathematics, these behaviors become rules that move a population of candidate solutions around a search space, gradually converging on the best answer.

What makes the CEED problem so stubborn is its mathematical character. The relationship between a generator’s output and its fuel consumption follows a cubic curve, and the emissions of carbon dioxide, nitrogen oxides and sulfur dioxide each follow their own cubic functions. The combined objective is nonlinear, nonconvex and discontinuous, riddled with local minima that can trap classical techniques such as lambda iteration or gradient-based methods. To make matters harder, the problem is genuinely multi-objective, so the researchers converted it into a single objective by applying price penalty factors that fold the cost of each pollutant into one total figure in dollars per hour.

Constraint handling is where many algorithms stumble. Traditional approaches add a fixed penalty whenever a solution violates the power balance requirement, the equation demanding that total generation minus transmission losses equals demand. But a fixed penalty is a blunt instrument. Set it too low and the algorithm happily returns infeasible schedules; set it too high and the search becomes so restricted that it cannot explore properly. Finding the right value through trial and simulation is time-consuming and inefficient. The team’s answer was to make the penalty itself dynamic, drawing on five different mathematical functions: sinusoidal, hyperbolic tangent, arctangent, linear and hyperbolic secant. Each function is raised to an order that changes during the run, so the penalty strength adapts to both the size of the violation and the stage of the optimization. Early on, mild penalties preserve exploration; later, penalties tighten to push solutions toward feasibility.

The second major innovation targets the heart of the MRFO algorithm. In the original design, a parameter called C governs the switch between exploration, casting the net wide across the search space, and exploitation, homing in on the best solution found so far. This parameter rises monotonically from zero to one, a predictable march that can leave the algorithm poorly balanced at certain stages. The researchers replaced it with a chaotic switching factor, drawn from nine well-known chaotic maps including Chebyshev, logistic, tent, sine, Gauss and singer maps, combined with an exponential adjustment function. The result is a parameter that still trends upward but oscillates along the way, periodically reopening the search space and preventing the population from settling prematurely into a local optimum.

To find out which chaotic map worked best, the team ran all nine enhanced variants, named C1MRFO through C9MRFO, through the CEC-2022 benchmark suite of twelve test functions, using a population of thirty and five hundred iterations averaged over thirty independent runs. The Chebyshev-based variant, C1MRFO, emerged as the overall winner, ranking first among the nine and beating the original MRFO, which finished last. Wilcoxon rank-sum tests confirmed that most of the improvements were statistically significant, with p-values below the conventional 0.05 threshold. When C1MRFO was pitted against a roster of state-of-the-art metaheuristics, including the dhole optimization algorithm, evolutionary mating algorithm, dung beetle optimizer, snow geese algorithm, particle swarm optimization, grey wolf optimizer, sparrow search algorithm and marine predators algorithm, it achieved the best final rank, taking first place on many of the test functions while keeping computational times competitive.

The real test, however, came when the algorithm was applied to actual dispatch problems. The researchers simulated a six-unit power system under four load demands of 150, 175, 200 and 225 megawatts. Across all four scenarios, C1MRFO produced the lowest comprehensive total cost, which combines fuel cost with the weighted penalties for carbon dioxide, nitrogen oxide and sulfur dioxide emissions. At 150 megawatts, for example, it reached a total cost of about 10,137 dollars per hour, edging out the grey wolf optimizer and clearly beating algorithms such as the arithmetic optimization algorithm, Harris hawks optimization, whale optimization algorithm, sine cosine algorithm, particle swarm optimization and the lambda iteration method. The pattern held at every load level, with the margin over the weakest classical method growing to thousands of dollars per hour at 225 megawatts.

The dynamic penalty functions delivered a further, striking gain. When the sinusoidal penalty function was applied with the parameter order set to dy5, the total cost at 150 megawatts dropped from roughly 10,137 to 9,065.71 dollars per hour, and similar reductions appeared at the higher loads, reaching 14,621.45 dollars per hour at 225 megawatts. Comparing the five penalty shapes, the sinusoidal function consistently produced the lowest minimum costs, beating the hyperbolic tangent, arctangent, linear and hyperbolic secant alternatives by margins ranging from hundreds to nearly two thousand dollars per hour depending on the load. Notably, when the penalty order was fixed rather than dynamic, performance was worst across the board, confirming that the adaptivity of the penalty, not merely its presence, is what matters.

To verify that each improvement contributes on its own, the team conducted an ablation study on a much larger twenty-unit system with a total demand of 2,500 megawatts. The results formed a clear hierarchy: the original MRFO achieved 61,835.16 dollars per hour, the chaos-enhanced C1MRFO reached 61,826.11, and the fully improved version combining chaos with the optimal sinusoidal dynamic penalty, sin(dy5)MRFO, achieved the best figure of 61,807.36 dollars per hour. All three outperformed the phasor particle swarm optimization benchmark, which came in at 61,895.57 dollars per hour. The stepwise improvement demonstrates that the two enhancements are not redundant but superimpose their benefits, and that the framework scales to systems far larger than the six-unit case.

The authors are careful to note where the work goes next. The current approach converts the multi-objective problem into a single objective through price penalty factors, but future versions could adopt Pareto-based multi-objective optimization to generate full sets of non-dominated trade-off solutions. The test cases also assume deterministic demand, so extending the algorithm to stochastic dispatch with wind and solar integration is another goal, as is combining the dynamic penalty with feasibility-first constraint-handling strategies. For now, the study offers a compelling demonstration that borrowing the foraging logic of a filter-feeding fish, perturbing it with chaos, and teaching it to calibrate its own punishments can squeeze meaningful savings and emission reductions out of one of power engineering’s hardest scheduling problems.

Subject of Research: An improved manta ray foraging optimization algorithm with chaotic switching factors and dynamic penalty functions for solving the combined economic emission dispatch problem in power systems.

Article Title: Manta ray foraging optimization algorithm with dynamic penalty function and chaos switching factor to solve combined economic emission dispatch problem

Article References: Sun, Y.-F., Wang, J.-S., Zhang, X.-Y., Li, Y.-X., Guan, X.-Y., & Liu, X. (2026). Manta ray foraging optimization algorithm with dynamic penalty function and chaos switching factor to solve combined economic emission dispatch problem. Results in Engineering, 32, Article 113383. https://doi.org/10.1016/j.rineng.2026.113383

Image Credits: AI Generated

DOI: Not provided

Keywords: manta ray foraging optimization, combined economic emission dispatch, metaheuristic algorithms, chaotic maps, dynamic penalty function, power system optimization, emission reduction, fuel cost minimization, CEC-2022 benchmarks, swarm intelligence, constraint handling, Results in Engineering

When a tunnel catches fire, temperatures inside the concrete lining can soar to 800 or 900 degrees Celsius, and the emergency response almost always involves dousing the structure with water. A new study suggests that this familiar firefighting image—flames followed by a torrent of cooling water—may be doing far more damage to the concrete than engineers have typically accounted for. Researchers investigating basalt fiber-reinforced concrete have found that how a fire-damaged structure is cooled can be as important as how hot it got in the first place, with water quenching consistently inflicting greater strength loss than slow, natural cooling across every temperature they tested.

The research, published in Case Studies in Construction Materials, systematically examined how heating temperature, cooling regime, fiber content, and impact loading combine to determine whether concrete survives a fire with its load-bearing capacity intact. The team, led by Yang Liu and colleagues including Delin Li, Jinyue Dong, Yingda Zhang, Xianliang Zhou, and Qiong Zhou, chose basalt fiber-reinforced concrete, or BFRC, as their test material because of the fiber’s unusual combination of properties. Basalt fibers, drawn from volcanic rock, offer high tensile strength, thermal stability up to melting points between 1350 and 1450 degrees Celsius, and a cost profile far more attractive than premium alternatives such as carbon fiber. This makes them an appealing candidate for tunnel linings, where fire resistance and economy must be balanced.

The experimental design was exhaustive. The researchers prepared a C35-grade concrete mix—the grade commonly recommended for tunnel design—using ordinary Portland cement, medium sand, and continuously graded gravel between 5 and 10 millimeters. Into this base mix they incorporated chopped basalt fibers, 12 millimeters long with a monofilament diameter of 17 micrometers, at volume fractions of 0, 0.1, and 0.2 percent. The fibers were added in small batches during dry mixing to prevent agglomeration, and the fresh concrete was cast into cylindrical specimens measuring 50 millimeters in diameter and 25 millimeters in height, then cured under standard conditions at 20 degrees Celsius and 95 percent relative humidity.

The thermal assault began in a high-temperature electric oven capable of reaching 1200 degrees Celsius. Specimens were heated at a controlled rate of 5 degrees Celsius per minute to target temperatures of 200, 400, 600, and 800 degrees Celsius, then held at temperature for a full hour to ensure the heat penetrated the interior. Half the specimens were then removed and allowed to cool naturally in ambient air, mimicking a structure left to cool on its own after a fire. The other half were plunged immediately into a tank of tap water, replicating the thermal shock of firefighting operations. After cooling, the water-immersed specimens were conditioned in air for seven days to return them to a natural moisture state before testing.

The dynamic mechanical properties were then probed using a split Hopkinson pressure bar, the workhorse instrument for measuring how materials behave under the extreme strain rates produced by impacts and blasts. The apparatus used bars 75 millimeters in diameter, made of high-strength steel with an elastic modulus of 210 gigapascals, with a striker bar propelled by compressed gas striking an incident bar that transmitted stress waves through the specimen. Three impact velocities—3.5, 4.4, and 5.2 meters per second—were tested, achieved by varying the operating gas pressure, with three specimens per condition. A thin pulse shaper was pasted at the impact end of the incident bar to smooth out high-frequency oscillations, and stress equilibrium between the two faces of each specimen was verified from strain-gauge signals, confirming the reliability of the data. Fragments from each impact were sieved through meshes ranging from 4.75 millimeters down to 75 micrometers so that a fractal dimension could be calculated, providing a quantitative measure of how thoroughly each specimen shattered.

The visual evidence of thermal damage followed a familiar but telling progression. At 200 degrees Celsius, the concrete surfaces turned bluish-gray, with no cracking. By 400 degrees, the color shifted to yellowish-gray, still structurally sound. At 600 degrees, fine cracks appeared along with slight surface peeling, and by 800 degrees the specimens showed severe cracking, spalling, and edge breakage. Water-cooled specimens generally exhibited more pronounced surface damage than their naturally cooled counterparts, and at 800 degrees the concave shape at the start of the stress-strain curve—evidence of pre-existing microcracks being closed under compression—was especially marked in the water-quenched samples.

The strength data told a stark quantitative story. For naturally cooled specimens containing 0.2 percent basalt fiber, the dynamic peak stress at an impact velocity of 3.5 meters per second fell from 51.1 megapascals at room temperature to 18.6 megapascals at 800 degrees—a reduction of 63.6 percent. Water cooling made matters worse across the board: at the same impact velocity, the reductions reached 23.9 percent at just 200 degrees and climbed to 64.9 percent at 800 degrees. Averaged across the temperature range, each 100-degree increment cost naturally cooled specimens roughly 6.7 to 8.0 percent of their strength, while water-cooled specimens lost about 7.5 to 8.1 percent per 100 degrees. The gap between the two cooling regimes widened at higher impact velocities: at 5.2 meters per second and 800 degrees, the water-cooled specimens suffered 6.4 percent greater strength loss than the naturally cooled ones, whereas at the lowest velocity the difference was a mere 1.3 percent. The explanation lies in competing timescales—at slow loading rates, cracks have time to propagate and dissipate energy, blunting the effect of water-induced damage, but at high rates the loading outpaces the material’s energy dissipation, and the microcracks created by thermal shock propagate catastrophically.

Against this backdrop of degradation, the basalt fibers demonstrated genuine protective power. At every temperature from 200 to 800 degrees, specimens with 0.2 percent fiber content outperformed both plain concrete and the 0.1 percent mix. At 400 degrees and an impact velocity of 3.5 meters per second, the naturally cooled fiber-reinforced specimens reached peak stresses of 40.1 megapascals, some 16.9 percent above plain concrete, while at 200 degrees under water cooling the 0.2 percent mix achieved 38.8 megapascals—a remarkable 32.4 percent improvement over unreinforced concrete. The mechanism, the authors explain, is architectural: randomly distributed fibers form a three-dimensional network that binds tightly with the cementitious matrix, reduces internal pores and cracks, and bridges microcracks before they can coalesce into fracture surfaces. But this protection has a ceiling. Above roughly 600 degrees, the fibers begin to degrade and lose their reinforcing capability, and the fiber-matrix interface—already weakened by thermal mismatch between fiber and paste—becomes a liability during water quenching.

Fractal analysis of the impact debris added another layer of insight. The fractal dimension of the fragments rose steadily with temperature, from around 1.82 at 200 degrees to above 2.3 at 800 degrees under the highest loading rate, indicating progressively finer and more complex fragmentation. Water-cooled specimens consistently shattered into finer pieces than naturally cooled ones, particularly at 600 and 800 degrees, confirming that the steep thermal gradients during quenching seed dense populations of microcracks that later explode into fracture networks under impact.

Microstructural investigation grounded these observations in chemistry and physics. X-ray diffraction of specimens heated to 800 degrees showed the disappearance of calcium hydroxide peaks and weakened calcium carbonate signals, evidence that these phases had thermally decomposed—reactions that release vapor, increase porosity, and weaken the matrix. Scanning electron microscopy revealed the deterioration in vivid detail: an unheated specimen showed a dense microstructure with abundant calcium-silicate-hydrate gel and a compact interfacial transition zone around each aggregate, while specimens heated to 800 degrees contained no detectable C-S-H gel at all, their structure riddled with pores and cracks and hosting instead a dicalcium silicate-like phase. Comparing cooling regimes at 400 degrees, the naturally cooled sample remained relatively dense and rich in hydration products, whereas its water-quenched twin displayed obvious cracks at the interfacial transition zone and a loosened surrounding matrix—the fingerprint of thermal shock stress racing through the cooling material.

The practical implications reach well beyond the laboratory. Tunnel linings that survive a fire are routinely assessed and, if necessary, repaired, but this study shows that the assessment itself must account for the cooling history. A lining quenched by firefighting water may be significantly weaker than its peak temperature alone would suggest, particularly if it subsequently faces dynamic loads from traffic vibration, rockfall impacts, or secondary events. Conversely, the results give engineers a positive tool: specifying basalt fiber reinforcement at around 0.2 percent by volume can meaningfully bolster the impact resistance of fire-exposed concrete below 600 degrees, at a modest material cost. The authors frame their work as a theoretical foundation for post-fire safety assessment of tunnel linings and for the wider engineering adoption of basalt fiber-reinforced concrete—a reminder that in structural fire engineering, what happens in the minutes after the flames matter almost as much as the fire itself.

Optical frequency combs have transformed the way scientists measure light and time, earning the Nobel Prize in Physics in 2005 and underpinning technologies ranging from atomic clocks to exoplanet detection. Yet the miniature on-chip versions of these combs, known as soliton microcombs, have long suffered from a frustrating paradox: the most stable and useful operating regimes have also been the most fragile and difficult to sustain. A newly published study in Light: Science & Applications reports a significant advance toward resolving this tension. Researchers led by J. Tang, J. Yang and E. Yan describe a multimodal locking-enabled soliton microcomb that combines robust, self-starting operation with an exceptionally low repetition rate on the scale of a day of stable measurement time, opening a path toward practical high-precision metrology systems that can run unattended for extended periods.

At the heart of the advance is a technique the authors call multimodal locking, a scheme that stabilizes not just a single resonance of the optical microresonator but multiple modes simultaneously, locking the comb dynamics to a regime in which dissipative Kerr solitons can form and persist without delicate manual tuning. In conventional microcomb systems, soliton formation typically requires a precise and rapid traversal of a “chaotic” modulation instability regime, a step that is notoriously sensitive to thermal transients in the resonator chip. Laser frequency must be swept with exquisite control while the soliton’s own heat dissipation shifts the resonator’s resonances, creating a feedback loop that can destroy the soliton state within microseconds. The multimodal locking approach circumvents this fragility by engaging several mode families of the resonator at once, so that the interplay among them pins the comb to a stable soliton state even as environmental and thermal conditions fluctuate.

The significance of achieving a low repetition rate in the same device cannot be overstated. A frequency comb’s repetition rate is the spacing between its individual comb teeth, and for many applications, coarser spacing is a decisive advantage. In laser ranging and distance metrology, for example, the ambiguity range of an interferometric measurement is set by the inverse of the comb’s free spectral range. Microcombs, generated in whispering-gallery-mode or ring resonators only a few millimeters across, usually exhibit repetition rates of tens to hundreds of gigahertz, which severely limits the unambiguous measurement range available to a comb-tooth-resolved distance measurement. By engineering a low-repetition-rate soliton state, the research team extends this ambiguity range dramatically, making the microcomb directly useful for absolute distance measurement, multi-wavelength interferometry and coherent communications, where each comb line must be spectrally resolvable and addressable.

Achieving both low repetition rate and robust soliton generation in a single microresonator is a formidable engineering challenge. Low repetition rates demand large resonator diameters, which in turn produce dense spectra of resonances and increased susceptibility to mode crossings with higher-order transverse and polarization mode families. Such crossings historically destabilize soliton formation, causing comb states to collapse or preventing single-soliton operation altogether. The multimodal locking strategy turns this liability into an asset: rather than engineering mode crossings away, the researchers harness the coupling among mode families to lock the comb dynamics into a favorable basin of attraction. The result is a soliton microcomb whose low repetition rate is preserved over day-scale timeframes, a duration that places the device among the most operationally stable microcomb demonstrations reported to date.

Day-scale stability matters because it transforms the microcomb from a laboratory demonstration into an instrument. Many high-precision experiments, from astronomical spectrograph calibration to the accumulation of statistical confidence in atomic spectroscopy, require comb sources that remain locked and stable across many hours or even days of continuous operation. Conventional fiber-laser-based frequency combs achieve this, but at the cost of table-top size, high power consumption and substantial maintenance overhead. Microcombs promised to shrink the technology onto a photonic chip, yet until now their operational fragility has confined most demonstrations to short observation windows. The robust low-repetition-rate soliton microcomb described in the new work demonstrates that chip-scale frequency combs can, with appropriate locking architecture, sustain the kind of long-duration, hands-off operation that real metrology campaigns demand.

The physics underlying soliton microcombs rests on the Kerr nonlinearity of the resonator material, typically silicon nitride, magnesium fluoride or fused silica. When continuous-wave pump light resonantly builds up inside the ring, the Kerr effect produces a four-wave-mixing process that seeds sidebands spaced by the resonator’s free spectral range. Under the right conditions of pump detuning and intracavity power, these sidebands evolve into a coherent train of ultrashort optical pulses, dissipative Kerr solitons, that circulate around the resonator and emit a broadband comb spectrum with tooth spacing equal to the repetition rate. The comb’s output frequencies are determined by two radio-frequency quantities, the pump laser frequency and the repetition rate, which is why the combs can serve as rulers for measuring optical frequencies against a reference. Any drift in either quantity propagates directly into measurement error, and the multimodal locking mechanism suppresses precisely this drift by anchoring the comb’s underlying mode structure to the resonator’s intrinsic stability.

Applications for such a device span an unusually broad swath of science and engineering. In high-precision laser ranging, low-repetition-rate combs with tooth spacings in the low-gigahertz or even megahertz regime allow interferometric distance measurements over ranges of meters to kilometers without the need for auxiliary synthetic-wavelength techniques. In astronomical spectroscopy, comb calibration sources with well-controlled line spacing enable the detection of Earth-like exoplanets through radial velocity shifts measured in centimeters per second. In telecommunications, low-repetition-rate combs align naturally with standard channel grids used in wavelength-division multiplexed fiber networks, making them attractive as multi-carrier sources. In spectroscopy, dual-comb techniques benefit from repetition rates that are low enough to permit resolution of closely spaced molecular absorption features. A microcomb that is simultaneously robust, low-repetition-rate and day-scale stable therefore functions as a versatile enabling component across these domains rather than a single-purpose device.

The demonstration also speaks to a broader trend in photonics: the migration from delicate, expert-operated laboratory setups to self-contained, self-stabilizing systems. The self-starting character implied by multimodal locking is particularly important for this transition. A comb that requires a skilled operator to initiate soliton formation, sweeping the laser with practiced timing past the chaotic regime, cannot be deployed in an observatory dome, on a satellite platform or inside an industrial instrument. A comb that engages its soliton state reliably on demand and then holds it through thermal cycling, laser drift and mechanical perturbation removes the last practical barrier to integration. The authors’ demonstration of day-scale continuous operation is, in effect, a stress test passed by a technology that researchers have been attempting to harden for more than a decade.

The implications extend to metrology standards and timing as well. Optical atomic clocks, which now achieve fractional frequency uncertainties below one part in 10^18, rely on frequency combs to divide down their optical oscillations to countable microwave rates. Miniaturizing the comb stage of such clocks is a prerequisite for portable optical timekeeping, which would benefit navigation in GPS-denied environments, geodesy that monitors sea level and crustal motion with millimeter precision, and fundamental physics tests that track the constancy of natural constants over months and years. A robust chip-scale soliton microcomb with day-scale stability contributes a key subsystem to this vision, and the multimodal locking principle may generalize to other resonator platforms, materials and spectral regions, including the mid-infrared where molecular fingerprinting awaits compact comb sources.

Looking forward, the study suggests several directions that the photonics community is likely to pursue. Combining multimodal locking with on-chip pump lasers and integrated electronics would complete the transformation toward fully self-contained comb modules with no free-space optics. Extending the approach to lower repetition rates still, into the megahertz regime where fiber-laser combs operate, would bring microcombs into direct competition with their bulkier predecessors for the most demanding ranging and spectroscopy tasks. Improvements in resonator fabrication, reducing mode-volume disorder and engineering dispersion with greater precision, will further expand the parameter space in which robust soliton states can be locked. The present work demonstrates that the fragility once considered intrinsic to soliton microcombs is not a fundamental limit but an engineering problem, and one that can be solved by embracing, rather than avoiding, the multimodal structure of real optical resonators.

The research, published in Light: Science & Applications, marks a milestone in the maturation of microcomb technology. By uniting three properties that had previously been achieved only separately, self-starting robustness, low repetition rate and day-scale operational stability, the work delivers a soliton microcomb genuinely suited to high-precision metrology. As the demand for compact frequency references accelerates across navigation, astronomy, communications and fundamental science, the multimodal locking strategy reported by Tang, Yang, Yan and colleagues offers a clear and technically grounded route from chip-scale optical chaos to dependable, precision measurement light.

Wildfire scientists have unveiled a new rapid risk-assessment tool that shows substantial wildfire danger to homes and buildings in three major Arctic communities, revealing that between 14% and 26% of structures in Anchorage, Fairbanks, and Whitehorse face high or very high risk from approaching flames. The study, led by Jennifer I. Schmidt of the University of Alaska Anchorage together with Robert H. Ziel, Monika P. Calef, Anna Varvak, and Julio C. Postigo, was published in the International Journal of Disaster Risk Science. Rather than relying on complex fire-behavior simulations that demand specialist analysts and expensive software, the approach combines datasets that most local governments already possess: building footprints, zoning maps, and readily available vegetation layers. The result is an updatable, affordable portrait of which buildings face the greatest threat and where mitigation dollars can do the most good.

The urgency behind the work is hard to overstate. Wildfire costs and damages in the United States alone are now estimated at between USD 394 billion and 893 billion annually, and the 2019 to 2020 Australian fire disasters cost nearly USD 100 billion. The Arctic, often imagined as too cold and wet to burn, is not immune: since 1980, wildfires have cost Alaska between USD 2 billion and 5 billion, with projections suggesting costs could rise 68% by the end of the century. Canadian suppression expenses are expected to climb by 60 to 110% relative to the 1980 to 2009 baseline. Meanwhile, structure loss from wildfire has doubled over the past two decades in the western United States, and western Canada has suffered major structure-loss events in four of the last eight fire seasons. Nearly half the global population now lives within the wildland-urban interface, the zone where people and forests meet or intermingle, making tools that translate fire science into actionable local maps increasingly vital.

The researchers’ central premise is conceptually simple but technically consequential: for a structure to be at high risk, a flammable hazard must exist nearby, within roughly 500 meters, and the surrounding landscape must exhibit high wildfire exposure. The method begins with a flammability hazard layer at 30-meter resolution, built from a modified version of NASA’s ABoVE Landcover dataset in which generic evergreen classes were refined to reflect locally flammable species such as black spruce, white spruce, lodgepole pine, and subalpine fir, each assigned hazard ratings from 0 to 100. Critically, the team treated structures themselves as fuel. Each building contributes three fuel units to the surrounding flammability score, a value chosen after systematic testing of one through ten fuel units; three proved optimal, nudging densely built environments up a hazard category without causing wholesale reclassification. Structure density was computed at 90-meter resolution, multiplied by three, and added to the vegetation-based hazard to produce an enhanced flammability hazard layer.

Wildfire exposure, the second pillar of the assessment, captures the potential for hazard to travel. It is calculated as the average flammability hazard within a 100-meter or 500-meter radius, representing the likelihood that fire can reach and burn a given location even if the local fuels themselves are sparse. Exposure values, also ranging from 0 to 100, are classified into five tiers from very low to extreme, and concentrated pockets of flammable vegetation drive scores upward. The third ingredient is structure vulnerability, produced by multiplying the enhanced flammability hazard by a structure importance factor derived from zoning data: residential buildings and critical infrastructure such as hospitals, schools, and nursing homes receive the highest weight, followed by commercial and other structures. Where zoning was unavailable, as on military installations like Joint Base Elmendorf-Richardson in Anchorage and Fort Wainwright and Eielson Air Force Base in Fairbanks, the team collaborated with base personnel or interpreted aerial imagery to classify neighborhoods. The final risk score synthesizes hazard, exposure, and vulnerability into six categories, from extremely low to very high.

The results paint a nuanced picture of Arctic fire danger. Fairbanks showed the greatest proportion of structures in the top two risk categories at 26%, followed by Whitehorse at 22% and Anchorage at 14%. Anchorage, paradoxically, had the highest share of buildings in the top vulnerability tiers, 42%, yet the lowest risk, a consequence of its large urban footprint where dense housing raises vulnerability but fragmented, less-flammable vegetation and abundant developed ground suppress exposure. Fairbanks and Whitehorse, with structures woven into spruce- and pine-dominated boreal landscapes, showed closer alignment between vulnerability and risk. Anchorage was also the only study area where more than 1% of cells jumped two or more flammability categories when buildings were counted as fuel, at 7%, reflecting its 27% share of structure area containing seven or more buildings per cell. The analysis also confirmed the primacy of exposure: when a structure’s wildfire exposure was extreme, at least 87% of such buildings landed in the highest risk category across all three communities.

A key innovation lies in how the team redrew the wildland-urban interface itself. Traditional WUI maps distinguish urban, interface, and intermix zones; the researchers added a fourth category, the transmission zone, marking areas roughly 2.4 kilometers beyond the WUI boundary where fire could move into a community and where fuel breaks might be sited. Within the 500-meter buffer around structures, urban areas were defined by exposure values of 20 or less, interface zones by 21 to 60, and intermix by values above 60, cross-referenced with Alaska and Yukon wildfire management designations that prioritize suppression from critical down to unplanned or wilderness levels. Fairbanks claimed the largest overall WUI footprint, and the transmission zone proved the largest category in every study area, spanning 23,347 square kilometers around Anchorage, 6,180 square kilometers around Fairbanks, and 665 square kilometers around Whitehorse. Across all sites, most structures fell within interface zones, 52 to 61%, followed by intermix at 20 to 35%. Thousands of buildings sit in the high-risk interface and intermix zones where mitigation should concentrate: 15,537 structures in Anchorage, 16,212 in Fairbanks, and 3,734 in Whitehorse.

The method faced its most demanding test against a real disaster. On 17 August 2019, during southcentral Alaska’s hottest and driest summer on record, the likely human-caused McKinley wildfire was driven by high winds across 13.3 square kilometers in three days, destroying 52 residential structures, 84 outbuildings, and three commercial buildings, affecting 1,028 properties and forcing roughly 350 to 400 people to evacuate. By comparing pre-fire aerial imagery from 2017 with post-fire imagery from 2020 at 0.3-meter resolution and linking vanished buildings to the Matanuska-Susitna Borough structure database, the team identified exactly which structures burned. The verdict was statistically compelling: burned structures had significantly higher wildfire exposure than those that survived, and the proportion of burned structures with exposure values above 60 was significantly larger than among unburned ones. Structure risk scores above 60 were also significantly associated with burning. Notably, the test held even under conditions so extreme that normally resistant, low-flammability fuels were igniting.

The study’s design was shaped by sustained community engagement. Through workshops in 2019, 2021, and 2023, along with 15 to 20 meetings organized by Anchorage community councils and landowner associations and public events ranging from wildfire awareness days to Arbor Day, residents and agency staff reshaped everything from color schemes to model parameters. One consequential example: Anchorage residents and fire managers insisted that the maps reflect the spruce beetle outbreak first detected in 2016 and peaking around 2020, which killed and removed many spruce trees. The team responded with rules derived from aerial detection surveys, reclassifying white spruce cells in areas of beetle intensity of 11% or greater to mixed forest, and shifting mixed forest cells to open shrub, while a succession model updated vegetation on areas burned within the last 30 years. This iterative dialogue ensured outputs that both wildfire practitioners and the general public can interpret, a prerequisite the authors see as central to building fire-adapted communities.

The practical implications for planning are direct. Because mitigation funding is often tied to officially designated WUI areas, the ability to rapidly recalculate boundaries from updated building layers expands funding opportunities and keeps pace with community growth. The framework also points to different interventions in different zones: structures in intermix areas, surrounded by vegetation, benefit most from fuel treatments, while those in interface zones, where flammable vegetation is scarcer, are better served by home hardening measures such as fire-resistant roofs, siding, decks, screens, and unburnable buffers. Given that wildfire activity tends to occur where exposure exceeds 60, the authors recommend concentrating resources on the top two risk categories. Where national products like wildfirerisk.org rely on generalized fuel layers that can overestimate unburnable area in WUI zones and are updated too slowly for annually burning landscapes, this approach trades complexity for transparency, using local data that governments already collect. The underlying wildfire exposure toolbox is publicly available on GitHub, and the team argues that empowered residents, equipped with accurate, current, neighborhood-scale information, are crucial to reducing wildfire risk, particularly where communities may prefer action without government intervention.