Stability and Resilience—A Systematic Approach

: Stability and resilience are two crucial concepts to the proper functioning and understanding of the behavior of both natural and man-made systems exposed to perturbations and change. However, although the two have covered a similar territory within dynamic systems, the terminology and applications differ signiﬁcantly. This paper presents a critical analysis of the two concepts by ﬁrst collating the wealth of modern stability concept literature within dynamics systems and then linking it to resilience thinking, deﬁned as adaptation where the system has the ability to respond perturbations and change through passive and active feedback structures. A lumped mass and simple pendulum, two simple linear and nonlinear dynamic systems following a state-space approach from modern control systems theory, are used to support the analysis and application. The research ﬁndings reveal that the two overarching categories of engineering resilience and socio-ecological resilience (extended ecological resilience) are in fact a reinvention of a closed-loop system dynamic stability with different types of active feedback mechanisms. Additionally, structural stability describes some vital aspects of social–ecological resilience such as critical thresholds where, under change, a system loses the ability to return to the starting form or move to another suitable form through active feedback mechanisms or direct management actions.


Introduction
Stability and resilience are two crucial concepts to the proper functioning and understanding of the behavior of both natural and man-made systems exposed to perturbations and change.However, although they cover a similar territory within dynamic systems, the terminology and applications differ significantly.The concept of stability, initially referred to as classical stability and significantly influenced by celestial mechanics, was first employed in 1749 by the Swiss mathematician Leonhard Euler in the context engineering mechanics-statics to describe the equilibrium of columns as rigid bodies under critical buckling load [1][2][3].In 1892, the Russian mathematician Aleksandr Mikhailovich Lyapunov produced the seminal PhD thesis: A general task about the stability of motion introducing the classical stability strand through a precise mathematical formulation.This is considered to be the beginning of the modern stability concept, which has profoundly contributed to the development of modern control systems theory.Since then, the stability concept has evolved greatly, with applications to many other disciplines including economics, numerical analysis, quantum mechanics, nuclear physics, and control systems theory through a unified, consistent, and systematic approach and terminology [1].A general definition for modern stability might either fall under a system's resistance to change (structural stability) or a system's state tendency to return to its initial state in response to a perturbation (dynamic stability) [4].
In addition, there is another strand of the stability concept, first introduced by Odum [5] as an ecological stability concept in the context of ecology and relevant sci-ences, which has benefited less from systematic and empirical work and its conceptual integrity remains fragmented with no consensus on its definition [6].Grimm and Wissel [7] give an exhaustive review of the ecological stability terminology and conclude that ecological stability within a certain ecological context can be fully or partially described by six properties of constancy, resilience (in a narrow sense of returning to the equilibrium after a temporary disturbance), persistence, resistance, elasticity and the domain of attraction.These, on their own, are vague and conflicting notions.
On the other hand, the word resilience originates from the Latin word 'resiliere', which translates as 'bounce back'.This term was first used in 1807 by the English physicist Thomas Young to explain elastic deformation within the context of material sciences [8,9].As an ecological concept resilience first appeared in academic parlance with Holling and was subsequently divided into two categories: ecological resilience and engineering resilience [10,11].The former, which is wider in scope, is defined by the dynamics far from any equilibrium steady-state and is measured by the magnitude of perturbation (change) bearable by a system before it flips to another stability domain (flipping to an irreversible or hard to reverse stability domain).The new system structure is achieved through changes in the variables and processes that govern the system.The latter, engineering resilience, is more limited in scope and is defined by the dynamics close to an equilibrium steady-state.Engineering resilience is measured by the resistance to perturbation, and the speed of return to an equilibrium steady-state after a disruption event takes place [11].Although there has been substantial evolution to the concept of resilience over the past five decades and its extension to applications beyond ecological systems, it is still mostly approached randomly and has benefited less from systematic research [12].
In contrast to the historical precedence and academic prevalence of the term stability (the classical stability strand, currently branded as the modern stability concept) it covers a similar sphere to the term resilience.Although there is a well-established and systematic body of knowledge regarding stability, this has not been applied to resilience thinking, which largely remains scattered and is less ubiquitous compared to stability.Figure 1 shows the publications concerned with stability, with engineering and related sciences on the top of the list.The journals with the most publications are the Proceedings of The IEEE Conference on Decision and Control and IEEE Transactions on Automatic Control, while ecological-related disciplines remain at a considerably lower level.On the contrary, publications related to resilience (Figure 2) are dominant in the social sciences journals, with the greatest number of articles appearing in Sustainability Switzerland and Ecology and Society [13].
Buildings 2022, 12, x FOR PEER REVIEW 2 of 25 In addition, there is another strand of the stability concept, first introduced by Odum [5] as an ecological stability concept in the context of ecology and relevant sciences, which has benefited less from systematic and empirical work and its conceptual integrity remains fragmented with no consensus on its definition [6].Grimm and Wissel [7] give an exhaustive review of the ecological stability terminology and conclude that ecological stability within a certain ecological context can be fully or partially described by six properties of constancy, resilience (in a narrow sense of returning to the equilibrium after a temporary disturbance), persistence, resistance, elasticity and the domain of attraction.These, on their own, are vague and conflicting notions.
On the other hand, the word resilience originates from the Latin word 'resiliere', which translates as 'bounce back'.This term was first used in 1807 by the English physicist Thomas Young to explain elastic deformation within the context of material sciences [8,9].As an ecological concept resilience first appeared in academic parlance with Holling and was subsequently divided into two categories: ecological resilience and engineering resilience [10,11].The former, which is wider in scope, is defined by the dynamics far from any equilibrium steady-state and is measured by the magnitude of perturbation (change) bearable by a system before it flips to another stability domain (flipping to an irreversible or hard to reverse stability domain).The new system structure is achieved through changes in the variables and processes that govern the system.The latter, engineering resilience, is more limited in scope and is defined by the dynamics close to an equilibrium steady-state.Engineering resilience is measured by the resistance to perturbation, and the speed of return to an equilibrium steady-state after a disruption event takes place [11].Although there has been substantial evolution to the concept of resilience over the past five decades and its extension to applications beyond ecological systems, it is still mostly approached randomly and has benefited less from systematic research [12].
In contrast to the historical precedence and academic prevalence of the term stability (the classical stability strand, currently branded as the modern stability concept) it covers a similar sphere to the term resilience.Although there is a well-established and systematic body of knowledge regarding stability, this has not been applied to resilience thinking, which largely remains scattered and is less ubiquitous compared to stability.Figure 1 shows the publications concerned with stability, with engineering and related sciences on the top of the list.The journals with the most publications are the Proceedings of The IEEE Conference on Decision and Control and IEEE Transactions on Automatic Control, while ecological-related disciplines remain at a considerably lower level.On the contrary, publications related to resilience (Figure 2) are dominant in the social sciences journals, with the greatest number of articles appearing in Sustainability Switzerland and Ecology and Society [13].Some researchers have tried to apply the concept of ecological stability to resilience thinking, where both concepts often lack a unified and systematic conceptual foundation.Figure 3 lists the publication counts of articles that mention both the word "stability" and "resilience" in their titles; here, the environmental and social science domains are the dominant disciplines and Ecological Indicators is the top journal [13].Van Meerbeek et al. [14] conducted a systematic review of stability versus resilience in the ecological context and argue that ecological stability is the overarching concept while resilience along with recovery, tolerance, and latitude are the constituent concepts.This view is also shared in other similar theoretical studies, either in full [15,16] or in part [17].However, the main issue in these studies is the narrow and diluted treatment of resilience, equating it only with a system's fixed rate of return to equilibrium without considering the fact that resilience can be achieved through both the inherent system characteristics and management that can be thought of as adaptation.Additionally, the treatment of stability in these studies mostly omits the modern strand of stability and tend to concentrate on the ecological strand of stability, where there is no consensus on its definition.On the other hand, some researchers who have tried to apply the modern strand of stability to the resilience thinking have either covered a sub-class of stability [10,18] or have considered resilience in a diluted and incomplete fashion [18].Some researchers have tried to apply the concept of ecological stability to resilience thinking, where both concepts often lack a unified and systematic conceptual foundation.Figure 3 lists the publication counts of articles that mention both the word "stability" and "resilience" in their titles; here, the environmental and social science domains are the dominant disciplines and Ecological Indicators is the top journal [13].Van Meerbeek et al. [14] conducted a systematic review of stability versus resilience in the ecological context and argue that ecological stability is the overarching concept while resilience along with recovery, tolerance, and latitude are the constituent concepts.This view is also shared in other similar theoretical studies, either in full [15,16] or in part [17].However, the main issue in these studies is the narrow and diluted treatment of resilience, equating it only with a system's fixed rate of return to equilibrium without considering the fact that resilience can be achieved through both the inherent system characteristics and management that can be thought of as adaptation.Additionally, the treatment of stability in these studies mostly omits the modern strand of stability and tend to concentrate on the ecological strand of stability, where there is no consensus on its definition.On the other hand, some researchers who have tried to apply the modern strand of stability to the resilience thinking have either covered a sub-class of stability [10,18] or have considered resilience in a diluted and incomplete fashion [18].
This paper attempts to apply the modern strand of the stability concept, with its established and systematic terminology from modern control systems theory to resilience thinking in the form of adaptation [19].Such a methodology will provide a much needed and original contribution to the body of the knowledge.The study is structured as follows.Section 2 introduces the methodological framework and case examples for supporting the arguments put forward in the subsequent sections.Section 3 provides a state-ofthe-art review of concepts of stability in dynamic systems from a control systems theory perspective.In Section 4 stability is categorized into two broad divisions of dynamic and structural stabilities.A comprehensive conceptual framework is introduced in Section 5 that connects the concepts of modern stability and their application to resilience thinking, supported by illustrations of two simple linear or nonlinear dynamic systems.Lastly, Section 6 presents a summary and discussion.This paper attempts to apply the modern strand of the stability concept, with its tablished and systematic terminology from modern control systems theory to resilien thinking in the form of adaptation [19].Such a methodology will provide a much need and original contribution to the body of the knowledge.The study is structured as follow Section 2 introduces the methodological framework and case examples for supporting t arguments put forward in the subsequent sections.Section 3 provides a state-of-thereview of concepts of stability in dynamic systems from a control systems theory persp tive.In Section 4 stability is categorized into two broad divisions of dynamic and stru tural stabilities.A comprehensive conceptual framework is introduced in Section 5 th connects the concepts of modern stability and their application to resilience thinking, su ported by illustrations of two simple linear or nonlinear dynamic systems.Lastly, Secti 6 presents a summary and discussion.

