Next Article in Journal
CW Direct Detection Lidar with a Large Dynamic Range of Wind Speed Sensing in a Remote and Spatially Confined Volume
Next Article in Special Issue
Lifting of International Sanctions and the Shadow Economy in Iran—A View from Outer Space
Previous Article in Journal
AMS-Net: An Attention-Based Multi-Scale Network for Classification of 3D Terracotta Warrior Fragments
Previous Article in Special Issue
Erratum: Liu et al. Nightlight as a Proxy of Economic Indicators: Fine-Grained GDP Inference around Chinese Mainland via Attention-Augmented CNN from Daytime Satellite Imagery. Remote Sens. 2021, 13, 2067
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Delineating Functional Urban Areas Using a Multi-Step Analysis of Artificial Light-at-Night Data

by
Nataliya Rybnikova
1,2,3,*,
Boris A. Portnov
2,
Igal Charney
3 and
Sviatoslav Rybnikov
4,5
1
Department of Mathematics, University of Leicester, Leicester LE1 7RH, UK
2
Department of Natural Resources and Environmental Management, University of Haifa, Haifa 3498838, Israel
3
Department of Geography and Environmental Studies, University of Haifa, Haifa 3498838, Israel
4
Institute of Evolution, University of Haifa, Haifa 3498838, Israel
5
Department of Evolutionary and Environmental Biology, University of Haifa, Haifa 3498838, Israel
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(18), 3714; https://0-doi-org.brum.beds.ac.uk/10.3390/rs13183714
Submission received: 10 August 2021 / Revised: 9 September 2021 / Accepted: 10 September 2021 / Published: 17 September 2021
(This article belongs to the Special Issue Nighttime Lights as a Proxy for Economic Performance of Regions)

Abstract

:
A functional urban area (FUA) is a geographic entity that consists of a densely inhabited city and a less densely populated commuting zone, both highly integrated through labor markets. The delineation of FUAs is important for comparative urban studies and it is commonly performed using census data and data on commuting flows. However, at the national scale, censuses and commuting surveys are performed at low frequency, and, on the global scale, consistent and comparable data are difficult to obtain overall. In this paper, we suggest and test a novel approach based on artificial light at night (ALAN) satellite data to delineate FUAs. As ALAN is emitted by illumination of thoroughfare roads, frequented by commuters, and by buildings surrounding roads, ALAN data can be used, as we hypothesize, for the identification of FUAs. However, as individual FUAs differ by their ALAN emissions, different ALAN thresholds are needed to delineate different FUAs, even those in the same country. To determine such differential thresholds, we use a multi-step approach. First, we analyze the ALAN flux distribution and determine the most frequent ALAN value observed in each FUA. Next, we adjust this value for the FUA’s compactness, and run regressions, in which the estimated ALAN threshold is the dependent variable. In these models, we use several readily available, or easy-to-calculate, characteristics of FUA cores, such as latitude, proximity to the nearest major city, population density, and population density gradient, as predictors. At the next step, we use the estimated models to define optimal ALAN thresholds for individual FUAs, and then compare the boundaries of FUAs, estimated by modelling, with commuting-based delineations. To measure the degree of correspondence between the commuting-based and model-predicted FUAs’ boundaries, we use the Jaccard index, which compares the size of the intersection with the size of the union of each pair of delineations. We apply the proposed approach to two European countries—France and Spain—which host 82 and 72 FUAs, respectively. As our analysis shows, ALAN thresholds, estimated by modelling, fit FUAs’ commuting boundaries with an accuracy of up to 75–100%, being, on the average, higher for large and densely-populated FUAs, than for small, low-density ones. We validate the estimated models by applying them to another European country—Austria—which demonstrates the prediction accuracy of 47–57%, depending on the model type used.

1. Introduction

More than 50% of the world’s population currently resides in urban areas, and this share is expected to increase to 70% by 2050 [1]. Due to a significant concentration of production factors, urban areas produce approximately 80% of the global GDP [2]. This makes spatial dynamics of urban areas to be important for policy-makers and researchers alike. Decision-makers can devise informed development policies, while in the research community, this information can be used to monitor the process of urban growth and the forces behind it [3,4,5], to assess the impact of urbanization on agriculture and natural landscapes [6], on biodiversity [7], on land surface temperature [8], and other socioeconomic and physical phenomena.
Urban growth is characterized by two distinctive components—physical growth and functional change. The former group of attributes reflect changes in impervious surfaces and built-up characteristics, such as building density, building volumes [9,10], as well as population size and density of individual urban settlements [11,12,13]. Concurrently, functional attributes of urban growth reflect factor mobility, associated with various economic activities, such as commuting, commerce, industrial production and services [14]. Such exchanges are especially intense between urban cores, where a large share of production factors is concentrated, and their surrounding areas. Functionally-integrated clusters, representing geographic entities that consist of a densely inhabited city and a less densely populated commuting zone, both highly integrated through labor markets, are commonly referred to as functional urban areas or FUAs [15]. A FUA is conceptually different from an urban agglomeration, which is commonly defined as a major city surrounded by an adjacent hinterland [16]. The major difference between the two is commuting, which is crucial for delineating FUAs, but is not a prime consideration for the definition of urban agglomerations.
According to the mainstream approach adopted by the European Union (EU) and the Organization for Economic Co-operation and Development (OECD), the boundaries of FUAs are defined in three consecutive steps. First, urban cores are identified as contiguities of high-density grid cells with population density of at least 1500 residents per km2 and the total population in the contiguous cells of at least 50,000 residents. Second, local administrative units (LAUs) with at least 50% of their residents living inside the urban core are identified. At the final step, the commuting zone, comprising LAUs, which have at least 15% of their residents employed in the core city, is determined. Together with the central city, these administrative units are assumed to form a single FUA [17].
However, commuting data, needed to perform such delineations, are laborious to collect and are infrequent and sporadic even in developed countries [18]. In addition, different countries and regions report communing data with different frequencies, and sometime collect them using different definitions and methodologies [18]. As a result, comparable cross-country estimates of FUA boundaries cannot always be obtained.
As artificial light-at-night (ALAN) data are freely available globally and provide a seamless global coverage, the idea of using them for the identification of human activities was investigated in several studies (see inter alia [19,20,21]). In previous studies, ALAN data were used in health geography [22,23,24,25,26,27], for the analysis of economic performance of countries and regions [28,29,30,31], and in population density research [20,32,33,34,35,36]. The use of such data in the studies of light pollution and its ecological effects is also common [24,37,38,39,40].
In recent years there have been attempts to use ALAN data for the identification of urban areas [20,21,41,42,43,44,45,46,47]. In one such study, Imhoff et al. [41] examined frequency-based ALAN thresholds for three large metropolitan areas in the U.S.—Miami, Chicago and Sacramento. After the authors analyzed the frequencies of differently lit pixels in the ALAN images, they determined that pixels present with 85%, 89% and 94% frequencies, occupy the areas of approximately same size, such as those reported in the Census for the corresponding metropolitan entities.
In another study, Sutton et al. [20] investigated 2000 cities across the globe, and compared their actual boundaries with those produced by three different frequency-based ALAN thresholds—40%, 80% and 90%. As the study revealed, pixels in the ALAN image, observed with a frequency of 80% or more, correspond to the actual municipal boundaries best, reaching a correlation level of about 68%.
In a separate study, Henderson et al. [21] examined frequency- and intensity-based ALAN thresholds that match the boundaries of San Francisco, Beijing and Lhasa. As the authors of this study have found, the optimal ALAN frequency-based thresholds that produce the total lit area comparable in size to the Landsat data-derived urban delineations, reach 88% for Lhasa, 97% for Beijing, and 92% for San Francisco, with the corresponding ALAN flux being equal to 19, 30 and 51 digital numbers (DN), respectively. However, the spatial correspondence between metropolitan boundaries, determined using ALAN thresholds, and actual metropolitan delineations was found to be relatively low, not exceeding 8–44%.
It should also be noted that the aforementioned studies focus on the identification of built-up urban contiguities, while, to the best of our knowledge, only one study by Bosker et al. [18] analyzed functional urban delineations based on commuting flows. The authors of this analysis compared varying percentiles of ALAN intensities, reported by the VIIRS/DNB satellite’s sensor for 2015, with commuting delineations in Malaysia. As this study revealed, the best fit of ~40% is observed when 7% commuting frequency delineations are compared with delineations based on the 25th percentile of ALAN intensities.
A possible reason for such a low fit of less than 40% is that FUAs even in the same country differ by the amount of ALAN they emit. As a result, different ALAN thresholds must be used for the delineation of FUA boundaries in different parts of the urban system. In Figure 1, we illustrate this point using two FUAs in France, as an example. As evidenced by this figure, the ALAN threshold of 0.71 nW/cm2/sr fits reasonably well the boundaries of the Paris FUA, but the same threshold fits rather poorly the much smaller Chateauroux FUA, ALAN flux at which boundary does not exceed 0.15 nW/cm2/sr.
Considering that ALAN emissions from different FUAs vary substantially, it is thus important to establish varying ALAN thresholds, which would fit individual FUAs. This task can be performed for each FUA separately. However, in order to be practical, the approach needs to be sufficiently general, to enable its application to different FUAs, both for countries and regions with well-established commuting data and for other locations with unavailable or sparsely available commuting information. In this paper, we develop such an approach and test it against actual FUA delineations.

2. Materials and Methods

2.1. Study Phases

The proposed approach is implemented in several steps, as detailed in Figure 2. The data sources and analysis stages are described in the subsections below.

2.2. Data Sources

