Understanding the Role of Optimized Land Use / Land Cover Components in Mitigating Summertime Intra-Surface Urban Heat Island E ﬀ ect: A Study on Downtown Shanghai, China

: In this study, 167 land parcels of downtown Shanghai, China, were used to investigate the relationship between parcel-level land use / land cover (LULC) components and associated summertime intra-surface urban heat island (SUHI) e ﬀ ect, and further analyze the potential of mitigating summertime intra-SUHI e ﬀ ect through the optimized LULC components, by integrating a thermal sharpening method combining the Landsat-8 thermal band 10 data and high-resolution Quickbird image, statistical analysis, and nonlinear programming with constraints. The results show the remarkable variations in intra-surface urban heat island (SUHI) e ﬀ ect, which was measured with the mean parcel-level blackbody sensible heat ﬂux in kW per ha (Mean_pc_BBF). Through measuring the relative importance of each speciﬁc predictor in terms of their contributions to changing Mean_pc_BBF, the inﬂuence of parcel-level LULC components on excess surface ﬂux of heat energy to the atmosphere was estimated using the partial least square regression (PLSR) model. Analysis of the present and optimized parcel-level LULC components and their contribution to the associated Mean_pc_BBF were comparable between land parcels with varying sizes. Furthermore, focusing on the gap between the present and ideally optimized area proportions of parcel-level LULC components towards minimizing the Mean_pc_BBF, the uncertainties arising from the datasets and methods, as well as the implications for sustainable land development and mitigating the UHI e ﬀ ect were discussed.

Since the 1830s, the meteorological station recorded air temperature (AT) has been widely used to monitor the UHI effect at canopy height, due to the advantage of long time series data for historical

Location of the Study Area
Downtown Shanghai, the economic center of China, is located between latitudes 31 • 32 N-31 • 27 N and longitudes 120 • 52 E-121 • 45 E. This region has a northern subtropical monsoon climate. The average annual temperature approximates 15 • C, with temperatures averaging 28 • C in the summer and 4 • C in the winter. The average annual precipitation ranges between 1000 and 1200 mm, with approximately 60% of the rainfall being received during spring and autumn. It covers an area of 4000 km 2 , housing approximately 84.03% of a total of 24.18 million permanent residents within the boundary of greater Shanghai [41]. In this study, concerning previous studies of urban function zones (UFZs) in China [42][43][44][45], four typical UFZs of downtown Shanghai were selected (see Figure 1 and Table 1).
Energies 2020, 13, x FOR PEER REVIEW 3 of 18 and 4 °C in the winter. The average annual precipitation ranges between 1000 and 1200 mm, with approximately 60% of the rainfall being received during spring and autumn. It covers an area of 4000 km 2 , housing approximately 84.03% of a total of 24.18 million permanent residents within the boundary of greater Shanghai [41]. In this study, concerning previous studies of urban function zones (UFZs) in China [42][43][44][45], four typical UFZs of downtown Shanghai were selected (see Figure 1 and Table 1).  The heart of downtown Shanghai, with a cluster of municipal administrative services, banking, headquarter economy, commercial center, and historical and cultural resorts.

Xujiahui 2.58
One of the sub-centers of downtown Shanghai, featuring a cluster of a commercial center, historical and cultural resorts, high education, advanced medical care, and innovative enterprises.

Materials
The data sources include the Landsat-8 images, high-resolution Quickbird image, and the auxiliary datasets. Due to the influence of cloud contamination, the available satellite images were limited. In this study, two cloud-free Landsat-8 satellite images (path/row:118/38) acquired on two summertime dates (August 13, 2013, and August 3, 2015), were used for retrieval of LST. A high-   [46] were used as the auxiliary datasets.

Methods
To better illustrate our research aims and related arrangement of the sections, an overall technical flowchart containing the basic and key procedures employed for this study were shown in Figure 2.

