Phase Transitions in Collective Damage of Civil Structures under Natural Hazards
Pith reviewed 2026-05-22 11:27 UTC · model grok-4.3
The pith
As hazard intensity increases, collective structural damage in cities shifts abruptly from safe to damaged states, analogous to a first-order phase transition.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
As hazard intensity increases, the system can shift abruptly from a largely safe to a largely damaged state, analogous to a first-order phase transition in statistical physics. Higher diversity in the building portfolio smooths this transition, but multiscale damage clustering traps the system in an extended critical-like regime, suppressing the emergence of a more predictable disordered phase. These patterns are captured by a random-field Ising model in which the external field, disorder strength, and temperature map to effective hazard demand, structural diversity, and modeling uncertainty. Application to real inventories reveals that widely used engineering modeling practices can shift a.
What carries the argument
The random-field Ising model, with hazard demand mapped to the external field, building diversity to disorder strength, and uncertainty to temperature, that captures collective interactions among damaged and undamaged structures.
If this is right
- Damage can switch from mostly safe to mostly damaged in an abrupt jump rather than a gradual ramp as hazard intensity rises.
- Greater variety among buildings reduces the sharpness of the transition and the size of the jump.
- Multiscale clustering of damage keeps the system inside a critical-like regime for a wider range of hazard levels.
- Common engineering modeling choices can move a city between synchronized and volatile damage regimes.
- Exceedance-based risk metrics can be biased by up to 50 percent under moderate earthquakes, producing large gaps in estimated repair costs.
Where Pith is reading between the lines
- The same phase-transition framing could be tested on other networked infrastructure such as transportation or power systems.
- Risk-assessment codes might need to add diagnostics for critical regimes when damage clusters span multiple scales.
- Urban planning that increases building diversity could be evaluated for its effect on reducing abrupt damage thresholds.
- Observational datasets from past events could be reanalyzed to check whether predicted critical regimes appear at the city scale.
Load-bearing premise
The collective process of structural damage across many buildings can be represented accurately by a random-field Ising model using the stated mappings of hazard, diversity, and uncertainty.
What would settle it
Post-event damage surveys from a moderate earthquake in a real city that show smooth, gradual transitions without evidence of an extended critical regime driven by multiscale clustering.
Figures
read the original abstract
The fate of cities under natural hazards depends not only on hazard intensity but also on the coupling of structural damage, a collective process that remains poorly understood. Here we show that urban structural damage exhibits phase-transition phenomena. As hazard intensity increases, the system can shift abruptly from a largely safe to a largely damaged state, analogous to a first-order phase transition in statistical physics. Higher diversity in the building portfolio smooths this transition, but multiscale damage clustering traps the system in an extended critical-like regime (analogous to a Griffiths phase), suppressing the emergence of a more predictable disordered (Gaussian) phase. These phenomenological patterns are characterized by a random-field Ising model, with the external field, disorder strength, and temperature interpreted as the effective hazard demand, structural diversity, and modeling uncertainty, respectively. Applying this framework to real urban inventories reveals that widely used engineering modeling practices can shift urban damage patterns between synchronized and volatile regimes, systematically biasing exceedance-based risk metrics by up to 50% under moderate earthquakes ($M_w \approx 5.5$--$6.0$), equivalent to a several-fold gap in repair costs. This phase-aware description turns the collective behavior of civil infrastructure damage into actionable diagnostics for urban risk assessment and planning.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript claims that collective structural damage in urban areas under natural hazards exhibits phase-transition behavior, with abrupt shifts from safe to largely damaged states as hazard intensity increases (analogous to a first-order transition). Higher building-portfolio diversity smooths the transition, while multiscale damage clustering extends a critical-like regime (Griffiths phase) that suppresses a predictable disordered phase. These patterns are characterized via a random-field Ising model in which the external field, disorder strength, and temperature are interpreted as effective hazard demand, structural diversity, and modeling uncertainty, respectively. Application to real urban inventories indicates that common engineering modeling choices can shift damage patterns between synchronized and volatile regimes, biasing exceedance-based risk metrics by up to 50% for moderate earthquakes (Mw ≈ 5.5–6.0).
Significance. If the RFIM mapping is shown to follow from mechanical principles rather than post-hoc fitting and the reported bias is reproducible on independent inventories, the work would supply a statistically grounded diagnostic for collective damage that could refine urban risk metrics and planning. The explicit link between modeling practices and systematic bias in repair-cost estimates is a concrete, actionable contribution.
major comments (2)
- [Model formulation and parameter interpretation] The central claim that abrupt transitions and Griffiths-phase trapping emerge from the physics of coupled structures rests on the RFIM Hamiltonian faithfully representing load redistribution, shared foundations, or correlated ground motion. The manuscript interprets the external field, disorder, and temperature directly as hazard, diversity, and uncertainty without a bottom-up coarse-graining derivation from finite-element failure propagation or similar mechanics; this mapping is load-bearing for the assertion that the observed phase behavior is not a model artifact.
- [Application to real urban inventories and risk-metric bias] The 50% bias result and the demonstration that engineering practices shift the system between synchronized and volatile regimes are obtained by applying the fitted RFIM to real inventories. If the same damage data are used both to calibrate the effective parameters and to identify the phase-transition signatures, the risk-metric bias claim risks circularity; an independent validation set or out-of-sample test is required to establish that the bias is a genuine prediction rather than a fitted description.
minor comments (2)
- [Methods] Notation for the effective temperature (modeling uncertainty) and its relation to the disorder strength should be clarified with an explicit equation or table entry, as the current interpretation leaves open whether temperature is held fixed or varied with hazard intensity.
- [Introduction] The abstract states that the patterns are 'phenomenological'; the main text should explicitly state whether any interaction terms in the RFIM are derived from mechanics or postulated to reproduce the desired phase diagram.
Simulated Author's Rebuttal
We thank the referee for the constructive report and for identifying two key areas where the manuscript's claims require additional clarification. We address each major comment below and indicate the revisions we will make.
read point-by-point responses
-
Referee: [Model formulation and parameter interpretation] The central claim that abrupt transitions and Griffiths-phase trapping emerge from the physics of coupled structures rests on the RFIM Hamiltonian faithfully representing load redistribution, shared foundations, or correlated ground motion. The manuscript interprets the external field, disorder, and temperature directly as hazard, diversity, and uncertainty without a bottom-up coarse-graining derivation from finite-element failure propagation or similar mechanics; this mapping is load-bearing for the assertion that the observed phase behavior is not a model artifact.
Authors: We agree that the RFIM is employed as an effective, phenomenological description rather than a direct coarse-graining from detailed mechanical models. The manuscript already qualifies the patterns as 'phenomenological' and motivates the parameter mapping through physical analogy to structural diversity and uncertainty. However, the referee is correct that a stronger defense would benefit from explicit discussion of this limitation. We will revise the model-formulation section to state clearly that the RFIM serves as a minimal effective model capturing collective statistics, not a first-principles derivation, and we will add a paragraph outlining possible routes for future bottom-up validation against finite-element simulations of small building clusters. revision: partial
-
Referee: [Application to real urban inventories and risk-metric bias] The 50% bias result and the demonstration that engineering practices shift the system between synchronized and volatile regimes are obtained by applying the fitted RFIM to real inventories. If the same damage data are used both to calibrate the effective parameters and to identify the phase-transition signatures, the risk-metric bias claim risks circularity; an independent validation set or out-of-sample test is required to establish that the bias is a genuine prediction rather than a fitted description.
Authors: The referee correctly identifies a potential circularity concern. In the current analysis the effective parameters were obtained from a combination of literature values for building diversity and uncertainty together with a calibration subset of the inventory, while the phase-transition diagnostics and bias quantification were performed on the remaining data. Nevertheless, to eliminate any ambiguity we will add an explicit out-of-sample test using a held-out portion of the inventory and, where possible, a second independent urban dataset. We will also report the sensitivity of the 50% bias figure to the choice of calibration subset. revision: yes
Circularity Check
No significant circularity; RFIM serves as interpretive characterization rather than self-referential derivation.
full rationale
The paper applies the random-field Ising model to characterize observed damage patterns in urban inventories, interpreting its parameters phenomenologically as effective hazard demand, structural diversity, and modeling uncertainty. The abstract explicitly frames the phase-transition and Griffiths-phase claims as patterns first identified in the data and then described by the model. No equations or steps reduce the central results to fitted inputs renamed as predictions, self-citations that bear the load of uniqueness, or ansatzes smuggled from prior author work. The derivation remains self-contained: the model is an external analogy from statistical physics used to organize empirical findings, without the target phase behavior being presupposed in the parameter mappings or data calibration in a circular manner.
Axiom & Free-Parameter Ledger
free parameters (3)
- external field
- disorder strength
- temperature
axioms (1)
- domain assumption Collective damage in civil structures obeys the statistics of a random-field Ising model.
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/AbsoluteFloorClosure.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We cast spatial heterogeneity among cities in the random-field Ising model (RFIM), whose Hamiltonian is H({s}) = −∑_i (H + h_i) s_i − ∑_{i<j} J_{ij} s_i s_j
-
IndisputableMonolith/Cost/FunctionalEquation.leanJ_uniquely_calibrated_via_higher_derivative unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
the external field, disorder strength, and temperature interpreted as the effective hazard demand, structural diversity, and modeling uncertainty, respectively
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Norman A. Abrahamson, Walter J. Silva, and Ronnie Kamai. Summary of the ASK14 Ground Mo- tion Relation for Active Crustal Regions.Earthquake Spectra, 30(3):1025–1055, August 2014. ISSN 8755-2930, 1944-8201. doi: 10.1193/070913EQS198M. URLhttps://journals.sagepub.com/doi/ 10.1193/070913EQS198M
-
[2]
Timothy D. Ancheta, Robert B. Darragh, Jonathan P. Stewart, Emel Seyhan, Walter J. Silva, Brian S.- J. Chiou, Katie E. Wooddell, Robert W. Graves, Albert R. Kottke, David M. Boore, Tadahiro Kishida, and Jennifer L. Donahue. NGA-West2 Database.Earthquake Spectra, 30(3):989–1005, August 2014. ISSN 8755-2930, 1944-8201. doi: 10.1193/070913EQS197M. URLhttps:/...
-
[3]
Seismic performance assessment of buildings, volume 1 – methodology, second edition
Applied Technology Council. Seismic performance assessment of buildings, volume 1 – methodology, second edition. Technical Report FEMA P-58-1, Federal Emergency Management Agency (FEMA), Washington, D.C., December 2018. URLhttps://www.atcouncil.org/. Prepared by the Applied TechnologyCouncilfortheFederalEmergencyManagementAgency.ProjectOfficers: MichaelMa...
work page 2018
-
[4]
Geoff Boeing. Modeling and Analyzing Urban Networks and Amenities With OSMnx.Geographical Analysis, 57(4):567–577, October 2025. ISSN 0016-7363, 1538-4632. doi: 10.1111/gean.70009. URL https://onlinelibrary.wiley.com/doi/10.1111/gean.70009
-
[5]
David M. Boore, Jonathan P. Stewart, Emel Seyhan, and Gail M. Atkinson. NGA-West2 Equations for Predicting PGA, PGV, and 5% Damped PSA for Shallow Crustal Earthquakes.Earthquake Spectra, 30(3):1057–1085, August 2014. ISSN 8755-2930, 1944-8201. doi: 10.1193/070113EQS184M. URL https://journals.sagepub.com/doi/10.1193/070113EQS184M
-
[6]
Abrahamson, Linda Al Atik, Timothy D
Yousef Bozorgnia, Norman A. Abrahamson, Linda Al Atik, Timothy D. Ancheta, Gail M. Atkinson, Jack W. Baker, Annemarie Baltay, David M. Boore, Kenneth W. Campbell, Brian S.J. Chiou, Robert Darragh, Steve Day, Jennifer Donahue, Robert W. Graves, Nick Gregor, Thomas Hanks, I. M. Idriss, Ronnie Kamai, Tadahiro Kishida, Albert Kottke, Stephen A. Mahin, Sanaz R...
-
[7]
PCA Meets RG.Journal of Statistical Physics, 167(3-4):462–475, May 2017
Serena Bradde and William Bialek. PCA Meets RG.Journal of Statistical Physics, 167(3-4):462–475, May 2017. ISSN 0022-4715, 1572-9613. doi: 10.1007/s10955-017-1770-6. URLhttp://link.springer. com/10.1007/s10955-017-1770-6
-
[8]
Kenneth W. Campbell and Yousef Bozorgnia. NGA-West2 Ground Motion Model for the Average Hori- zontal Components of PGA, PGV, and 5% Damped Linear Acceleration Response Spectra.Earthquake Spectra, 30(3):1087–1115, August 2014. ISSN 8755-2930, 1944-8201. doi: 10.1193/062913EQS175M. URLhttps://journals.sagepub.com/doi/10.1193/062913EQS175M
-
[9]
Xu-Yang Cao, De-Cheng Feng, and Michael Beer. A KDE-based non-parametric cloud approach for efficient seismic fragility estimation of structures under non-stationary excitation.Mechanical Systems and Signal Processing, 205:110873, December 2023. ISSN 08883270. doi: 10.1016/j.ymssp.2023.110873. URLhttps://linkinghub.elsevier.com/retrieve/pii/S0888327023007811
-
[10]
Brian S.-J. Chiou and Robert R. Youngs. Update of the Chiou and Youngs NGA Model for the Average Horizontal Component of Peak Ground Motion and Response Spectra.Earthquake Spectra, 30(3):1117–1153, August 2014. ISSN 8755-2930, 1944-8201. doi: 10.1193/072813EQS219M. URL https://journals.sagepub.com/doi/10.1193/072813EQS219M
-
[11]
Analysis neighborhoods.https://data.sfgov.org/dataset/ Analysis-Neighborhoods/p5b7-5n3h, 2024
City and County of San Francisco. Analysis neighborhoods.https://data.sfgov.org/dataset/ Analysis-Neighborhoods/p5b7-5n3h, 2024. Accessed October 21, 2025
work page 2024
-
[12]
Hazus Earthquake Model Technical Manual: Hazus 6.1
Federal Emergency Management Agency (FEMA). Hazus Earthquake Model Technical Manual: Hazus 6.1. Technical report, Federal Emergency Management Agency, Washington, D.C., July 2024. URL https://www.fema.gov/flood-maps/products-tools/hazus
work page 2024
-
[13]
Pablo Heresi and Eduardo Miranda. Structure-to-structure damage correlation for scenario-based re- gional seismic risk assessment.Structural Safety, 95, March 2022. ISSN 01674730. doi: 10.1016/j. strusafe.2021.102155. Publisher: Elsevier B.V
work page doi:10.1016/j 2022
-
[14]
Wenjian Hu, Rajiv R. P. Singh, and Richard T. Scalettar. Discovering phases, phase transitions, and crossovers through unsupervised machine learning: A critical examination.Physical Review E, 95(6): 062122, June 2017. ISSN 2470-0045, 2470-0053. doi: 10.1103/PhysRevE.95.062122. URLhttps: //link.aps.org/doi/10.1103/PhysRevE.95.062122. 78
-
[15]
Nirmal Jayaram, Jack W Baker, N Jayaram, and J W Baker. Correlation model for spatially dis- tributed ground-motion intensities.Earthquake Engineering & Structural Dynamics, 38(15):1687–1708, 12 2009. ISSN 1096-9845. doi: 10.1002/EQE.922. URLhttps://onlinelibrary.wiley.com/ doi/full/10.1002/eqe.922https://onlinelibrary.wiley.com/doi/abs/10.1002/eqe.922htt...
-
[16]
Chulyoung Kang, Oh Sung Kwon, and Junho Song. Evaluation of correlation between engineering demand parameters of structures for seismic system reliability analysis.Structural Safety, 93, 11 2021. ISSN 01674730. doi: 10.1016/j.strusafe.2021.102133
-
[17]
Bahrampouri, and Renmin Pretell
Albert Kottke, Greg Lavrentiadis, Nasser Marafi, M. Bahrampouri, and Renmin Pretell. arkot- tke/pygmm: v0.8.0 (2025-07-24), 2025. URLhttps://doi.org/10.5281/zenodo.16422451
-
[18]
Xinzheng Lu, Frank McKenna, Qingle Cheng, Zhen Xu, Xiang Zeng, and Stephen A Mahin. An open- source framework for regional earthquake loss estimation using the city-scale nonlinear time history analysis.Earthquake Spectra, 36(2):806–831, May 2020. ISSN 8755-2930, 1944-8201. doi: 10.1177/ 8755293019891724. URLhttps://journals.sagepub.com/doi/10.1177/875529...
-
[19]
Chu Mai, Katerina Konakli, and Bruno Sudret. Seismic fragility curves for structures using non- parametric representations.Frontiers of Structural and Civil Engineering, 11(2):169–186, June 2017. ISSN 2095-2430, 2095-2449. doi: 10.1007/s11709-017-0385-y. URLhttp://link.springer.com/10. 1007/s11709-017-0385-y
-
[20]
Frank McKenna. OpenSees: A Framework for Earthquake Engineering Simulation.Computing in Science & Engineering, 13(4):58–66, 2011
work page 2011
-
[21]
Frank McKenna, Stevan Gavrilovic, Jinyan Zhao, Kuanshi Zhong, Adam Zsarnoczay, Barbaros Cetiner, Sina Naeimi, Sang ri Yi, Akash Bangalore Satish, Amin Packzad, Pedro Arduino, and Wael Elhaddad. NHERI-SimCenter/R2DTool: Version 5.5.0 (v5.5.0), 2025. URLhttps://doi.org/10.5281/zenodo. 16884221. NHERI SimCenter, Rapid Regional Resilience Determination Tool (R2DTool)
-
[22]
Long-range Ising model for regional-scale seismic risk analysis
Sebin Oh, Sangri Yi, and Ziqi Wang. Long-range Ising model for regional-scale seismic risk analysis. Earthquake Engineering & Structural Dynamics, 53(12):3904–3923, 2024
work page 2024
-
[23]
Planet dump retrieved from https://planet.osm.org .https://www
OpenStreetMap contributors. Planet dump retrieved from https://planet.osm.org .https://www. openstreetmap.org, 2017. 79
work page 2017
-
[24]
Eric M. Thompson. An updated Vs30 map for california with geologic and topographic constraints (ver. 2.0, july 2022), 2018. URLhttps://doi.org/10.5066/F7JQ108S. U.S. Geological Survey data release
-
[25]
Dimitrios Vamvatsikos and C. Allin Cornell. Incremental dynamic analysis.Earthquake Engineering & Structural Dynamics, 31(3):491–514, March 2002. ISSN 0098-8847, 1096-9845. doi: 10.1002/eqe.141. URLhttps://onlinelibrary.wiley.com/doi/10.1002/eqe.141
-
[26]
Lei Wang. Discovering phase transitions with unsupervised learning.Physical Review B, 94(19):195105, November 2016. ISSN 2469-9950, 2469-9969. doi: 10.1103/PhysRevB.94.195105. URLhttps://link. aps.org/doi/10.1103/PhysRevB.94.195105
-
[27]
Donald L. Wells and Kevin J. Coppersmith. New empirical relationships among magnitude, rup- ture length, rupture width, rupture area, and surface displacement.Bulletin of the Seismologi- cal Society of America, 84(4):974–1002, August 1994. ISSN 1943-3573, 0037-1106. doi: 10.1785/ BSSA0840040974. URLhttps://pubs.geoscienceworld.org/bssa/article/84/4/974/11...
work page 1994
-
[28]
MengjieXiang, JiaxuShen, ZekunXu, andJunChen. Structure-to-structureseismicdamagecorrelation model.Earthquake Engineering & Structural Dynamics, 53(10):3205–3229, 2024. ISSN 1096-9845. doi: 10.1002/eqe.4172. URLhttps://onlinelibrary.wiley.com/doi/abs/10.1002/eqe.4172. _eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/eqe.4172
-
[29]
Tian You and Solomon Tesfamariam. Spatial correlation in building seismic performance for regional resilience assessment.Resilient Cities and Structures, 3(2):57–65, June 2024. ISSN 27727416. doi: 10.1016/j.rcns.2024.06.004. URLhttps://linkinghub.elsevier.com/retrieve/pii/ S2772741624000267
-
[30]
Jian Zhong, Sien Zhou, Hao Wang, and Huimin Hu. Regional seismic fragility of bridge network derived by covariance matrix model of bridge portfolios.Engineering Structures, 309, June 2024. ISSN 18737323. doi: 10.1016/j.engstruct.2024.118035. Publisher: Elsevier Ltd
-
[31]
Minjie Zhu, Frank McKenna, and Michael H. Scott. OpenSeesPy: Python library for the OpenSees finite element framework.SoftwareX, 7:6–11, January 2018. ISSN 23527110. doi: 10.1016/j.softx.2017.10.009. URLhttps://linkinghub.elsevier.com/retrieve/pii/S2352711017300584. 80
-
[32]
SimCenter Earthquake Testbed: San Francisco, CA, 2023
Adam Zsarnóczay, Wael Elhaddad, Barbaros Cetiner, Kuanshi Zhong, Frank McKenna, and Gre- gory Deierlein. SimCenter Earthquake Testbed: San Francisco, CA, 2023. URLhttps://www. designsafe-ci.org/data/browser/public/designsafe.storage.published/PRJ-3880. 81
work page 2023
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.