Data for the present study were drawn from the following four main sources:
(1)
The ALAN raster maps of France and Spain (see Figure 3), used in the study for model training and validation, and ALAN raster for Austria, used for additional validation of the models’ performance, were clipped from 2015 radiance-calibrated ALAN image downloaded from the VIIRS/SNPP website [48]. The ALAN data used in the study are free of background noise, solar and lunar contamination, and also free from data degraded by cloud cover, and features unrelated to electric lighting (e.g., fires, flares, volcanoes) [49]. In addition, the data underwent an outlier removal procedure, applied to abnormally high radiance pixels that occur infrequently over a year [49]. The image in question is the closest temporal match for other data sources used in the analysis, specifically for the FUA delineations, available for 2011 only (Figure 4). Although ALAN images are available today from the VIIRS-SNPP website on a monthly basis, and, since 2018, as daily composites [50], we opted to use an annual composite image, so as to minimize disturbances resulting from ALAN seasonal fluctuations and weather conditions, such as, e.g., cloud cover, which are often present in monthly and daily composites [50]. The subject image is of a ~500 × 500 m spatial resolution and reports the summarized intensity of nighttime light in nW/cm2/sr for different wavelengths in the 500–900 ηm diapason [46]. In the image, ALAN levels vary from 0 to 4187 nW/cm2/sr for France, and from 0 to 550 nW/cm2/sr for Spain (see Figure 3 and Table 1).
(2)
Boundaries of FUAs and their cores (see Figure 4) were obtained as shapefiles from the OECD website [51]. These shapefiles are generated using GeoStat grids, based on 2011 commuting data reported in national censuses [51].
(3)
The latitudes of the FUA cores’ centroids and distances to the closest major city, used to explain the variance of the optimal ALAN thresholds, were calculated using the above FUA cores’ shapefiles by applying ArcGIS−10.x software tools.
(4)
Population density of the FUA cores, and population densities of their 5–15–25 km buffers, also used as explanatory variables for the estimation of the optimal ALAN thresholds, were calculated using 1 × 1 km population grids obtained from the LandScan database for 2011 [52].

2.3. Initial Determination of the ALAN Thresholds

For the sake of simplicity, let’s assume that the nighttime light source of highest intensity is located at the center of a FUA, and light intensities drop monotonically and uniformly towards the FUA’s periphery (see Figure 5).
Such an assumption might be fully plausible for compact and monocentric urban areas (Figure 6). Under these conditions, the territorial footprint of the FUA’s ALAN emissions follows a perfect circle, and the most frequently observed (i.e., modal) ALAN values are found at the FUA’s outer boundary (Figure 5a). These modal values are also the dimmest ones, and, as such, they effectively define the FUA’s outer boundary (Figure 5b).
If the above assumptions are upheld, the analysis of the frequency distribution of the observed ALAN values can help to identify the ALAN level, which coincides best with the FUA’s boundary. In particular, the researcher needs to choose the modal ALAN value, for which ALAN intensity is expected to be close to zero (Figure 5b).

2.4. Correction for Compactness

The above assumption of monotonic and concentric distribution of ALAN emissions (Figure 5) is upheld only if the boundaries of FUAs that are circularly shaped. However, if a FUA’s shape is not circular, using the modal ALAN value as a delineation threshold would underestimate the actual area of the FUA. Figure A1 in Appendix A, which reports different FUAs’ footprints, helps to illustrate this point. As this figure shows, the more distant the shape of a FUA from a perfect circle, the brighter ALAN values emerge as the most frequent. For such non-circular FUAs, it is thus necessary to correct for compactness, so as to account for a FUA’s shape deviation from a perfect circle.
To perform such a correction, we first estimate the FUA’s compactness (c), calculating it as the ratio between the area of a FUA and the area of its bounding circle [53,54]:
c F U A = S F U A S B C
where SFUA = area of a FUA; SBC = area of the bounding circle, calculated using the Minimum Bounding Geometry tool in the ArcGIS software.
Next, to represent FUAs, which deviate from circular shapes, we model them as ellipses of the same compactness:
c E l = S E l S B C = π a b π a 2 = b a = c F U A
where SEl = area of an ellipse with semi-axes a and b (a > b).
At the next step, to correct the initially estimated ALAN threshold (see Section 2.3) for a FUA’s compactness, we calculate the radius of the circle, r, which has the maximal intersection with ellipse, CEl. As shown in Box A1 in Appendix A, this radius is equal to:
r = a b
Lastly, we estimate the percentile of the ALAN value distribution, p*, corrected for compactness (see Box A1 in the Appendix A for the justification):
p * = 2 π a r c s i n 1 c 1 + c
According to (4), for compact shapes, which are close to a circle, i.e., for which c → 1, the optimal ALAN threshold percentile (p*) tends to the dimmest ALAN value (p* → 0), while for prolongated shapes with c → 0, p* → 100, that is, the optimal ALAN threshold will tend to the highest ALAN percentile (see Figure 7).

2.5. Regression Modelling

After the optimal ALAN threshold is identified for each FUA by determining the modal ALAN value (see Section 2.3), and corrected for compactness (Section 2.4), we link the estimated threshold values to several explanatory variables, characterizing the FUA cores, so as to determine these variables’ load on the optimal ALAN threshold value. To model these relationships, the following generic regression equation is used:
A L A N i = b 0 + b 1 * L a t i + b 2 * D i + b 3 * P D i + b 4 * P D D i + ε i
where ALANi is the optimal ALAN threshold for FUA i (nW/cm2/sr); Lati is latitude of the FUA core’s centroid (decimal degrees, dd); Di is distance to the nearest major city, calculated between a given FUA core’s centroid and the centroid of the nearest FUA with more than 1.5M residents (dd); PD is population density of the FUA core (persons per km2); PDD is population density decline gradient, calculated as the ratio between the FUA core’s population density and population density in the FUA core’s buffer with a 5 km width for small FUAs (under 100,000 residents), a 15 km width for mid-sized FUAs (100,000–250,000 residents), and a 25 km width for large FUAs (over 250,000 residents); b0..b4—regression coefficients, and ε is a random error term.
The predictors used in the model are expected to contribute to the ALAN threshold’s variance due to varying reasons. In particular, population density is known to be closely associated with ALAN flux (see inter alia [33,35,55]). Concurrently, population density gradient might capture changes in the pattern of population density around the FUA core. Concurrently, distance to the nearest major city is likely to show how local development patterns are modulated by proximity to major urban concentrations [55]. In addition, as population concentrations in high latitudes often require more artificial illumination, especially during long winters [40], FUA’s latitude is also included into the model as a potential predictor.
In the analysis, we tested different functional forms of the models, and determined that the logarithmic transformation of the PD and PDD variables provides the best results, by improving the regression fit substantially (p < 0.05). The initial analysis was performed in the IBM SPSSv.25 software using its multiple regression module. To ensure the normality distribution of the dependent variable, ALANi, we applied Box-Cox transformation procedure, to redefine the ALAN thresholds [56].
In addition to ordinary least square regressions (OLS), we also tested “random forest” regressions. Such regressions imply building an ensemble of “decision trees”, each of which “voting” for a certain level of the dependent variable, with subsequent averaging of the estimates across all the decision trees [57]. In the present analysis, we implemented a standard realization of the “random forest” regression (the TreeBagger module) in the MATLAB v.R2020x software [58]. During the estimation procedure, two parameters were a matter of choice—the number of independent variables used for the individual decision tree construction and the number of decision trees that comprise the forest. To ensure the comparability of the results, we used all independent variables, covered by the analysis, for the decision trees’ construction, and defined number of trees to be equal to 100, which is usually considered to be a reasonable number for reaching a generalization error convergence (see for example [57,59]). Each decision tree was built for 80% of randomly selected observations.

2.6. Adjustment for Contiguity

When the analysis is performed, any given ALAN threshold level might identify several clusters of identically lit pixels, some of which might be related to a given FUA, while other pixels might be located elsewhere. Therefore, to identify the ALAN pixels relevant to a given FUA, the following analytical procedure was implemented. First, for each FUA, we identified pixels that overlap the FUA’s core area, considering the core boundary information as an initial input (see Section 2.2: Data Sources). Next, for each pixel selected thereby, we analyzed all the pixels in its surroundings. If the ALAN value of a neighboring pixel was lower or equal to that of the pixel under analysis but greater than the ALAN threshold identified for the FUA (see Section 2.3 and Section 2.4), the pixel in question was considered to be a part of the FUA analyzed. We have continued this procedure as long as all the pixels, which satisfy the above criteria, maintained a spatial contiguity. Then, for each FUA, we selected local administrative areas (LAUs), most of which area (that is, >50%) is occupied by the pixels identified thereby. These LAUs were considered to be a part of a given FUA (the MATLAB code for contiguity adjustment can be obtained from the authors upon request).

2.7. Initial Validation

To assess the performance of the estimated models (see Section 2.5), we analyzed the degree of correspondence between the empirically determined (see Section 2.4) and model-predicted ALAN thresholds adjusted for contiguity (see Section 2.6). To this end, the model estimated for France was used to predict the ALAN thresholds for individual FUAs in Spain and vice versa. In order to assess the extent to which the empirically determined and model-predicted ALAN thresholds coincide, we used different metrics, including Pearson correlation coefficients, standard error of the estimates (SEE), and weighted mean squared errors (WMSE).
Next, we compared the FUAs’ delineations, either empirically determined using the modal ALAN values (see Section 2.4) and adjusted for contiguity (see Section 2.6), or model-predicted (see Section 2.5), with commuting-based FUAs’ delineations (see Section 2.2). To perform such a comparison, we used the Jaccard Index (JI), which estimates the share of intersection within the union of the two sets relative to these shapes’ union [60]:
J I F U A C , F U A T = F U A C F U A T F U A C F U A T
where FUAC = the set of local autonomous units (LAUs) forming a FUA defined by commuting, and FUAT = set of LAUs within either an empirically determined or model-predicted FUA boundary. The value of the index in question ranges from zero, when no intersection between the two sets is present, to one, when the two sets completely coincide and their intersection is equal to their union [60].