Methodology and Analysis Tools
This paper utilizes a state-space approach from the modern control systems theo where every system has inputs (controls), states, and outputs (responses).Such an a proach is shaped by the interaction of inputs and outputs across the system bounda (external environment) through a system state that is not directly measurable and obse able.The inputs are those that could potentially be influenced by the system operato (excluding disturbance) while the outputs and states represent the system performan [20].To support and consolidate stability classification and its application to resilien thinking throughout the paper, two simple linear and nonlinear dynamic systems lumped mass (a single degree of freedom structure subjected to a lateral loading) a simple pendulum, are utilized (Table 1; Figure 4).

Methodology and Analysis Tools
This paper utilizes a state-space approach from the modern control systems theory where every system has inputs (controls), states, and outputs (responses).Such an approach is shaped by the interaction of inputs and outputs across the system boundary (external environment) through a system state that is not directly measurable and observable.The inputs are those that could potentially be influenced by the system operators (excluding disturbance) while the outputs and states represent the system performance [20].To support and consolidate stability classification and its application to resilience thinking throughout the paper, two simple linear and nonlinear dynamic systems, a lumped mass (a single degree of freedom structure subjected to a lateral loading) and simple pendulum, are utilized (Table 1; Figure 4).
Alternative frequency form: . y(rate of displacement) u = .. y 0 (seismic acceleration as the input variable) y = x 1 (displacement as the system output) m = 1, mass of the structure; C = damping coefficient; K = 4, elastic stiffness of the structural materials; y 0 = ground displacement due to earthquake; y = displacement of the mass relative to the ground due to earthquake.System natural frequency: System damping ratio expressed as a ratio of the damping coefficient: Eigenvalues/poles: [21] Table 1.Cont.