Methods
To better illustrate our research aims and related arrangement of the sections, an overall technical flowchart containing the basic and key procedures employed for this study were shown in Figure 2.

Classification of LULC Components and Delimitation of Land Parcels
By using a land surface classification system with nine LULC types in downtown Shanghai [31], the classification of LULC components was performed with an object-oriented classification (OOC) method. This OOC method could make full use of contextual data-such as spectral information, texture, spatial neighborhood properties, and fractal dimensions-from very high-resolution images and delineate the objects of interest [47] with much higher accuracy than traditional per-pixel classification methods [48][49][50]. The overall accuracy is approximately 80.03%. The post-classification product was further manually corrected by overlapping it with Google Earth and Baidu Map layers and then validated in a field survey, with an overall correction accuracy of 91.1%. Subsequently, based on our prior knowledge of land use and the developmental intensity of the UFZs, 167 regular polygons were drawn to delimit the land parcels along the directional and geometric features of traffic roads, bounding walls, and urban creeks enclosing the clustered land units.
Herein, we want to emphasize that the classification of LULC components was used for the estimation of land surface emissivity. Given the mixing-pixel effect of the land surface showing the same/similar LST, and many LULC types may complicate the analysis between LULC components and the associated thermal effect, for this study, the nine LULC classes were aggregated into four

Classification of LULC Components and Delimitation of Land Parcels
By using a land surface classification system with nine LULC types in downtown Shanghai [31], the classification of LULC components was performed with an object-oriented classification (OOC) method. This OOC method could make full use of contextual data-such as spectral information, texture, spatial neighborhood properties, and fractal dimensions-from very high-resolution images and delineate the objects of interest [47] with much higher accuracy than traditional per-pixel classification methods [48][49][50]. The overall accuracy is approximately 80.03%. The post-classification product was further manually corrected by overlapping it with Google Earth and Baidu Map layers and then validated in a field survey, with an overall correction accuracy of 91.1%. Subsequently, based on our prior knowledge of land use and the developmental intensity of the UFZs, 167 regular polygons were drawn to delimit the land parcels along the directional and geometric features of traffic roads, bounding walls, and urban creeks enclosing the clustered land units.
Herein, we want to emphasize that the classification of LULC components was used for the estimation of land surface emissivity. Given the mixing-pixel effect of the land surface showing the

Retrieval and Validation of High-Resolution Thermally Sharpened LST
In this study, retrieval and validation of high-resolution thermally sharpened LST include the following steps. Firstly, for Landsat 8 TIR bands, given band 11 is subject to higher uncertainty of telescope stray light disturbance [51,52], band 10 was used to generate the top of the atmosphere (TOA) radiance as follows, where L sensor,λ is the at-sensor radiance of thermal band pixels in W/(m 2 ster µm), the gain is the slope of the radiance/DN conversion function, and offset the slope of the radiance/DN conversion, respectively [53]. Secondly, a spatial interpolation-based method known as co-Kriging interpolation was employed to generate the high-resolution thermally sharpened LST, via a combination of the high-resolution land surface products and Landsat-8 TIR band 10 data [31]. Given the bias arising from the different resolutions between the high-resolution LULC components product and coarse TOA radiance data, both datasets must be scaled to the same resolution. Therefore, the high-resolution LULC components product were resampled with multiple resolutions (1-9 m) and set as the base maps to overlap and delimit the TOA radiance layer. We assumed that same/similar surfaces would have the same/similar radiance values since there is no available reference data for this study. Then, for each scene of the TOA radiance layer, we tried many times from hundreds of points per km 2 to 6000 points per km 2 to compare the co-Kriging interpolation results, and we found the remarkable decrease in the pairwise bias between the raw and simulated points in response to the increasing point number. It was found that the threshold of 3000 points was reasonable as the pairwise bias was inclined to get flat and changed little when the point number was greater than 3000 points [31]. Therefore, 3000 spatially random points per km 2 with the same interval as the multiple resolutions of the LULC component products were generated to extract the point-based radiance values falling within the polygon of each same/similar LULC category. The extracted points with radiance values were used to rebuild the TOA radiance layer, using the co-Kriging interpolation method. Besides, based on the multiple-resolution LULC layers, a surface emissivity correction for the same/similar LULC components was performed according to empirical studies [54,55] and laboratory testing [56]. Subsequently, based on a generalized single-channel method known as the range transfer equation (RTE) for retrieving LST [57], the emissivity corrected LST using was computed as follows, where B λ (T S ) is the black body radiance given by Plank's law, T S the black body LST in Kelvin (K), ε λ the corrected emissivity of the specific land surface (details see reference [31]). L ↓ atm, λ is the downwelling atmospheric radiance, L ↑ atm, λ the upwelling atmospheric radiance, | \textpipe \textpipevar \textvertline λ the total atmospheric transmissivity between the surface and sensor, all of which were retrieved from a web-based Atmospheric Correction Tool [58]. k 2 and k 1 are band 10 thermal conversion constants included in the metadata file of Landsat 8 data [53].
Finally, validation for the sharpened LST products was performed by comparing the pixel-based root-mean-square error (RMSE) between the original 30 m LST products and the resampled 30 m products from the sharpened products (Appendix A Table A1). The sharpened LST products at 1 m resolution were employed for further analysis due to their superior visual quality and minimal RMSE.

