Coupling Coordination Relationship between Urban Sprawl and Urbanization Quality in the West Taiwan Strait Urban Agglomeration, China: Observation and Analysis from DMSP / OLS Nighttime Light Imagery and Panel Data

: Urban sprawl is the most prominent characteristic of urbanization, and increasingly a ﬀ ects local and regional sustainable development. The observation and analysis of urban sprawl dynamics and their relationship with urbanization quality are essential for framing integrative urban planning. In this study, the urban areas of the West Taiwan Strait Urban Agglomeration (WTSUA) were extracted using nighttime light imagery from 1992 to 2013. The spatio-temporal characteristics and pattern of urban sprawl were quantitatively analyzed by combining an urban expansion rate index and a standard deviation ellipse model. The urbanization quality was assessed using an entropy weight model, and its relationship with urban sprawl was calculated by a coupling coordination degree model. The results showed that the urban area in the WTSUA experienced a signiﬁcant increase, i.e., 18,806.73 km 2 , during the period 1992–2013. The central cities grew by 11.08% and noncentral cities by 27.43%, with a general uneven city rank-size distribution. The urban sprawl showed a circular expansion pattern, accompanied by a gradual centroid migration of urban areas from the southeast coast to the central-western regions. The coupling coordination level between urban expansion and urbanization quality increased from serious incoordination in 1992 to basic coordination in 2013. Dual driving forces involving state-led policies and market-oriented land reform had a positive inﬂuence on the harmonious development of urban sprawl and urbanization quality of the WTSUA. This research o ﬀ ers an e ﬀ ective approach to monitor changes in urban sprawl and explore the coupling coordination relationship between urban sprawl and urbanization quality. The study provides important scientiﬁc references for the formulation of future policies and planning for sustainable development in urban agglomerations.


Introduction
Urban sprawl is a prevalent phenomenon worldwide because of accelerating population growth and global urbanization [1][2][3].Its manifestation through the evolution of regional urban systems has far-reaching impacts on urban economic, social, and eco-environmental development.The harmonization of urban sprawl with eco-environment and socio-economic interests is crucial to the sustainable development of urban systems [4].In recent decades, the urbanization process has been gradually transformed from the traditional urbanization mode to the new-type urbanization path, with a focus on the "quality" of urbanization, which has a greater emphasis of eco-environment conservation and human well-being improvement [5].Furthermore, when urban centers expand and become clustered around large cities [6], urban agglomerations form and dominate regional urban economic and social development [7,8].Then, the relationship between urban sprawl and urbanization quality becomes more complicated [9].Under such circumstances, studies on the coordination relationship between urban sprawl and urbanization quality are a top priority because they help to develop knowledge-based urban planning strategies and actualize effective urban sprawl paths, which are vital for social stability, economic prosperity, and ecological safety [10,11].
Urban sprawl in the world has followed different paths.In most western countries, market-oriented land use contributes to urban sprawl, showing a leapfrog (i.e., noncontiguous, or scattered), low density, and single-use urban development pattern [12,13].In contrast, urban sprawl in China is characterized by high-density housing projects, multifunctional community construction, and rapid expansion, driven by the coexistence of state-led, top-down planning policy and market-regulated, bottom-up demand [10,14].Many empirical studies on urban sprawl have been carried out to determine its trends, intensity, driving forces, evolution, and future scenarios [5,15,16].However, in contrast to the rich literature on developed regions and countries, urban sprawl in emerging urban agglomerations has not been fully investigated.As such, it is essential for urban agglomerations to analyze the range and intensity of urban sprawl systematically, as well as to identify the main influencing factors to facilitate effective strategic decision making.
An enormous number of studies on urbanization quality have been conducted by qualitative and quantitative methods [17][18][19].For example, Zhou et al. [9] developed a conceptual comprehensive index system of urbanization quality assessment from the perspective of the structure and function of an urban system.Lin et al. [20] proposed a comprehensive assessment system of urbanization quality including demographic urbanization, spatial urbanization, economic urbanization, and social urbanization to analyze the effects of urbanization on the eco-environment.Wang et al. [18] and Martínez-Zarzoso [21] evaluated urbanization quality and examined the relationship between it and CO 2 emissions in China and several developing countries, respectively.However, existing research has seldom focused on the interaction mechanisms between urbanization quality and urban sprawl.
Remote sensing (RS) and geographic information systems (GIS) are common tools for analyzing urban sprawl at different spatial scales [22,23].Remote sensing allows the change and distribution of urban areas to be delineated more effectively [24], while GIS enables valuable spatial information to be analyzed [25].For example, QuickBird and Worldview-2 images were used to detect the changing trend and spatial characteristics of urban sprawl in one Indian city [24].Landsat TM/ETM+ images were used to detect and analyze urban change in a rapidly urbanizing region of China and in Pune metropolis, India [23,26].The spatio-temporal evolution of urban agglomerations of bay areas in the US, China and Japan were detected and compared using Landsat TM/OLI [27].Because of the low classification efficiency of traditional optical image data, time-discontinuous images are often used to extract urban areas, leading to the omission of detailed changes over long periods [28].In contrast, nighttime light images well represent the spatial extent of urban areas and have advantages in studying urban dynamics on large spatial scales or with high temporal frequencies [29][30][31][32].For example, Castells-Quintana [33] and Li et al. [34] used nighttime light images to analyze the urban pattern and economic growth in sub-Saharan Africa and in the countries along "The Belt and Road", respectively.Huang et al. [35] detected changes in the city size dynamics in China from 1992 to 2008 and discovered that urban areas extracted from nighttime light imagery were more accurate than census data.Stokes and Seto divided cities into different categories using nighttime light data to explore whether infrastructure was synchronized with land and demographic transitions [36].Generally, nighttime light images are valuable for characterizing urban expansion at regional, national, and global scales [28,37], and it is worth exploring the potential of this data type in the study of urban sprawl dynamics.However, few studies have systematically explored urbanization quality based on nighttime light radiance variation, which has seriously impeded the application of nighttime light data in urban development [38].
The West Taiwan Strait Urban Agglomeration (WTSUA) is situated on the southeast coast, across the Taiwan Strait from Taiwan Province, China, and plays a vital role in strengthening the communication and cooperation between Mainland China and Taiwan Province [39].Given the advantages of nighttime light imagery in reflecting urban dynamics, as well as the importance of revealing the interaction between urban sprawl and urbanization quality, the objectives of this research were to: (1) extract and map urban areas of the WTSUA; (2) characterize the dynamics of urban sprawl in terms of intensity, rank-size, and expansion pattern; (3) assess the urbanization quality; and (4) determine the coupling relationship between urban sprawl and urbanization quality over time.This analysis will provide valuable information for urban policymakers and planners to formulate urban development strategies.To fulfill these objectives, time-series nighttime light images were used to extract urban areas.The urban expansion rate index (UERI), standard deviation ellipse (SDE), and rank-size rule (RSR) model were used to calculate the urban sprawl intensity, sprawl pattern, and city-size distribution dynamics, respectively.The entropy weight (EW) model was used to evaluate the urbanization quality of the WTSUA.A coupling coordination degree (CCD) model was applied to quantify the relationship between urban sprawl and the urbanization quality of urban agglomeration.