System System Model System States, Input, and Output System Parameters and Equilibrium States Source
Case example 2: Simple pendulum system-a nonlinear time-invariant second-order dynamic system.

A Review of Stability Concept Terminology
This section introduces the main concepts and tools available under the realm of the modern stability concepts within the context of control systems theory.Table 2 introduces each of these well-established major concepts in system terms.

Stability Concept Definition
A system is an assemblage of functionally related components forming a unity whole to achieve a certain purpose [20,[24][25][26].A system is described by the overarching fundamental variables of state (x), input (u), and output (y).Equations (13-16) describe the linear time invariant and nonlinear time variant continuous systems, where A is the system state or dynamic matrix, B is the input or control matrix, C is the output matrix, and D is

A Review of Stability Concept Terminology
This section introduces the main concepts and tools available under the realm of the modern stability concepts within the context of control systems theory.Table 2 introduces each of these well-established major concepts in system terms.

System
A system is an assemblage of functionally related components forming a unity whole to achieve a certain purpose [20,[24][25][26].A system is described by the overarching fundamental variables of state (x), input (u), and output (y).Equations ( 13)-( 16) describe the linear time invariant and nonlinear time variant continuous systems, where A is the system state or dynamic matrix, B is the input or control matrix, C is the output matrix, and D is the direct transfer or feedforward matrix.f i and g j are scalar arguments of the state, input, and time (t) vectors [27,28]. . . .
Perturbations and change Perturbations and changes cause a system state or a system form to deviate from its initial state or initial form [29,30].
Equilibrium Equilibrium is a system state value at the system state space where the system lies at rest (with zero rate of change).It has either a stable or unstable region around it [31].

• Stable equilibrium
Equilibrium is said to be dynamically stable (attractor) when a system state perturbed by a bounded external perturbation from its equilibrium state remains bounded, including a return to the equilibrium [32].It can also be dynamically stable when the system state matrix/Jacobian eigenvalues/poles lie within the left half-plane (LHP) [33].

• Unstable equilibrium
Equilibrium is said to be dynamically unstable (repellor) when a system state perturbed by a bounded external perturbation from its equilibrium state remains unbounded [32], or when the system state matrix/Jacobian eigenvalues/poles lie within the right half-plane (RHP) [33].

Local stability
Common in nonlinear systems where the system domain of attraction covers only a certain area in the phase space.If the system state is perturbed within the boundary of this domain of attraction, it will asymptotically return to the equilibrium state [34].

Global stability
Common in linear systems.Here, the system domain of attraction covers its entire phase space.If the system state is perturbed, it will converge asymptotically to the equilibrium state [34].

Lyapunov stability
Also known as internal stability, it is applied to autonomous systems.A system is stable about an equilibrium state in the sense of Lyapunov, if all initial values of states starting near the equilibrium state, stay near the equilibrium state [27,35].Lyapunov analysis includes the approaches.

• Lyapunov first (indirect) method
Better suited for smaller perturbations, it is based on the linearization of a nonlinear system around an equilibrium state and subsequently finding its eigenvalues, which need to be equal to or more than zero for the system to be stable [27,35].

Stability Concept Definition
• Lyapunov second (direct) method Provides invaluable insights into the qualitative behavior of nonlinear systems, including their domains of attraction, thresholds, and global stability.The method is directly applied to nonlinear systems without going through any linearization process.In this approach, using a Lyapunov function, if the total energy of a system is continuously dissipating/decreasing along its state trajectory, it will eventually reach a stable equilibrium state where it will remain at rest.However, both stable in the sense of Lyapunov, and asymptotically stable are mathematically stable systems; from an engineering perspective, they are not desirable as the time taken for the system state to return to the equilibrium state is infinite [27,32].

Asymptotic stability
If a system is perturbed from its equilibrium state and it ultimately returns to the equilibrium state, is called asymptotically stable.More precisely, a system is asymptotically stable if it is stable in the sense of Lyapunov and there exists a positive constant for which the system state deviations converge to zero as time goes to infinity [32].

Exponential stability
If a Lyapunov stable system returns to the equilibrium state with an exponential rate of decay, it is termed exponentially stable.From exponential stability, the time required for the perturbed system to return to the equilibrium state can be readily calculated and therefore this is the most desirable property in engineering systems [32,36].

Monotonic stability
As a special case of asymptotic stability, a monotonically stable system is one in which the perturbed state returns to the stable equilibrium state monotonically (through a monotonic decay of perturbations) [37].Monotonic here translates into a constant decreasing trend of perturbations over the entire domain of application.The energy function for a more visual representation of this trend would be a uniform decrease in perturbations without any oscillations (no imaginary part of the system state/Jacobian matrix eigenvalues).It is the absence of oscillations that renders monotonic stability a distinctive case of asymptotic stability.

Finite-time stability (FTS)
A dynamical system state is finite-time stable if its trajectory starts from an initial value within a prescribed bound of the equilibrium state and returns to zero in a finite time frame.Various methods are used for calculating finite time stability measures (such as the settling time function), including Lyapunov theory [38][39][40].A further practical application of FTS is in the analysis and control of input-output dynamical systems [41].It should be mentioned that finite-time stability and exponential stability are similar, closely related concepts.

Bounded-input bounded-output (BIBO) stability
Also referred to as external stability or forced response mode stability.In BIBO stability a system is a casual operator mapping bounded inputs into bounded outputs, which means that a stable system should render a bounded/limited/finite output for all time when a bounded/limited/finite input is exerted.This type of stability is well-suited for application in the context of linear systems [42,43].
Input-state stability (ISS).i.e., bounded-input bounded state (BIBS) A unified approach of internal and external stabilities that is well-suited for application in the nonlinear systems context.Here, for any bounded input there should be a bounded state [27] (Khalil and Grizzle, 2002).In a special case, if the output is equal to the system state, then the system is also BIBO stable (externally stable).

Structural stability
Structural stability is indicative of the robustness of the qualitative behavior of system equilibrium trajectories that are not affected by small internal perturbations (model/parameter uncertainties) or not changed radically (such as the emergence of a new domain of attraction, bifurcating branches, etc.) as a result of the small parameter changes [17,36].Depending on the system focus and context, parameter uncertainty can be studied either in the form of assemblage or certain/individual parameters of interest.

Stability Concepts Classification
It should first be mentioned that any system referred to in this paper is an open-loop system without zero control action (unless stated otherwise) where the flow of information is one way, contrary to the closed-loop control where the flow of information is circular (Figure 5).Given the perturbation event type, modern stability concepts can be divided into two broad categories: dynamic stability, related to perturbations in the form of either (i) system state initial values other than equilibrium or (ii) temporary input disturbance after it is discontinued, and structural stability, which is related to changes in the system in the form of (i) parameter uncertainty that affects system matrix A (Jacobian J) or (ii) parameter uncertainty that moves the equilibrium state including control action.Perturbation type (i) are usually the result of a direct management action, while perturbation type (ii) vary based on the system type and act through the system input port (e.g., the seismic ground acceleration in case example (1) or acceleration in case example (2)).Examples of the changes in the system form (type i) are usually the degradation/obsolescence in the system structure, which are generally countered by active control action (change type ii; case examples 1 and 2).

Stability Concepts Classification
It should first be mentioned that any system referred to in this paper is an open-loop system without zero control action (unless stated otherwise) where the flow of information is one way, contrary to the closed-loop control where the flow of information is circular (Figure 5).Given the perturbation event type, modern stability concepts can be divided into two broad categories: dynamic stability, related to perturbations in the form of either (i) system state initial values other than equilibrium or (ii) temporary input disturbance after it is discontinued, and structural stability, which is related to changes in the system in the form of (i) parameter uncertainty that affects system matrix A (Jacobian J) or (ii) parameter uncertainty that moves the equilibrium state including control action.Perturbation type (i) are usually the result of a direct management action, while perturbation type (ii) vary based on the system type and act through the system input port (e.g., the seismic ground acceleration in case example (1) or acceleration in case example (2)).Examples of the changes in the system form (type i) are usually the degradation/obsolescence in the system structure, which are generally countered by active control action (change type ii; case examples 1 and 2).

Dynamic Stability
Dynamic or dynamical stability entails all those conditions where the time component cannot be ignored and some sort of dynamic process is involved in the system, either through the input (i.e., temporary input disturbance right after it is discontinued), or the system state including the output value (i.e., state initial value other than equilibrium) [45,46].Therefore, dynamic stability is essentially the internal stability of the system and is a measure of the tendency for the system to return (smooth and monotonic or oscillatory) to its equilibrium (initial/original state) value overt time (finite, infinite, or never) after being perturbed (in the absence of any control action) [46].The tendency to return to equilibrium is an identical concept to the system passive control/dissipation/damping, which is further explored in Section 5.In control systems terms (generally for linear and linearized systems), dynamic stability is equivalent to the rate of return of the system's state (and/or settling time for a time-varying system) to equilibrium, which is measured by the dominant eigenvalue horizontal distance (real part) from the imaginary axis in the complex plane on a root locus diagram.A dominant eigenvalue is related to the slow-moving state of the system and is located closest to the imaginary axis which corresponds to the slowest and dominant decay/return rate to the equilibrium [47].Depending on the nature of the system's dominant eigenvalue and its linearity, various types of dynamic stability, explained in Section 3, exist and are outlined in Table 3 and Figure 6.Tools such as the root locus diagram for linear systems and Lyapunov functions, are generally used for nonlinear systems, where no finite rate of return time to the equilibrium state can be extracted and are an important means for assessing dynamic stability.Remark: BIBO and BIBS can be considered a part of the dynamic stability for the period immediately after the temporary input disturbance is discontinued.If a permanent input in the sense of an active control is considered, then it becomes a closed-loop dynamic stability (Figure 5)

Structural Stability
Structural stability indicates that a system model is stable for all model uncertainties until reaching a critical threshold (stability radius, defined as the minimum destabilizing effect of combined system changes/model uncertainty), which leads to instability [48].Loosely speaking, the former definition of structural stability (Section 3), where parameter changes can cause emergence or disappearance of new equilibrium domains is wellpositioned in the context of complex dynamic systems.The latter definition of structural stability is related to changes in the shape of the domain(s) around the existing equilibrium state(s) or simply changes to the dynamics within a single domain of equilibrium.Borrowing terminology from the resilience scenario, the two subject definitions for structural stability can be labeled under general structural stability and specified structural stability, which subsequently determine hard and soft thresholds (unstable points along the system state trajectory) for a dynamic system (Figure 6).
Specified structural stability, also referred to as persistence and endurance [49], inertia [50] (Murdoch, 1970), constancy [51], robustness and resistance [52], is normally used as a measure of the system's functional persistence under uncertainty parameters (inverse of sensitivity) and is generally related to the dynamics within a domain of equilibrium [53][54][55].
To avoid confusion and to treat stability systematically and with precision, mentioning a sub-class or adjective prior to the word 'stability', such as dynamic stability and structural stability, or their low tier sub-classes is a critical consideration.

Structural Stability
Structural stability indicates that a system model is stable for all model uncertainties until reaching a critical threshold (stability radius, defined as the minimum destabilizing effect of combined system changes/model uncertainty), which leads to instability [48].Loosely speaking, the former definition of structural stability (Section 3), where parameter changes can cause emergence or disappearance of new equilibrium domains is well-positioned in the context of complex dynamic systems.The latter definition of structural stability is related to changes in the shape of the domain(s) around the existing equilibrium state(s) or simply changes to the dynamics within a single domain of equilibrium.Borrowing terminology from the resilience scenario, the two subject definitions for structural stability can be labeled under general structural stability and specified structural stability, which subsequently determine hard and soft thresholds (unstable points along the system state trajectory) for a dynamic system (Figure 6).

Application of Stability to Resilience Thinking
Anderies et al. [56] argue that specified resilience (resilience to a known type of perturbation) is equivalent to robustness and the robust control analysis used by Csete and Doyle [57] in the context of biological systems with parameter uncertainties through the application of feedbacks.However, most of the available literature fails to provide a unified, comprehensive, and systematic treatment of the broad stability concepts in terms of perturbation, changes, dynamic and structural stabilities, and relevant control actions versus the two prevailing categories of engineering resilience and socio-ecological resilience, as well as the means to obtain resilience.
Resilience here is defined as adaptation proposed in the authors' earlier work [19] where the system has the ability to respond to perturbations and change through passive and active feedback structures; the system state or system form is returned to a starting position or transitions to another suitable state or form.Simulations from the lumped mass and simple pendulum linear and nonlinear dynamic systems are used to support the modern stability ideas to resilience thinking in the form of adaptation.This is further discussed under two main features of the system state and form return abilities through passive and active feedback mechanisms.The two overarching categories of resilience, engineering resilience, and socio-ecological resilience, are distinguished in the analysis process; the range of perturbations and change in engineering resilience is narrow and predictable (normally within a single domain of attraction), while socio-ecological resilience is broader and unpredictable.

Resilience as the Ability of the System to Return to Initial State
The ability of a dynamic system to respond to perturbations and return the system state to a starting position or another suitable state is a local, or within a domain of attraction, phenomenon that is mostly used in the engineering resilience category.Here, constancy, or a target system performance, usually in the form of a system equilibrium/steady state or output as part of it, is the main objective.Such an objective for engineering resilience is precisely translated in the form of a faster system state return rate as well as larger system resistance to input disturbance, which translates into smaller system state deviations from its equilibrium.Both of these elements can be achieved by various sub-classes of dynamic stability where system has the ability to return to its equilibrium state including the quality and form of the state return trajectory.The dissipative measures built within the system, also referred to as resistance to perturbation and robustness, determine the extent of the deviation from the system equilibrium state under input disturbances.These can be instantaneous, such as impulse/pulse input disturbances (e.g., shocks and impact loading) or permanent, such as step/press input disturbances (e.g., stressors such as climate change), which subsequently change the system form and move its equilibrium state.
The rate of return of the system to equilibrium state or the settling time is determined by the system's dominant eigenvalue or dynamic stability and the system resistance to perturbation is determined by the passive feedback/control that is built into the system.Tools such as the root locus diagram for linear systems (Figure 7) and Lyapunov functions, generally used for nonlinear systems where no finite rate of return time to equilibrium state can be extracted (Figures 8 and 9), are important methods for determining the system state rate of return to equilibrium or, alternatively, the settling time.Based on the nature of the system's dominant eigenvalue and its linearity, various sub-classes of dynamic stability and its applicability level to engineering resilience are summarized in Table 4.As indicated in Table 5 and Figure 10, two identical systems might have different amounts of passive control (damping) but as long as their dominant eigenvalues' real parts are equal, the resulting settling times will be equal irrespective of their behavior (monotonic or oscillatory).In the context of linear systems and a single domain of attraction, within a certain vicinity of the stable equilibrium, both passive control and rate of return to equilibrium increase until an optimum amount of resistance is reached (e.g., critical damping for second-order linear systems).Any further increase in passive control will negatively affect In the context of linear systems and a single domain of attraction, within a certain vicinity of the stable equilibrium, both passive control and rate of return to equilibrium increase until an optimum amount of resistance is reached (e.g., critical damping for secondorder linear systems).Any further increase in passive control will negatively affect the settling time.Thus, engineering resilience is directly proportional to dynamic stability (within an open-loop setting with zero active control action) for a certain interval of passive control (e.g., critical damping) only.Beyond the interval of passive control, any added amount of passive control might curtail the system's ability for innovation or might lock the system into domains with unproductive conditions such as economic recessions.The relationship between passive feedback and settling time, the two elements of engineering resilience, is demonstrated for two simple linear and nonlinear dynamic systems in Table 6 and Figures 11 and 12.A good measure for engineering resilience can be considered as the volume o system's 3-dimenstional state portrait, where the two horizontal axes indicate the sy state deviations under input disturbance (passive control indicator) while the vertical indicates the system settling time (system state rate of return to equilibrium).This c bines both of the engineering resilience elements and the smaller the volume the la the resilience becomes (Figure 13).For an equal value of settling time, monotonic e nential stability renders a slightly higher value of engineering resilience than that of o latory exponential stability given the slightly higher area of the system state trajec around equilibrium of the former compared to the latter (Figure 10).A good measure for engineering resilience can be considered as the volume of the system's 3-dimenstional state portrait, where the two horizontal axes indicate the system state deviations under input disturbance (passive control indicator) while the vertical axis indicates the system settling time (system state rate of return to equilibrium).This combines both of the engineering resilience elements and the smaller the volume the larger the resilience becomes (Figure 13).For an equal value of settling time, monotonic exponential stability renders a slightly higher value of engineering resilience than that of oscillatory exponential stability given the slightly higher area of the system state trajectory around equilibrium of the former compared to the latter (Figure 10).
A good measure for engineering resilience can be considered as the volume of the system's 3-dimenstional state portrait, where the two horizontal axes indicate the system state deviations under input disturbance (passive control indicator) while the vertical axis indicates the system settling time (system state rate of return to equilibrium).This combines both of the engineering resilience elements and the smaller the volume the larger the resilience becomes (Figure 13).For an equal value of settling time, monotonic exponential stability renders a slightly higher value of engineering resilience than that of oscillatory exponential stability given the slightly higher area of the system state trajectory around equilibrium of the former compared to the latter (Figure 10).However, for linear systems the dynamic stability remains constant within the domain of attraction.In the context of nonlinear systems, the dynamic stability critically decreases (lengthy settling times) in the vicinity of unstable equilibrium (instability threshold), which results in smaller values of engineering resilience.Figure 14 depicts how the system state response critically slows down near the unstable equilibrium threshold ( x 1 = 6.264 or x 2 = 3.14 Rad) as does the resistance to perturbation.However, for linear systems the dynamic stability remains constant within the domain of attraction.In the context of nonlinear systems, the dynamic stability critically decreases (lengthy settling times) in the vicinity of unstable equilibrium (instability threshold), which results in smaller values of engineering resilience.Figure 14 depicts how the system state response critically slows down near the unstable equilibrium threshold (  6.264   3.14  as does the resistance to perturbation.The ability of the system to return to its equilibrium state ceases once the passive feedback of the system reaches zero (open-loop system with zero active control) or, alternatively, when the system eigenvalues are on the imaginary axis of the complex plane.By crossing the zero limit (c = 0), which can referred to as a soft threshold, the system stable The ability of the system to return to its equilibrium state ceases once the passive feedback of the system reaches zero (open-loop system with zero active control) or, alternatively, when the system eigenvalues are on the imaginary axis of the complex plane.By crossing the zero limit (c = 0), which can referred to as a soft threshold, the system stable equilibrium state [0, 0] turns into an unstable one (Figure 15).

Resilience as System's Form Return Ability
The ability of a dynamic system to respond to change by returning the system starting position or another suitable form can be a local or global phenomenon a mostly discussed in the socio-ecological resilience category, where maintaining the sy identity and/or avoiding critical thresholds are the main objectives.Critical thresho the system form can be divided into soft and hard thresholds, where the former is and the latter is hard or impossible to reverse.The reversibility of the system form management action under the resilience scenario.It is thoroughly covered by the a feedback/control action under the modern stability concept, which also means th active feedback action changes the system form and subsequently alters the system' sive control measures that define open-loop dynamic stability (Figure 5).Since chan the system form under the engineering resilience definition are limited and known erally within a domain of attraction), the control law for the closed-loop system doe change and is labeled as a robust (active) control.For socio-ecological resilience o other hand, which has a wide and unknown range of system changes, the control l not constant and therefore is categorized as adaptive (active) control (Figure 5).Add ally, there is input disturbance of either an instantaneous (impulse input) or perm (step input) nature that can flip the system state into another stability domain; this, with the system's passive control, are the two main elements that are crucial for meas socio-ecological resilience.
In the context of linear systems and a single domain of attraction and global dyn stability (Case example 1), there is no limit on the input disturbance that leads to ins ity.Although the system is theoretically a globally stable system, performance/spec tion requirements can be imposed to limit the system domain of equilibrium to a c boundary, e.g., maximum story drift/displacement in Case example 1, which pract

Resilience as System's Form Return Ability
The ability of a dynamic system to respond to change by returning the system to a starting position or another suitable form can be a local or global phenomenon and is mostly discussed in the socio-ecological resilience category, where maintaining the system identity and/or avoiding critical thresholds are the main objectives.Critical thresholds in the system form can be divided into soft and hard thresholds, where the former is easy, and the latter is hard or impossible to reverse.The reversibility of the system form is a management action under the resilience scenario.It is thoroughly covered by the active feedback/control action under the modern stability concept, which also means that the active feedback action changes the system form and subsequently alters the system's passive control measures that define open-loop dynamic stability (Figure 5).Since changes to the system form under the engineering resilience definition are limited and known (generally within a domain of attraction), the control law for the closed-loop system does not change and is labeled as a robust (active) control.For socio-ecological resilience on the other hand, which has a wide and unknown range of system changes, the control law is not constant and therefore is categorized as adaptive (active) control (Figure 5).Additionally, there is input disturbance of either an instantaneous (impulse input) or permanent (step input) nature that can flip the system state into another stability domain; this, along with the system's passive control, are the two main elements that are crucial for measuring socio-ecological resilience.
In the context of linear systems and a single domain of attraction and global dynamic stability (Case example 1), there is no limit on the input disturbance that leads to instability.
Although the system is theoretically a globally stable system, performance/specification requirements can be imposed to limit the system domain of equilibrium to a certain boundary, e.g., maximum story drift/displacement in Case example 1, which practically renders local stability.Therefore, socio-ecological resilience tilts toward engineering resilience, which is mostly about the passive control built into the system (Table 7).In the context of nonlinear systems with multiple domains of attraction, e.g., Case example 2, or more complex nonlinear systems with hard thresholds described by fold bifurcations, a higher value of input disturbance that can flip the system into another stability domain along with a higher amount of passive control measures indicate a high level of socio-ecological resilience (Table 7 and Figure 16).resilience, which is mostly about the passive control built into the system (Table 7).In the context of nonlinear systems with multiple domains of attraction, e.g., Case example 2, or more complex nonlinear systems with hard thresholds described by fold bifurcations, a higher value of input disturbance that can flip the system into another stability domain along with a higher amount of passive control measures indicate a high level of socioecological resilience (Table 7 and Figure 16).Contrary to the simple nonlinear systems (Case example 2-Figure 17), for complex nonlinear dynamic systems the relationship between parameters affecting the system matrix A (or Jacobian) (dynamic stability and passive control) and those moving the equilibrium states (e.g., input or active control action) is an area for further exploration; since the subject parameters are often entangled and/or of a double-edged sword nature, there is increased complexity in the trade-offs between engineering resilience vs. socio-ecological resilience.
Buildings 2022, 12, x FOR PEER REVIEW 20 Contrary to the simple nonlinear systems (Case example 2-Figure 17), for com nonlinear dynamic systems the relationship between parameters affecting the system trix A (or Jacobian) (dynamic stability and passive control) and those moving the eq rium states (e.g., input or active control action) is an area for further exploration; sin subject parameters are often entangled and/or of a double-edged sword nature, th increased complexity in the trade-offs between engineering resilience vs. socio-ecol resilience.Bifurcation diagrams are useful tools from nonlinear dynamics and structural s ity that can be utilized for nonlinear dynamic systems with single parameter change can help define hard and soft thresholds for socio-ecological resilience.For a sing rameter system, changes can assess how far a critical threshold is located from the cu equilibrium state, as well as to indicate whether a threshold that has been crossed versible (Figure 18) or irreversible (Figure 19), which corresponds to a soft and har of resilience, respectively [58].A hard loss of resilience is known as hysteresis (pa pendence or path memory) or fold bifurcation in complex nonlinear systems and it s that the system trajectory not only depends on the parameter value but also depen which side the parameter value is approached (e.g., parameter R in the spruce budw model-Figure 19).This occurs where there are the multiple domains of attraction [ Bifurcation diagrams are useful tools from nonlinear dynamics and structural stability that can be utilized for nonlinear dynamic systems with single parameter changes and can help define hard and soft thresholds for socio-ecological resilience.For a single parameter system, changes can assess how far a critical threshold is located from the current equilibrium state, as well as to indicate whether a threshold that has been crossed is reversible (Figure 18) or irreversible (Figure 19), which corresponds to a soft and hard loss of resilience, respectively [58].A hard loss of resilience is known as hysteresis (path dependence or path memory) or fold bifurcation in complex nonlinear systems and it shows that the system trajectory not only depends on the parameter value but also depends on which side the parameter value is approached (e.g., parameter R in the spruce budworm model-Figure 19).This occurs where there are the multiple domains of attraction [27].When system resilience under passive control is not satisfactory, active control action (given that the system is controllable) is required to alter the system form by shifting the system eigenvalues to a desired location and adjusting the system behavior accordingly (Figure 20).For instance, the subject linear system given in Case example 1 is fully controllable, but its current behavior is oscillatory and not satisfactory.A full state active feedback tool can be used to bring the underdamped oscillatory state behavior to a critical damped monotonic behavior (or, in alternative scenarios, to speed up the system settling   When system resilience under passive control is not satisfactory, active control action (given that the system is controllable) is required to alter the system form by shifting the system eigenvalues to a desired location and adjusting the system behavior accordingly (Figure 20).For instance, the subject linear system given in Case example 1 is fully controllable, but its current behavior is oscillatory and not satisfactory.A full state active feedback tool can be used to bring the underdamped oscillatory state behavior to a critical damped monotonic behavior (or, in alternative scenarios, to speed up the system settling When system resilience under passive control is not satisfactory, active control action (given that the system is controllable) is required to alter the system form by shifting the system eigenvalues to a desired location and adjusting the system behavior accordingly (Figure 20).For instance, the subject linear system given in Case example 1 is fully controllable, but its current behavior is oscillatory and not satisfactory.A full state active feedback tool can be used to bring the underdamped oscillatory state behavior to a critical damped monotonic behavior (or, in alternative scenarios, to speed up the system settling time or stabilize the system in the case of instability).This will mean shifting the system poles from the underdamped case (λ 1 = −0.2000+ 1.9900 i, λ 2 = −0.2000− 1.9900 i) to a critically damped one (λ 1 = λ 2 = −2.000).The reference scale (r) is set to zero (zero input/equilibrium state) and instead, the disturbance input is introduced in the form of the initial condition of [20 0] (Figure 20).The control action bears a cost, which needs to be assessed versus the performance improvement it brings about.
Buildings 2022, 12, x FOR PEER REVIEW 22 of 25 time or stabilize the system in the case of instability).This will mean shifting the system poles from the underdamped case ( λ 0.2000 1.9900 i, λ 0.2000 1.9900 i to a critically damped one (λ λ −2.000).The reference scale (r) is set to zero (zero input/equilibrium state) and instead, the disturbance input is introduced in the form of the initial condition of [20 0] (Figure 20).The control action bears a cost, which needs to be assessed versus the performance improvement it brings about.

Discussion and Conclusions
Both dynamic and structural stabilities are principally dependent on the system parameters where the former is about the rate of return to equilibrium and the latter is about the current equilibrium distance from a fold bifurcation point, which means staying in the same domain of equilibrium.Engineering resilience, which focuses on maintaining the system performance under limited and known types of perturbations, is therefore partially covered by the open-loop dynamic stability with zero active control.Dynamic stability without active control action can be seen as a static return rate to equilibrium through the utilization of the passive control system after the perturbation stimulus is removed.Moreover, if the system return behavior and return rate are not satisfactory, the control action is the only means of rectifying the unsatisfactory situation.On the other hand, socio-ecological resilience, which focuses on maintaining the system identity and avoiding critical thresholds, is partially covered by structural stability (distance to soft or hard thresholds) by keeping a static distance from equilibrium.Here, control action is needed to adjust that distance as per the perturbation value and type, which are often unknown and with large variability, and requires a closed-loop dynamic stability with adaptive control as an alternative for socio-ecological resilience.Additionally, control systems tools, such as state controllability, are potential directions for developing tools that

Discussion and Conclusions
Both dynamic and structural stabilities are principally dependent on the system parameters where the former is about the rate of return to equilibrium and the latter is about the current equilibrium distance from a fold bifurcation point, which means staying in the same domain of equilibrium.Engineering resilience, which focuses on maintaining the system performance under limited and known types of perturbations, is therefore partially covered by the open-loop dynamic stability with zero active control.Dynamic stability without active control action can be seen as a static return rate to equilibrium through the utilization of the passive control system after the perturbation stimulus is removed.Moreover, if the system return behavior and return rate are not satisfactory, the control action is the only means of rectifying the unsatisfactory situation.On the other hand, socio-ecological resilience, which focuses on maintaining the system identity and avoiding critical thresholds, is partially covered by structural stability (distance to soft or hard thresholds) by keeping a static distance from equilibrium.Here, control action is needed to adjust that distance as per the perturbation value and type, which are often unknown and with large variability, and requires a closed-loop dynamic stability with adaptive control as an alternative for socio-ecological resilience.Additionally, control systems tools, such as state controllability, are potential directions for developing tools that can systematically describe the reversibility of equilibrium states in complex dynamic systems and the trade-offs between dynamic and structural stabilities.This study reveals that the two major categories of engineering resilience and extended ecological resilience are intrinsically the reinvention of a closed-loop system dynamic stability with different types of active feedback structures.Moreover, structural stability describes key aspects of social-ecological resilience-including critical thresholds where, under change, a system loses the ability to return to the starting form or move to another suitable form through active feedback mechanisms or direct management actions.

Figure 3 .
Figure 3. Publication counts listing both stability and resilience in their titles.Data source: [13].

Figure 3 .
Figure 3. Publication counts listing both stability and resilience in their titles.Data source: [13].

Figure 4 .
Figure 4. Two linear and nonlinear dynamic systems utilized in the study: (a) single degree of freedom structure subjected to a later loading.Source: adapted from Carmichael [21].(b) A schematic diagram of a simple pendulum.

Figure 4 .
Figure 4. Two linear and nonlinear dynamic systems utilized in the study: (a) single degree of freedom structure subjected to a later loading.Source: adapted from Carmichael [21].(b) A schematic diagram of a simple pendulum.

Figure 5 .
Figure 5. Flow of information for: (a) open-loop and (b) closed-loop representation of a dynamic system.Source: adapted from Carmichael [44].

Figure 5 .
Figure 5. Flow of information for: (a) open-loop and (b) closed-loop representation of a dynamic system.Source: adapted from Carmichael [44].

Figure 6 .
Figure 6.Linkages between dynamic stability, structural stability, and active control action-a unified and big picture.

Figure 6 .
Figure 6.Linkages between dynamic stability, structural stability, and active control action-a unified and big picture.

Figure 7 .
Figure 7. Root locus diagram of the lumped mass dynamic system (Case example 1) for four different arrangements of dominant eigenvalue/pole based on the system's internal resistance/damping.

Figure 7 . 25 Figure 8 .
Figure 7. Root locus diagram of the lumped mass dynamic system (Case example 1) for four different arrangements of dominant eigenvalue/pole based on the system's internal resistance/damping.Buildings 2022, 12, x FOR PEER REVIEW 13 of 25

Figure 8 .
Figure 8. Lyapunov function (V), indicative of global asymptotic stability for the lumped mass dynamic system (Case example 1) along with the system phase portraits under [2, 1.5] state initial value.

Figure 8 .
Figure 8. Lyapunov function (V), indicative of global asymptotic stability for the lumped mass dynamic system (Case example 1) along with the system phase portraits under [2, 1.5] state initial value.

Figure 9 .
Figure 9. Lyapunov function (V) indicative of local asymptotic stability for the simple pendulum system (Case example 2) along with the system phase portraits under [−1.5 rad, 4] state initial values.

Figure 9 .
Figure 9. Lyapunov function (V) indicative of local asymptotic stability for the simple pendulum system (Case example 2) along with the system phase portraits under [−1.5 rad, 4] state initial values.

Figure 10 .
Figure 10.System equal dominant eigenvalues rendering equal settling times for a unit impulse in the form of a temporary input disturbance (Case example 1).

Figure 10 .
Figure 10.System equal dominant eigenvalues rendering equal settling times for a unit impulse in the form of a temporary input disturbance (Case example 1).

Figure 13 .
Figure 13.System phase portrait (domain of attraction) 3D representation as an indicator of engineering resilience.The vertical axis indicates settling time for a unit impulse input and [0, 0] initial conditions (Case example 1).

Figure 13 .
Figure 13.System phase portrait (domain of attraction) 3D representation as an indicator of engineering resilience.The vertical axis indicates settling time for a unit impulse input and [0, 0] initial conditions (Case example 1).

Figure 14 .
Figure 14.System state response indicative of the critical slowing down near an unstable equilibrium threshold (Case example 2 with [0, 6.264] initial conditions).

Figure 14 .
Figure 14.System state response indicative of the critical slowing down near an unstable equilibrium threshold (Case example 2 with [0, 6.264] initial conditions).

Buildings 2022 , 18 Figure 15 .
Figure 15.An approximate 3D representation of the equilibrium domain changes due to the s passive control levels.Dominant eigenvalue/pole real part indicate the rate of return to the eq rium (Case example 1).

Figure 15 .
Figure 15.An approximate 3D representation of the equilibrium domain changes due to the system passive control levels.Dominant eigenvalue/pole real part indicate the rate of return to the equilibrium (Case example 1).

Figure 16 .
Figure 16.Phase portrait indicative of the system's soft thresholds disturbed by initial conditions (Case example 2).

Figure 16 .
Figure 16.Phase portrait indicative of the system's soft thresholds disturbed by initial conditions (Case example 2).

Figure 17 .
Figure 17.Permanent input (step input) acting as a parameter (input perturbation) that mov system equilibrium state and affects resilience (Case example 2).

Figure 17 .
Figure 17.Permanent input (step input) acting as a parameter (input perturbation) that moves the system equilibrium state and affects resilience (Case example 2).

Figure 19 .
Figure 19.Spruce Budworm model S-shaped bifurcation for parameter (R) indicating a hard loss of resilience.

Figure 19 .
Figure 19.Spruce Budworm model S-shaped bifurcation for parameter (R) indicating a hard loss of resilience.

Figure 19 .
Figure 19.Spruce Budworm model S-shaped bifurcation for parameter (R) indicating a hard loss of resilience.

Figure 20 .
Figure 20.System state response with and without full state feedback (closed-loop control) for a zero reference/set point (equilibrium state) and [20 0] initial condition (Case example 1).

Figure 20 .
Figure 20.System state response with and without full state feedback (closed-loop control) for a zero reference/set point (equilibrium state) and [20 0] initial condition (Case example 1).

Table 1 .
Dynamic system case examples for stability and resilience thinking.

Table 3 .
Classification of dynamic stability concepts with a focus on the system performance.

Table 4 .
Dynamic stability application to engineering resilience with a focus on the system performance.Infinite Least favorable Nonlinear systems (particularly time-varying), multiple equilibriums and local stabilities Perturbation/change: system changes/parameter uncertainty Lyapunov stability (second method) Infinite Less favorable Nonlinear systems (particularly time-varying), multiple equilibriums and local stabilities Perturbation/change: system changes/parameter uncertainty Lyapunov stability (second method) Infinite Less favorable Mainly nonlinear systems (particularly time-varying), multiple equilibriums and local stabilities Perturbation: state initial value and input temporary disturbance after it is disconnected Lyapunov stability (first method, linearization)

Table 5 .
System dominant eigenvalue's real part as a direct indication of the system settling time (Case example 1).

Table 6 .
Dynamic stability as a measure of engineering resilience for impulse input disturbance (δ(t)) (Case examples 1 and 2).

Table 7 .
Dynamic stability as a measure of engineering resilience for impulse input disturbance (δ(t)) (Case examples 1 and 2).

Table 7 .
Dynamic stability as a measure of engineering resilience for impulse input disturbance (δ(t)) (Case examples 1 and 2).