Estimation of Parcel-Level Sensible Heat Flux
Based on the high-resolution sharpened LST products, the parcel-level blackbody sensible heat flux (BBF) was estimated follows [32]: where Φ BBFD is the pixel-based BBFD (W·m −2 ), C 1 = 3.7404×10 8 (W·µ 4 ·m −2 ), C 2 = 14387, λ is TOA radiance, and T is LST in Kelvin (K), A i is area of the ith land parcel (m 2 ). However, due to the uneven size of the land parcels, there exists a remarkable variation of parcel-level BBF. Therefore, regardless of the size effect of land parcels, the parcel-level BBFs were converted to per ha BBF (pc_BBF), namely a cumulative number of the BBFD multiplying by the pixels within an idea parcel sized 1 ha, to make them comparable.

Statistical Analysis and Estimating the Optimized Lulc Components for Minimizing Mean_pc_BBF
As a regular step, the exploring data analysis procedure, including descriptive statistics, normality test, Box-Cox transformation for skewed data (if there were), and Pearson's product-moment correlation, was performed. In this study, Pearson's product-moment correlation analysis helped quantify the assumed relationship between the dependent and independent variables. However, the result of correlation analysis indicated the multicollinearity between the independent variables (for details see Section 4.1). Thus, the partial least square regression (PLSR) model that was developed for fixing the problem of the multilinearity was employed. Further, to avoid the over-fitting problem and determine a reasonable model with the appropriate number of components that has good predictive ability, the leave-one-out (LOO) method was used for cross-validation, by selecting the model with the highest average predicted R 2 and the lowest average prediction sum of squares (PRESS). The validated PLSR model was written as follows, where Mean_pc_BBF is the averaged pc_BBF of two summertime days since the overall parcel-level BBFD products on two dates were highly similar to each other (for details, see Appendix A Figure A1), so were the pc_BBF; α 1 is the intercept/constant item, β 1~β4 are the partial coefficients; X 1~X4 are the Box-Cox transformed area proportions of the paved surface, waterbody, vegetation, and building within each land parcel, respectively (for details see Appendix A Table A2). ε is the error term. Subsequently, based on the validated PLSR model, nonlinear programming with constraints towards minimizing Mean_pc_BBF through the optimized LULC components was written as follows, Min α 1 +β 1 ·X 1 +β 2 ·X 2 +β 3 ·X 3 +β 23 ·X 2 ·X 3 +β 4 ·X 4 (8) where the definitions of α 1 , β 1~β4 , X 1~X4 are the same as Equation (7). X ' 1~X ' 4 are the raw area proportions of the paved surface, waterbody, vegetation, and building within each land parcel. α 0 is the constant item of the cubic regression model, β' 1~β ' 3 are the coefficients. BBF_Threshold is the threshold value of pc_BBF in response to the changing area proportion of vegetation.
The statistical processes employed in this study were performed with R 3.6.2 [59] and the pls library for PLSR [60]. The Rsolnp library [61] was used to estimate the optimized LULC components for minimizing Mean_pc_BBF.

