Spatial Differentiation and Elements Influencing Urban Resilience in the Middle Reaches of the Yangtze River under the COVID-19 Pandemic

Since the outbreak of COVID-19 in Wuhan City, Hubei Province, in December 2019, the middle reaches of the Yangtze River became the key areas of the spread of the pandemic and association, and also as the urban economic recovery process after the pandemic eased, it provided an excellent opportunity to research urban resilience. From the viewpoint of urban social-ecological system resilience in public health emergencies, this study comprehensively applies the spatial econometrics, geodetector model, and other methods to investigate the urban resilience level, spatial differentiation, and dominant elements in the middle reaches of the Yangtze River under the impact of the pandemic. -is study would aid in providing a scientific basis for sustainable spatial planning and governance. -e results demonstrated that the urban resilience in the middle reaches of the Yangtze River had notable spatial agglomeration features, eight elements including tertiary industry proportion possessed a robust explanatory power to the spatial differentiation of urban resilience, and the explanatory power was markedly enhanced after the interaction between influential elements of economic and ecological subsystems. -us, to upgrade the system cycle mechanism and augment the endogenous power for urban development, we need to focus more on the flow of innovative elements in central cities, the optimization of ecological and safe spatial patterns in Hunan and Hubei Provinces, and the innovation of sustainable supply chain in the entire region.


Introduction
After the first case of the coronavirus disease 2019 (COVID- 19) was reported in December 2019 in Wuhan City, Hubei Province [1] (hereafter referred to as "pandemic"), the disease spread quickly to other regions, especially the middle reaches of the Yangtze River, which was the most severely affected area being centered on Wuhan City. e operation of socioeconomic orders in each city was seriously affected; the large-scale supply chain was interrupted, reflecting on spatial planning, production, supply modes, and urban governance.
e World Health Organization defined the pandemic as a public health emergency of international concern [2]. Currently, the prevention and control of the pandemic have entered into the normalized stage, and the recovery process of urban social economy has provided an opportunity for the research and data collection of urban resilience. With the background of growing uncertain perturbation factors, such as public health emergency represented by the pandemic, enhancing the adaptive circulating ability of city as a social-ecological system (SES) and conducting effective spatial planning and governance are of great significance for urban sustainable development. e research on the resilience of urban SES after public health emergency and other disaster events has always been a hot topic in public management, economic geography, and urban and rural planning, involving the comprehensive impact of various elements of the SES and primarily focusing on the role of government and market or the valuation of urban social economy recovery. Ma and Ren [3] and Olcott and Oliver [4] conducted relevant research and quantitative analyses on the allocation relation and dynamic stability of resource elements in complex systems. Regarding the COVID-19 pandemic, Wang and Zhao [5], Ma et al. [6], and Li et al. [7] examined the spatial differentiation features of a domestic epidemic from different temporal and spatial scales and revealed the correlation features of such elements as landscape patterns of ecological and construction lands, emergency response management and government regulation, and regional centrality with the pandemic. Besides, Yang et al. [8] proposed a new calculation and application method for economic complex systems that could be used in engineering applications, including population-scale control and economic and market activities. However, considering the data timeliness and accuracy, related studies have primarily explored the spatial correlation and macrocontrol between epidemiological transmission and socioeconomic elements, but the empirical studies on urban resilience in the case of the pandemic and the research into proper governance measures in light of local conditions remain unexplored.
In this study, empirical analysis is performed for the urban resilience in the middle reaches of the Yangtze River under the pandemic based on the existing literature. In addition, this study explores the spatial differentiation and dominant elements of urban resilience using a geodetector model.
is study also provides a scientific basis for the ecological and safe spatial pattern planning and governance and the sustainable innovation and development of the supply chain, as well as enhancing the urban resilience and sustainable development capacity.

