  Perturbed by natural hazards, community-level infrastructure networks operate
like many-body systems, with behaviors emerging from coupling individual
component dynamics with group correlations and interactions. It follows that we
can borrow methods from statistical physics to study the response of
infrastructure systems to natural disasters. This study aims to construct a
joint probability distribution model to describe the post-hazard state of
infrastructure networks and propose an efficient surrogate model of the joint
distribution for large-scale systems. Specifically, we present maximum entropy
modeling of the regional impact of natural hazards on civil infrastructures.
Provided with the current state of knowledge, the principle of maximum entropy
yields the ``most unbiased`` joint distribution model for the performances of
infrastructures. In the general form, the model can handle multivariate
performance states and higher-order correlations. In a particular yet typical
scenario of binary performance state variables with knowledge of their mean and
pairwise correlation, the joint distribution reduces to the Ising model in
statistical physics. In this context, we propose using a dichotomized Gaussian
model as an efficient surrogate for the maximum entropy model, facilitating the
application to large systems. Using the proposed method, we investigate the
seismic collective behavior of a large-scale road network (with 8,694 nodes and
26,964 links) in San Francisco, showcasing the non-trivial collective behaviors
of infrastructure systems.