Study Area
The WTSUA is located to the west of Taiwan Strait, China.The longitude ranges from 22.82 • N to 29.57• N and the latitude ranges from 113.86 • E to 121.23 • E, with a total area of 2.70 × 10 5 km 2 and a total population of about 9.42 × 10 7 in 2018.It consists of four provinces, Fujian, Zhejiang, Jiangxi, and Guangdong, including 20 prefecture-level administrative cities (Table 1 and Figure 1).The five central cities are in the eastern coastal region, and the fifteen noncentral cities are mostly situated in the central-western region.The average annual precipitation is 1249-2985 mm, and the annual mean temperature is 15-23 • C. The prevalent climate of this region is a subtropical monsoon climate, characterized by a hot and rainy summer and a warm and wet winter.More than 75% of the region is covered by natural forests, which provides suitable habitat for numerous wildlife species, and it possesses more ecological resources than other urban agglomerations in China.
In December 2009, the "Coordinated Development Planning for West Taiwan Strait Urban Agglomeration", a major national program for promoting the development of the WTSUA, was officially issued, marking the start of a more detailed and concrete development strategy and target for the WTSUA.By virtue of its convenient transport networks, geographical location, and support from the state, the WTSUA has become one of China's most important economic zones.Furthermore, the region is becoming a platform for exchange and cooperation with global influence.The GDP of the WTSUA reached 5.87 trillion Yuan in 2018, accounting for 6.52% of China's GDP.

Nighitime Light Imagery and Pre-Processing
To analyze the urban sprawl dynamics of the WTSUA, a time series dataset from the DMSP/OLS (Defense Meteorological Satellite Program's Operational Linescan System) was used with a spatial resolution is 1 km.The DMSP/OLS dataset provides yearly nighttime light images for the period of 1992-2013, and records nighttime light through a DN (digital number) value for each pixel, i.e., from 0 to 63.A value of '0' represents no light in the area, and a value of '63' indicates maximum brightness.
Compared with other publicly available nighttime light imagery, i.e., NPP-VIIRS (National Polar-orbiting Partnership Visible Infrared Imaging Radiometer Suite), Luojia 1-01, and Jilin 1-03b, the DMSP/OLS data eliminate the effects of fires, lighting, moonlight, clouds and auroras, reducing the error caused by background noise.Furthermore, each DMSP/OLS image is acquired from 8:30 to 9:30 pm locally, while the NPP-VIIRS satellite passes over at about 01:30 am.Nighttime lights captured earlier, especially at around 9:00 pm, are likely to include light sources from residents, roads, commercial areas, and industrial sites.When nighttime lights are measured around midnight, some of these light sources may no longer contribute, depending on lifestyles in different regions [40].On the other hand, although Luojia 1-01 (spatial resolution is 130 m) [41] and Jilin 1-03b (spatial resolution is 0.89 m) [42] nighttime light data are superior in spatial resolution, radiometric detection range, and onboard calibration, their available time spans are 2018 to the present and 2017 to the present, respectively; as such, it would be difficult to undertake long-term research using them [38,43,44].Therefore, DMSP/OLS data are more representative and suitable for long-term analyses of urban areas than other nighttime light data.
The DMSP/OLS time series data were acquired from six different DMSP satellites, i.e., F10 (1992-1994), F12 (1994)(1995)(1996)(1997)(1998)(1999), F14 (1997-2003), F15 (2000-2007), F16 (2004)(2005)(2006)(2007)(2008)(2009), and F18 (2010-2013).Datasets acquired by different satellites in the same year or by the same satellite over different years are not comparable until they are radiometrically calibrated.As referenced by Hu et al. [40], the DMSP/OLS data were calibrated from 1992 to 2013.After systematic calibration, the urban areas of the WTSUA from 1992 to 2013 were extracted using the thresholding technique [28].Due to considerable regional variation in socio-economic development status and physical environment, it was difficult to select a single threshold to distinguish the urban areas.Therefore, the WTSUA was divided into two regions, i.e., coastal region and inland region.The former includes WZ, ND, FZ, PT, QZ, XM, ZZ, CZ, JY, and ST.The latter includes LS, NP, SM, LY, MZ, QZa, SR, YT, FZa, and GZ.The threshold was determined automatically by the method of maximum interclass variance, also known as the OTSU algorithm [45].First, 10,000 random points in 1995, 2000, 2005, and 2010 were generated in the urban areas according to the auxiliary data which were the land cover data for 1995, 2000, 2005, and 2010, provided by the National Earth System Science Data Center, China.Then, the DN value of each point was extracted in the corresponding years.Next, the optimal thresholds for 1995, 2000, 2005, and 2010 were calculated by the OTSU method for each region.Finally, the determined thresholds were applied to different stages, i.e., the thresholds for 1995, 2000, 2005, and 2010 were applied to the periods 1992-1997, 1998-2002, 2003-2008, and 2009-2013, respectively.After the urban areas of each year had been extracted, the accuracy verification dataset of the extracted urban areas was produced.Specifically, 100 urban patches were chosen randomly from the 2010 land cover data.Then, the urban patches in 1995, 2000, and 2005 were also extracted based on land cover data for 1995, 2000, and 2005, which corresponded spatially to the 100 patches in 2010.Finally, the urban areas of the remaining years located in the 100 random patches were mapped by visual interpretation using Landsat TM/OLI images.All of these constituted the accuracy validation dataset, which was used to evaluate the accuracy of extracted urban areas by the DMSP/OLS data.As referenced by Ma et al. [46], a polygon-based validation method was used to validate the extraction accuracy of urban areas by DMSP/OLS data, showing that the overall accuracy for each year was higher than 0.89, and the average overall accuracy is 0.91 ± 0.02.

Urban Expansion Rate Index
To describe the speed of urban expansion objectively in spatial terms, an UERI was calculated, which may be expressed as follows [47]: where U a is the urban area at the beginning year of the study periods, U b is the end year of the study periods, and T is the number of years.Time intervals considered in this research were 1992-1999, 1999-2006, and 2006-2013.

Standard Deviation Ellipse
The morphological characteristics of the distribution of urban areas, e.g., range, direction, and centroid, make it possible to provide a more intuitive description of the trends of urban sprawl to determine the pattern and process of urban change.To achieve this, the SDE model was chosen to calculate the centrality, orientation, and shape of the urban distribution [48].The area of the ellipse represents the region in which urban areas are spatially dominant.The azimuth of the ellipse reflects the principal tendency of urban distribution.The centroid is the center of all urban areas in two-dimensional space.Changes in the urban centroid over time can reflect the overall evolutionary route of urban sprawl.The minor and major axes of the ellipse indicate the directions and scope of the urban distribution and the dispersion degree among different urban regions [49].A greater difference between the minor and major axes indicates a more obvious polarization phenomenon.The elements of SDE are expressed as follows [50]: where (x i , y i ) are the geographic coordinates of the ith urban area, w i is a weight, (x i , y i ) are the weighted means of the ith urban area.(X w , Y w ) is the centroid of the ellipse.SDE x is the length of the long axis of the ellipse, and SDE y is the short axis of the ellipse.θ represents the ellipse directional orientation, indicating the clockwise rotation degree from north to the ellipse long axis.

Rank-Size Rule Model
A comparison of the evenness of change in city size over time is vital to measure the intercity concentration in the urban agglomeration system [35].To quantify the distribution dynamics of city size for the WTSUA, the RSR model was applied to quantify the vertical dimension of the urban agglomeration system, which is calculated as follows [51]: where P 1 is the largest city's urban area in WTSUA, P i is the urban area of the city which is ranked ith by urban area.R i is the number of cities with city size bigger than P i .q is the Pareto coefficient that denotes the distribution evenness of city size.When q = 1, the city size is directly proportional to the largest city size and the reciprocal of its rank.When q < 1, the distribution of city size is less even, indicating that there are only a small number of large-sized cities.In contrast, when q > 1, the distribution of city size is more even.The temporal change of q can reflect the relationship between the centralization force and dispersion force of city-size distribution.If the value of q presents an increasing tendency, the intercity centralization force is greater than its dispersion force, and vice versa.

Establishment of an Assessment System for Urbanization Quality
Understanding urbanization quality is a complex and multidimensional task.To assess urbanization quality, a well-diversified evaluation system was developed [9,20].Given the data availability and comparability over the study period, as well as their relevance to the assessment objective, panel data of twenty-eight basic variables were selected to assess urbanization quality from 1992 to 2013.These indicators are described in Table 2.The original values for all evaluation variables were derived from the China City Statistical Yearbook and the China Statistical Yearbook.Per capita green areas (m 2 ) X 5 Treatment rate of domestic sewage (%) X 6 Landfill rate for urban waste (%) Social Well-being Level X 7 Proportion of built-up area to total area (%) X 8 Per capita urban road area in municipal district (m 2 ) X 9 Per capita public library holdings (volume) X 10 Proportion of science and technology expenditure to GDP (%) X 11 Proportion of educational expenditure to GDP (%) X 12 Number of teachers per 1000 students (unit) X 13 Per capita housing (m 2 /person) X 14 Number of medical beds per 1000 persons (unit) X 15 Registered urban unemployment rate (%) X 16 Number of buses per 1000 persons (unit)

Econo-demographic Development X 17
Proportion of population in municipal district to total population (%) X 18 Population density in municipal district (person/km 2 ) Proportion of secondary and tertiary industry employment to total employment (%) X 20 Per capita investment in fixed assets (Yuan/person) X 21 Per capita GDP (Yuan/person) X 22 Per capita GDP in the secondary and tertiary industry (Yuan/person) X 23 GDP (100 million yuan) X 24 Per capita disposable income of urban residents (Yuan) X 25 Total retail sales of consumer goods (100 million yuan) X 26 Proportion of secondary and tertiary industry GDP to total GDP (%) X 27 Per capita local government fiscal revenue (Yuan) X 28 Financial institution deposits at year-end (100 million Yuan)

Entropy Weight Model
The EW model is used to describe the relative importance of each indicator for the overall assessment objective.Several evaluation methods can be used for calculating weights, e.g., analytic hierarchy process [52], Delphi [53], principal component analysis [54], and artificial neural networks [55].The EW model can determine an index's weight based on the value of information entropy for the index itself, which avoids the influence of subjective bias from different decision-makers [56,57].Moreover, the approach is not only time-saving, but could also reduce the deviations caused by the discrepancies between different evaluated objects [58].Therefore, it can reflect the relationship among indices objectively and efficiently.In this study, urbanization quality was evaluated by the EW model.

Coupling Coordination Degree Model
A CCD model was used to measure the coordination relationship between urban sprawl and urbanization quality.The CCDM can be expressed as follows [59]: where g(y) and f (x) are the function of urban sprawl and urbanization quality, respectively.C is the coupling degree of urban sprawl and urbanization quality.T reflects the overall coordination level of the urban sprawl and urbanization quality.D is the CCD of urban sprawl and urbanization quality; the higher its value, the higher the coupling coordination level.Both α and β are the weights of urban sprawl and urbanization quality to be determined.Considering that urban sprawl plays a vital role in urbanization quality in the WTSUA, the weights are supposed to be equal, i.e., α = β = 0.5 [60].
As referenced by Xing et al. [59] as well as Wen and Wen [61], the coupling coordination relationship between urban sprawl and urbanization quality can be described by dividing the CCD into different categories.In this study, the coupling coordination level was compartmentalized into six types, as shown in Table 3.

Rate, Intensity, and Pattern of Urban Sprawl
After automatically extracting urban areas in each year based on the OTSU threshold method, the speed, intensity, and pattern of urban expansion were analyzed.Table 4 shows the expansion area and spatial expansion rate of urban areas for different city types.The results indicated that the urban areas of each city experienced a general increasing trend from 1992 to 2013, with a total growth rate of 895.56 km 2 /y for the WTSUA.During 1992-1999 and 1999-2006, the urban areas of

Rate, Intensity, and Pattern of Urban Sprawl
After automatically extracting urban areas in each year based on the OTSU threshold method, the speed, intensity, and pattern of urban expansion were analyzed.Table 4 shows the expansion area and spatial expansion rate of urban areas for different city types.The results indicated that the urban areas of each city experienced a general increasing trend from 1992 to 2013, with a total growth rate of 895.56 km 2 /y for the WTSUA.During 1992-1999 and 1999-2006, the urban areas of central cities expanded more than those of noncentral cities.In contrast, the urban areas of noncentral cities exhibited a dramatic increase from 2006 to 2013, which was more than 1.5 times the total expansion area of central cities in this stage.In terms of the sprawl intensity, the spatial expansion rate of noncentral cities in urban areas was always higher than that of central cities from 1992 to 2013, indicating a higher value of UERI.Except for Xiamen city, the urban expansion areas of central cities all exceeded those of noncentral cities over the 22 years.Zhangzhou city showed the largest expansion in noncentral cities, increasing by 1517.16 km 2 from 1992 to 2013.
Figure 3 illustrates the spatio-temporal dynamics of SDE elements for urban sprawl from 1992 to 2013.Despite slight fluctuations, the SDE centroid showed a process of shifting northwest throughout the study period (Figure 3a), implying a stronger expansion power for cities located in the central-western regions of the WTSUA.The spatial distribution area of SDE for urban areas experienced a persistent increasing trend (Figure 3b), suggesting an increasing expansion of the main urban distribution areas.It is noteworthy that the SDE azimuth changed from 35.53 • in 1992 to 44.67 • in 2013, presenting a trend of gradually rotating southeast (Figure 3c).These results suggest that the urban sprawl force driven by the central-western inland cities (belonging to noncentral cities) exceeded that from the eastern coastal cities (mostly belonging to central cities) over time.
The spatio-temporal sprawl pattern of urban areas in the WTSUA is displayed in Figure 4.A cycloid expansion model of urban areas was found for each city from 1992 to 2013, showing that urban areas expanded gradually in circular or semicircular patterns, with their initial area as the center.The urban sprawl of the total urban agglomeration demonstrated a multi-nuclei pattern as a whole, including three main cores, i.e., the Wenzhou, Fuzhou-Quanzhou-Xiamen, and Shantou regions.From 1992 to 2000, despite some small fluctuations, the Pareto coefficient increased, indicating that more nonuniformities had occurred during this stage.From 2001 to 2009, the Pareto coefficient experienced a general decrease.It can be inferred that during this period, WTSUA's cities experienced a "balancing trend", with large-sized cities growing more slowly and decreasing in rank, while small-sized cities grew more rapidly and increased in rank.From 2010 to 2013, the Pareto coefficient was basically stable with a lower value, suggesting a more balanced distribution of urban rank-size than the previous two stages.It can be inferred that during this period, WTSUA's cities experienced a "balancing trend", with large-sized cities growing more slowly and decreasing in rank, while small-sized cities grew more rapidly and increased in rank.From 2010 to 2013, the Pareto coefficient was basically stable with a lower value, suggesting a more balanced distribution of urban rank-size than the previous two stages.

Urbanization Quality of Urban Agglomeration
The weight of each indicator for the urbanization quality evaluation system was calculated by the EW model shown in Figure 6.It indicated that econo-demographic development accounted for the largest weight (42.34%), followed by social well-being level (36.18%).Eco-environment status was the smallest weight (20.58%) in the first-grade indices.In terms of the basic indices, their weight was roughly similar.
Figure 7 illustrates the urbanization quality dynamics of the WTSUA.A progressive increase trend was shown for the urbanization quality of each city from 1992-2013.Among different subsystems of urbanization quality, econo-demographic development changed most dramatically, with an increase of 1402.86% during the study period.The social well-being level and eco-environment status increased by 78.80% and 98.69% from 1992 to 2013, respectively.

Urbanization Quality of Urban Agglomeration
The weight of each indicator for the urbanization quality evaluation system was calculated by the EW model shown in Figure 6.It indicated that econo-demographic development accounted for the largest weight (42.34%), followed by social well-being level (36.18%).Eco-environment status was the smallest weight (20.58%) in the first-grade indices.In terms of the basic indices, their weight was roughly similar.
Figure 7 illustrates the urbanization quality dynamics of the WTSUA.A progressive increase trend was shown for the urbanization quality of each city from 1992-2013.Among different subsystems of urbanization quality, econo-demographic development changed most dramatically, with an increase of 1402.86% during the study period.The social well-being level and eco-environment status increased by 78.80% and 98.69% from 1992 to 2013, respectively.

Coupling Coordination Relationship between Urban Sprawl and Urbanization Quality
Figure 8 visualizes the spatio-temporal dynamics of the coupling coordination relationship between urbanization quality and urban sprawl for the WTSUA from 1992 to 2013.During the study period, the level of coupling coordination increased gradually for each city.Moreover, the CCD value of the central-western inland areas was obviously lower than that of the eastern coastal regions, showing a significant spatial variation (Figure 8a).As for the total WTSUA, the value of the CCD presented a gradual increase over time, suggesting the enhancement of improving coordination between urban sprawl and urbanization quality from 1992 to 2013 (Figure 8b).Specifically, the coupling coordination level of urban sprawl and urbanization quality experienced four stages:

Coupling Coordination Relationship between Urban Sprawl and Urbanization Quality
Figure 8 visualizes the spatio-temporal dynamics of the coupling coordination relationship between urbanization quality and urban sprawl for the WTSUA from 1992 to 2013.During the study period, the level of coupling coordination increased gradually for each city.Moreover, the CCD value of the central-western inland areas was obviously lower than that of the eastern coastal regions, showing a significant spatial variation (Figure 8a).As for the total WTSUA, the value of the CCD presented a gradual increase over time, suggesting the enhancement of improving coordination between urban sprawl and urbanization quality from 1992 to 2013 (Figure 8b).Specifically, the coupling coordination level of urban sprawl and urbanization quality experienced

Urban Development Associated with State-led Policies and Market-oriented Land Reform
Urban areas, considered as the core centers of social and economic development, are crucial for ensuring national stabilization and people's livelihoods.Since the implementation of 'Reform and Opening-up', China has experienced rapid urbanization [62,63].The ratio of urbanization was less than 18% in 1978 but up to 36% in 2000; it will be more than 80% by 2050 [8].In this context, the urban system has driven unprecedented development, accompanied by urban spatial expansion and improvements in the urban structure.In the WTSUA, our results suggest that the distribution and expansion speed of urban areas have increased, with a rise of 1602.86% from 1992 to 2013 (Figures 2 and 3, and Table 4).Previous studies in China concerning the urban sprawl dynamics in Beijing, Shanghai, Guangzhou, and Wuhan show that these cities present a similar trend of urban expansion [64][65][66][67].In the process of rapid urban expansion, a series of issues associated with the eco-environment and society have arisen [68][69][70].Facing the challenges of urban sprawl, several plans and policies have been promulgated to promote coordinated urban development.Figure 9 illustrates the characteristics, emphasis, and pivotal policies and events in the different urban sprawl stages of China.It shows that the guiding principles of national urban development were issued and implemented in the early 1980s, which aimed to regulate orderly urban development by actively developing small-sized cities, slowly developing medium-sized cities, and controlling the scale of large-sized cities [62].Driven by these policies, small-and medium-sized cities have developed rapidly.This probably explains why the expansion speed of the noncentral cities of the WTSUA was higher than that of the central cities, and why the centroid of urban areas shifted towards noncentral cities from 1992 to 2013.Through the government's macro-control of urban development, more people have migrated to small-and medium-sized cities, which alleviates, to a certain extent, the eco-environmental and social pressures caused by excessive population growth in large-sized cities.Meanwhile, owing to the gradual elimination of the rural and urban household registration system since 1997 [71], people migrating from rural to urban areas can obtain more and better social welfare and benefits than before.Our research found that econo-demographic development elevated dramatically, with an increase of 1402.86%, and the social well-being level increased by 78.80% in the WTSUA from 1992 to 2013.This has contributed to the promotion of social stability in urban areas while stimulating urban sprawl.In 2006, the Chinese government issued a new General Land Use Plan (2006-2020) and strengthened its control over the scale of urban areas in large-sized cities [72].Moreover, as eco-environmental issues have become more widely recognized, a variety of environmental conservation policies and plans have been

Urban Development Associated with State-Led Policies and Market-Oriented Land Reform
Urban areas, considered as the core centers of social and economic development, are crucial for ensuring national stabilization and people's livelihoods.Since the implementation of 'Reform and Opening-up', China has experienced rapid urbanization [62,63].The ratio of urbanization was less than 18% in 1978 but up to 36% in 2000; it will be more than 80% by 2050 [8].In this context, the urban system has driven unprecedented development, accompanied by urban spatial expansion and improvements in the urban structure.In the WTSUA, our results suggest that the distribution and expansion speed of urban areas have increased, with a rise of 1602.86% from 1992 to 2013 (Figures 2 and 3, and Table 4).Previous studies in China concerning the urban sprawl dynamics in Beijing, Shanghai, Guangzhou, and Wuhan show that these cities present a similar trend of urban expansion [64][65][66][67].In the process of rapid urban expansion, a series of issues associated with the eco-environment and society have arisen [68][69][70].Facing the challenges of urban sprawl, several plans and policies have been promulgated to promote coordinated urban development.Figure 9 illustrates the characteristics, emphasis, and pivotal policies and events in the different urban sprawl stages of China.It shows that the guiding principles of national urban development were issued and implemented in the early 1980s, which aimed to regulate orderly urban development by actively developing small-sized cities, slowly developing medium-sized cities, and controlling the scale of large-sized cities [62].Driven by these policies, small-and medium-sized cities have developed rapidly.This probably explains why the expansion speed of the noncentral cities of the WTSUA was higher than that of the central cities, and why the centroid of urban areas shifted towards noncentral cities from 1992 to 2013.Through the government's macro-control of urban development, more people have migrated to small-and medium-sized cities, which alleviates, to a certain extent, the eco-environmental and social pressures caused by excessive population growth in large-sized cities.Meanwhile, owing to the gradual elimination of the rural and urban household registration system since 1997 [71], people migrating from rural to urban areas can obtain more and better social welfare and benefits than before.Our research found that econo-demographic development elevated dramatically, with an increase of 1402.86%, and the social well-being level increased by 78.80% in the WTSUA from 1992 to 2013.This has contributed to the promotion of social stability in urban areas while stimulating urban sprawl.In 2006, the Chinese government issued a new General Land Use Plan (2006-2020) and strengthened its control over the scale of urban areas in large-sized cities [72].Moreover, as eco-environmental issues have become more widely recognized, a variety of environmental conservation policies and plans have been introduced to improve the eco-environment of urban systems [73].As observed from our findings, the eco-environment status increased by 98.69%, the urbanization quality exhibited an upward trend, and the coupling coordination level between urbanization quality and urban sprawl increased gradually from 1990 to 2013 (Figures 7 and 8).A previous study also stated that high-levels of economy, eco-environment, and social well-being would be conducive to achieving the desired effects of urban sprawl and ameliorating the quality of urbanization [4], which is consistent with our research.
In addition to national policies, several regional urban planning schemes are also indirectly or directly responsible for the urban expansion and urbanization quality of the WTSUA.In 1992, Fujian Province's government formulated a development strategy aiming at accelerating the economic progress of the southeast coastal zone [74].Thus, the central cities situated in the coastal zone were given priority for development driven by this policy, resulting in a gradual increase in the differentiation of rank-size distribution between central cities and noncentral cities.In 2001, the separation system of urban and rural household registration was abolished, which led to a further migration wave of rural people into different sized cities [75].With the rapid development of small-and medium-sized cities, the disparities in rank-size distribution between central cities and noncentral cities declined.In 2009, the WTSUA was formally established, which allowed systematic planning to better regulate urban sprawl; this also led to a more stable distribution of the city rank-size.Our findings showed that the change trend of the urban rank-size distribution of the WTSUA's cities could be divided into three stages, i.e., 1992-2000, 2001-2009, and 2010-2013 (Figure (Figure 5b), which was consistent with the time of regional policy implementation.
Market-oriented land reform is another common underlying force in the regulation of urban sprawl and urbanization quality; it is characterized by spontaneous institutional change induced by market-forces and initiated by private/nonstate sectors such as rural communities and township-village-enterprises [76].On one hand, feasible land reform regulated by the market can effectively perfect the redistribution of land development power.Generally, the larger the city, the higher the commuting cost and land price, which can restrain the city's expansion and reduce urban sprawl.Notably, some empirical studies on large metropolitan areas proved that the top ten most compact cities are almost mega cities [77,78].Under these circumstances, small-and medium-sized cities are more likely to have opportunities to expand, thereby optimizing the distribution of human, economic, and social resources.As referenced by Gao et al. [1], small-and medium-sized cities in China showed a more significant sprawl trend than large-sized cities from 1990 to 2010.On the other hand, compared with government intervention, market regulation is more rapid and flexible, and can respond to the process of urban expansion in a timely manner according to the current urban development situation.The market-driven mechanism can, thus, contribute to driving the urban sprawl pattern typified by the circle-layer structure, which seems to be a highly efficient method of urban development [12].From 1992 to 2013, each city in the WTSUA was found to expand in a cycloid pattern.
As mentioned above, state-led planning and market-oriented land reform have played critical roles in regulating the urban space to accommodate the increasing demand for economic and population growth while promoting the social stability of the urban system.Consequently, this dual-driving force has resulted in a relatively high degree of urbanization quality.

Lessons and Suggestions for WTSUA's Further Urban Planning
Although the coupling relationship between urban sprawl and urbanization quality of the WTSUA became more coordinated during the study period, there are still some issues that need to be addressed.With regard to central cities and noncentral cities, there were distinct spatial differences between them in urban expansion intensity and the coupling coordination level between urban sprawl and urbanization quality (Table 4, Figures 3, 4 and 8).For the central cities, the urban expansion rate was low but the coupling coordination level was high, while for the noncentral cities, the opposite pattern was observed, implying that the urban development in the eastern coastal and central-western areas has not been balanced.Furthermore, with continuing urban expansion, cities, especially large-sized cities, can be disrupted by many issues, such as population dilatation, resource depletion, environmental deterioration, traffic congestion, and housing shortages.In contrast, some smalland medium-sized cities have undergone slow economic development, population loss, and limited urban functionality [87].Our findings showed that the eco-environment status and social well-being level count for about 60% of the urbanization quality for the WTSUA (Figure 6), indicating that the improvement of urbanization quality should pay more attention to eco-environment and social well-being.In addition, some land-use regulations and land reform arbitrarily segregate "rural" and "urban" land, which is not conducive to the integrated and coordinated development of rural and urban areas [88].
According to the results of this study, some suggestions can be proposed to guide the formulation and implementation of sustainable urban development planning.First, the government should propose more feasible policies and measures to strictly limit the expansion of large-sized cities and encourage the development of small-and medium-sized cities. Decision-makers and scientists should develop more effective plans to solve city problems and stop the uncontrolled expansion of urban areas.Second, it is necessary to establish an effective mechanism for the government to intervene reasonably in the market economy, which would avoid or reduce the further deterioration of city functions.Third, the "Desakota" regions, referred to as areas with an intense mixture of nonagricultural and agricultural activities that often stretch along corridors between large-sized city centers, should be the focus of attention, because the development processes of these regions will directly or indirectly influence urban sprawl [89].Fourth, eco-environmental conservation and social welfare should be continued to be reinforced and improved by promulgating related plans and regulations, with the aim of achieving regional sustainable development.Fifth, existing effective urban sprawl patterns should be used as a reference in future urban planning.

Advantages and Limitations of the Study
Previous studies have focused on the relationship between urban sprawl and some attributes involved in urban development, such as revenue, pollution, land-use change, and population migration.Nevertheless, urban areas form large-scale complex systems, which not only have their own, inherent elements, but are also closely associated with other open ecosystems.Therefore, it may be difficult to discern the elements of urban development thoroughly by focusing only on one or a few correlative factors of urban sprawl.In our study, urbanization quality was assessed comprehensively, and the coupling relationship between urban sprawl and urbanization quality was quantified systemically.This will improve our understanding of the causal relationship between these parameters.
The efficient and accurate determination of indicator weight is crucial to the evaluation process and result.Due to the differentiation of diverse regions in terms of resource endowment, development policy, and technology level, the weight of each indicator in the evaluation system would vary from region to region for the urbanization quality.Hence, based on the subjective weight methods, it could be difficult to obtain unbiased results if the expert lacks a systemic understanding of the study area.However, the EM model used in the research, as an objective weight approach, not only eliminates the influence of the subjectivity of experts and the differences of evaluated objects, but also reduces the time spent establishing an expert panel [9,18].Therefore, the weight approach adopted in the study has a certain applicability and advantages in assessing urbanization quality in the WTSUA and even in other regions.
Many studies involving urban area extraction have used traditional optical remote sensing images as the basic data sources.However, the processing of information based on these data are more often than not time-consuming, especially for some large-scale study areas.Moreover, in previous studies, the thresholds of urban area extraction based on nighttime light data were usually determined by subjective judgment or with the aid of auxiliary data, which would be not accurate enough and omit some urban area information [28].In contrast, the long-term urban patches in our study were automatically identified by a thresholding approach, i.e., the OTSU method, which was more effective, labor-saving, and precise.Furthermore, the frequency of DMSP/OLS nighttime light observations is greater, making it possible to detect the dynamic trends of urban sprawl and undertake in-depth analyses of the driving mechanisms of urban sprawl at a regional scale.
Image resolution is a pivotal prerequisite for geo-object recognition, as it largely determines the ability to identify different objects.However, because of its low spatial resolution, DMSP/OLS nighttime light imagery may not reflect urban information well at a small-scale.As a result, a loss of precision might occur when extracting details of small urban areas.Thus, it is imperative to improve the resolution of the DMSP/OLS nighttime light data by using high-resolution optical images in future research.In addition, other limitations existed in the time-series DMSP/OLS nighttime light images.The nighttime light information recorded on the DMSP/OLS dataset only became available in 1992, and has been unavailable since 2014, leading to a lack of data for the study of the urban dynamics of previous and more recent years.Therefore, it is necessary that a longer time series dataset be formed by fusing multisource night light data; this would be conducive to more scientific and reasonable analyses of the dynamics of urban expansion.

Conclusions
Analyses and assessments of urban sprawl and urbanization quality for urban agglomerations are essential for defining the urban sprawl dynamics, understanding the reciprocal relationship between urban sprawl and urbanization quality, and underpinning knowledge-based urban development planning.In the study, taking the DMSP/OLS dataset as the basic data source, the spatio-temporal dynamics of urban areas were determined using the UERI, SDE model, and a RSR model for the WTSUA from 1992 to 2013.Based on an EW model and a CCD model, the urbanization quality and the coupling coordination level between urban sprawl and urbanization quality were evaluated and quantified, respectively.Our results showed that the urban areas experienced a significant increase in a cycloid expansion pattern, and that their rank-size distribution was uneven and changed as a result of the implementation of different policies.As urbanization quality increased, the coupling relationship between urban sprawl and urbanization quality gradually increased from serious incoordination to basic coordination.Two main driving forces, i.e., state-led policies and market-oriented land reform, played crucial roles in promoting the coupling coordination development between urbanization quality and urban sprawl in the WTSUA.These results could improve our understanding of urban sprawl dynamics and their coupling relationship with urbanization quality.These conclusions can be used as a guide to effectively facilitate orderly urban expansion and urbanization quality improvement.

3 . 1 .
Spatio-Temporal Characteristics of Urban Sprawl 3.1.1.Dynamic Change of Total Nighttime Light Inventory By calibrating the DMSP/OLS data for each year systemically, we could compare and analyze all the nighttime light images of the WTSUA during the period 1992-2013 under a unified standard.

Figure 3 .
Figure 3. Spatio-temporal dynamics of standard deviation ellipse (SDE) elements for urban areas from 1992 to 2013.(a) spatial distribution and centroid migration route of SDE from 1992 to 2013; (b) areal change trend of SDE spatial distribution from 1992 to 2013; (c) azimuth change of SDE from 1992 to 2013.

Figure 3 .
Figure 3. Spatio-temporal dynamics of standard deviation ellipse (SDE) elements for urban areas from 1992 to 2013.(a) spatial distribution and centroid migration route of SDE from 1992 to 2013; (b) areal change trend of SDE spatial distribution from 1992 to 2013; (c) azimuth change of SDE from 1992 to 2013.

Figure 4 .
Figure 4. Spatio-temporal sprawl pattern of the WTSUA from 1992 to 2013.(a) spatio-temporal distribution of urban areas of the WTSUA from 1992 to 2013; (b) expansion dynamics of urban areas in Wenzhou city; (c) expansion dynamics of urban areas in Fuzhou city; (d) expansion dynamics of urban areas in Quanzhou city; (e) expansion dynamics of urban areas in Xiamen city; (f) expansion dynamics of urban areas in Shantou city.

Figure 4 .
Figure 4. Spatio-temporal sprawl pattern of the WTSUA from 1992 to 2013.(a) spatio-temporal distribution of urban areas of the WTSUA from 1992 to 2013; (b) expansion dynamics of urban areas in Wenzhou city; (c) expansion dynamics of urban areas in Fuzhou city; (d) expansion dynamics of urban areas in Quanzhou city; (e) expansion dynamics of urban areas in Xiamen city; (f) expansion dynamics of urban areas in Shantou city.3.1.3.City Rank-Size Distribution Change of Urban Agglomeration System The urban rank-size distribution dynamics are shown in Figure 5.The value of the estimated Pareto coefficient ranged from 0.58 to 0.78, indicating an uneven city size distribution in the WTSUA, namely, there were a small number of large-sized cities.According to the temporal change of the Pareto coefficient, three different stages were defined: 1992-2000, 2001-2009, and 2009-2013.From 1992 to 2000, despite some small fluctuations, the Pareto coefficient increased, indicating that more nonuniformities had occurred during this stage.From 2001 to 2009, the Pareto coefficient experienced a general decrease.It can be inferred that during this period, WTSUA's cities experienced a "balancing trend", with large-sized cities growing more slowly and decreasing in rank, while small-sized cities grew more rapidly and increased in rank.From 2010 to 2013, the Pareto coefficient was basically stable with a lower value, suggesting a more balanced distribution of urban rank-size than the previous two stages.

Figure 5 .
Figure 5. Temporal changes of the urban rank-size distribution and Pareto coefficient (q) in the WTSUA from 1992 to 2013.(a) change trend of the urban rank-size distribution; (b) change trend of Pareto coefficient (q).

Figure 5 .
Figure 5. Temporal changes of the urban rank-size distribution and Pareto coefficient (q) in the WTSUA from 1992 to 2013.(a) change trend of the urban rank-size distribution; (b) change trend of Pareto coefficient (q).

Figure 6 .
Figure 6.Weight of each indicator for the urbanization quality evaluation system ((X1: Per capita water resources, X2: Per capita land area, X3: Per capita gas supply volume, X4: Per capita green areas, X5: Domestic sewage treatment rate, X6: Landfill rate for urban waste, X7: Proportion of built-up area and total area, X8: Per capita urban road area in municipal district, X9: Per capita public library holdings, X10: Proportion of educational expenditure to GDP, X11: Proportion of science and technology expenditure to GDP, X12: Number of teachers per 10,000 students, X13: Number of medical beds per 10,000 persons, X14: Number of buses per 10,000 persons, X15: Per capita housing, X16: Registered urban unemployment rate, X17: Proportion of population in municipal district to total population, X18: Population density in municipal district, X19: Proportion of secondary and tertiary industry employment to total employment, X20: Per capita investment in fixed assets, X21: Per capita GDP, X22: GDP, X23: Per capita GDP in the secondary and tertiary industry, X24: Per capita disposable income of urban residents, X25: Per capita local government fiscal revenue, X26: Proportion of secondary and tertiary industry GDP to total GDP, X27: Total retail sales of consumer goods, X28: Deposits of financial institutions at year-end).

Figure 6 .Figure 7 .
Figure 6.Weight of each indicator for the urbanization quality evaluation system ((X 1 : Per capita water resources, X 2 : Per capita land area, X 3 : Per capita gas supply volume, X 4 : Per capita green areas, X 5 : Domestic sewage treatment rate, X 6 : Landfill rate for urban waste, X 7 : Proportion of built-up area and total area, X 8 : Per capita urban road area in municipal district, X 9 : Per capita public library holdings, X 10 : Proportion of educational expenditure to GDP, X 11 : Proportion of science and technology expenditure to GDP, X 12 : Number of teachers per 10,000 students, X 13 : Number of medical beds per 10,000 persons, X 14 : Number of buses per 10,000 persons, X 15 : Per capita housing, X 16 : Registered urban unemployment rate, X 17 : Proportion of population in municipal district to total population, X 18 : Population density in municipal district, X 19 : Proportion of secondary and tertiary industry employment to total employment, X 20 : Per capita investment in fixed assets, X 21 : Per capita GDP, X 22 : GDP, X 23 : Per capita GDP in the secondary and tertiary industry, X 24 : Per capita disposable income of urban residents, X 25 : Per capita local government fiscal revenue, X 26 : Proportion of secondary and tertiary industry GDP to total GDP, X 27 : Total retail sales of consumer goods, X 28 : Deposits of financial institutions at year-end).

Figure 7 .
Figure 7. Urbanization quality evolution of the WTSUA from 1992 to 2013.Figure 7. Urbanization quality evolution of the WTSUA from 1992 to 2013.
Figure8visualizes the spatio-temporal dynamics of the coupling coordination relationship between urbanization quality and urban sprawl for the WTSUA from 1992 to 2013.During the study period, the level of coupling coordination increased gradually for each city.Moreover, the CCD value of the central-western inland areas was obviously lower than that of the eastern coastal regions, showing a significant spatial variation (Figure8a).As for the total WTSUA, the value of the CCD presented a gradual increase over time, suggesting the enhancement of improving coordination between urban sprawl and urbanization quality from 1992 to 2013 (Figure8b).Specifically, the coupling coordination level of urban sprawl and urbanization quality experienced four stages: a serious incoordination phase (low coupling state) from 1992-1994, a moderate incoordination phase (low coupling state) from 1995-2005, a slight incoordination phase (antagonism state) from 2006-2009, and a basic coordination phase (running-in state) from 2010-2013.It is noteworthy that although the overall coupling coordination level of the study area was basic coordination by 2013, the coupling coordination relationship between urban sprawl and urbanization quality had developed into moderate coordination or high coordination stage in some eastern coastal cities.

Figure 8 Figure 8 .
Figure8visualizes the spatio-temporal dynamics of the coupling coordination relationship between urbanization quality and urban sprawl for the WTSUA from 1992 to 2013.During the study period, the level of coupling coordination increased gradually for each city.Moreover, the CCD value of the central-western inland areas was obviously lower than that of the eastern coastal regions, showing a significant spatial variation (Figure8a).As for the total WTSUA, the value of the CCD presented a gradual increase over time, suggesting the enhancement of improving coordination between urban sprawl and urbanization quality from 1992 to 2013 (Figure8b).Specifically, the coupling coordination level of urban sprawl and urbanization quality experienced four stages: a serious incoordination phase (low coupling state) from 1992-1994, a moderate incoordination phase (low coupling state) from 1995-2005, a slight incoordination phase (antagonism state) from 2006-2009, and a basic coordination phase (running-in state) from 2010-2013.It is noteworthy that although the overall coupling coordination level of the study area was basic coordination by 2013, the coupling coordination relationship between urban sprawl and urbanization quality had developed into moderate coordination or high coordination stage in some eastern coastal cities.

Figure 8 .
Figure 8. Spatio-temporal dynamics of coupling coordination level between urbanization quality and urban sprawl for the WTSUA from 1992 to 2013.(a) spatial variations of coupling coordination level between urbanization quality and urban sprawl for the WTSUA in different years; (b) change trend of CCD between urbanization quality and urban sprawl for the entire WTSUA from 1992 to 2013.

Table 1 .
Basic information description of each city in WTSUA.

Table 2 .
Evaluation indicator system of urbanization quality.

Table 3 .
Classification type of coupling coordination degree between urban sprawl and urbanization quality.

Table 4 .
Changes in area and expansion rate of urban areas for different city types in the WTSUA.