Synoptic Analysis of Parcel-level LULC Components and Mean_pc_BBF
Overall, Figures 3 and 4 show the highly similar patterns of parcel-level LULC components and associated intra-SUHI effect measured with Mean_pc_BBF at four UFZs. As shown in Figure 3, there Energies 2020, 13, 1678 8 of 17 were significant variations of parcel-level LULC components measured by their area proportions. The buildings occupied the highest area proportion of the land parcels (averaged 47.72%), followed by the paved surfaces (averaged 24.59%) and vegetation (averaged 23.98%), whereas waterbody occupied the lowest area proportion (averaged 1.25%). In response to the fragmented vegetation and waterbodies within the total landscape, Figure 4 shows the overall unevenly distributed pattern of Mean_pc_BBF. As shown, the clusters of pixels with a paved surface and dense buildings exhibited higher BBFD than the spatially sparse pixels with waterbodies and vegetation. By contrast, the aggregated pixels with dominant vegetation and waterbodies, such as the university campus, big parks, boulevards, and a few well-vegetated residential communities, exhibited significantly lower Mean_pc_BBF (

Interpretation of Relationship between Parcel-Level LULC Components and Mean_pc_BBF
In the sense of qualitative statistics, Table 2 shows the statistically significant relationships between the Box-Cox transformed parcel-level LULC components and the intra-SUHI effect measured with Mean_pc_BBF, as evidenced with the correlation coefficients indicating the multicollinearity that may mislead our understanding of the relationships as mentioned above. However, it can be seen that except for the positive relationship between X2 and X3, there were the

Interpretation of Relationship between Parcel-Level LULC Components and Mean_pc_BBF
In the sense of qualitative statistics, Table 2 shows the statistically significant relationships between the Box-Cox transformed parcel-level LULC components and the intra-SUHI effect measured with Mean_pc_BBF, as evidenced with the correlation coefficients indicating the multicollinearity that may mislead our understanding of the relationships as mentioned above. However, it can be seen that except for the positive relationship between X2 and X3, there were the negative relationships between the other LULC components, indicating the overall exclusive competition of land use purpose for urban land development and associated influence on Mean_pc_BBF. Figure 5 visually addresses the bivariate relationships between fine-scale parcel-level LULC components (e.g., building and vegetation) and Mean_pc_BBF. The observed positive/negative bivariate relationships with remarkable variations in Mean_pc_BBF are wholly in agreement with the heterogeneity of parcel-level LULC components and associated Mean_pc_BBF shown in Figures 3 and 4. Based on Figure 5b, the BBF_Threshold in equation 10 was estimated.     Table 2) that causes the multicollinearity, only two of them were selected to draw the scatterplots.
Further, beyond addressing the bivariate relationships shown in Table 2 and Figure 5, Table 3 comprehensively shows the quantitative relationship between the Box-Cox transformed parcel-level LULC components and their roles in controlling Mean_pc_BBF, explaining 56.1% variation of the Mean_pc_BBF associated with independent variables. As shown, along with the waterbody, the interaction of waterbody and vegetation exhibited the statistically significant negative coefficients, indicating their capacity for suppressing Mean_pc_BBF when fixing the other variables. In contrast, paved surface, vegetation, and building exhibited statistically significant positive coefficients, indicating their capacity for increasing Mean_pc_BBF. However, herein, the unstandardized coefficients (Coefs) of the Box-Cox transformed variables result in puzzling understanding as they cannot be directly comparable for indicating each specific positive/negative relationship with Mean_pc_BBF.   Table 2) that causes the multicollinearity, only two of them were selected to draw the scatterplots.
Further, beyond addressing the bivariate relationships shown in Table 2 and Figure 5, Table 3 comprehensively shows the quantitative relationship between the Box-Cox transformed parcel-level LULC components and their roles in controlling Mean_pc_BBF, explaining 56.1% variation of the Mean_pc_BBF associated with independent variables. As shown, along with the waterbody, the interaction of waterbody and vegetation exhibited the statistically significant negative coefficients, indicating their capacity for suppressing Mean_pc_BBF when fixing the other variables. In contrast, paved surface, vegetation, and building exhibited statistically significant positive coefficients, indicating their capacity for increasing Mean_pc_BBF. However, herein, the unstandardized coefficients (Coefs) of the Box-Cox transformed variables result in puzzling understanding as they cannot be directly comparable for indicating each specific positive/negative relationship with Mean_pc_BBF. Alternatively, by measuring the changing standard deviation of Mean_pc_BBF in response to each one standard deviation increase in a given Box-Cox transformed variable, the standardized coefficients (S-Coefs) can better indicate the relative strength of the relationship between the Box-Cox transformed variables and Mean_pc_BBF. As indicated, the waterbody was found to be the strongest negative determinant of Mean_pc_BBF when fixing the other variables, as with each one standard deviation increase in waterbody would approximately decrease 0.291 standard deviations of Mean_pc_BBF. The interaction of waterbody and vegetation was found to be the second strongest negative determinant of Mean_pc_BBF when fixing the other variables, as with each one standard deviation increase in the interaction of waterbody and vegetation would approximately decrease 0.196 standard deviations of Mean_pc_BBF. Besides, vegetation was found to be the smallest positive determinant of Mean_pc_BBF when fixing the other variables, as with each one standard deviation increases in vegetation would approximately increase 0.085 standard deviations of Mean_pc_BBF. In contrast, when fixing the other variables, paved surface and building were found to be the mediate and highest positive determinants of Mean_pc_BBF, as with the relatively remarkable increase in standard deviations of Mean_pc_BBF, which respond to each one standard deviation increase in paved surface and building, respectively. The use of the S-Coefs made sense for the results. Table 4 shows the remarkable changes in summertime Mean_pc_BBFs associated with present and optimized parcel-level LULC components. As shown, except for the slight change of proportion of the other impervious surface in the present and optimized parcel-level LULC components, the changes in area proportions of the other LULC components were remarkable. As a whole, compared with present LULC components, the optimized parcel-level LULC components with the thresholding area proportions would cause an overall 51.80% decrease in summertime Mean_pc_BBFs.