2.8. Second-Step Validation

For an additional validation, we applied the models estimated for France and Spain to FUAs in another European country—Austria (Figure 8).
Although Austria differs from the two other countries under analysis in terms of size, urbanization level, topography, and FUAs’ location, it was chosen for an additional model validation, to demonstrate that the estimated models perform reasonably well even in this specific case. As all FUAs in this country are located apart from each other (see Figure 8), this country is considered particularly suitable for the intended validation.
The validation procedure was carried out in the following four steps. First, we determined the optimal ALAN thresholds for each FUA empirically (see Section 2.3 and Section 2.4). Second, we used the ALAN-threshold identification models, estimated for France and Spain (see Section 2.5), to predict optimal ALAN thresholds for FUAs in Austria, using relevant input variables (see Section 2.2 and Section 2.5), and, then, adjusted these estimates for contiguity (see Section 2.6). Considering that FUAs in Austria are located in close proximity to international borders, the input information was not limited to the areas inside Austria only. For instance, population density-decline gradient and distance to the closest major city were calculated regardless of the state borders. Third, we assessed the correspondence between the empirically determined and models-predicted ALAN thresholds using Pearson correlation coefficients, SEE, and WMSE. Finally, we compared delineations, based on model-predicted ALAN thresholds, with commuting-based delineations, while expanding the study area by a 50-km buffer around the Austrian border, to cover the parts of FUAs located outside Austria and potentially extending into neighboring countries. As in the previous stage of the analysis (Section 2.7), the comparison of the shapes was performed using JI.

3. Results

3.1. Optimal ALAN Thresholds

The descriptive statistics of the ALAN thresholds, estimated by the multi-step approach described in Section 2.1, are reported in Table 2, separately for France and Spain, both as ALAN percentiles and actual ALAN levels in nW/cm2/sr. As evidenced by this table, the optimal ALAN thresholds identified for individual FUAs appear to vary widely, ranging from 0.15 to 9.91 nW/cm2/sr for France, and from 0.13 to 8.23 nW/cm2/sr for Spain.
The most frequent (i.e., modal) ALAN values for all FUAs under analysis are reported in Figure 9, separately for France and Spain. In both countries, the modal ALAN values are not identical to the dimmest ones, thus pointing out that the “circularity” assumption (see Section 2.3) is violated. The bottom sub-figures report ALAN thresholds, corrected for compactness using the approach described in Section 2.4. As it can be seen from the comparison of the upper and bottom diagrams, modal ALAN thresholds, corrected for compactness are closer to the dimmest ALAN values than before the correction (especially for France), albeit differences in distributions are still valid.

3.2. Explaining the Variance of the Observed ALAN Thresholds

In Table 3, we report the results of OLS analysis, linking individually determined optimal ALAN thresholds with geographic and socio-economic attributes of the FUAs’ core areas. As evidenced by Table 3, the predictors used in the analysis help to explain ~74% of the ALAN threshold variance (R2 = 0.739–0.740). Characteristically, in both models, significant predictors are nearly identical and exhibit the same signs: population density (+); population density gradient (−); latitude (+), and distance to the nearest major city () (p < 0.01).
As random forest regressions do not provide explicit estimates of the explanatory variables’ coefficients, we do not report these models here, but should remark that these estimates in terms of correlation with the actual ALAN threshold levels are similar to the OLS estimates reported in Table 3 (r = 0.856 for France and r = 0.883 for Spain, as opposed to r = 0.866 for France and r = 0.812 for Spain in the OLS models), while in terms of SEE they are poorer (SEE = 0.913 for France and SEE = 0.817 for Spain in comparison to SEE = 0.533 for France and SEE = 0.629 for Spain; see Table 3).
However, in terms of WMSE, random forest regressions are much better (WMSE = 0.945 for France and WMSE = 0.545 for Spain in comparison to WMSE = 4.521 for France and WMSE = 2.718 for Spain in the OLS models; see Table 3). Considering this result, we use the ALAN threshold estimates, produced by the random forest regressions, in the following analysis.

3.3. Model Cross-Validation

In Figure 10, we report the correspondence between the empirically determined and model-predicted ALAN thresholds. For this analysis, the model estimated for France (see Table 3) is applied to FUAs in Spain and vice versa. As evidenced by this figure, the estimates are fairly congruent, with r > 0.819.

3.4. Model-Estimated vs. Commuting-Based FUAs’ Delineations

Figure 11 shows several most successful examples of FUAs’ delineations, generated by the proposed approach. Concurrently, in Figure 12, we report actual FUA delineations and model estimates for all FUAs in continental France and Spain. In addition, in Table 4, we report the degree of correspondence between the model-estimated and commuting-based delineations, assessed using JI.
As can be seen in Table 4, the calculated JI values range between 0.30 and 0.64, being higher for large FUAs (JI = 0.499–0.507) than for small FUAs (JI = 0.33–0.34). For densely populated FUAs, the match between the commuting-based and ALAN-based delineations is especially high, reaching 0.557–0.638, or 56–64% (see Table 4). [The JI values for all the French and Spanish FUAs are reported in Figure A2 in Appendix A].

3.5. Second-Step Validation

In Table 5, we report ALAN threshold values for FUAs in Austria, calculated using the ‘French’ and ‘Spanish’ models (Table 3), and compared to individually fitted ALAN thresholds. As evidenced by this table, the ALAN thresholds, estimated using the French and Spanish models, correspond to the individually fitted ALAN thresholds quite well, with r > 0.77 and SEE < 0.82. Yet, in terms of WMSE, the French model performs poorer in comparison to the Spanish model (WMSE = 0.711 vs. WMSE = 10.102, respectively). In Figure 13, we report FUAs’ delineations obtained by averaging the estimates obtained using the French and Spanish models (see Table 3).

4. Discussion and Conclusions

The delineation of geographic boundaries of FUAs is important for comparative urban studies. However, using commuting data for this task is not always feasible due to difficulties in data collection. In the present study, we suggested and tested an approach, based on the analysis of ALAN data. As ALAN is emitted from roads, frequented by commuters, and by buildings surrounding roads, ALAN emissions can be used, as we hypothesize, for the identification of FUAs.
We verify this hypothesis using data on commuting-based delineations available for France and Spain, applying a multi-step approach. First, we fit the ALAN threshold for each individual FUA, using the modal value of the ALAN frequency distribution. Next, we explain this threshold by a multiple regression analysis, using several characteristics of the FUAs’ cores, such as latitude of the core’s centroid, distance to the closest major city, population density, and density decline gradient. Although the boundaries of the FUA core areas used as an initial input are not generated by the analysis per se, such boundaries, if not a priori available, can be identified easily using Global Human Settlement [61] or LandScan [52] as contiguities of densely populated grids. Lastly, we cross-validate the obtained models for three European countries.
As our analysis indicates, the degree of correspondence between the individually fitted and model-predicted ALAN thresholds is relatively high (r > 0.819), with Jaccard Index values reaching up to 75% for France and up to 100% for Spain.
Our results are more robust than those obtained by Bosker and colleagues [18] for FUAs in Indonesia, according to which the correspondence between ALAN-based and 15% commuting-based FUA delineations did not exceed 28%. We explain the improvement, obtained in the present study, by the use of individually-fitted ALAN thresholds, based on the analysis of modal values, corrected for compactness.
To the best of our knowledge, this study is the first that estimates the optimal ALAN thresholds that approximate the boundaries of individual FUAs, using readily available, or easy-to-compute, characteristics of the FUAs’ cores, such as latitude of the core’s centroid, distance to the closest major city, population density and population density decline gradient, combined with ALAN flux data.
The proposed modelling approach might be useful for FUA delineations in countries and regions, for which commuting data are unavailable, as well as for places, in which commuting data are not updated on a regular basis, and for a comparative analysis of countries and regions, which use different commuting-assessment procedures. Using our modeling approach, FUAs’ boundaries can be determined in the following steps. First, the boundaries of FUAs’ cores should be identified. If such boundaries are not readily available, they can be determined as contiguities of high-density grid cells, using input sources, such as Global Human Settlement [61] or LandScan [52] grids. The procedure might follow the algorithm described in Dijkstra et al. (2019): Grid cells with population density of at least 1500 residents per km2 are identified. Afterward, the grid cells identified thereby are grouped into contiguous area with a total population of at least 50,000 residents. For such areas, the development and locational characteristics are identified next, including the latitude of the contiguity’s centroid, distance to the closest major city, population density and population density decline gradient (see Section 2.5). These characteristics of the core areas are then used as predictors in the ALAN-threshold estimation models, reported in Section 2.5 for either France or Spain, or both, to obtain the optimal ALAN threshold estimates for each individual FUA. Finally, a VIIRS-DNB raster is used, to select pixels, corresponding to the estimated ALAN threshold, and to identify LAUs associated with such pixels’ contiguities, as detailed in Section 2.6.
The present study has several limitations. While for some FUAs, our estimates are quite accurate, reaching the levels of accuracy of 74–100%, for other, typically smaller FUAs, our estimations are less accurate. We assume that the reason might be that commuting-based boundaries rely mainly on work-related commuting, while omitting other human flows, such as travels for leisure, services, and social activities. In contrast, the suggested ALAN approach captures human activities at large. In addition, the ALAN-based approach might omit areas occupied by functions that operate mainly at daytime and emit much less light at night. For smaller FUAs, this source of error might by more pronounced than for large FUAs, where many functions operate around the clock. Another possible reason for a relatively low correspondence between some commuting- and ALAN-based delineations might be due to the fact that many FUAs are not monocentric, or might have a shape which is far from circular or elliptic, which we considered for modelling. For such cases, further studies might be needed to reflect more complex situations, in which FUA is either polycentric, or adjacent to other FUAs and their boundaries overlap or merge.
It should also be noted that in this study, we investigated the performance of the proposed method by applying it to three well-developed countries in Europe—France, Spain, and Austria. Yet, question remains about the models’ applicability to countries outside Europe and to countries in mid-latitudes, and, especially, to less-developed countries. We expect that applying the models to such countries might result in the overestimation of the optimal ALAN thresholds and thus in the underestimation of the commuting extent (the evidence for this conclusion is provided in [21]). Therefore, a follow-up investigation of the applicability of the proposed models to less developed countries and regions might be needed. Additionally, we need to acknowledge that a temporal mismatch between ALAN and actual FUAs’ delineation exists. Whether it might affect the results of the analysis should be clarified in future studies, after newer commuting data become available.