Literature Review
Resilience, first originating from Latin word "resilio," implies "returning to the original state" [9], and the concept connotation has experienced development in multidisciplinary context and social transformation and is progressively applied to studies on urban problems. e standpoints of the study evolve in three stages: engineering resilience, ecological resilience, and social-ecological resilience. A recent study defined urban resilience as the ability of the socialecological network of the urban system to sustain or rapidly restore the good operation status in case of imbalance [10]. Under the theoretical framework of SES adaptive governance, it is essential to coordinate the interrelationship between environment, economy, and society to ascertain the resilience management strategy and regulate the state of the complex adaptive system, thereby augmenting the urban resilience and further coping with nonlinear change [11].
Regarding dynamic adaptive governance, Ma implemented a series of analytical and applied studies on management, economy, and dynamics systems based on the system theory [3,[12][13][14][15][16]. In particular, researching the complexity of the closed-loop supply chain system for resource renewability is conducive to elucidating the supplydemand correlation between government and market in the context of urban resilience and sustainable development [17,18]. Regarding the construction of an urban resilience evaluation model and the selection of indicators, the Resilience Alliance proposed the strategic management target and research framework for urban resilience under the SES. Fox-Lent et al. [19] established the Resilience Matrix Framework for four stages of coping with perturbation factors (prepare, absorb, recover, and adapt) from the physical, information, cognitive, and social perspectives. Orencio et al. built the Localized Disaster-Resilience Index that focused on environmental and natural resource management, financial instruments, social protection, technical measures, and planning systems [20]. e modeling method for urban resilience measurement has different focuses owing to specific research objectives and objects, and there is lack of a set of authentic measurement index systems for index selection. Nevertheless, as the system resilience is affected by multiple elements and changes dynamically, the research on specific events in a specified region should be progressive and dynamic, and the guiding function of the research conclusion on spatial coordinated development and adaptive governance should be emphasized.
is pandemic could serve as an efficient empirical test for related studies on urban resilience assessment, which is advantageous to further explore and comprehend the influential elements and interrelation of urban resilience, thus optimizing the supply and demand from the source and actively increasing the resilience of urban SES.

Model Construction
3.1. Research Areas. Wang and Zhao [5] illustrated that the distribution of the number of COVID-19 infections and infection growth rate exhibited notable spatial agglomeration features in the entire country, which is in line with the range of urban agglomerations in the middle reaches of the Yangtze River and the three major coastal city agglomerations, and the aggregation is slowly weakened from the urban agglomerations in the middle reaches of the Yangtze River to other urban agglomerations. us, we selected 28 cities in the middle reaches of the Yangtze River, which were markedly affected by the perturbation factors of the pandemic, as the research areas. Considering the consistency of the research units in terms of the administrative and geographical-economic elements, as well as the data statistics, we excluded three counties directly under the provincial government, namely, Xiantao City, Qianjiang City, and Tianmen City, from the scope of this study.
Based on the Development Plan of the Urban Agglomeration in the Middle Reaches of the Yangtze River approved by the State Council in 2015, the urban agglomeration in the middle reaches of the Yangtze River implies the national megaurban agglomeration, with Wuhan City as the center and Wuhan metropolitan area, Changsha-Zhuzhou-Xiangtan agglomeration, and urban agglomeration around the Poyang Lake as the host components. In 2017, the economic aggregate, population, and area of the urban agglomeration accounted for 9.6%, 11.1%, and 3.3% of the corresponding national indicators, respectively. e research scope ( Figure 1) covered 10 cities of Hubei Province (Wuhan City, Huangshi City, Ezhou City, Huanggang City,

Model Assumption.
In this study, we adopted three empirical elements of urban resilience level (Table 1) to measure the specific numerical value and spatial differentiation of the urban resilience level. In addition, the three empirical elements (Table 1) and the elements in the index system of urban resilience level ( Table 2) were detected and analyzed to identify the dominant driving elements and the interactions between elements.
Martin et al. claimed that the resilience could be better measured on the basis of recovery extent and speed of an area from external disturbance [28]. According to the specific situation of the pandemic and considering the limitation of single index element during characterization and the subjectivity of the weighted index element system, the three empirical elements (the gross domestic product (GDP) year-on-year growth rate (Y1), the year-on-year growth rate of total electricity consumption (Y2), and the in-city travel intensity recovery index (Y3)) were selected to jointly characterize the resilience level of urban social economy within the research areas. In particular, the GDP year-on-year growth rate is a vital macroeconomic observation index, the total electricity consumption can be considered a barometer of social economy operation, and the in-city travel intensity often acts as a crucial measurement index for the socioeconomic activity in the field of urban big data analysis. Considering the differences in characterization features and statistical caliber of different empirical elements, the three indexes were defined on the basis of varying time ranges (detailed in Table 1 and Section 3.3).
e index system for elements influencing the urban resilience level was constructed in combination with the connotation of SES resilience. Based on the measurement model of urban SES resilience level of Cutter et al. [21], together with the attention of the public health to issues like the health system, public organization, culture contact, and vulnerable groups [24], we proposed a hypothesis about the index system for influential elements of the urban resilience level regarding the problems exposed in the pandemic ( Table 2). e system comprised 4 subsystems, 10 indexes, and 20 elements, involving the resistance ability, adaptive ability, and learning ability that should be possessed by the SES. e four subsystems were social resilience, economic resilience, environmental resilience, and engineering resilience, which interacted with each other and were connected internally and externally, thereby jointly and dynamically adjusting the system status.

Data Sources.
Based on the epidemiological analyses on the prevalent features and trends of the pandemic within the research areas [29,30], the pandemic in the research areas can be categorized into four main stages, namely, outbreak period (from January 10 to February 3), peak period (from February 4 to February 20), plateau period (from February 21 to March 8), and involuting period (from March 9 to April 26). On February 21, 2020, the Central Committee held a Political Bureau meeting to deploy the work for pandemic prevention and control, as well as socioeconomic development, so that the pandemic prevention and control could fully enter the second stage. Hence, February 21 was set as a crucial time point in this study. In consideration of the pandemic analysis in the research areas and different characterization features and calculation methods of the three empirical elements, the selection and calculation for different time ranges were accomplished. Table 3 provides a detailed explanation of the major data source.
It is worth noting that the data of the in-city travel intensity (Y3) were obtained from the Baidu Migration Big Data Platform, which maps the trajectory of population flow through mobile phone location information and is characterized by daily continuity and immediacy. In this study, the data of 2019 and 2020 were crawled and cleaned. However, only the data from January 1 to March 15, 2019, were disclosed, and the mean of the numerical values of 2 consecutive weeks (from March 2 to March 15) was used as the comparative data of the in-city travel intensity in 2019 to exclude the impact of the passenger transport during Spring Festival. In 2020, considering that the pandemic entered a plateau period and a new stage of prevention and control since February 21, with February 24 (Monday) being the start point and special values of the May Day holiday being eliminated, the data for a total of 9 weeks from February 24 to April 26 were obtained after crawling and processing of daily data.

Spatial Differentiation Analysis Method: Spatial
Autocorrelation. Spatial autocorrelation reflects the correlativity of a certain geographic phenomenon or attribute value in an area unit with that in adjacent area units [31]. In this study, we used the Global Moran's I and Local Indicators of Spatial Association (LISA) index to measure the degrees of global autocorrelation and local autocorrelation of urban resilience, respectively, as per the following formula: where n denotes the total number of urban units; x i and x j indicate the observed values of variable x in city i and city j, respectively; x denotes the mean of the observed values; and W ij implies the spatial weight matrix value. e range of values is [−1, 1]; the numerical value tending to 1 suggests that the distribution of similar units is more concentrated, that is, spatial agglomeration; the value tending to −1 Discrete Dynamics in Nature and Society 3 suggests the high degree of spatial dispersion; and the value close to 0 indicates that the degree of global autocorrelation is low.
where x i ′ and x j ′ signify standardized values of variable x ′ in city i and city j, respectively, and W ij denotes the spatial weight matrix value. Local autocorrelation can be categorized into five types of agglomeration: HH type, LL type, HL type, LH type, and "not significant" type. Of these, HH type and LL type are positive correlation, suggesting a high degree of spatial aggregation; HL type and LH type are negative correlation, suggesting a high degree of spatial dispersion; and "not significant" implies no obvious agglomeration and dispersion characteristics.

Dominant Element Analysis Method: Geodetector.
Geodetector is a set of statistical methods to detect spatial differentiation and related driving elements [32]. We used differentiation and factor detection in modeling tools and interaction detection to identify the dominant driving elements and the interactions between elements. In the differentiation and factor detection, q value was applied to measure the explanatory power of element X to the spatial differentiation of the urban resilience attribute Y, with a range [0, 1]; the larger the value was, the stronger the explanatory power to the spatial differentiation was. It was calculated as follows: where N indicates the total number of urban units and σ 2 is the variance of Y value. In the interaction detection, there were five types of interactions: nonlinear diminution (q(X1∩X2) < min(q(X1), q(X2))); single-factor nonlinear diminution (min(q(X1), q(X2))<q(X1∩X2) < max(q(X1), q(X2))); two-factor enhancement (q(X1∩X2) > max(q(X1), q(X2))); independent enhancement (q(X1∩X2) � q(X1) + q(X2)); and nonlinear enhancement (q(X1∩X2) > q(X1) + q(X2)).

Urban Spatial Pattern Analysis Method: Landscape
Metrics. According to the land-cover data, we divided landscape types into cultivated land, ecological land, construction land, and unoccupied land. Among them, we specifically measured and calculated the ecological land and construction land using FRAGSTATS software platform and selected six indexes (Table 4) to characterize the urban spatial patterns in the middle reaches of the Yangtze River. e formulas for the abovementioned landscape metrics are as follows: where n i denotes the total area of landscape element i and A denotes the total area of all landscape types, namely, the number of patches per unit area (number of patches/100 hectares).
e AI equals the ratio of the number of self-adjacent public edges to the most probable number of public edges of the landscape element network. In formula (5), g ii signifies the number of similar adjacent patches of the landscape type, which is calculated according to the length of the common boundary between pixels of the same type of patches. e AI is the highest when the length of the common boundary of all pixels of a specific type of patches reaches the maximum value; otherwise, it is the lowest.
where m denotes the total number of patch types in the landscape; e ik means the edge length of patch i adjacent to the patch of landscape element k; and the numerical value equal to 100 illustrates that the edge length of neighboring patches is equal, that is, the probability of vicinity between patches is equal.
where P * ij represents the perimeter of the patch (cell count); a * ij denotes the patch area (cell count); and Z implies the number of cells in the landscape.
where m denotes the total number of patch types in the landscape; P i implies the percentage of area of patch type i; and g ik implies the number of adjacent patch type i and patch type k.
where p i signifies the proportion of patch type i in the landscape.

Study on Spatial Differentiation of the Urban Resilience
Level. According to Table 6, the urban resilience level represented by three elements revealed a significant spatial aggregation effect, which was stronger than the pandemic disaster severity. We conducted LISA clustering to further interpret local autocorrelation ( Figure 3). Precisely, when the urban resilience level was represented by the GDP yearon-year growth rate, a Low-Low Cluster (low-efficiency type) was formed in the northern wing of the urban agglomeration in the middle reaches of the Yangtze River, which comprised Wuhan City, Ezhou City, Huangshi City, Huanggang City, Xiaogan City, Jingmen City, and Xiangyang City, suggesting that the resilience of these areas and surrounding areas was at a low value. A High-High Cluster (high-efficiency type) of urban resilience was formed in the Discrete Dynamics in Nature and Society southeast corner when the urban resilience level was represented by the year-on-year growth rate of total electricity consumption, which comprised Zhuzhou City, Fuzhou City, Yichun City, Ji'an City, and Pingxiang City, and a highresilience cluster area was developed by these areas and the surrounding areas. A Low-Low Cluster was formed in Wuhan City and its surrounding areas. When the urban resilience level was represnted by the in-city travel intensity recovery index, a Low-Low Cluster of urban resilience was formed, comprising Wuhan City, Ezhou City, and Huangshi City. Hence, the resilience level of cities around Wuhan City and in the north of the middle reaches of the Yangtze River was relatively low. For the urban resilience with the total electricity consumption as a single indicator, the cities in the south of the middle reaches of the Yangtze River exhibited a better recovery trend. Regarding the situation in each province, the resilience of Xinyu City was relatively low in Jiangxi Province.
We concretely analyzed the in-city travel intensity recovery index by time periods (Table 7, Figure 4); its spatial relationship underwent the process of "significant aggregation," "random, "short-term discretization," and "long-term random." It had been roughly in the dynamic change of lowvalue cluster in the surrounding areas of Wuhan City and high-value cluster in the southeast corner since the beginning of the study, and there existed a short-term situation of a high-value cluster in some cities of Hubei Province. Considering the particularity of time nodes, we believe that it had a specific relationship with the local policies on pandemic prevention and control. As of April 26, the travel intensity in all cities (except Wuhan) was basically the same : GDP year-on-year growth rate GDP F(j) g dp � (G j1 − G j2 )/G j2 , where G j1 denotes the constant-price GDP of city j in the first three quarters of 2020 and G j2 denotes the constant-price GDP of city j in the first three quarters of 2019.
Y2: year-on-year growth rate of total electricity consumption Tec where E j1 denotes the total electricity consumption with constant energy consumption of city j in July 2020 and E j2 denotes the total electricity consumption with constant energy consumption of city j in July 2019.
Y3: in-city travel intensity recovery index Td where T j means the in-city travel intensity of city j, namely, the indexation result of the ratio of the number of travelers to the resident population in the city; T j1 denotes the daily data of the in-city travel intensity of city j from February 24 to April 26, 2020 (a total of 9 weeks); and T j2 denotes the mean of the daily data of the in-city travel intensity from March 2 to 15, 2019.
Y0: pandemic disaster severity Ds where pop t denotes the total number of confirmed cases in city j as of February 21, 2020. Discrete Dynamics in Nature and Society as that in 2019, and the spatial agglomeration effect of this attribute gradually weakened during the study period. We used the three empirical elements to measure the urban resilience level. e results revealed that, in the event of the COVID-19 pandemic, urban resilience in the middle reaches of the Yangtze River displayed significant characteristics of spatial agglomeration, with a roughly spatial trend of decreasing from south to north. As of the end of April, the urban travel intensity in all cities (except Wuhan) was the same as in 2019. e cities along the southeast wing had a good recovery trend in the year-on-year growth rate of total electricity consumption and in-city travel intensity recovery index. e overall economic recovery and development of the cities in Hunan Province were relatively

Landscape metrics Ecological significance
Patch density (PD) It reflects the degree of landscape fragmentation; a higher numerical value indicates a smaller patch and a higher degree of fragmentation.
Aggregation index (AI) It reflects the aggregation degree of the same type of patches; a higher numerical value indicates a greater aggregation degree of similar landscape types.

Interspersion juxtaposition index (IJI)
It reflects the overall interspersion juxtaposition of the same patch; a low numerical value implies that this patch is only adjacent to a few other types of patches.

Patch cohesion index (COHESION)
It reflects the natural connectivity of patches; the higher the numerical value, the better the connectivity of the landscape type.

Contagion index (CONTAG)
It reflects the agglomeration degree or extension trend of different landscape types, including the spatial information; a high numerical value implies that good connectivity is formed for a specific dominant landscape type. Shannon's diversity index (SHDI) It reflects the heterogeneity and complexity of landscape types; a higher numerical value indicates more abundant land use, which tends to be evenly distributed.
Note. Refer to the help file of FRAGSTATS v4.2.1 for details.  Table 8 reports the top 10 indicators by the q value. Considering the ranking of the three elements, the proportion of tertiary industry (E3), GDP growth rate (E2), registered unemployment rate in urban areas (S2), AQI (N5), NDVI (N4), green coverage rate in built-up areas (N2), population density (N3), and proportion of Internet users (P3) had a high explanatory power for the spatial differentiation of urban resilience. is indicates that the innovative economic development, stable employment environment, reasonable population size, safe and ecological landscape pattern, and efficient information transmission in the research areas under the disturbance element of the pandemic markedly affect the short-term urban social and economic recovery.

Interaction Detection.
According to the interaction detection matrix (Table 9), the driving force of all detection elements after interaction on the urban resilience in the research areas was more significant than the explanatory   power of a single element. e interaction type of elements included double-element enhancement and nonlinear enhancement, with no independent elements present. Some elements with relatively weak explanatory power in singleelement detection demonstrated strong explanatory power in interaction detection, suggesting that the effect of elements like the number of beds per 10,000 people on the urban resilience in this event could be fully shown when interacting with other elements. When the urban resilience was characterized by the GDP year-on-year growth rate, the explanatory power of the proportion of the tertiary industry, registered unemployment rate in urban areas, and air quality was markedly enhanced when interacting with other elements. Among them, the air quality ∩ proportion of tertiary industry had the strongest explanatory power (0.975). From the standpoint of the year-on-year growth rate of total electricity consumption, the explanatory power of GDP growth rate significantly increased when interacting with other elements. e explanatory power of the proportion of urban construction land, which had a relatively mediocre performance in single-element detection, was as high as 0.793 after  Discrete Dynamics in Nature and Society  10 Discrete Dynamics in Nature and Society interacting with the proportion of tertiary industry. From the viewpoint of the in-city travel intensity recovery index, the explanatory power of the number of schools per 10,000 people ∩ the registered unemployment rate in urban areas was the highest (0.910), and the explanatory power of the proportion of tertiary industry, proportion of Internet users, NDVI, population density, and green coverage rate in builtup areas was high when interacting with other elements.
e above-described research revealed that good industrial pattern and employment environment, excellent natural ecology and pollution control, harmonious city scale and form, and complete information infrastructure construction were of utmost significance for the short-term urban socioeconomic recovery. e explanatory power was increased after the interaction of various influential elements, including social, economic, ecological, and engineering subsystems. Among them, the interaction of various elements of the ecological and economic subsystems was more significant, acting together on the dynamic cycle of the urban social ecosystem. Hence, it is essential to enhance the endogenous development momentum of cities by building economic markets and improving space quality.

Recommendations on Adaptive Governance Based on
Influential Elements

Ecological and Safe Landscape Pattern Control.
Close attention was paid to the significant effects of the economic and ecological subsystems. We further measured the specific landscape metrics to analyze the urban spatial pattern (Figures 5-7). Comparatively, the cities of Jiangxi Province as a whole displayed better ecological connectivity and homogeneity. Hunan Province had fragmented ecological land. Wuhan and its surrounding areas and the cities of Hunan Province revealed a relatively large population and land scale, and relatively concentrated construction land. e SHDI data formed a "double T" high-value distribution pattern with provincial capital cities as nodes in the spatial pattern, with rich landscape patterns and high fragmentation and uncertainty. Coordinating urban construction by integrating existing ecological advantage resources, sorting out urban landscape and construction land structure, enhancing the malleability and connectivity of ecological land, and articulating urban policies combining scale and form are crucial to partially eliminate the pressure on urban scale and density and enhance urban resilience.  show the economic subsystem data and the spatial autocorrelation of the urban innovation and entrepreneurship index and indicate that the Wuhan-Changsha main axis developed vigorously under the connection of the Beijing-Guangzhou line in the early years and exhibited the trend of conversion from development along the axis to contiguous distribution. Wuhan City had strong population and economic agglomeration effects; however, its surrounding cities displayed slow economic development and hindered transmission of central functions. Notably, Jiangxi Province exhibited a strong economic situation in recent years, and an innovation economy radiation circle with Nanchang City as the core had been developed. e advantages of each cluster and the economic belt along the river should be connected by constructing transportation and other infrastructure to break down faults and barriers and promote the outward flow of innovative elements in central cities. Furthermore, new ways of ecological priority and green development should be actively investigated to build an area regional economic development community with complementary advantages.

Multiple Efforts, Adaptation to Local Conditions, and
Interval Coordination. Adaptive governance measures should be actively investigated to respond to uncertain disturbances similar to pandemic events to enhance the risk adaptability and sustainable development capability of cities. Besides, measures like multiple efforts, adaptation to local conditions, and interval coordination should be considered to upgrade the system cycle mechanism. Based on the local spatial autocorrelation of eight dominant elements in Figure 10 and the spatial distribution of 20 indicators, it can be inferred that the three provincial capital cities had unique location advantages and policy advantages; however, more space was available to develop the radiation and driving effect of economy and informatization. Comparatively, Jiangxi Province had significant advantages in the environmental resilience subsystem, while Wuhan City and Ezhou City had great pressure on the urban human settlement environment and low vegetation index and green coverage in built-up areas. e northern cities in the middle reaches of the Yangtze River were relatively weak in withstanding and adaptation to risks at the economic level. It was a must to focus on the     economic growth momentum and the degree of advanced industrial structure of cities. Overall, there remains scope for improvement in economic innovation and green development in Hubei Province; thus, the absorption capacity of Wuhan's surrounding cities to receive spillovers from central cities must be strengthened. Jiangxi Province had significant advantages in resource utilization and sustainable environmental development; therefore, attention must be paid to the efficiency of economic output and the cultivation of high-quality talents. Hunan Province had good economic, scientific, and educational development trends; however, issues like balanced economic development and environmental governance must be settled urgently to avoid extensive and inefficient use of resources and excessive concentration of chemical industry. It is imperative to further play the economic radiation role of central cities and augment the absorption capacity of smalland medium-sized cities to create the multicenter development pattern with Wuhan, Changsha, and Nanchang as centers and Ezhou, Zhuzhou, and Jiujiang as regional subcenters. Besides, it is crucial to further explore the advantageous and characteristic industries of subgroups of

Conclusions
In the event of the pandemic, urban resilience in the middle reaches of the Yangtze River showed significant characteristics of spatial agglomeration, with a roughly spatial trend of decreasing from south to north. e cities along the southeast wing had a good recovery trend in the year-onyear growth rate of total electricity consumption and in-city travel intensity recovery index. e overall economic recovery and development of the cities of Hunan Province were relatively balanced. e spatial differentiation of urban resilience in the middle reaches of the Yangtze River is a comprehensive outcome of the interaction of various elements. e independent and interactive forces of ecological and economic subsystems are significant. e proportion of the tertiary industry, GDP growth rate, registered unemployment rate in urban areas, AQI, NDVI, green coverage rate in built-up areas, population density, and proportion of Internet users play key roles in short-term social and economic recovery of cities.
Currently, attaining the transformation and upgrading of traditional industries by increasing the ability of cities around the provincial capital to receive spillover resource, optimizing the spatial pattern of ecological land and construction land in Hunan and Hubei Provinces, and promoting the overall sustainable supply chain innovation in the region constitute the major breakthrough promoting the resilience and sustainable development of the region. In this study, the level of urban resilience in the research areas was jointly characterized by three elements (GDP year-onyear growth rate, year-on-year growth rate of total electricity consumption, and in-city travel intensity recovery index), which have different characteristics and limitations, and a more scientific description of urban resilience warrants further investigation. e selection of elements has a clearer event orientation and geographical orientation, and its external validity also warrants verification. e socioeconomic recovery status of each city under the pandemic closely correlated with administrative power, which was not extensively discussed in this study owing to the limited length of the research. e urban SES is a dynamic structure under the interaction of multiple subjects and elements, which is challenging to analyze. Bias is inevitable in the process of scientific argumentation in this study. e ability to comprehend the urban system should be improved in future research and practice, and dynamic governance strategies should be adopted, which is of utmost significance for the construction of a resilient and sustainable urban development spatial coordination mechanism.
Data Availability e data source is specifically described in Section 3.3, and the data used to support the findings of this study are included within the article.

Additional Points
(1) Only some counties (districts) in Fuzhou City and Ji'an City are included in the urban agglomeration in the middle reaches of the Yangtze River, but the entire cities of Ji'an and Fuzhou are included in this study to facilitate data acquisition and maintain the integrity of the research areas. (2) e Statistical Yearbook 2018 used in this study publishes statistical data of 2017.

Conflicts of Interest
e authors declare that they have no conflicts of interest.