Comparison of the Present and Optimized Parcel-level LULC Components and the Associated Mean_pc_BBF
On the other hand, if fixing the area portion of the other impervious surface constant, it is an excellent way to enhance the parcel-level area proportions of the waterbody and vegetation and simultaneously decrease area proportion of Building. However, due to the remarkable variation of parcel-level LULC components, there were only 7.79% of the land parcels (big parks and memorial cemetery), ultimately meeting the optimized thresholding area proportions of the parcel-level LULC components. In contrast, the land parcels with dominant impervious surface approximated 44.31% of total land parcels. Then, focusing on the gap between the present and ideally optimized area proportions of parcel-level LULC components, there would be higher uncertainties of minimizing the Mean_pc_BBF through adjusting urban land development and optimizing the parcel-level LULC components. How to seek practical solutions will be further emphasized in the discussion section.

On the Data, Methods, and Findings of this Study
The occurrence of the intra-UHI effect depends on the complicated human-nature processes in the built environment. In addition to the macro-scale climate system, the micro-scale artificial modification of urban climate due to varying thermal properties of LULC is one of the key drivers of the intra-SUHI effect. Most of the datasets used in this study are publicly accessible, and statistical analysis methods and optimization algorithm adopted are well known among the researchers. Together with the spatial interpolation-based thermal sharpening method, which was proven to be capable to capture the varying thermal properties of fine-scale LULC components [31], the PLSR analysis quantitatively examined the relative importance of parcel-level LULC components in determining the Mean_pc_BBF. Furthermore, results of the present and optimized parcel-level LULC components and their contribution to the associated Mean_pc_BBF were comparable across the land parcels, regardless of their varying size.
In the existing studies for investigating SUHI effect via satellite thermal remote sensing, most of which only used the MODIS, Landsat TM/ETM+/8, and ASTER thermal band products [62,63], or additionally combined with high-resolution aerial photos/commercial satellite satellites for added interpretation [17,64,65], the results based on the unsharpened thermal band data were too coarse to examine the relationship between the observed LST and fine-scale LULC components. In contrast, our findings, though exhibiting the uncertainties of the Mean_pc_BBF between the present and ideally optimized area proportions of parcel-level LULC components, can provide more detailed information on the fine-scale intra-SUHI effect and optional choices for decision-making towards mitigating UHI effect. Thus, the application of the accessible public datasets and the methods in this study can be referenced to the case studies elsewhere.

Implications for Practical Solutions to Optimized Parcel-Level LULC Components towards Mitigating Intra-SUHI Effect
In the modern history of downtown Shanghai, this city was in the track of international trading and industrialization manipulated by the western colonialist authorities. This city's compact form had been criticized for deteriorating the natural watery landscape and lacking green space during its growth [66]. In the period of 'urban industrialization ', this city became China's biggest industrial center, and most of the land patches were encouraged for intensive development of industrial facilities, traffic roads, and residential quarters at the loss of the creeks and farmlands within the former urban-rural mosaics [67][68][69]. The resultant remnants of waterbodies and vegetation within urban areas accounted for tiny area proportions of downtown Shanghai, making this city a typically hot city during the summer.
Recently, with the growing awareness of the better environment, the local government has been dedicated to urban regeneration and improve this city's green infrastructure via land replacement for new greenspace. The consequent increase in greenspace from 8278 ha in 1998 to 127,332 ha in 2015 was rewarded with a remarkable decrease in previous UHI hotspots [70][71][72]. However, as addressed in Section 4, the overall competitive parcel-level LULC components embodied the heterogeneity of urban land development intention and associated SUHI effect. As evidence, the cumulative area proportions of waterbodies and vegetation are still less than the recommended minimum 30% baseline of the study area, which can better interpret the weakness of present parcel-level waterbodies and vegetation in modifying the SUHI effect. To achieve the goal of minimizing Mean_pc_BBF, it is essential to trade-off the conflicts rooting in the facets of urban land development intention, then further reduce and remove the significant mismatching in area proportions between the present and optimized parcel-level LULC components. To do so, regeneration of the land parcels with dominant buildings and hard-top pavements well as architectural renovation should be officially given the top priority. For example, practice in ecological resilience, such as the well-managed vertical planting and roof garden, has exemplified its remarkable cooling capacity in microclimatic UHI mitigation. In addition to increasing the green ratio, it is crucial to shape the openness of urban morphology and alleviate the blocking of wind corridors that reduces convective heat removal and transfer by the wind [73]. We recommend that future land parcel design should guide the moderate density clusters of middle-and high-rises with acceptable building distance (e.g., 30 m or more, rather than merely the sunshine spacing).

Limitation of this Study and Future Researching Tasks
There are several noticeable shortcomings of this study. Firstly, due to the 16-day revisiting interval of Landsat 8 satellite and cloud contamination, only two cloud-free images are available for this study. Thus, the estimated parcel-level BBF, which was generated from the thermally sharpened LST products, can only capture the instantaneous field of view rather than the continuous scenes of the study area. Secondly, until the present, except for a few five-year intervals of aerial photography performed in earlier years, there have no available airborne sensors like the ATLAS at 10-m for investigating the BBF in downtown Shanghai. Alongside with lacking the airborne sensed data, the absence of in-situ measured parcel-level BBF datasets made the cross-validation of the estimated parcel-level BBF impossible. Thirdly, the three-dimensional features of urban morphology and their influences on micro-climate and human thermal comfort were not considered in this study. Consequently, in the sense of data demand for near real-time, temporal continuity, spatial representativeness, the direct linkage between the retrieved parcel-level BBF, human thermal comfort, and energy requirement for cooling was insufficiently illustrated. To fill these emphasized knowledge gaps, in future research design, the in-situ measurement of micro-climatic parameters, high-resolution land cover types, the three-dimensional features of urban morphology, should be incorporated into the computational fluid dynamics (CFD) simulation platforms (e.g., ENVI-met and WindperfectDX) to cross-validate the satellite observation and in-situ measurements, and then be further refined to generate the tempo-spatially continuous datasets.
In short, we argue that focusing on these above-mentioned shortcomings, to find the good questions and attempt to comprehensively put the conceived ideas into practice, then we can get close to the truth and find the practical ways to mitigate the adverse UHI effect and further enhance the adaption capacity to climate change.

Conclusions
Selecting downtown Shanghai as a case, this study quantitatively examined the influences of parcel-level LULC components on summertime intra-SUHI effect and further analyzed the potential of mitigating summertime intra-SUHI effect through the optimized LULC components. The major findings were summarized as follows, (1) Overall, there exists a statistically significant relationship between fine-scale parcel-level LULC components and the intra-SUHI effect measured with Mean_pc_BBF. The observed remarkable variations in Mean_pc_BBF can be attributed to the heterogeneity of parcel-level LULC components.
(2) Focusing on the influence of parcel-level LULC components on intra-SUHI effect, the relative importance of waterbodies and vegetation in terms of their contributions to decreasing associated Mean_pc_BBF urgently require increasing their area proportions. Theoretically, the optimized parcel-level LULC components would cause an estimated 51.798% decrease in summertime Mean_pc_BBFs. However, how to balance the conflicts between the present and ideally optimized parcel-level LULC components towards minimizing the Mean_pc_BBF should be carefully considered. It is still a long and formidable challenge for decision-making of sustainable land development and mitigating the UHI effect.
(3) In contrast to the relative coarse scales of region and city levels, our study demonstrated a practical approach linking fine-scale parcel-level LULC components and intra-SUHI effect, using the integration of satellite-based thermal sharpening, statistical analysis, and nonlinear programming with constraints. The methodology and findings presented beneficial insights for guiding sustainable urban land development towards enhancing the city's adaption capacity to climatic change for megacities like Shanghai.

Acknowledgments:
The authors thank the Geospatial Data Cloud site, Computer Network Information Center, Chinese Academy of Sciences (http://www.gscloud.cn) for providing the free Landsat TM/OLI images. The authors are indebted to the R Foundation for Statistical Computing and Beijing Piesat Information Technology Co., Ltd for free usage of the PIE 6.0 Remote Sensing image processing system.

Conflicts of Interest:
The authors declare no conflict of interest. Table A1. Pixel-based root-mean-square error (RMSE) between the original 30 m and the resampled 30 m LST products from the sharpened products (1-9 m).  Chinese Academy of Sciences (http://www.gscloud.cn) for providing the free Landsat TM/OLI images. The authors are indebted to the R Foundation for Statistical Computing and Beijing Piesat Information Technology Co., Ltd for free usage of the PIE 6.0 Remote Sensing image processing system.

Conflicts of Interest:
The authors declare no conflict of interest.