Author Contributions

Conceptualization, B.A.P. and I.C.; methodology, N.R., S.R. and B.A.P.; software, S.R. and N.R.; formal analysis, N.R.; data curation, N.R.; writing—original draft preparation, N.R.; writing—review and editing, B.A.P., I.C. and N.R.; overall study supervision, B.A.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Council for Higher Education of Israel (postdoctoral scholarship of N.R.).

Data Availability Statement

Initial data and processing codes are available from N.R. upon request.

Acknowledgments

We express our gratitude to three anonymous reviewers for the highly valuable comments.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

Table A1. FUAs’ identification (ID number = number in axis X in Figure 9 in the main manuscript).
Table A1. FUAs’ identification (ID number = number in axis X in Figure 9 in the main manuscript).
Upper Subfigures
ID NumberFUA Code FUA NameID NumberFUA Code FUA Name
FranceSpain
1FR067Quimper1ES053Ciudad Real
2FR061Niort2ES538Avila
3FR021Poitiers3ES546Merida
4FR053La Rochelle4ES016Toledo
5FR051Troyes5ES057Ponferrada
6FR086Evreux6ES527Jaen
7FR066Saint-Brieuc7ES050Manresa
8FR077Roanne8ES011Santiago de Compostela
9FR035Tours9ES040Talavera de la Reina
10FR059Chalon-sur-Saone10ES540Chiclana de la Frontera
11FR093Brive-la-Gaillarde11ES545Lorca
12FR104Chalons-en-Champagne12ES528Lleida
13FR076Belfort13ES529Ourense
14FR025Besancon14ES059Zamora
15FR073Tarbes15ES523Leon
16FR505Charleville-Mezieres16ES031Lugo
17FR038Le Mans17ES043Ferrol
18FR037Brest18ES034Caceres
19FR068Vannes19ES515Burgos
20FR096Albi20ES519Albacete
21FR050Montbeliard21ES014Pamplona
22FR074Compiegne22ES041Palencia
23FR022Clermont-Ferrand23ES542Basin
24FR023Caen24ES017Badajoz
25FR506Colmar25ES544Linares
26FR036Angers26ES510Donostia-San Sebastian
27FR019Orleans27ES062Sanlucar de Barrameda
28FR049Lorient28ES033Girona
29FR058Chambery29ES013Oviedo
30FR069Cherbourg30ES009Valladolid
31FR018Reims31ES022Vigo
32FR090Chateauroux32ES054Benidorm
33FR056Angouleme33ES501Granada
34FR063Beziers34ES044Pontevedra
35FR020Dijon35ES537Alcoy
36FR064Arras36ES514Almeria
37FR057Boulogne-sur-Mer37ES552Igualada
38FR016Nancy38ES547Sagunto
39FR014Amiens39ES026Coruna (A)
40FR048Annecy40ES007Murcia
41FR079Saint-Quentin41ES037Puerto de Santa Maria, El
42FR045Pau42ES516Salamanca
43FR006Strasbourg43ES021Alicante
44FR215Rouen44ES073Elda
45FR082Beauvais45ES048Guadalajara
46FR304Melun46ES046Gandia
47FR011Saint-Etienne47ES020Cordoba
48FR084Creil48ES023Gijon
49FR214Valence49ES533Marbella
50FR046Bayonne50ES532Algeciras
51FR026Grenoble51ES004Seville
52FR065Bourges52ES035Torrevieja
53FR060Chartres53ES018Logrono
54FR099Frejus54ES028Reus
55FR039Avignon55ES522Cadiz
56FR024Limoges56ES005Saragossa
57FR205Nice57ES006Malaga
58FR034Valenciennes58ES508Jerez de la Frontera
59FR008Nantes59ES065Linea de la Concepcion, La
60FR010Montpellier60ES012Vitoria
61FR040Mulhouse61ES505Elche/Elx
62FR047Annemasse62ES521Huelva
63FR007Bordeaux63ES001Madrid
64FR004Toulouse64ES070Irun
65FR043Perpignan65ES520Castellon de la Plana/Castello de la Plana
66FR044Nimes66ES015Santander
67FR052Saint-Nazaire67ES525Tarragona
68FR017Metz68ES002Barcelona
69FR009Lille69ES003Valencia
70FR003Lyon70ES019Bilbao
71FR012Le Havre71ES506Cartagena
72FR519Cannes72ES039Aviles
73FR207Lens - Lievin
74FR209Douai
75FR032Toulon
76FR062Calais
77FR001Paris
78FR203Marseille
79FR042Dunkerque
80FR324Martigues
81FR208Henin - Carvin
82FR013Rennes
Bottom Subfigures
ID NumberFUA CodeFUA NameID NumberFUA CodeFUA Name
FranceSpain
1FR324Martigues1ES013Oviedo
2FR047Annemasse2ES034Caceres
3FR039Avignon3ES012Vitoria
4FR040Mulhouse4ES014Pamplona
5FR048Annecy5ES021Alicante
6FR065Bourges6ES023Gijon
7FR082Beauvais7ES065Linea de la Concepcion, La
8FR208Henin - Carvin8ES041Palencia
9FR304Melun9ES547Sagunto
10FR505Charleville-Mezieres10ES059Zamora
11FR049Lorient11ES050Manresa
12FR067Quimper12ES035Torrevieja
13FR066Saint-Brieuc13ES037Puerto de Santa Maria, El
14FR214Valence14ES062Sanlucar de Barrameda
15FR209Douai15ES540Chiclana de la Frontera
16FR068Vannes16ES070Irun
17FR207Lens - Lievin17ES514Almeria
18FR053La Rochelle18ES046Gandia
19FR506Colmar19ES053Ciudad Real
20FR064Arras20ES057Ponferrada
21FR084Creil21ES528Lleida
22FR050Montbeliard22ES532Algeciras
23FR077Roanne23ES028Reus
24FR056Angouleme24ES039Aviles
25FR079Saint-Quentin25ES054Benidorm
26FR069Cherbourg26ES521Huelva
27FR012Le Havre27ES519Albacete
28FR086Evreux28ES537Alcoy
29FR519Cannes29ES527Jaen
30FR063Beziers30ES033Girona
31FR058Chambery31ES525Tarragona
32FR090Chateauroux32ES011Santiago de Compostela
33FR096Albi33ES040Talavera de la Reina
34FR052Saint-Nazaire34ES031Lugo
35FR057Boulogne-sur-Mer35ES505Elche/Elx
36FR059Chalon-sur-Saone36ES026Coruna (A)
37FR060Chartres37ES522Cadiz
38FR061Niort38ES016Toledo
39FR022Clermont-Ferrand39ES544Linares
40FR010Montpellier40ES018Logrono
41FR020Dijon41ES533Marbella
42FR019Orleans42ES510Donostia-San Sebastian
43FR026Grenoble43ES552Igualada
44FR025Besancon44ES073Elda
45FR076Belfort45ES019Bilbao
46FR023Caen46ES015Santander
47FR038Le Mans47ES003Valencia
48FR045Pau48ES529Ourense
49FR007Bordeaux49ES520Castellon de la Plana/Castello de la Plana
50FR021Poitiers50ES004Seville
51FR034Valenciennes51ES516Salamanca
52FR044Nimes52ES001Madrid
53FR035Tours53ES002Barcelona
54FR008Nantes54ES005Saragossa
55FR006Strasbourg55ES006Malaga
56FR004Toulouse56ES007Murcia
57FR042Dunkerque57ES009Valladolid
58FR017Metz58ES044Pontevedra
59FR003Lyon59ES017Badajoz
60FR104Chalons-en-Champagne60ES020Cordoba
61FR036Angers61ES022Vigo
62FR016Nancy62ES043Ferrol
63FR073Tarbes63ES048Guadalajara
64FR046Bayonne64ES506Cartagena
65FR093Brive-la-Gaillarde65ES501Granada
66FR099Frejus66ES508Jerez de la Frontera
67FR203Marseille67ES523Leon
68FR205Nice68ES542Basin
69FR001Paris69ES538Avila
70FR051Troyes70ES546Merida
71FR009Lille71ES545Lorca
72FR011Saint-Etienne72ES515Burgos
73FR013Rennes
74FR014Amiens
75FR018Reims
76FR024Limoges
77FR032Toulon
78FR037Brest
79FR062Calais
80FR043Perpignan
81FR215Rouen
82FR074Compiegne
Figure A1. Light emission distribution from the center of a monocentric FUA, modelled by different geometric shapes (left panel) and distributions of ALAN in corresponding FUAs (right panel).
Figure A1. Light emission distribution from the center of a monocentric FUA, modelled by different geometric shapes (left panel) and distributions of ALAN in corresponding FUAs (right panel).
Remotesensing 13 03714 g0a1aRemotesensing 13 03714 g0a1b
Figure A2. Jaccard Index for the estimated delineations, derived from the compactness-based (a,b) and model-based (c,d) ALAN thresholds: FUAs in France (a,c) and Spain (b,d). Note: The column numbering refers to FUA numbers listed in Table A2 below. In the graphs, FUAs are sorted in descending order according to their JI values.
Figure A2. Jaccard Index for the estimated delineations, derived from the compactness-based (a,b) and model-based (c,d) ALAN thresholds: FUAs in France (a,c) and Spain (b,d). Note: The column numbering refers to FUA numbers listed in Table A2 below. In the graphs, FUAs are sorted in descending order according to their JI values.
Remotesensing 13 03714 g0a2aRemotesensing 13 03714 g0a2b
Table A2. FUAs’ identification (ID number = number on the X-axis in Figure A2 above).
Table A2. FUAs’ identification (ID number = number on the X-axis in Figure A2 above).
ID NumberSubfigure (a)Subfigure (b)Subfigure (c)Subfigure (d)
FUA CodeFUA NameFUA CodeFUA NameFUA CodeFUA NameFUA CodeFUA Name
1FR006StrasbourgES065Linea de la Concepcion, LaFR006StrasbourgES065Linea de la Concepcion, La
2FR037BrestES015SantanderFR049LorientES501Granada
3FR047AnnemasseES001MadridFR047AnnemasseES001Madrid
4FR024LimogesES540Chiclana de la FronteraFR039AvignonES540Chiclana de la Frontera
5FR069CherbourgES002BarcelonaFR003LyonES515Burgos
6FR039AvignonES501GranadaFR203MarseilleES057Ponferrada
7FR042DunkerqueES520Castellon de la PlanaFR066Saint-BrieucES004Seville
8FR043PerpignanES506CartagenaFR043PerpignanES041Palencia
9FR001ParisES525TarragonaFR007BordeauxES538Avila
10FR023CaenES514AlmeriaFR025BesanconES009Valladolid
11FR062CalaisES018LogronoFR062CalaisES506Cartagena
12FR052Saint-NazaireES041PalenciaFR001ParisES516Salamanca
13FR067QuimperES516SalamancaFR068VannesES533Marbella
14FR022Clermont-FerrandES522CadizFR063BeziersES523Leon
15FR505Charleville-MezieresES009ValladolidFR008NantesES053Ciudad Real
16FR205NiceES054BenidormFR042DunkerqueES532Algeciras
17FR009LilleES022VigoFR052Saint-NazaireES022Vigo
18FR034ValenciennesES019BilbaoFR046BayonneES529Ourense
19FR073TarbesES533MarbellaFR093Brive-la-GaillardeES013Oviedo
20FR046BayonneES004SevilleFR061NiortES522Cadiz
21FR066Saint-BrieucES521HuelvaFR012Le HavreES014Pamplona
22FR049LorientES529OurenseFR053La RochelleES521Huelva
23FR050MontbeliardES013OviedoFR010MontpellierES039Aviles
24FR057Boulogne-sur-MerES026Coruna (A)FR214ValenceES034Caceres
25FR003LyonES039AvilesFR045PauES037Puerto de Santa Maria, El
26FR008NantesES003ValenciaFR205NiceES015Santander
27FR032ToulonES021AlicanteFR009LilleES546Merida
28FR506ColmarES037Puerto de Santa Maria, ElFR096AlbiES520Castellon de la Plana
29FR040MulhouseES053Ciudad RealFR505Charleville-MezieresES002Barcelona
30FR012Le HavreES014PamplonaFR024LimogesES542Basin
31FR519CannesES532AlgecirasFR506ColmarES003Valencia
32FR004ToulouseES062Sanlucar de BarramedaFR017MetzES062Sanlucar de Barrameda
33FR084CreilES528LleidaFR048AnnecyES019Bilbao
34FR065BourgesES057PonferradaFR022Clermont-FerrandES026Coruna (A)
35FR203MarseilleES033GironaFR057Boulogne-sur-MerES043Ferrol
36FR208Henin - CarvinES043FerrolFR215RouenES059Zamora
37FR036AngersES046GandiaFR004ToulouseES046Gandia
38FR068VannesES035TorreviejaFR036AngersES021Alicante
39FR010MontpellierES044PontevedraFR032ToulonES544Linares
40FR053La RochelleES011Santiago de CompostelaFR065BourgesES545Lorca
41FR074CompiegneES059ZamoraFR040MulhouseES044Pontevedra
42FR044NimesES510Donostia-San SebastianFR519CannesES017Badajoz
43FR017MetzES508Jerez de la FronteraFR016NancyES005Saragossa
44FR082BeauvaisES515BurgosFR073TarbesES020Cordoba
45FR007BordeauxES006MalagaFR019OrleansES035Torrevieja
46FR063BeziersES546MeridaFR034ValenciennesES510Donostia-San Sebastian
47FR076BelfortES050ManresaFR023CaenES537Alcoy
48FR038Le MansES523LeonFR026GrenobleES006Malaga
49FR093Brive-la-GaillardeES005SaragossaFR013RennesES054Benidorm
50FR214ValenceES527JaenFR044NimesES050Manresa
51FR016NancyES023GijonFR059Chalon-sur-SaoneES525Tarragona
52FR096AlbiES020CordobaFR021PoitiersES011Santiago de Compostela
53FR324MartiguesES070IrunFR076BelfortES527Jaen
54FR021PoitiersES544LinaresFR084CreilES007Murcia
55FR207Lens - LievinES007MurciaFR077RoanneES040Talavera de la Reina
56FR077RoanneES031LugoFR035ToursES519Albacete
57FR209DouaiES552IgualadaFR099FrejusES508Jerez de la Frontera
58FR045PauES048GuadalajaraFR082BeauvaisES023Gijon
59FR025BesanconES016ToledoFR067QuimperES514Almeria
60FR056AngoulemeES017BadajozFR209DouaiES528Lleida
61FR013RennesES073EldaFR324MartiguesES031Lugo
62FR026GrenobleES028ReusFR037BrestES028Reus
63FR048AnnecyES040Talavera de la ReinaFR014AmiensES016Toledo
64FR019OrleansES537AlcoyFR011Saint-EtienneES070Irun
65FR018ReimsES519AlbaceteFR074CompiegneES073Elda
66FR099FrejusES547SaguntoFR018ReimsES552Igualada
67FR304MelunES034CaceresFR104Chalons-en-ChampagneES033Girona
68FR064ArrasES538AvilaFR050MontbeliardES547Sagunto
69FR011Saint-EtienneES542BasinFR064ArrasES018Logrono
70FR060ChartresES505Elche/ElxFR069CherbourgES048Guadalajara
71FR059Chalon-sur-SaoneES545LorcaFR038Le MansES505Elche/Elx
72FR035ToursES012VitoriaFR079Saint-QuentinES012Vitoria
73FR086EvreuxFR207Lens - Lievin
74FR079Saint-QuentinFR058Chambery
75FR014AmiensFR090Chateauroux
76FR215RouenFR020Dijon
77FR058ChamberyFR051Troyes
78FR051TroyesFR056Angouleme
79FR061NiortFR060Chartres
80FR020DijonFR304Melun
81FR104Chalons-en-ChampagneFR086Evreux
82FR090ChateaurouxFR208Henin - Carvin
Box A1. Estimation of the compactness-based ALAN threshold (derivation).
The figure below illustrates our assumptions: We model actual FUAs’ shapes by ellipses with axes a and b (a > b).
Remotesensing 13 03714 i001
A. Optimal Radius Calculation
To calculate the radius of the circle (r*), ensuring maximal intersection with the ellipse we should define and maximize Jaccard index: J I = S C S E S C S E m a x , where Sc = area of the circle, and SE = area of the ellipse. For this sake, we should calculate and differentiate the following function:
J I r = 0 i y e l l i p s e x d x + i r y c i r c l e x d x 0 i y c i r c l e x d x + i a y e l l i p s e x d x
where y c i r c l e = r 2 x 2 is equation of circle and y ellipse = b a a 2 x 2 is equation of ellipse.
Limit of integration i is defined as x coordinate of intersection ycircle and yellipse:
r 2 x 2 = b a a 2 x 2 »   x = a r 2 b 2 a 2 b 2
Both integrals in Equation (A1) are of the same type, which are calculated in the same way:
k 2 x 2 d x = x = k s i n y d x = k c o s y d y = k 2 k 2 s i n 2 y * k c o s y d y = k 2 c o s 2 y d y = k 2 2 cos 2 y + 1 d x = k 2 2 sin 2 y 2 + y =   k 2 2 cos y sin y + y = k 2 2 x k 1 x k 2 + a s i n x k
Proceeding from the equations of circle and ellipse, limit of integration i (formula (A2)), and the integral calculation (formula (A3)), let’s consequentially compute the integrals in JI(r). Thus, the first integral in nominator will look like the following:
0 i y e l l i p s e x d x =   b a a 2 2 x a 1 x a 2 + a s i n x a a r 2 b 2 a 2 b 2 0 = =   a b 2 r 2 b 2 a 2 r 2 a 2 b 2 + a s i n r 2 b 2 a 2 b 2
The second one will be equal to
i r y c i r c l e x d x =   r 2 2 x r 1 x r 2 + a s i n x r r a r 2 b 2 a 2 b 2 = = π r 2 4 a b r 2 b 2 a 2 r 2 2 a 2 b 2 r 2 2 a s i n a r r 2 b 2 a 2 b 2
The nominator of JI(r) will equal to (A4) + (A5):
0 i y e l l i p s e x d x + i r y c i r c l e x d x =   a b 2 a s i n r 2 b 2 a 2 b 2 r 2 2 a s i n a r r 2 b 2 a 2 b 2 + π r 2 4
Actually, the denominator of JI(r), representing the union of SC and SE, is equal to the sum of the quarter of corresponding areas diminished by the intersection of SC and SE, calculated in formula (A6):
0 i y c i r c l e x d x + i a y e l l i p s e x d x = =   1 4 π r 2 + π a b a b 2 a s i n r 2 b 2 a 2 b 2 r 2 2 a s i n a r r 2 b 2 a 2 b 2 + π r 2 4 =   r 2 2 a s i n a r r 2 b 2 a 2 b 2 a b 2 a s i n r 2 b 2 a 2 b 2 + π a b 4
Thus, JI(r) is equal to
J I r = a b * a s i n r 2 b 2 a 2 b 2 r 2 * a s i n a r r 2 b 2 a 2 b 2 + π r 2 2 a b * a s i n r 2 b 2 a 2 b 2 + r 2 * a s i n a r r 2 b 2 a 2 b 2 + π a b 2 = = a b * a s i n r 2 b 2 a 2 b 2 r 2 * a s i n a r r 2 b 2 a 2 b 2 = z =   z + π r 2 2 z + π a b 2
Derivative of JI(r) will be equal to
d J I d r = z a b + r 2 2 z r + π r a b = = z =   a b * a s i n r 2 b 2 a 2 b 2 r 2 * a s i n a r r 2 b 2 a 2 b 2 = =   2 r * a s i n a r r 2 b 2 a 2 b 2 a b + r 2 2 r a b * a s i n r 2 b 2 a 2 b 2 r 2 * a s i n a r r 2 b 2 a 2 b 2 + π r a b = = 2 r a b a s i n a r r 2 b 2 a 2 b 2 + a s i n r 2 b 2 a 2 b 2 + π r a b
To find the maximum of the function JI(r), let’s put equal to zero its derivative and define r:
2 r a b a s i n a r r 2 b 2 a 2 b 2 + a s i n r 2 b 2 a 2 b 2 + π r a b = 0 ; a s i n a r r 2 b 2 a 2 b 2 + a s i n r 2 b 2 a 2 b 2 =   π 2 ; s i n a s i n a r r 2 b 2 a 2 b 2 + a s i n r 2 b 2 a 2 b 2 = =   s i n α + β = sin α 1 s i n 2 β + 1 s i n 2 α sin β = =   r 2 b 2 a 2 r 2 r a b = s i n π 2 = 1 r 2 b 2 * a 2 r 2 = r a b   r 2 b 2 * a 2 r 2 =   r 2 a b 2 r 2 a b 2 = 0 r 2 = a b r = a b
B. Optimal Percentile Calculation
Optimal percentile p* will equal to the share of area (2) (see figure above) of the area of ellipse:
p * = i a y e l l i p s e x d x i r y c i r c l e x d x π a b 4
Under defined optimal radius r = a b , limit of integration i is equal to:
a r 2 b 2 a 2 b 2 = a b a + b
Thus, proceeding from the equations of circle and ellipse, limit of integration i (formula (A12)), and the integral resolution (formula (A3)), p* will equal to:
p * = 1 π a b 4 ( b a a 2 2 x a 1 x a 2 + a s i n x a | b a + b a a b a a 2 2 x a b 1 x a b 2 + a s i n x a b b a + b a a b ) = 1 π a b 4 ( π a b 4 a b 2 a b a + b + a s i n b a + b π a b 4 a b 2 a b a + b + a s i n a a + b ) =   2 π a s i n a a + b + a s i n b a + b
Since s i n α + β = sin α 1 s i n 2 β + 1 s i n 2 α sin β ,
s i n a s i n a a + b + a s i n b a + b =   a b a + b , and then
p * = 2 π a s i n a b a + b
Finally, putting compactness c of the ellipse with axes a and b (a > b) to be equal to their ratio between the ellipse’s area and the area of the bonding circle ( c = A r e a   o f E l l i p s e A r e a o f B o n d i n g C i r c l e = π a b π a 2 = b a , optimal percentile p* will be equal to
p * = 2 π a r c s i n a b a + b = 2 π a r c s i n a a b a a a + b a = 2 π a r c s i n 1 c 1 + c

References

  1. 68% of the World Population Projected to Live in Urban Areas by 2050, Says UN|UN DESA|United Nations Department of Economic and Social Affairs. Available online: https://www.un.org/development/desa/en/news/population/2018-revision-of-world-urbanization-prospects.html (accessed on 8 March 2021).
  2. Urban Development Overview. Available online: https://www.worldbank.org/en/topic/urbandevelopment/overview (accessed on 8 March 2021).
  3. Dos Santos, S.; Adams, E.A.; Neville, G.; Wada, Y.; de Sherbinin, A.; Mullin Bernhardt, E.; Adamo, S.B. Urban growth and water access in sub-Saharan Africa: Progress, challenges, and emerging research directions. Sci. Total Environ. 2017, 607–608, 497–508. [Google Scholar] [CrossRef] [PubMed]
  4. Li, C.; Li, J.; Wu, J. What drives urban growth in China? A multi-scale comparative analysis. Appl. Geogr. 2018, 98, 43–51. [Google Scholar] [CrossRef]
  5. Wolff, M.; Wiechmann, T. Urban growth and decline: Europe’s shrinking cities in a comparative perspective 1990–2010. Eur. Urban Reg. Stud. 2018, 25, 122–139. [Google Scholar] [CrossRef]
  6. Martellozzo, F.; Amato, F.; Murgante, B.; Clarke, K.C. Modelling the impact of urban growth on agriculture and natural land in Italy to 2030. Appl. Geogr. 2018, 91, 156–167. [Google Scholar] [CrossRef] [Green Version]
  7. McDonald, R.I.; Mansur, A.V.; Ascensão, F.; Colbert, M.; Crossman, K.; Elmqvist, T.; Gonzalez, A.; Güneralp, B.; Haase, D.; Hamann, M.; et al. Research gaps in knowledge of the impact of urban growth on biodiversity. Nat. Sustain. 2020, 3, 16–24. [Google Scholar] [CrossRef]
  8. Zullo, F.; Fazio, G.; Romano, B.; Marucci, A.; Fiorini, L. Effects of urban growth spatial pattern (UGSP) on the land surface temperature (LST): A study in the Po Valley (Italy). Sci. Total Environ. 2019, 650, 1740–1751. [Google Scholar] [CrossRef]
  9. Van de Voorde, T.; Jacquet, W.; Canters, F. Mapping form and function in urban areas: An approach based on urban metrics and continuous impervious surface data. Landsc. Urban Plan. 2011, 102, 143–155. [Google Scholar] [CrossRef]
  10. Yoshida, H.; Omae, M. An approach for analysis of urban morphology: Methods to derive morphological properties of city blocks by using an urban landscape model and their interpretations. Comput. Environ. Urban Syst. 2005, 29, 223–247. [Google Scholar] [CrossRef]
  11. Ng, E. Policies and technical guidelines for urban planning of high-density cities - air ventilation assessment (AVA) of Hong Kong. Build. Environ. 2009, 44, 1478–1488. [Google Scholar] [CrossRef]
  12. Ng, E.; Yuan, C.; Chen, L.; Ren, C.; Fung, J.C.H. Improving the wind environment in high-density cities by understanding urban morphology and surface roughness: A study in Hong Kong. Landsc. Urban Plan. 2011, 101, 59–74. [Google Scholar] [CrossRef] [PubMed]
  13. Yuan, C.; Ng, E.; Norford, L.K. Improving air quality in high-density cities by understanding the relationship between air pollutant dispersion and urban morphologies. Build. Environ. 2014, 71, 245–258. [Google Scholar] [CrossRef] [PubMed]
  14. Vanderhaegen, S.; Canters, F. Mapping urban form and function at city block level using spatial metrics. Landsc. Urban Plan. 2017, 167, 399–409. [Google Scholar] [CrossRef]
  15. Glossary: Functional Urban Area—Statistics Explained. Available online: https://ec.europa.eu/eurostat/statistics-explained/index.php/Glossary:Functional_urban_area (accessed on 9 March 2021).
  16. Portnov, B.A.; Schwartz, M. URBAN CLUSTERS AS GROWTH FOCI*. J. Reg. Sci. 2009, 49, 287–310. [Google Scholar] [CrossRef]
  17. Dijkstra, L.; Poelman, H.; Veneri, P. The EU-OECD definition of a functional urban area. OECD Reg. Dev. Work. Pap. 2019, 11. [Google Scholar] [CrossRef]
  18. Bosker, M.; Park, J.; Roberts, M. Definition Matters: Metropolitan Areas and Agglomeration Economies in a Large Developing Country. J. Urban Econ. 2020, 7, 103275. [Google Scholar] [CrossRef]
  19. Zhou, Y.; Smith, S.J.; Elvidge, C.D.; Zhao, K.; Thomson, A.; Imhoff, M. A cluster-based method to map urban area from DMSP/OLS nightlights. Remote Sens. Environ. 2014, 147, 173–185. [Google Scholar] [CrossRef]
  20. Sutton, P.; Roberts, D.; Elvidge, C.; Baugh, K. Census from Heaven: An estimate of the global human population using night-time satellite imagery. Int. J. Remote Sens. 2001, 22, 3061–3076. [Google Scholar] [CrossRef]
  21. Henderson, M.; Yeh, E.T.; Gong, P.; Elvidge, C.; Baugh, K. Validation of urban boundaries derived from global night-time satellite imagery. Int. J. Remote Sens. 2003, 24, 595–609. [Google Scholar] [CrossRef]
  22. Kloog, I.; Haim, A.; Stevens, R.G.; Portnov, B.A. Global co-distribution of light at night (LAN) and cancers of prostate, colon, and lung in men. Chronobiol. Int. 2009, 26, 108–125. [Google Scholar] [CrossRef]
  23. Kloog, I.; Stevens, R.G.; Haim, A.; Portnov, B.A. Nighttime light level co-distributes with breast cancer incidence worldwide. Cancer Causes Control 2010, 21, 2059–2068. [Google Scholar] [CrossRef]
  24. Haim, A.; Portnov, B.A. Light Pollution as a New Risk Factor for Human Breast and Prostate Cancers; Springer: Dordrecht, The Netherlands, 2013; ISBN 9789400762206. [Google Scholar]
  25. Rybnikova, N.A.; Portnov, B.A. GDP per capita and obesity prevalence worldwide: An ambiguity of effects modification. Int. J. Obes. 2017, 41, 352. [Google Scholar] [CrossRef]
  26. Garcia-Saenz, A.; Sánchez de Miguel, A.; Espinosa, A.; Valentin, A.; Aragonés, N.; Llorca, J.; Amiano, P.; Martín Sánchez, V.; Guevara, M.; Capelo, R.; et al. Evaluating the Association between Artificial Light-at-Night Exposure and Breast and Prostate Cancer Risk in Spain (MCC-Spain Study). Environ. Health Perspect. 2018, 126, 047011. [Google Scholar] [CrossRef]
  27. Garcia-Saenz, A.; De Miguel, A.S.; Espinosa, A.; Costas, L.; Aragonés, N.; Tonne, C.; Moreno, V.; Pérez-Gómez, B.; Valentin, A.; Pollán, M.; et al. Association between outdoor light-at-night exposure and colorectal cancer in Spain. Epidemiology 2020, 31, 718–727. [Google Scholar] [CrossRef]
  28. Doll, C.H.; Muller, J.-P.; Elvidge, C.D. Night-time Imagery as a Tool for Global Mapping of Socioeconomic Parameters and Greenhouse Gas Emissions. AMBIO J. Hum. Environ. 2000, 29, 157–162. [Google Scholar] [CrossRef]
  29. Ebener, S.; Murray, C.; Tandon, A.; Elvidge, C.C. From wealth to health: Modelling the distribution of income per capita at the sub-national level using night-time light imagery. Int. J. Health Geogr. 2005, 4, 5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Ghosh, T.; L Powell, R.; D Elvidge, C.; E Baugh, K.; C Sutton, P.; Anderson, S. Shedding light on the global distribution of economic activity. Open Geogr. J. 2010, 3, 147–160. [Google Scholar]
  31. Wu, R.; Yang, D.; Dong, J.; Zhang, L.; Xia, F. Regional Inequality in China Based on NPP-VIIRS Night-Time Light Imagery. Remote Sens. 2018, 10, 240. [Google Scholar] [CrossRef] [Green Version]
  32. Amaral, S.; Monteiro, A.M.; Camara, G.; Quintanilha, J.A. DMSP/OLS night-time light imagery for urban population estimates in the Brazilian Amazon. Int. J. Remote Sens. 2006, 27, 855–870. [Google Scholar] [CrossRef]
  33. Zhuo, L.; Ichinose, T.; Zheng, J.; Chen, J.; Shi, P.J.; Li, X. Modelling the population density of China at the pixel level based on DMSP/OLS non-radiance-calibrated night-time light images. Int. J. Remote Sens. 2009, 30, 1003–1018. [Google Scholar] [CrossRef]
  34. Anderson, S.J.; Tuttle, B.T.; Powell, R.L.; Sutton, P.C. Characterizing relationships between population density and nighttime imagery for Denver, Colorado: Issues of scale and representation. Int. J. Remote Sens. 2010, 31, 5733–5746. [Google Scholar] [CrossRef]
  35. Mellander, C.; Lobo, J.; Stolarick, K.; Matheson, Z. Night-time light data: A good proxy measure for economic activity? PLoS ONE 2015, 10, e0139779. [Google Scholar] [CrossRef] [Green Version]
  36. Hopkins, G.R.; Gaston, K.J.; Visser, M.E.; Elgar, M.A.; Jones, T.M. Artificial light at night as a driver of evolution across urban-rural landscapes. Front. Ecol. Environ. 2018, 16, 472–479. [Google Scholar] [CrossRef] [Green Version]
  37. Hölker, F.; Wolter, C.; Perkin, E.K.; Tockner, K. Light pollution as a biodiversity threat. Trends Ecol. Evol. 2010, 25, 681–682. [Google Scholar] [CrossRef] [PubMed]
  38. Narisada, K.; Schreuder, D. Light Pollution Handbook; Astrophysics and Space Science Library; Springer: Dordrecht, The Netherlands, 2004. [Google Scholar]
  39. Bennie, J.; Duffy, J.; Davies, T.; Correa-Cano, M.; Gaston, K. Global Trends in Exposure to Light Pollution in Natural Terrestrial Ecosystems. Remote Sens. 2015, 7, 2715–2730. [Google Scholar] [CrossRef] [Green Version]
  40. Falchi, F.; Furgoni, R.; Gallaway, T.A.A.; Rybnikova, N.A.A.; Portnov, B.A.A.; Baugh, K.; Cinzano, P.; Elvidge, C.D.D. Light pollution in USA and Europe: The good, the bad and the ugly. J. Environ. Manag. 2019, 248, 109227. [Google Scholar] [CrossRef]
  41. Imhoff, M.L.; Lawrence, W.T.; Stutzer, D.C.; Elvidge, C.D. A technique for using composite DMSP/OLS “city lights” satellite data to map urban area. Remote Sens. Environ. 1997, 61, 361–370. [Google Scholar] [CrossRef]
  42. Shi, K.; Huang, C.; Yu, B.; Yin, B.; Huang, Y.; Wu, J. Evaluation of NPP-VIIRS night-time light composite data for extracting built-up urban areas. Remote Sens. Lett. 2014, 5, 358–366. [Google Scholar] [CrossRef]
  43. He, C.; Shi, P.; Li, J.; Chen, J.; Pan, Y.; Li, J.; Zhuo, L.; Toshiaki, I. Restoring urbanization process in China in the 1990s by using non-radiance-calibrated DMSP/OLS nighttime light imagery and statistical data. Chin. Sci. Bull. 2006, 51, 1614–1620. [Google Scholar] [CrossRef]
  44. Liu, Z.; He, C.; Zhang, Q.; Huang, Q.; Yang, Y. Extracting the dynamics of urban expansion in China using DMSP-OLS nighttime light data from 1992 to 2008. Landsc. Urban Plan. 2012, 106, 62–72. [Google Scholar] [CrossRef]
  45. Yu, B.; Shu, S.; Liu, H.; Song, W.; Wu, J.; Wang, L.; Chen, Z. Object-based spatial cluster analysis of urban landscape pattern using nighttime light satellite images: A case study of China. Int. J. Geogr. Inf. Sci. 2014, 28, 2328–2355. [Google Scholar] [CrossRef]
  46. Elvidge, C.D.; Baugh, K.E.; Zhizhin, M.; Hsu, F.-C. Why VIIRS data are superior to DMSP for mapping nighttime lights. Proc. Asia-Pac. Adv. Netw. 2013, 35, 62. [Google Scholar] [CrossRef] [Green Version]
  47. Dou, Y.; Liu, Z.; He, C.; Yue, H. Urban Land Extraction Using VIIRS Nighttime Light Data: An Evaluation of Three Popular Methods. Remote Sens. 2017, 9, 175. [Google Scholar] [CrossRef] [Green Version]
  48. Earth Observation Goup. Available online: https://eogdata.mines.edu/products/vnl/ (accessed on 8 August 2021).
  49. Elvidge, C.D.; Baugh, K.; Zhizhin, M.; Hsu, F.C.; Ghosh, T. VIIRS night-time lights. Int. J. Remote Sens. 2017, 38, 5860–5879. [Google Scholar] [CrossRef]
  50. Román, M.O.; Wang, Z.; Sun, Q.; Kalb, V.; Miller, S.D.; Molthan, A.; Schultz, L.; Bell, J.; Stokes, E.C.; Pandey, B.; et al. NASA’s Black Marble nighttime lights product suite. Remote Sens. Environ. 2018, 210, 113–143. [Google Scholar] [CrossRef]
  51. Functional Urban Areas by Country—OECD. Available online: https://www.oecd.org/cfe/regional-policy/functionalurbanareasbycountry.htm (accessed on 3 June 2020).
  52. LandScan Datasets|LandScan™. Available online: https://landscan.ornl.gov/landscan-datasets (accessed on 17 March 2020).
  53. Measuring Compactness. Available online: https://fisherzachary.github.io/public/r-output.html (accessed on 24 January 2021).
  54. Li, W.; Goodchild, M.F.; Church, R.L. An Efficient Measure of Compactness for 2D Shapes and its Application in Regionalization Problems. Int. J. Geogr. Inf. Sci. 2013, 27, 1227–1250. [Google Scholar] [CrossRef] [Green Version]
  55. Brody, S. The characteristics, causes, and consequences of sprawling development patterns in the United States. Nat. Educ. Knowl. 2013, 4, 2. [Google Scholar]
  56. SPSS Library: MANOVA and GLM. Available online: https://stats.idre.ucla.edu/spss/library/spss-librarymanova-and-glm-2/ (accessed on 4 April 2021).
  57. Breiman, L. Random forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  58. Create Bag of Decision Trees—MATLAB. Available online: https://www.mathworks.com/help/stats/treebagger.html (accessed on 5 January 2021).
  59. Bernard, S.; Adam, S.; Heutte, L. Dynamic Random Forests. Pattern Recognit. Lett. 2012, 33, 1580–1586. [Google Scholar] [CrossRef] [Green Version]
  60. Chung, N.C.; Miasojedow, B.; Startek, M.; Gambin, A. Jaccard/Tanimoto similarity test and estimation methods. BMC Bioinform. 2019, 20, 644. [Google Scholar] [CrossRef]
  61. Global Human Settlement—GHS POPULATION GRID—European Commission. Available online: https://ghsl.jrc.ec.europa.eu/ghs_pop.php (accessed on 20 June 2021).
Figure 1. Commuting-based boundaries (black lines) of the Paris (a) and Chateauroux (b) FUAs vs. the ALAN contours (blue lines), representing the 0.71 nW/cm2/sr threshold level.
Figure 1. Commuting-based boundaries (black lines) of the Paris (a) and Chateauroux (b) FUAs vs. the ALAN contours (blue lines), representing the 0.71 nW/cm2/sr threshold level.
Remotesensing 13 03714 g001aRemotesensing 13 03714 g001b
Figure 2. Flowchart of study stages.
Figure 2. Flowchart of study stages.
Remotesensing 13 03714 g002
Figure 3. ALAN maps for continental France (a) and Spain (b). Note: Areas located outside the national borders are marked in blue.
Figure 3. ALAN maps for continental France (a) and Spain (b). Note: Areas located outside the national borders are marked in blue.
Remotesensing 13 03714 g003
Figure 4. FUAs and their cores in continental France (a) and Spain (b).
Figure 4. FUAs and their cores in continental France (a) and Spain (b).
Remotesensing 13 03714 g004
Figure 5. A simplified distribution of ALAN emissions (a) and the associated frequency distribution of ALAN values (b).
Figure 5. A simplified distribution of ALAN emissions (a) and the associated frequency distribution of ALAN values (b).
Remotesensing 13 03714 g005
Figure 6. Examples of compact monocentric FUAs, which territorial footprints are close to a circular shape: Le Mans (a) and Limoges (b) in France. Note: Thin grey lines mark FUAs’ boundaries.
Figure 6. Examples of compact monocentric FUAs, which territorial footprints are close to a circular shape: Le Mans (a) and Limoges (b) in France. Note: Thin grey lines mark FUAs’ boundaries.
Remotesensing 13 03714 g006
Figure 7. Relationship between a FUA’s compactness (c) and the optimal ALAN percentile (p*). Note: Shapes deviating from a perfect circle are assumed to be elliptical; see text for explanations.
Figure 7. Relationship between a FUA’s compactness (c) and the optimal ALAN percentile (p*). Note: Shapes deviating from a perfect circle are assumed to be elliptical; see text for explanations.
Remotesensing 13 03714 g007
Figure 8. FUAs in Austria used for the models’ validation.
Figure 8. FUAs in Austria used for the models’ validation.
Remotesensing 13 03714 g008
Figure 9. The modal and dimmest ALAN values estimated for individual FUAs before correcting for compactness ((a) = France; (b) = Spain) and after correcting for compactness ((c) = France; (d) = Spain). Notes: The column numbering (axis X) refers to FUA numbers listed in Table A1 of the Appendix. FUAs are sorted in an ascending order according to their modal ALAN values (upper diagrams) or according to compactness-based ALAN thresholds (bottom diagrams).
Figure 9. The modal and dimmest ALAN values estimated for individual FUAs before correcting for compactness ((a) = France; (b) = Spain) and after correcting for compactness ((c) = France; (d) = Spain). Notes: The column numbering (axis X) refers to FUA numbers listed in Table A1 of the Appendix. FUAs are sorted in an ascending order according to their modal ALAN values (upper diagrams) or according to compactness-based ALAN thresholds (bottom diagrams).
Remotesensing 13 03714 g009
Figure 10. Models cross-validation results for France (a) and Spain (b).
Figure 10. Models cross-validation results for France (a) and Spain (b).
Remotesensing 13 03714 g010
Figure 11. Examples of FUAs featuring compactness-based boundaries (blue lines), model-based boundaries (green lines) and commuting-based boundaries (black lines): Paris (a) and Madrid (b) (see text for explanations).
Figure 11. Examples of FUAs featuring compactness-based boundaries (blue lines), model-based boundaries (green lines) and commuting-based boundaries (black lines): Paris (a) and Madrid (b) (see text for explanations).
Remotesensing 13 03714 g011
Figure 12. Commuting-based (a) vs. model-estimated (b) delineations of FUAs in France and Spain.
Figure 12. Commuting-based (a) vs. model-estimated (b) delineations of FUAs in France and Spain.
Remotesensing 13 03714 g012
Figure 13. Commuting-based vs. models-estimated delineations of FUAs in Austria (see text for explanations).
Figure 13. Commuting-based vs. models-estimated delineations of FUAs in Austria (see text for explanations).
Remotesensing 13 03714 g013
Table 1. Descriptive statistics of the research variables.
Table 1. Descriptive statistics of the research variables.
VariableMinimumMaximumMeanSD
France (82 FUAs)
Latitude of the FUA core’s centroid (dd)42.75751.00147.1372.366
Population density of the FUA core (persons per km2)89.5292586.130470.764382.596
Population density decline gradient a1.11015.4093.9862.224
Distance to the nearest major city (dd) a0.0005.9001.7811.359
Average ALAN level (nW/cm2/sr)0.64920.4113.6403.614
Spain (72 FUAs)
Latitude of the FUA core centroid (dd)36.11043.56140.0932.381
Population density of the FUA core (persons per km2)20.5543485.360831.839794.144
Population density decline gradient b0.92616.2324.8633.724
Distance to the nearest major city (dd)0.0005.1202.2021.542
ALAN averaged level (nW/cm2/sr)0.74523.1296.4424.583
Notes: a Calculated as straight line distance between a FUA core’s centroid and centroid of the closest FUA with 1.5M+ residents; b Calculated as the ratio between the population density of the FUA core and that of the core’s buffer with a 5 km width for small FUAs (less than 100,000 residents), a 15 km buffer for medium-size FUAs (100,000–250,000 residents), and a 25 km buffer for large FUAs (over 250,000 residents).
Table 2. Descriptive statistics of the identified ALAN thresholds.
Table 2. Descriptive statistics of the identified ALAN thresholds.
Country/VariableMinimumMaximumMeanSD
France (Number of FUAs = 82)
• ALAN percentile (0–100)13.61045.29025.6236.533
• ALAN threshold (nW/cm2/sr)0.1509.9100.6641.218
Spain (Number of FUAs = 72)
• ALAN percentile (0–100)12.66046.46029.3557.279
• ALAN threshold (nW/cm2/sr)0.1308.2301.0261.518
Table 3. Factors affecting ALAN threshold values estimated for individual FUAs (Method—OLS; Dependent variable—ALAN optimal threshold level, Box-Cox transformed with α = −0.55).
Table 3. Factors affecting ALAN threshold values estimated for individual FUAs (Method—OLS; Dependent variable—ALAN optimal threshold level, Box-Cox transformed with α = −0.55).
PredictorModel 1 (France)Model 2 (Spain)
B aBeta bt cB aBeta bt c
(Constant)−11.192-−8.466 *−11.846-−8.090 *
Latitude (dd)0.1060.2464.209 *0.1400.2783.909 *
Population density of the FUA core, persons per km2 (ln)1.1510.75111.143 *1.1581.11413.124 *
Population density gradient (ln)−1.369−0.685−10.349 *−1.209−0.724−7.977 *
Distance to the nearest major city (dd)−0.137−0.183−3.035 *−0.147−0.190−2.945 *
N of obs.8272
R20.7390.740
r0.8660.812
SEE0.5330.629
WMSE4.5212.718
F54.43 *46.304 *
Notes: a unstandardized regression coefficients; b standardized regression coefficients; c t-statistic and its significance level; SEE = standard error of the estimates; WMSE = weighted mean squared error; F = F-statistics; * 0.01 significance level.
Table 4. Values of the Jaccard index (JI) calculated for FUAs of different size and population density.
Table 4. Values of the Jaccard index (JI) calculated for FUAs of different size and population density.
FUA TypeN. of Obs.Delineations Derived from Compactness-Based ALAN ThresholdsDelineations Derived from Model-Based ALAN Thresholds
MeanSDMeanSD
All FUAs under analysis1540.3420.1580.3510.150
FUAs in:
• France820.3040.1160.3260.124
• Spain720.3850.1860.3780.171
FUAs by class:
• 1&2 (Smallest)930.3270.1540.3350.152
• 3 (Medium)550.3510.1540.3610.138
• 4 (Largest)60.4990.1870.5070.134
Population density in the FUA core, people per km2 (ln)
• ≤5180.2350.0970.3490.139
• >51360.3560.1590.3510.151
• >6.5480.4400.1770.4240.174
• >7.5120.5510.1900.4870.204
Population density in the core’s buffer zone, people per km2 (ln)
• ≤51010.3130.1370.3400.135
• >5530.3980.1800.3720.174
• >6110.5250.2200.4690.235
• >6.340.6380.2770.5570.299
Ratio between population density in the core and the core’s buffer zone
• ≤2250.3030.1510.3190.157
• >21290.3500.1590.3570.148
• >5430.3630.1530.3570.155
• >10110.4670.1840.4580.197
Table 5. Individually fitted vs. model-estimated ALAN thresholds for FUAs in Austria.
Table 5. Individually fitted vs. model-estimated ALAN thresholds for FUAs in Austria.
FUA ALAN Threshold, nW/cm2/sr
Individually FittedEstimated Using the “French” ModelEstimated Using the “Spanish” Model
Vienna0.340.651.56
Graz0.230.481.05
Linz0.240.491.04
Salzburg0.270.491.08
Innsbruck0.260.300.78
Klagenfurt0.170.300.72
Performance indicators (in relation to individually fitted ALAN threshold)
r-0.7830.771
SEE-0.2160.819
WMSE-0.71110.102
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rybnikova, N.; Portnov, B.A.; Charney, I.; Rybnikov, S. Delineating Functional Urban Areas Using a Multi-Step Analysis of Artificial Light-at-Night Data. Remote Sens. 2021, 13, 3714. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13183714

AMA Style

Rybnikova N, Portnov BA, Charney I, Rybnikov S. Delineating Functional Urban Areas Using a Multi-Step Analysis of Artificial Light-at-Night Data. Remote Sensing. 2021; 13(18):3714. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13183714

Chicago/Turabian Style

Rybnikova, Nataliya, Boris A. Portnov, Igal Charney, and Sviatoslav Rybnikov. 2021. "Delineating Functional Urban Areas Using a Multi-Step Analysis of Artificial Light-at-Night Data" Remote Sensing 13, no. 18: 3714. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13183714

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop