Access the full text.
Sign up today, get DeepDyve free for 14 days.
References for this paper are not available at this time. We will be adding them shortly, thank you for your patience.
Article A Preliminary Study on the Impact of Landscape Pattern Changes Due to Urbanization: Case Study of Jakarta, Indonesia 1,2,3, 4 1,2 Dikman Maheng *, Assela Pathirana and Chris Zevenbergen Faculty of Civil Engineering and Geosciences, Delft University of Technology, Stevinweg 1, 2628CN Delft, The Netherlands; firstname.lastname@example.org Department of Coastal and Urban Risk and Resilience, IHE Delft, Institute for Water Education, WestVest 7, 2611AX Delft, The Netherlands Department of Environmental Engineering, Universitas Muhammadiyah Kendari, Jl. Ahmad Dahlan 10, Kendari 93117, Indonesia UNDP Maldives Country Office, United Nations Development Programme (UNDP), 4th Floor, H. Aaage (Bank of Ceylon Building), Boduthakurufaanu Magu, Malé 20026, Maldives; email@example.com * Correspondence: firstname.lastname@example.org Abstract: Urbanization is changing land use–land cover (LULC) transforming green spaces (GS) and bodies of water into built-up areas. LULC change is affecting ecosystem services (ES) in urban areas, such as by decreasing of the water retention capacity, the urban temperature regulation capacity and the carbon sequestration. The relation between LULC change and ES is still poorly examined and quantified using actual field data. In most ES studies, GS is perceived as lumped areas instead of distributed areas, implicitly ignoring landscape patterns (LP), such as connectivity and aggrega- tion. This preliminary study is one of the first to provide quantitative evidence of the influence of landscape pattern changes on a selection of urban ecosystem services in a megacity as Jakarta, In- donesia. The impact of urbanization on the spatiotemporal changes of ES has been identified by Citation: Maheng, D.; Pathirana, A.; considering connectivity and aggregation of GS. It reveals that LP changes have significantly de- Zevenbergen, C. A Preliminary creased carbon sequestration, temperature regulation, and runoff regulation by 10.4, 12.4, and Study on the Impact of Landscape Pattern Changes Due to 11.5%, respectively. This indicates that the impact of GS on ES is not only determined by its area, Urbanization: Case Study of Jakarta, but also by its LP. Further detailed studies will be needed to validate these results. Indonesia. Land 2021, 10, 218. https://doi.org/10.3390/land10020218 Keywords: urbanization; land use–land cover change; landscape pattern; green space area; ecosys- tem services; Jakarta Received: 12 December 2020 Accepted: 19 February 2021 Published: 21 February 2021 1. Introduction Publisher’s Note: MDPI stays neu- The urban population had increased to 3.9 billion (54% of the total population) by tral with regard to jurisdictional claims in published maps and insti- 2014 and will increase to around 6.3 billion (67% of the total population) by 2050 . Pop- tutional affiliations. ulation increase drives land use–land cover (LULC) change transforming green spaces as well as bodies of water into built-up areas resulting in decrease and loss of ecosystem services (ES) in urban areas [2–5]. Ecosystem services are the varied benefits to humans from the natural environment and healthy ecosystems. They can be categorized as provi- Copyright: © 2021 by the authors. Li- sioning services, regulating services, cultural services, and supporting services . A city censee MDPI, Basel, Switzerland. or an urban area can be seen as a single ecosystem where the urban population directly This article is an open access article benefits from ecosystem services generated by urban ecosystems, such as green and blue distributed under the terms and con- infrastructure . Ecosystem services in urban areas include air quality improvement, ditions of the Creative Commons At- noise reduction, controlling air temperature, carbon storage and sequestration, water reg- tribution (CC BY) license (http://crea- ulation, and recreation [7–13] and therefore enhance citizens well-being [14,15]. tivecommons.org/licenses/by/4.0/). Many studies show that urbanization and LULC change are key anthropogenic driv- ers affecting urban ecosystem’s functions and services [11,16–19]. Some examples of the impact of LULC change on ES due to urbanization are a decrease in the capacity to retain Land 2021, 10, 218. https://doi.org/10.3390/land10020218 www.mdpi.com/journal/land Land 2021, 10, 218 2 of 27 stormwater leading to urban flooding , and an increase in the land surface tempera- ture (LST)  in cities. Yuan et al. showed that in Nanjing, China, there was a decrease in the provision and delivery of supply services, the supporting services, and the regulat- ing services by 17, 0.9, and 4.3%, respectively, due to urbanization between 2000 and 2015. In Beijing, urban growth between 1985 and 2015 decreased carbon storage and water yield by 3.3 and 4%, respectively, while sediment export increased by 15.5% (Sun et al., 2018b). Another study from Atlanta Metropolitan area, US, showed that urban expansion be- tween 1985 and 2012 reduced carbon storage by 23% and water purification capacity by 28 and 49%, for nitrogen and phosphorus, respectively, but increased sediment transport by 17% (Sun et al., 2018). Jaligot et al.  showed that the cultural services in Yaoundé, Cameroon, decreased by more than 90% due to urbanization between 2000 and 2018. LULC change from 1990 to 2017 in the West Bengal, India, resulted in some of the regu- lating services increasing while others decreased. Water regulation, water supply, and waste treatment increased by 57, 22, and 8%, respectively, while climate regulation and nutrient cycling decreased by 31.82 and 31.25%, respectively . The relationship between urbanization, LULC change, and ES has been studied from various perspectives. Panagopoulos et al.  discussed the importance of urban ES and landscape planning to sustainable urban development. Das and Das  investigated the impact of LULC change on ES by quantifying the total ecosystem service value (ESV) based on the area of the ecosystem, and a value coefficient (VC) . The extension of built-up areas in West Bengal decreased agricultural and vegetation covers. As a result, the total ESV of landscape decreased by 24.30% from 1990 to 2017. Wang et al.  inves- tigated the impact of urbanization from 2000 to 2010 on multiple ES at hotspots and urban megaregion scale of Beijing–Tianjin–Hebei (BTH), the capital economic zone of China, by considering three components of urbanization: the built-up area proportion, the popula- tion density, and the gross domestic product (GDP) density. The ES in their study was calculated using the area of the ecosystem, for instance, the total food production was based on the area of cropland. It was revealed that providing services of food supply had an “inverse U” shape relation with the population and the GDP density, while urban ex- pansion decreased regulating services. The impact of urbanization on ES was not only influenced by an increase in urban areas or a decrease in green spaces but also by land- scape pattern changes, such as a number of patches (NP), aggregation index (AI), and “clumpiness” index, as discussed by Asadolahi et al. . The ecosystem services in their study were based on the area of ecosystems, including water body, agriculture, rangeland, forest, and urban. Moreover, the influence of landscape patterns to ES provision was an- alyzed by comparing the landscape metric changes and ES. They saw that increasing food supply had a positive correlation with the good connectivity and connectedness of agri- cultural land. Haas et al.  discussed the relationship between a decrease in ES with an increase in built-up areas and landscape pattern changes due to urbanization in the period from 1989 to 2001 in Stockholm and Shanghai. They observed that an increase in built-up areas along with landscape pattern changes, such as a decrease in connectivity, an increase in fragmentation, and more complex shapes of the natural landscape, had a negative im- pact on ES. The ES in their study was based on the area of ecosystems, such as wetland, water, forest, agriculture, urban green space, high-density building, and low-density building. Moreover, the relationship between ES and landscape pattern changes was dis- cussed by comparing the changes of the landscape metrics and ES. For urban flood miti- gation, Bai et al.  conducted a field survey study in Luohe, China, which observed that inundated areas had a low proportion of green spaces and bare land. Moreover, landscape metrics analysis showed that the mean patch size of green spaces of flooded areas was smaller compared with that of roofs. Meanwhile, the connectivity of green spaces had no significant contribution to mitigating urban flooding. The importance of landscape pat- terns on surface runoff control was also observed by . They directly integrated large patch index (LPI) and aggregation index (AI) of green spaces with runoff coefficients to estimate runoff volume controlled in Beijing, China. It was found that increasing runoff Land 2021, 10, 218 3 of 27 was related to a decrease in green spaces, as well as increasing the fragmentation and disaggregation of green spaces. The landscape patterns can also exhibit a relation with an increase in LST in the Olympic forest park of Beijing, China, as investigated by using the Pearson correlation analysis . It was revealed that an increase in LST was related to increasing urbanization in the period 2000 to 2015. Moreover, there was a negative corre- lation between LST and the landscape metrics of large patch index (LPI) and aggregation index (AI) of green spaces. A growing number of studies have investigated the impact of urbanization and LULC change on ES. These studies reveal that a decrease in green spaces affects ES, feed- ing the perception that there is a positive relationship between the surface area of green spaces and urban ecosystem functioning [4,32]. However, these studies generally ignore the impact of spatial heterogeneity of the land surfaces as they lump together the area of each land-use type. Only a few of them acknowledged the differences in landscape pat- terns but still did not consider landscape patterns in calculating the impact of green spaces. As a result, an area of green spaces, as shown in Figure 1, could still be perceived as a lumped area instead of a distributed area, where the green space units have different landscape patterns, such as connectivity, aggregation, mean patch size, and large patch index, which could affect ecosystem services [33,34]. Figure 1. Schematic illustration of the loss of information due to lumping green areas. (a) Real landscape pattern; (b) the lumped area used in many previous ecosystem services (ES) studies. In many existing studies, the impact of urbanization and LULC change on ES is gener- ally discussed using the lumped area of green spaces. Some studies have discussed the im- pact by comparing the ES changes and the changes of landscape patterns, which show that landscape patterns can affect ES. However, there is a lack of studies that directly examine the impact of landscape pattern changes on the area of green spaces to ES calculation. This preliminary study aims to identify the impact of landscape pattern changes on urban ecosystem services using two spatial configurations of green spaces, as shown in Figure 1. Firstly, the area of green spaces is assumed as a lumped area, as depicted in Figure 1b. Secondly, the area of green spaces is assumed as a distributed area, as shown in Figure 1a, which is represented by directly combining the area and landscape patterns defined by the aggregation and the connectivity of the green spaces and bodies of water. As a preliminary study, this study is designed to focus on the investigation of the land- scape pattern changes and their influence on urban ES. It is expected that the results of this study can provide insight into the influence of landscape patterns on ES calculation. Moreover, the validation of the results is not a part of this study since there is a lack of field survey data for all municipalities in Jakarta. Hence, the validation would be consid- ered as a part of further studies. Land 2021, 10, 218 4 of 27 This study is divided into three main parts: (1) assessing the impact of urbanization on LULC change; (2) analyzing the impact of LULC change on the landscape patterns; and (3) analyzing the impact of the landscape pattern changes on ecosystem services. 2. Materials and Methods This study was carried out in Jakarta, Indonesia, and it was conducted in six steps. First, LULC maps for 1995 and 2014 were generated using Landsat images. The LULC classification used six LULC classes, which were urban, suburban, grassland, cropland, trees, and water (sea and surface water). Second, the results of LULC classification were imported into the FRAGSTATS, the spatial pattern analysis program, for landscape metric analysis. The analysis considered four landscape metrics and the relationship between them: Class Area (CA), Land Proportion (PLAND), Aggregation Index (AI), and connec- tivity (COHESION). The landscape metrics analysis was carried out for five municipalities of Jakarta. Third, the impact of urbanization on LULC change was analyzed, which was based on the CA of LULC classes. Fourth, the impact of LULC change on the landscape patterns was analyzed, which was based on PLAND, AI, and COHESION. Fifth, ecosys- tem services were calculated using two methods: area-based estimation and landscape metric area (LMA). The area-based estimation followed commonly ecosystem services area calculation that uses only the area of green spaces, while the landscape metric area directly incorporated the landscape metric representing the spatial characteristic/config- uration of green spaces to ecosystem services area calculation. Sixth, ecosystem services changes between 1995 and 2014 from two different methods were then analyzed using ecosystem services index (ESI). 2.1. Study Area Jakarta is the capital of Indonesia, which has experienced rapid urbanization and the extensive transformation of vegetative green spaces into impermeable urban surfaces [35– 38]. Urbanization in Jakarta is driven by several factors, including its strategic political and economic position as the capital city where important activities like the functions of the national government, education, manufacturing, and commerce are taking place . Jakarta has six municipalities, which are Central Jakarta, North Jakarta, West Jakarta, South Jakarta, and East Jakarta, as shown in Figure 2, and one regency, which is the Thou- sand Islands (not shown in the figure). This study focuses on the municipalities that con- sist of 42 subdistricts, which are home of approximately 10 million people, with the pop- ulation growth rate of 1.02% per year in 2015 . East Jakarta is the most populated area, whereas West Jakarta is the most densely populated area with a population density of 19,018 per km . Located on the northwest coast of Java, Jakarta’s climate is a tropical mon- soon climate according to the Köppen climate classification system and it has two seasons, which are the rainy season from October through May and the dry season from June to September. During the rainy season, the average monthly rainfall is 300 mm, while it is 43 mm in August during the dry season . Land 2021, 10, 218 5 of 27 Figure 2. (a) Map of Indonesia; (b) Jakarta is located on Java island, which is the densest island in Indonesia; (c) Jakarta consists of five municipalities. 2.2. Data Sources The impact of urbanization on LULC change was analyzed using land use–land cover maps of Jakarta for the year 1995 and 2014. The maps were prepared using Landsat images selected for the area of Jakarta and its surrounding. For the year 1995, the Landsat 5 TM L1TP from 24 August 1995, was used, while the Landsat 8 OLI TIRS L1TP obtained on 13 September 2014, was selected for the land use–land cover map of the year 2014. Both da- tasets had 30 m resolution. The land cover classification was based on the land cover classes used to investigate the relationship between green spaces and increasing urban temperature . We used the semiautomatic classification in QGIS , a free and open-source geo- graphic information system, to generate six land cover classes: urban, suburban, grassland, cropland, trees, and water. An urban area is an area where there is a dense mix of compact high-rise, compact midrise, or low-rise buildings, while a suburban area is an open arrange- ment of low-rise buildings with pervious land covers, such as grass, low plants, and scat- tered trees. A grassland area is defined as an area where the vegetation is dominated by grasses, and there are few or no other plants. A cropland area is an area mainly used for agricultural and farming purposes, where rice, corn, vegetables, and fruits are grown. An area of trees is an area where trees are dominant. Water represents both sea and surface water, such as rivers and lakes. Moreover, the LULC change analysis was based on the CA from five land cover classes: urban, trees, grassland, cropland, and bodies of water. The boundary area of Jakarta and its municipalities was obtained online from the Geospatial Information Agency of Indonesia (BIG) database, www.tanahair.indone- sia.go.id (accessed on 20 September, 2019). 2.3. Landscape Metrics Landscape metrics are quantified characteristics of the spatial pattern at different lev- els (i.e., patch level, class level, and landscape-level) used to describe landscape structures and LULC change [43–45]. Information about landscape structures is essential to analyze ecosystem structures, functions, and services  and to evaluate their spatial effects at the patch–class–landscape level [47,48]. Patch level metrics are important for defining in- dividual patches and for characterizing the spatial context of the patches, which are used Land 2021, 10, 218 6 of 27 as the computational basis for other landscape metrics. At the class level, the metrics are calculated by integrating all patches from a given class. The class level calculation can be done by a simple aggregation of the area of patches to reflect the total area of a given class. Moreover, landscape-level metrics are the integration of all classes over the landscape . This study uses four landscape metrics to capture the spatial features which affect urban ecosystem processes: class area (CA), connectivity (COHESION), land proportion (PLAND), and aggregation index (AI). CA is a measure of landscape composition; specif- ically, to indicate how much of the landscape is comprised of a patch type. The class area is the primary component used in the ecosystem service quantification . The connec- tivity of a class is measured by the patch cohesion index (COHESION) that measures the physical connectedness of the corresponding patch type. Patch cohesion increases as the patch type become more clumped or aggregated in its distribution; hence, more physically connected . The decrease in patch connectivity in a landscape could be a negative sign for ecosystem services . The land proportion (PLAND) is the percentage of a land- scape, which quantifies the proportional abundance of each patch type in the landscape, while the AI is used to understand the spatial aggregation of a patch in a landscape [31,45,51]. The landscape metrics in this study were calculated using FRAGSTATS 4.2, the spatial pattern analysis program, which is developed by McGarigal and Cushman of Uni- versity of Massachusetts in the United States of America . The landscape metrics used in this study were selected based on the premise that there is a correlation between ES and the area, the aggregation, and the connectivity of ecosystem [31,45,52]. In general, many ES studies focus on the relationship between eco- system services and the area of green spaces revealing a positive correlation between the size of the green spaces and the weight of the ecosystem services provided [53–55]. How- ever, some studies showed that ecosystem services were also influenced by landscape pat- terns [29,30,55,56]. In this study, ES was analyzed by considering the “effective area” of green spaces and bodies of water by taking into account the influence of aggregation and connectivity on the area of green spaces and bodies of water. The landscape metric com- binations for this study are selected based on the FRAGSTATS 4.2 guidelines and the ex- isting research on landscape metric and ES. There are some indices that can be used to analyze the relationship between landscape metrics and ES, such as Number of Patch (NP), Patch Density (PD), Large Patch Index (LPI), patch cohesion index (COHESION), land proportion (PLAND), CONTAGION, and Aggregation Index (AI). Among them, the AI was selected because it can give more accurate results compared to that of other aggre- gation indices, such as CONTAGION , while the COHESION was selected for connec- tivity analysis because it can represent the connectivity of a corresponding patch type of a particular class [28,50]. This resulted in the “Landscape Metric Area” (LMA) calculated by multiplying the metrics class area (CA), aggregation index (AI), and cohesion index (COHESION) for each land cover class, as follows: (1) 𝐿𝑀𝐴 = 𝐶𝐴 ×𝐴𝐼 ×𝐶𝑂𝐻𝐼𝑂𝑁𝐸𝑆 LMA was calculated at the municipality level following the Indonesian regulation for spatial planning. The standard focuses on the city level and defines the spatial config- uration of different land use–land cover types, including the minimum proportion of open green spaces . 2.4. Urban Ecosystem Services Many studies discussed ES in urban areas [7,9,28,45,58,59]. Back in the 1990s, one of the pioneering urban ecosystem services studies considered six ecosystem services: air filtering, microclimate regulation, noise reduction, rainwater drainage, sewage treatment, and recreation . Derkzen et al.  studied ecosystem services bundles in Rotterdam, the Netherlands, considering six ecosystem services: air purification, carbon storage, noise reduction, runoff regulation, cooling, and recreation. The importance of urban ecosystem services inspires local and national governments to incorporate ES in their urban planning Land 2021, 10, 218 7 of 27 process and documents [58,60–62]. For example, the government of Indonesia requires a minimum of 30% open green spaces in each city and urban area [57,63]. In this study, ES changes analysis was limited to carbon sequestration, temperature regulation, and runoff regulation, due to the limitation of the medium resolution satellite images to capture several urban details, such as street corridors and small city parks. Fur- thermore, in the ES selection increasing carbon emissions [64,65], and recent climate-re- lated hazard events in Jakarta, such as urban flooding in 2013, 2014, 2015, and 2020, have also been considered [66,67], as well as the potential increase in the urban temperature [37,68]. The ES estimation adopted a weighting factor used in , where the green space with high-temperature regulation capacity was given a weight of 1.0, whereas 0.5 was given for low capacity. In this study, the weighting factor was implemented for all ES estimations for which the factor was given on the basis of the importance of a particular green space to a particular ES. For example, grassland was given the weight of 0.5 for temperature regulation, but it was given 1.0 of the weight factors for runoff regulation. 2.4.1. Carbon Sequestration Green spaces have an important role in the carbon cycle, capturing CO2 from the at- mosphere during photosynthesis, and releasing it during respiration. In an urban area, carbon sequestration mainly takes place in green spaces [69–71]. Each green space has a different carbon sequestration rate, and urban trees generally have a higher rate compared to that of other green spaces because of their higher leaf area index (LAI), which is the ratio of the area of trees’ leaves per area of the land surface . For example, the urban forest of Punggol forest in Singapore had the annual carbon sequestration rate of 1.6 ton/year/ha, where trees with high LAI such as Delonix regia and Artocarpus heterophyllus gave a significant contribution . Studies from Indonesian cities showed that trees had a carbon sequestration rate higher than that of other vegetation, including grass and her- baceous. For instance, the carbon sequestration of tree-covered areas in a part of Jakarta was 129.92 kg/ha/hour, much higher than 2.74 kg/ha/hour for grass areas . Moreover, trees in the city of Singaraja on Bali island had 112.751 tons/ha of carbon sequestration higher than 4.845 tons/ha of carbon sequestration of herbaceous . In this study, carbon sequestration was based on the LMA of trees. The LMA of trees was given a weight of 1.0 considering its carbon sequestration, which is higher than that of other vegetation. Other types of green spaces were given a weight of 0, assuming no significant contribution to capturing carbon in the atmosphere. 2.4.2. Temperature Regulation Increasing LST and air temperature in urban areas can be influenced by several factors, including a decrease in green spaces and bodies of water. The contribution of green spaces to the control of urban temperatures can be achieved through ecosystem functions such as evapotranspiration, trees shading, and modifying air movement, while the greatest contrib- uting factor to the cooling effect is shading and evapotranspiration . The effectivity of green spaces in providing the ecosystem functions is influenced by various factors, includ- ing the area and the spatial distribution of green spaces , the configuration and the shape of green spaces [52,77], and the connectivity of green spaces . Among green spaces, trees have high cooling potential through shading and evapotranspiration, while other short vegetative covers such as grass and crops might have a lower capacity [9,78]. The temperature in urban areas can also be regulated through evaporation and heat ab- sorption from bodies of water . The cooling effect of bodies of water is influenced by the area of the bodies of water , and the effect decreases with the increase in distance from the water body, which is similar to the cooling effect from green spaces [81,82]. Temperature regulation in this study was estimated using the combination of the LMA from trees, grassland, cropland, and bodies of water. The LMA of trees was given a weight of 1.0, while a weight of 0.5 was given for the LMA of grassland, cropland, and bodies of water. Land 2021, 10, 218 8 of 27 2.4.3. Runoff Regulation Due to urbanization, cities and urban areas are more vulnerable to urban flooding. The increasing vulnerability is influenced by some factors such as a decrease in green spaces and bodies of water. As a result, many urban areas have an insufficient capacity for infiltration and interception, as well as an increase in the runoff coefficient. Urban flooding-related studies showed that the area of green spaces could be considered as one of the important factors that determine runoff regulation [30,83,84]. For instance, increasing tree-cover by nearly 11% in Beijing could increase the total runoff reduction volume by more than 30% . Moreover, the capacity of green spaces in controlling surface runoff is also influenced by several factors, including the interception of rainfall. Among many types of green spaces, trees could provide the high interception capacity, particularly in the short– moderate rainfall or in the earlier periods of precipitation, and it decreases in the longer duration or during high rainfall intensity . A field experiment from the city of Uruaçu, Goiás, Brazil, showed that interception rate during short, low-intensity precipitation was about 40%, while the rate was about 3.6% in the high intensity and long duration rainfall . Besides interception, green spaces also provide runoff regulation through infiltration, which is influenced by numerous factors including rainfall intensity, green space types, soil characteristics, and green space coverage [19,88–90]. In addition, green spaces con- tribute to runoff coefficients by increasing flow resistance, which can be associated with its density . Furthermore, bodies of water like rivers, lakes, and wetlands also have the potential to control surface runoff through storage functions . The estimation of runoff regulation service in this study considered the combination of the LMA from trees, cropland, grass, and bodies of water. All green spaces were given a weight of 1.0 based on the assumption that green spaces and bodies of water had a sim- ilar contribution to runoff regulation from one or the combination of infiltration, intercep- tion, flow resistance, and storage functions. 2.5. Ecosystem Services Index (ESI) This study used ecosystem services index (ESI) to assess the ES changes. The ESI was based on the ES calculated using the combination of the LMA from green spaces and bod- ies of water. The ES calculation was divided into two steps. First, the ES was only based on the area of green spaces and bodies of water called area-based estimation. Second, the ES was based on the LMA calculated using the following equation. = 𝑎 𝐿𝑀𝐴 +𝑏 𝐿𝑀𝐴 +𝑐 𝐿𝑀𝐴 +𝑑 𝐿𝑀𝐴 (2) 𝑟𝑒𝑒𝑡 𝑑𝑛𝑠𝑙𝑎𝑔𝑟𝑎𝑠 𝑛𝑑𝑝𝑙𝑎𝑐𝑟𝑜 𝑟𝑤𝑎𝑡𝑒 where a, b, c, and d were the weighting factors for each green space and bodies of water. The weighting factor was different for each green space in each ES estimation, as given in Table 1. Table 1. The weighting factors (modified from Derkzen et al. ). Weighting Factors Carbon Sequestration Temperature Regulation Runoff Regulation a (Trees) 1 1 1 b (Grassland) 0 0.5 1 c (Cropland) 0 0.5 1 d (Water) 0 1 1 𝐸𝑆 Land 2021, 10, 218 9 of 27 The ESI was then calculated by normalizing the ES values using the following for- mula: 𝑖𝑛𝑚 = (3) 𝑎𝑥𝑚 𝑖𝑛𝑚 where is the index of ES, 𝑖 for the municipality 𝑗 , is the value from ES 𝑖 for the municipality 𝑗 , and and are the minimum and the maximum value of 𝑖𝑛𝑚 𝑎𝑥𝑚 ES 𝑖 , respectively. 3. Results 3.1. Land Use–Land Cover (LULC) Change Analysis The high rate of land transformation marked urbanization in Jakarta from 1995 to 2014, as shown in Figure 3. It shows that there has been a substantial change in the spatial distribution of green spaces due to increasing urban areas. The growing urban areas, about 44.3%, lost about 58.1% of the areas with trees, 36.1% of grassland, 32.3% of cropland, 15.9% of suburban areas, and 6.1% of bodies of water (lakes and rivers). Figure 3. Land use–land cover change from 1995 to 2014. The LULC change in the study was analyzed by focusing on the spatial extension of the urban areas and the decrease in vegetative cover and bodies of water from 1995 to 2014, as shown in Figure 4. In 1995, the spatial distribution of urban areas was concen- 𝐸𝑆 𝐸𝑆 𝑖𝑗 𝑖𝑗 𝐸𝑆 𝐸𝑆𝐼 𝐸𝑆 𝐸𝑆 𝑖𝑗 𝐸𝑆𝐼 𝑖𝑗 𝐸𝑆 𝐸𝑆 Land 2021, 10, 218 10 of 27 trated in three municipalities, which were East Jakarta (65.3 km ), followed by North Ja- 2 2 karta (63.9 km ) and South Jakarta (58.5 km ). On the other hand, the smallest urban area was about 35.4 km , which was observed in Central Jakarta. After 20 years, urban areas increased in all municipalities. The largest urban area in 2014 was about 97.8 km in West Jakarta, which increased by about 93.7% from that in 1995. Meanwhile, urban areas in East Jakarta and North Jakarta increased by almost a half from that in 1995. On the other hand, the smallest urban area, by approximately 38.4 km , was identified in Central Jakarta, which increased by 8.2%. The increase in built-up areas was also observed in suburban areas from 1995 to 2014, as shown in Figure 4. In 1995, the largest suburban area was detected in East Jakarta (36.4 2 2 2 km ), followed by South Jakarta (31.6 km ) and West Jakarta (26.9 km ). Central Jakarta had the smallest suburban area by about 4.30 km . The spatial distribution of suburban areas in 2014 was a bit different from urban areas, which increased in all municipalities. The changes of suburban areas were marked by a decrease in West Jakarta (−64.0%), fol- lowed by North Jakarta (−51.0%) and Central Jakarta (−19.6%), but they increased in East Jakarta and South Jakarta. The largest suburban area of 38.37 km was identified in East Jakarta, which increased by 5.3% from that in 1995. This increase, however, was lower than 17.7% in South Jakarta. On the other hand, the smallest suburban area of 3.4 km was detected in Central Jakarta. Figure 4. The area of land use–land cover of Jakarta in 1995 and 2014. (a) The area of urban class; (b) the area of suburban class; (c) the area of trees; (d) the area of grassland; (e) the area of cropland; (f) the area of bodies of water. Land 2021, 10, 218 11 of 27 The dynamic changes of the built-up areas affected the area of green spaces, where trees were by far the most affected, followed by grassland and cropland, as shown in Fig- ures 4 and 5. In 1995, trees were mainly distributed in East Jakarta and South Jakarta, where the area of trees was more than two times of that in the other municipalities, such as in North Jakarta and West Jakarta. The largest area of trees (39.4 km ) was observed in East Jakarta, while the smallest area of 5.0 km was in Central Jakarta. In 2014, the area of trees decreased above the average in all municipalities except in South Jakarta and Central Jakarta. A significant decrease of 65.8% was observed in North Jakarta, followed by West Jakarta (−62.2%) and East Jakarta (−61.1%). However, the area of trees in that area was still higher than that in other municipalities. For instance, East Jakarta had the largest area of 2 2 trees (15.4 km ), which was higher than 14.5 km in South Jakarta. On the other hand, the smallest area of trees (3.2 km ) was detected in Central Jakarta. Besides trees, areas of grassland were also impacted by an increase in the built-up areas, as depicted in Figures 4 and 5. In 1995, the widespread occurrence of areas of grass- land was mainly detected in East Jakarta, West Jakarta, and South Jakarta. The largest grassland area (18.3 km ) was observed in East Jakarta, which was similar to that of West Jakarta. On the other hand, Central Jakarta had the smallest area of about 2.3 km . In 2014, the grassland area remained almost the same in East Jakarta, while a significant decrease of 84.5% was observed in West Jakarta, followed by North Jakarta (−54.5%). The smallest area of grassland was observed in Central Jakarta (1.6 km ), followed by West Jakarta (2.8 2 2 km ) and North Jakarta (4.4 km ). Figure 5. The spatial and temporal distribution of green spaces and bodies of water in 1995 and 2014. A decrease in the area of green spaces was also found in the area of cropland during urbanization from the year 1995 to 2014, as shown in Figures 4 and 5. The spatial distri- bution of cropland in 1995 was dominant in East Jakarta and North Jakarta, where the 2 2 areas of cropland were 23.4 km and 22.1 km , respectively. In 2014, there was a decrease in the area of cropland. The highest decrease was identified in West Jakarta, where the area of cropland decreased by 47.1%, followed by East Jakarta (−35.9%) and North Jakarta (−25.2%). Moreover, the spatial distribution of cropland in 2014 was similar to that in 1995, Land 2021, 10, 218 12 of 27 where the dominant distribution of cropland was observed in North Jakarta (16.5 km ), followed by East Jakarta (14.9 km ). The existence of bodies of water between 1995 and 2014 was identified by a small decrease by −6.1%, as shown in Figures 4 and 5. The spatial distribution of the bodies of water was mainly observed in North Jakarta, representing lakes, fishponds, and wetlands. As a result, North Jakarta had the largest area of bodies of water between 1995 and 2014, where small patches of bodies of water, including rivers, creeks and small wetlands, were evenly distributed in all municipalities. 3.2. Landscape Metrics Analysis The rapid urbanization in Jakarta not only changed the spatial distribution of green spaces but also altered the landscape patterns indicated by changes in the proportion, the aggregation, and the connectivity of green spaces, as well as of bodies of water. In 1995, landscape metrics analysis showed the highest proportion of trees at the municipality level of about 20% in East Jakarta and South Jakarta, while other municipalities had about 10% of trees in their area. Meanwhile, trees showed a dispersed distribution in Central Jakarta, North Jakarta, and West Jakarta, where the AI was lower than 50%. However, trees had good connectivity, as shown in Figure 6. It shows that trees had the COHESION higher than 60% in all municipalities, with the highest values being observed in South Jakarta. Land 2021, 10, 218 13 of 27 Figure 6. Landscape metrics of trees. (a) Land proportion; (b) aggregation index; (c) cohesion. In 2014, the proportion of trees decreased in all municipalities, where the highest de- crease of 65.8% was observed in North Jakarta, followed by West Jakarta and East Jakarta, where the proportion decreased by 62.2 and 61.1%, respectively. Meanwhile, the highest trees proportion of 9.9% was identified in South Jakarta, followed by East Jakarta (8.3%) and Central Jakarta (6.6%), while a similar proportion was observed in West Jakarta and North Jakarta. The decrease in the proportion of trees was followed by a decrease in the aggregation and the connectivity of trees except in Central Jakarta. It was found that trees in Central Jakarta were more aggregated and connected in a small proportion. The aggre- gation and the connectivity of trees in Central Jakarta increased by 12.1 and 6.7%, respec- tively. In other municipalities, trees were more disaggregated, as indicated by a decrease in the AI values. For instance, the AI in East Jakarta and South Jakarta decreased by 27.4 and 26.3%, respectively. Furthermore, the connectivity of trees also decreased in East Ja- karta and South Jakarta, where the connectivity decreased by 19.7 and 23.8%, respectively. However, the connectivity of trees was still high in all municipalities since the COHESION was higher than 50%. Meanwhile, the dominant proportion of the urban areas also influenced the land- scape patterns of the grassland, as shown in Figure 7. In 1995, the highest proportion of grassland was about 14.5% in West Jakarta, followed by South Jakarta (10.9%) and East Jakarta (9.9%). The proportion of grassland in West Jakarta was found higher than the proportion of trees in the same year, which showed a different green space pattern in a municipality. Meanwhile, the low AI percentage indicated that grassland had dispersed distribution in all municipalities. The connectivity of grassland in 1995 was almost similar in four municipalities where the connectivity varied from 71.4% to 72.4%. The highest connectivity of grassland was observed in West Jakarta, where COHESION was 84.9%. In 2014, all municipalities saw a decrease in the proportion of grassland. The highest decrease was observed in West Jakarta, where the proportion decreased from 14.5% in 1995 to 2.2% in 2014. In addition, the location of the highest proportion shifted from West Jakarta to East Jakarta. The highest grassland proportion of 9.9% was identified in East Jakarta, which was slightly higher than 9.6% in South Jakarta. On the other hand, a simi- larly low proportion was observed in Central Jakarta, North Jakarta, and West Jakarta, where the proportion of grassland was 3.3, 3.2, and 2.2%, respectively. Furthermore, the area of grassland underwent a small amount of disaggregation in all municipalities except in East Jakarta, where the area of grassland was aggregated by 12.7%. The connectivity of grassland in 2014 was also affected by urbanization as it decreased in four municipalities, but it increased in East Jakarta. The highest decrease was 29.3%: this was observed in West Jakarta, where the connectivity reduced from 84.9% in 1995 to 60.0% in 2014. Moreover, a small decrease was found in Central Jakarta, South Jakarta, and North Jakarta, where the connectivity decreased by 5.9, 5.3, and 3.5%, respectively. Land 2021, 10, 218 14 of 27 Figure 7. Landscape metrics of grassland. (a) Land proportion; (b) aggregation index; (c) cohesion. Moreover, the landscape patterns of cropland were also subject to dynamic changes due to urbanization, as shown in Figure 8. In 1995, the highest proportion of cropland was mainly observed in North Jakarta, East Jakarta, and West Jakarta. The highest cropland proportion of 15.9% was observed in North Jakarta, while the lowest cropland proportion of 2.3% was observed in Central Jakarta. Meanwhile, the cropland distribution in the four municipalities exhibited a scattered pattern, as indicated by the AI percentage. For in- stance, cropland in Central Jakarta, East Jakarta, and South Jakarta was below 50%. More- over, the connectivity of crops was high in North Jakarta, West Jakarta, and East Jakarta, where the COHESION was 95.3, 87.9, and 84.6%, respectively. Land 2021, 10, 218 15 of 27 Figure 8. Landscape metrics of cropland. (a) Land proportion; (b) aggregation index; (c) cohesion. In 2014, the proportion of cropland decreased sharply. The highest decrease of 47.1% was identified in West Jakarta, followed by East Jakarta and North Jakarta, where cropland decreased by 35.9 and 25.2%, respectively. However, North Jakarta still had the highest cropland proportion of 11.9%, while the smallest cropland proportion of 2.6% was observed in Central Jakarta. The decrease in the proportion of cropland affected its aggre- gation in some municipalities in 2014. For instance, cropland in West Jakarta was dis- aggregated by 13.9%, followed by North Jakarta and Central Jakarta, where cropland was disaggregated by 11.3 and 13.8%, respectively. Furthermore, cropland had a small de- crease in the connectivity, but the percentage was still high. For instance, the highest de- crease of 8.4% was observed in West Jakarta, where the connectivity of crops was 87.9%. In contrast with the landscape patterns of green spaces, urbanization did not have a significant impact on the landscape patterns of bodies of water, as depicted in Figure 9. The proportion of bodies of water in 2014 was similar to that in 1995. The dynamic changes Land 2021, 10, 218 16 of 27 of the landscape patterns of bodies of water were marked by a small increase in the pro- portion of bodies of water in East Jakarta and South Jakarta, while this proportion de- creased in North Jakarta and West Jakarta. However, the aggregation and the connectivity of bodies of water increased in all municipalities in 2014, as shown in Figure 9. Figure 9. Landscape metrics of water. (a) Land proportion; (b) aggregation index; (c) cohesion. 3.3. Spatial and Temporal Distribution of ES Changes The spatial and temporal distribution of ES changes were analyzed using Equations (2) and (3). The results showed that the landscape pattern changes had a significant impact on the decrease in the carbon sequestration, temperature regulation, and runoff regulation. The spatiotemporal changes of the carbon sequestration are given in Figure 10. The picture shows that the area-based estimation and the LMA-based estimation give an al- most similar pattern of the carbon sequestration, but the decrease rates are different. The area-based estimation that used a lumped area of trees gave a higher carbon sequestration Land 2021, 10, 218 17 of 27 capacity compared to that from the LMA-based estimation, which used an effective area. The area-based estimation showed that East Jakarta had the highest carbon sequestration in 2014, followed by South Jakarta, North Jakarta, West Jakarta, and Central Jakarta. The highest decrease in the carbon sequestration was identified in North Jakarta (−65.8%), fol- lowed by West Jakarta (−62.2%), East Jakarta (−61.1%), South Jakarta (−51.8%), and Central Jakarta (−35.5%). On the other hand, the landscape pattern changes in the LMA-based es- timation gave a significant change in the carbon sequestration capacity. For instance, the LMA-based estimation suggested that North Jakarta had the highest decrease in the car- bon sequestration (−78.4%), followed by East Jakarta (−77.3%), West Jakarta (−77.1%), South Jakarta (72.9%), and Central Jakarta (−22.8%). Furthermore, East Jakarta had the highest carbon sequestration in 2014, followed by South Jakarta, Central Jakarta, North Jakarta, and West Jakarta. Figure 10. ESI for carbon sequestration in the year 1995 and 2014. (a) ESI in 1995 without LMA; (b) ESI in 2014 without LMA; (c) ESI in 1995 with LMA; and (d) ESI in 2014 with LMA. The different spatial patterns of the temperature regulation indicate the influence of the landscape patterns of green spaces and bodies of water, as shown in Figure 11. It shows that the temperature regulation from the LMA-based estimation had a higher de- crease rate compared to that from the area-based estimation. The impact of the landscape pattern changes on the temperature regulation was reflected in the LMA-based estimation that suggested the highest decrease rate in West Jakarta (−74.6%), followed by South Ja- karta (−64.1%), East Jakarta (−61.1%), North Jakarta (−50.1%), and Central Jakarta (−23.7%). Moreover, the highest temperature regulation in 2014 was in East Jakarta, followed by North Jakarta, South Jakarta, West Jakarta, and Central Jakarta. On the other hand, the area-based estimation showed a similar spatial distribution of the temperature regulation, but the decrease rate was lower than that in the LMA-based estimation. For instance, the decrease rate of the temperature regulation with the area-based estimation in East Jakarta was −44.7% lower than that from the LMA-based estimation. This decreased pattern was also observed in all municipalities. Land 2021, 10, 218 18 of 27 Figure 11. Ecosystem services index (ESI) for temperature regulation in the year 1995 and 2014. (a) ESI in 1995 without landscape metric area (LMA); (b) ESI in 2014 without LMA; (c) ESI in 1995 with LMA; and (d) ESI in 2014 with LMA. Furthermore, the dynamic changes of the landscape patterns affect the spatial distri- bution of the runoff regulation, as shown in Figure 12. It shows that the area-based esti- mation indicated a decrease of the runoff regulation between 1995 and 2014 in all munic- ipalities where the highest decrease was found in West Jakarta (−64.6%), followed by North Jakarta (−40.5%), East Jakarta (−38.4%), South Jakarta (−33.6%), and Central Jakarta (−25.8%). Meanwhile, East Jakarta had the highest runoff regulation in 2014, followed by South Jakarta and North Jakarta, while the lowest runoff regulation was in Central Jakarta. The LMA-based estimation, on the other hand, showed a higher decrease rate of the runoff regulation. By considering landscape metrics in the effective area of green spaces, the LMA-based estimation suggested that the highest decrease rate was in West Jakarta (−74.6%), followed by South Jakarta (−57.9%), East Jakarta (−51.8%), North Jakarta (−46.3%), and Central Jakarta (−25.0%). Furthermore, the highest runoff regulation was in East Jakarta, followed by North, South Jakarta, West Jakarta, and Central Jakarta. Land 2021, 10, 218 19 of 27 Figure 12. ESI for runoff regulation in the year 1995 and 2014. (a) ESI in 1995 without LMA; (b) ESI in 2014 without LMA; (c) ESI in 1995 with LMA; and (d) ESI in 2014 with LMA. 4. Discussion 4.1. Landscape Pattern Changes and Ecosystem Services The LULC change analysis shows a considerable change in the spatial distribution and the landscape patterns of green spaces between 1995 and 2014, which in turn affect ecosystem services, as given by the LMA-based estimation. An increase in built-up areas results in a decrease in the size and the proportion of green spaces, as well as in the dis- aggregation and disconnection of green spaces. A change in the landscape patterns affects the capacity of green spaces to provide ecosystem services [8,45]. This study indicates that the decrease in the carbon sequestration between 1995 and 2014 was influenced by a decrease in the area of green spaces, as well as a decrease in the aggregation and the connectivity of trees. For instance, Central Jakarta was covered with about 3.2 km of trees in 2014, which is smaller than the area of trees in North Jakarta (namely 5.2 km ). However, the carbon sequestration capacity in Central Jakarta is higher than that in North Jakarta, which can be explained by differences in the landscape pat- terns. In Central Jakarta, the trees are more aggregated, as indicated by 49.6% of the AI index, while the AI of the trees in North Jakarta is 36.1%. Moreover, the trees in Central Jakarta are more connected, as seen from 78.9% of the COHESION, which is higher than 59.9% in North Jakarta. However, the influence of the connectivity and the aggregation of the trees is not always significant compared to its area. For instance, Central Jakarta has a smaller area and higher connectivity of trees, but carbon sequestration in that area is lower than that in East Jakarta, which has a higher area of trees. Furthermore, the aggregation of trees in East Jakarta is slightly lower than that in Central Jakarta and South Jakarta. The relationship between the carbon sequestration with the area and the landscape patterns of the green spaces is consistent with a previous study in Jakarta, which showed that a higher carbon sequestration rate was observed in the location of a large area of trees . Land 2021, 10, 218 20 of 27 The essential influence of the area of trees in the carbon sequestration was also observed in tropical cities such as Singapore  and in cities at high latitudes such as in Espoo, Southern Finland . Moreover, it is important to underline that the carbon sequestration discussed in this study and the examples from other cities are limited to a city or an urban scale. The urban carbon sequestration may not give a significant contribution to the global carbon sequestration process because the total carbon uptake from urban vegetation is smaller than carbon emissions in urban areas, as reported in some cities, including Mexico City and Singapore , Beijing , and Cebu City, Philippines . The significant influence of the landscape patterns on ecosystem services is also observed when comparing the temperature regulation estimates based on the area and the landscape patterns of green spaces, and bodies of water. Between 1995 and 2014, tem- perature regulation was found to decrease in all municipalities where the contribution of green spaces is more significant than that of bodies of water, which have small areas and proportions as a result of the land cover classification. As a result, the bodies of water in this study only provide a small contribution to temperature regulation as well as to other ES. Moreover, temperature regulation is more influenced by the area of green spaces than by their connectivity. For instance, the lowest temperature regulation is identified in Cen- tral Jakarta and West Jakarta, where the area of green spaces was small in 2014, while the connectivity is high. On the other hand, the highest temperature regulation in 2014 was observed in East Jakarta, where the large area and the high connectivity of green spaces are observed. In addition, the proportion of green spaces in this municipality was rela- tively higher than that in other municipalities except for West Jakarta. A positive relation- ship between the increasing urban temperature and the dynamic changes of the landscape patterns has been found in this study which is consistent with the spatial distribution of LST and air temperature in Jakarta discussed in Ramdhoni and Rushayati . In their study, the high LST (higher than 30 °C) was distributed in all municipalities. At the same time, the high air temperature was identified in Central Jakarta, East Jakarta, and North Jakarta, mainly where the small and the fragmented green spaces are observed in our study. The relationship between LST and the landscape patterns of green spaces was also observed in three megacities in Southeast Asia: Jakarta, Bangkok, and Manila . It revealed that Jakarta had the highest mean LST of the three cities considered. It was concluded that this high LST was caused by the relatively large, complex, and more aggregated impervious built-up areas of Jakarta. Furthermore, the impact of the decreasing aggregation and con- nectivity of green spaces was also identified in the urban center of the Olympic Forest Park in Beijing, China. The observed increase in LST correlated with an increase in the green space fragmentation and a decrease in the green space connectivity . Urbanization in Jakarta between 1995 and 2014 decreased runoff regulation in all municipalities, which could be related to a decrease in green spaces and bodies of water, and the landscape pattern changes. The changes in the runoff regulation are mainly influ- enced by the dynamic changes of the area and the landscape patterns of green spaces, while the influence of bodies of water is relatively small. The most significant decrease in runoff regulation is observed in West Jakarta, where there has been a high decrease in the area and the proportion of green spaces, which mean West Jakarta has the lowest runoff regulation. Furthermore, in all municipalities the connectivity of green spaces is relatively high while the proportion of green spaces is relatively low and more disaggregated, as shown by a decrease in the aggregation index. The low surface runoff regulation in North Jakarta and West Jakarta is consistent with the area inundated by the 2013 flood . Moreover, the capacity of green spaces to control surface runoff can also be influenced by local physical characteristics such as the slope of the land. However, these characteristics have not been addressed in this study. Increasing the area of green spaces for surface run- off control can be considered in the gentlest sloping area with the high topographic wet- ness index (TWI), while this is not an appropriate solution in an area with the steep slope and the low TWI . The capacity of green spaces to minimize surface runoff is also as- sociated with LPI and AI. In Beijing, for instance, a decrease of LPI and AI reduced the Land 2021, 10, 218 21 of 27 runoff reduction capacity from 23% in 2000 to 17% in 2010 . The low connectivity and aggregation could be associated with a decrease in the capacity of green spaces to mini- mize surface runoff. A study from Texas, US, showed that large area, less fragmentation, and high connectivity of green spaces could play an important role in anticipating the mean annual peak runoff. On the other hand, urban development could increase the built- up land cover class, which has more edges, more complex shapes, and is more connected, which facilitates high runoff . 4.2. Limitations of This Study and Future Perspectives The spatial analysis in this study used a medium satellite image resolution of 30 m from Landsat 5 and 8 to analyze three ecosystem services: carbon sequestration, temper- ature regulation, and runoff regulation. The Landsat products provide a cost-effective op- tion for spatial analysis purposes, such as the LULC change analysis and the landscape metrics analysis. However, some specific urban details such as city parks, small wetland areas, and lakes, as well as street corridors, could not be captured properly. Meanwhile, this study’s results have not yet been validated, which could limit the applicability of the proposed method. Even though it indicates a more realistic ES calculation, the proposed method could potentially bias the results due to different spatial scales used in this study and field survey data. In future studies, the use of high spatial resolution satellite images for land use–land cover classification is recommended. This increases the spatial analysis quality, and it al- lows researchers to include more aspects of ecosystem services analysis, such as air puri- fication from trees on street corridors or tourism and educational ecosystem services from the city and the neighborhood parks. The proposed ES estimation could be further im- proved in the future by incorporating empirical values based on runoff coefficients de- rived from land cover classes. For validation purposes, future studies should consider methods which apply remote sensing techniques using a range of spatial resolutions com- bined with actual land management information to identify carbon sequestration rates, land surface temperatures, and flood-prone areas. 5. Conclusions The present study shows that the landscape pattern changes in the period 1995–2014 have a significant impact on urban ES in Jakarta. The observed landscape pattern changes reduce the carbon sequestration, the temperature regulation, and the runoff regulation with an average value of 66, 55, and 51%, respectively, during that period. The impact of the landscape pattern changes on a decrease of ES is higher compared to the area-based calculation in which the average decrease of the carbon sequestration is 55%, followed by the temperature regulation (44%) and the runoff regulation (41%). This indicates that the ability of green spaces to provide ES is strongly influenced by the spatial characteristics or the landscape patterns of green spaces. From 1995 to 2014, it is observed that urbanization has not only reduced the number and total surface area of green spaces but also has changed its spatial characteristics. After 20 years, an increase in built-up areas has affected the spatial characteristics of green spaces marked by a scattered distribution, and a decrease in the aggregation and the con- nectivity of green spaces, making green spaces more fragmented. The average decrease in aggregation is about 8.7 and 9.6% for its connectivity, which can be associated with a de- crease in urban ES. The reduced number and area of green spaces along with the land- scape pattern changes and a decrease in urban ES indicate that there is a relationship be- tween those factors, which need to be considered as part of the ES calculation. This study proposes a new method for ES calculation that explicitly acknowledges the landscape patterns of green spaces resulting in the landscape metric area (LMA) of an urban area. The LMA reflects the area which is defined as green space that provides an effective contribution to ES. The spatial characteristics of green spaces are not static, but Land 2021, 10, 218 22 of 27 change over time. This study has shown that urbanization is a factor that impacts the sepa- ration and disaggregation as well as the disconnection of green spaces—and thus to a large extent the spatial distribution of green spaces. Ignoring the influence of the landscape pat- terns may result in an overestimation of the ES of an urban area, as shown in this study. The importance of the area of green spaces to ES provision has been discussed in many previous studies, which feed the general perception that there is a positive relation- ship between the area of green spaces and ES. The findings of our study contribute to the existing body of knowledge in this field and suggest that ES calculations should consider not only the area of green spaces but also the landscape patterns, which represent the spatial configuration of green spaces. By considering the landscape patterns in this study, the proposed method indicates a more realistic ES calculation compared to commonly area-based approaches in which green spaces can be perceived as a lumped area. Our study is one of the first to provide quantified and convincing evidence of the influence of the landscape pattern changes on urban ecosystem services in a megacity such as Jakarta. We acknowledge that further detailed studies will be needed and that the present study will provide incentives and directions to initiate further research. It is recommended to further improve this new method by using land cover classification from a high satellite image resolution, considering more landscape metric combinations, incorporating empir- ical values such as runoff coefficients or albedo values based on LULC, using field survey data or using model simulations to verify the ES calculation. Author Contributions: Conceptualization, D.M.; methodology, D.M.; formal analysis, D.M.; inves- tigation, D.M.; writing—original draft preparation, D.M., A.P. and C.Z.; writing—review and edit- ing, D.M., A.P. and C.Z.; visualization, D.M. All authors have read and agreed to the published version of the manuscript. Funding: This research was financially supported by the Indonesia Endowment Fund for Education (BUDI-LPDP) for the doctoral scholarship for Dikman Maheng. The Laminga fund of the Delft Uni- versity of Technology for providing the travel grant to present a part of this research at the IWA Water and Development Congress and Exhibition 2019. Institutional Review Board Statement: Not applicable. Informed Consent Statement: Not Applicable. Data Availability Statement: The data presented in this study are available on request from the corresponding author. Acknowledgments: The authors acknowledge the use of data from Geospatial Information Agency of Indonesia (BIG), and Landsat 5 and Landsat 8 images from U.S. Geological Survey. The authors appreciate all anonymous reviewers who provide valuable and constructive comments which have improved the quality of the paper. The authors also thank Mohanasundar Radakrishnan for the English correction of the first draft of the paper. Conflicts of Interest: The authors declare no conflict of interest. References 1. United Nations Trends in Urbanization. In World Urbanization Prospects; United Nations: New York, 2014; Vol. 12, pp. 7–12 ISBN 9789211515176. 2. Dou, Y.; Kuang, W. A comparative analysis of urban impervious surface and green space and their dynamics among 318 different size cities in China in the past 25 years. Sci. Total Environ. 2020, 706, 135828, doi:10.1016/j.scitotenv.2019.135828. 3. Wu, Y.; Tao, Y.; Yang, G.; Ou, W.; Pueppke, S.; Sun, X.; Chen, G.; Tao, Q. Impact of land use change on multiple ecosystem services in the rapidly urbanizing Kunshan City of China: Past trajectories and future projections. Land use policy 2019, 85, 419–427, doi:10.1016/j.landusepol.2019.04.022. 4. Xu, Q.; Dong, Y. xiang; Yang, R. Influence of land urbanization on carbon sequestration of urban vegetation: A temporal cooperativity analysis in Guangzhou as an example. Sci. Total Environ. 2018, 635, 26–34, doi:10.1016/j.scitotenv.2018.04.057. Land 2021, 10, 218 23 of 27 5. Ye, Y.; Bryan, B.A.; Zhang, J.; Connor, J.D.; Chen, L.; Qin, Z.; He, M. Changes in land-use and ecosystem services in the Guangzhou-Foshan Metropolitan Area, China from 1990 to 2010: Implications for sustainability under rapid urbanization. Ecol. Indic. 2018, 93, 930–941, doi:10.1016/j.ecolind.2018.05.031. 6. Alcamo, J.; Ash, N.J.; Butler, C.D.; Callicot, J.B.; Capistrano, D.; Carpenter, S.R.; Castilla, J.C.; Chambers, R.; Chopta, K.; Cropper, A.; et al. Ecosystem and Their Services. In Ecosystems and human well-being: A framework for assessment; Reid, W. V., Ed.; Island Press: Washington, D.C., USA, 2003; pp. 49–70 ISBN 1559634022. 7. Bolund, P.; Hunhammar, S. Ecosystem services in urban areas. Ecol. Econ. 1999, 29, 293–301, doi:10.1016/S0921-8009(99)00013- 8. Burkhard, B.; Kroll, F.; Nedkov, S.; Mü ller, F. Mapping ecosystem service supply, demand and budgets. Ecol. Indic. 2012, 21, 17–29, doi:10.1016/j.ecolind.2011.06.019. 9. Derkzen, M.L.; van Teeffelen, A.J.A.; Verburg, P.H. Quantifying urban ecosystem services based on high-resolution data of urban green space: An assessment for Rotterdam, the Netherlands. J. Appl. Ecol. 2015, 52, 1020–1032, doi:10.1111/1365- 2664.12469. 10. Gkatsopoulos, P. A methodology for calculating cooling from vegetation evapotranspiration for use in urban space microclimate simulations. Procedia Environ. Sci. 2017, 38, 477–484, doi:10.1016/j.proenv.2017.03.139. 11. Gómez-Baggethun, E.; Barton, D.N. Classifying and valuing ecosystem services for urban planning. Ecol. Econ. 2013, 86, 235– 245, doi:10.1016/j.ecolecon.2012.08.019. 12. Tratalos, J.; Fuller, R.A.; Warren, P.H.; Davies, R.G.; Gaston, K.J. Urban form, biodiversity potential and ecosystem services. Landsc. Urban Plan. 2007, 83, 308–317, doi:10.1016/j.landurbplan.2007.05.003. 13. Van Oudenhoven, A.P.E.; Petz, K.; Alkemade, R.; Hein, L.; De Groot, R.S. Framework for systematic indicator selection to assess effects of land management on ecosystem services. Ecol. Indic. 2012, 21, 110–122, doi:10.1016/j.ecolind.2012.01.012. 14. Panagopoulos, T.; González Duque, J.A.; Bostenaru Dan, M. Urban planning with respect to environmental quality and human well-being. Environ. Pollut. 2016, doi:10.1016/j.envpol.2015.07.038. 15. Loures, L.; Santos, R.; Panagopoulos, T. Urban parks and sustainable city planning - The case of Portimão, Portugal. WSEAS Trans. Environ. Dev. 2007. 16. Clerici, N.; Cote-Navarro, F.; Escobedo, F.J.; Rubiano, K.; Villegas, J.C. Spatio-temporal and cumulative effects of land use- land cover and climate change on two ecosystem services in the Colombian Andes. Sci. Total Environ. 2019, 685, 1181–1192, doi:10.1016/j.scitotenv.2019.06.275. 17. Depietri, Y.; Renaud, F.G.; Kallis, G. Heat waves and floods in urban areas: A policy-oriented review of ecosystem services. Sustain. Sci. 2012, 7, 95–107, doi:10.1007/s11625-011-0142-4. 18. Dupras, J.; Marull, J.; Parcerisas, L.; Coll, F.; Gonzalez, A.; Girard, M.; Tello, E. The impacts of urban sprawl on ecological connectivity in the Montreal Metropolitan Region. Environ. Sci. Policy 2016, 58, 61–73, doi:10.1016/j.envsci.2016.01.005. 19. Fox, D.M.; Witz, E.; Blanc, V.; Soulié, C.; Penalver-Navarro, M.; Dervieux, A. A case study of land cover change (1950-2003) and runoff in a Mediterranean catchment. Appl. Geogr. 2012, 32, 810–821, doi:10.1016/j.apgeog.2011.07.007. 20. Farrugia, S.; Hudson, M.D.; McCulloch, L. An evaluation of flood control and urban cooling ecosystem services delivered by urban green infrastructure. Int. J. Biodivers. Sci. Ecosyst. Serv. Manag. 2013, 9, 136–145, doi:10.1080/21513732.2013.782342. 21. Marando, F.; Salvatori, E.; Sebastiani, A.; Fusaro, L.; Manes, F. Regulating ecosystem services and green infrastructure: Assessment of urban heat island effect mitigation in the municipality of Rome, Italy. Ecol. Modell. 2019, 392, 92–102, doi:10.1016/j.ecolmodel.2018.11.011. 22. Yuan, Y.; Wu, S.; Yu, Y.; Tong, G.; Mo, L.; Yan, D.; Li, F. Spatiotemporal interaction between ecosystem services and urbanization: Case study of Nanjing City, China. Ecol. Indic. 2018, 95, 917–929, doi:10.1016/j.ecolind.2018.07.056. 23. Jaligot, R.; Kemajou, A.; Chenal, J. Cultural ecosystem services provision in response to urbanization in Cameroon. Land use Land 2021, 10, 218 24 of 27 policy 2018, 79, 641–649, doi:10.1016/j.landusepol.2018.09.013. 24. Das, M.; Das, A. Dynamics of urbanization and its impact on urban ecosystem services (UESs): A study of a medium size town of West Bengal, Eastern India. J. Urban Manag. 2019, 1–15, doi:10.1016/j.jum.2019.03.002. 25. Costanza, R.; D’Arge, R. The value of the world’s ecosystem services and natural capital. LK - https://royalroads.on.worldcat.org/oclc/4592801201. Nat. TA - TT - 1997, 387, 253–260. 26. Wang, J.; Zhou, W.; Pickett, S.T.A.; Yu, W.; Li, W. A multiscale analysis of urbanization effects on ecosystem services supply in an urban megaregion. Sci. Total Environ. 2019, 662, 824–833, doi:10.1016/j.scitotenv.2019.01.260. 27. Asadolahi, Z.; Salmanmahiny, A.; Sakieh, Y.; Mirkarimi, S.H.; Baral, H.; Azimi, M. Dynamic trade-off analysis of multiple ecosystem services under land-use change scenarios: Towards putting ecosystem services into planning in Iran. Ecol. Complex. 2018, 36, 250–260, doi:10.1016/j.ecocom.2018.09.003. 28. Haas, J.; Furberg, D.; Ban, Y. Satellite monitoring of urbanization and environmental impacts—A comparison of Stockholm and Shanghai. Int. J. Appl. Earth Obs. Geoinf. 2015, 38, 138–149, doi:10.1016/j.jag.2014.12.008. 29. Bai, T.; Mayer, A.L.; Shuster, W.D.; Tian, G. The hydrologic role of urban green space in mitigating flooding (Luohe, China). Sustain. 2018, 10, 1–13, doi:10.3390/su10103584. 30. Zhang, B.; Xie, G. di; Li, N.; Wang, S. Effect of urban green space changes on the role of rainwater runoff reduction in Beijing, China. Landsc. Urban Plan. 2015, 140, 8–16, doi:10.1016/j.landurbplan.2015.03.014. 31. Amani-Beni, M.; Zhang, B.; Xie, G. Di; Shi, Y. Impacts of urban green landscape patterns on land surface temperature: Evidence from the adjacent area of Olympic Forest Park of Beijing, China. Sustain. 2019, 11, 1–16, doi:10.3390/su11020513. 32. Jaganmohan, M.; Knapp, S.; Buchmann, C.M.; Schwarz, N. The bigger, the better? The influence of urban green space design on cooling effects for residential areas. J. Environ. Qual. 2016, 45, 134–145, doi:10.2134/jeq2015.01.0062. 33. Mitchell, M.G.E.; Suarez-Castro, A.F.; Martinez-Harms, M.; Maron, M.; McAlpine, C.; Gaston, K.J.; Johansen, K.; Rhodes, J.R. Reframing landscape fragmentation’s effects on ecosystem services. Trends Ecol. Evol. 2015, 30, 190–198, doi:10.1016/j.tree.2015.01.011. 34. Mitchell, M.G.E.; Bennett, E.M.; Gonzalez, A. Linking landscape connectivity and ecosystem service provision: Current knowledge and research gaps. Ecosystems 2013, 16, 894–908, doi:10.1007/s10021-013-9647-2. 35. Nagasawa, R.; Fukushima, A.; Yayusman, L.F.; Novresiandi, D.A. Urban expansion and its influences on the suburban land use change in Jakarta metropolitan region (JABODETABEK). Urban Plan. Des. Res. 2015, 3, 7–16, doi:10.14355/updr.2015.03.002. 36. Pravitasari, A.E. Study on impact of urbanization and rapid urban expansion in Java and Jabodetabek, megacity in Indonesia, Kyoto University, Japan, 2015. 37. Ramdhoni, S.; Rushayati, S.B. Open green space development priority based on distribution of air temperature change in capital city of Indonesia, Jakarta. Procedia Environ. Sci. 2016, 33, 204–213, doi:10.1016/j.proenv.2016.03.071. 38. Rustiadi, E.; Zain, A.M.; Trisasongko, B.H.; Carolita, I. Land cover change in Jabotabek region; Himiyama, Y., Mather, A., Bicik, I., Milanova, E. V., Eds.; II.; International Geographical Union Commission on Land Use/Cover Change: Hokkaido, 2002; 39. Firman, T. The continuity and change in mega-urbanization in Indonesia: A survey of Jakarta-Bandung Region (JBR) development. Habitat Int. 2009, 33, 327–339, doi:10.1016/j.habitatint.2008.08.005. 40. BPS DKI Jakarta Jakarta in Figures 2016; BPS Provinsi DKI Jakarta/Jakarta Statistics Bureau: Jakarta, 2016; 41. De Ridder, K.; Adamec, V.; Bañuelos, A.; Bruse, M.; Bü rger, M.; Damsgaard, O.; Dufek, J.; Hirsch, J.; Lefebre, F.; Pérez- Lacorzana, J.M.; et al. An integrated methodology to assess the benefits of urban green space. Sci. Total Environ. 2004, 334– 335, 489–497, doi:10.1016/j.scitotenv.2004.04.054. 42. Leroux, L.; Congedo, L.; Bellón, B.; Gaetano, R.; Bégué, A. Land cover mapping using Sentinel-2 images and the semi- automatic classification plugin: A Northern Burkina Faso case study. In QGIS and applications in agriculture and forest; 2018 Land 2021, 10, 218 25 of 27 ISBN 9781119457107. 43. Arowolo, A.O.; Deng, X.; Olatunji, O.A.; Obayelu, A.E. Assessing changes in the value of ecosystem services in response to land-use/land-cover dynamics in Nigeria. Sci. Total Environ. 2018, 636, 597–609, doi:10.1016/j.scitotenv.2018.04.277. 44. Gökyer, E. Understanding landscape structure using landscape metrics. In Advances in Landscape Architecture; Ö zyavuz, M., Ed.; IntechOpen: London, 2013. 45. Mugiraneza, T.; Ban, Y.; Haas, J. Urban land cover dynamics and their impact on ecosystem services in Kigali, Rwanda using multi-temporal Landsat data. Remote Sens. Appl. Soc. Environ. 2019, 13, 234–246, doi:10.1016/j.rsase.2018.11.001. 46. Daniels, B.; Zaunbrecher, B.S.; Paas, B.; Ottermanns, R.; Ziefle, M.; Roß-Nickoll, M. Assessment of urban green space structures and their quality from a multidimensional perspective. Sci. Total Environ. 2018, 615, 1364–1378, doi:10.1016/j.scitotenv.2017.09.167. 47. Estoque, R.C.; Murayama, Y. Landscape pattern and ecosystem service value changes: Implications for environmental sustainability planning for the rapidly urbanizing summer capital of the Philippines. Landsc. Urban Plan. 2013, 116, 60–72, doi:10.1016/j.landurbplan.2013.04.008. 48. McGarigal, K.; Cushman, S.; Ene., E. FRAGSTATS v4: Spatial pattern analysis program for categorical and continuous maps. Computer software program, University of Massachusetts, Amherst, MA, USA Available online: http://www.umass.edu/landeco/research/fragstats/fragstats.html (accessed on Sep 30, 2019). 49. De Smith, M.J.; Goodchild, M.F.; Longley, P. Geospatial analysis: a comprehensive guide to principles, techniques and software tools . Troubador Publishing Ltd.; Troubador Publishing Ltd: Leicester, 2007; ISBN 1-905886-60-8. 50. Kong, F.; Yin, H.; Nakagoshi, N.; Zong, Y. Urban green space network development for biodiversity conservation: Identification based on graph theory and gravity modeling. Landsc. Urban Plan. 2010, 95, 16–27, doi:10.1016/j.landurbplan.2009.11.001. 51. He, H.S.; DeZonia, B.E.; Mladenoff, D.J. An aggregation index (AI) to quantify spatial patterns of landscapes. Landsc. Ecol. 2000, doi:10.1023/A:1008102521322. 52. Bao, T.; Li, X.; Zhang, J.; Zhang, Y.; Tian, S. Assessing the distribution of urban green spaces and its anisotropic cooling distance on urban heat island pattern in Baotou, China. ISPRS Int. J. Geo-Information 2016, 5, 12, doi:10.3390/ijgi5020012. 53. Chen, J.; Goh, J. Carbon Accounting in Local-Scale Land Use and Land Cover Change. Cons. J. Sustaianble Dev. 2017, 17, 46– 54. Estoque, R.C.; Murayama, Y.; Myint, S.W. Effects of landscape composition and pattern on land surface temperature: An urban heat island study in the megacities of Southeast Asia. Sci. Total Environ. 2017, 577, 349–359, doi:10.1016/j.scitotenv.2016.10.195. 55. Li, X.; Zhou, W.; Ouyang, Z.; Xu, W.; Zheng, H. Spatial pattern of greenspace affects land surface temperature: Evidence from the heavily urbanized Beijing metropolitan area, China. Landsc. Ecol. 2012, 27, 887–898, doi:10.1007/s10980-012-9731-6. 56. Ramesh, R.; Chen, Z.; Cummins, V.; Day, J.; D’Elia, C.; Dennison, B.; Forbes, D.L.; Glaeser, B.; Glaser, M.; Glavovic, B.; et al. Land-Ocean Interactions in the Coastal Zone: Past, present & future. Anthropocene 2015, 12, 85–98, doi:10.1016/j.ancene.2016.01.005. 57. Werner, C. Green open spaces in Indonesian cities : Schisms between law and practice. Pacific Geogr. 2014, 41, 26–31. 58. Cortinovis, C.; Geneletti, D. Ecosystem services in urban plans: What is there, and what is still needed for better decisions. Land use policy 2018, 70, 298–312, doi:10.1016/j.landusepol.2017.10.017. 59. Vincent, S.U.; Radhakrishnan, M.; Hayde, L.; Pathirana, A. Enhancing the economic value of large investments in Sustainable Drainage Systems (SuDS) through inclusion of ecosystems services benefits. Water (Switzerland) 2017, 9, doi:10.3390/w9110841. 60. Rotterdam, C. of Rotterdam programme on sustainability and climate change 2015-2018; City of Rotterdam: Rotterdam, the Land 2021, 10, 218 26 of 27 Netherlands, 2015; 61. Oslo, C. of Urban ecology programme 2011-2026; City of Oslo, Department of Environmental Affairs and Transport: Oslo, 2011; 62. Sydney, C. of Environmental action 2016 – 2021: Strategy and action plan; City of Sydney: Sydney, 2017; 63. Arifin, H.S.; Nakagoshi, N. Landscape ecology and urban biodiversity in tropical Indonesian cities. Landsc. Ecol. Eng. 2011, 7, 33–43, doi:10.1007/s11355-010-0145-9. 64. Rusiawan, W.; Tjiptoherijanto, P.; Suganda, E.; Darmajanti, L. System dynamics modeling for urban economic growth and CO2 emission: A case study of Jakarta, Indonesia. Procedia Environ. Sci. 2015, 28, 330–340, doi:10.1016/j.proenv.2015.07.042. 65. Surahman, U.; Kubota, T.; Wijaya, A. Life cycle assessment of energy and CO2 emissions for residential buildings in Jakarta, Indonesia. IOP Conf. Ser. Mater. Sci. Eng. 2016, 128, doi:10.1088/1757-899X/128/1/012002. 66. Siswanto; Van Oldenborgh, G.J.; Van Der Schrier, G.; Lenderink, G.; Van Den Hurk, B. Trends in high-daily precipitation events in Jakarta and the flooding of January 2014. Bull. Am. Meteorol. Soc. 2015, 96, S131–S135, doi:10.1175/BAMS-D-15- 00128.1. 67. Wijayanti, P.; Zhu, X.; Hellegers, P.; Budiyono, Y.; van Ierland, E.C. Estimation of river flood damages in Jakarta, Indonesia. Nat. Hazards 2017, 86, 1059–1079, doi:10.1007/s11069-016-2730-1. 68. Darmanto, N.S.; Varquez, A.C.G.; Kawano, N.; Kanda, M. Future urban climate projection in a tropical megacity based on global climate change and local urbanization scenarios. Urban Clim. 2019, 29, 100482, doi:10.1016/j.uclim.2019.100482. 69. Churkina, G. Carbon cycle of urban ecosystems. In Carbon Sequestration in Urban Ecosystems; Lal, R., Augustin, B., Eds.; Springer Science+Business Media B.V., 2011; p. 315 ISBN 978-94-007-2365-8. 70. Kuittinen, M.; Moinel, C.; Adalgeirsdottir, K. Carbon sequestration through urban ecosystem services: A case study from Finland. Sci. Total Environ. 2016, 563–564, 623–632, doi:10.1016/j.scitotenv.2016.03.168. 71. Nero, B.F.; Callo-Concha, D.; Anning, A.; Denich, M. Urban green spaces enhance climate change mitigation in cities of the Global South: The case of Kumasi, Ghana. Procedia Eng. 2017, 198, 69–83, doi:10.1016/j.proeng.2017.07.074. 72. Forman, R.T.T. Urban Ecology; Cambridge University Press: Cambridge, 2013; ISBN 9781139030472. 73. Azaria, L.; Wibowo, A.; Putut Ash Shidiq, I.; Rokhmatuloh Carbon sequestration capability analysis of urban green space using geospatial data. E3S Web Conf. 2018, 73, 4–8, doi:10.1051/e3sconf/20187303009. 74. Oviantari, M. V.; Gunamantha, I.M.; Ristiati, N.P.; Santiasa, I.M.P.A.; Astariani, P.P.Y. Carbon sequestration by above-ground biomass in urban green spaces in Singaraja city. IOP Conf. Ser. Earth Environ. Sci. 2018, 200, 0–6, doi:10.1088/1755- 1315/200/1/012030. 75. Oke, T.R.; Mills, G.; Christen, A.; Voogt, J.A. Urban Climates; Cambridge University Press: Cambridge, 2017; ISBN 76. Maheng, D.; Ducton, I.; Lauwaet, D.; Zevenbergen, C.; Pathirana, A. The sensitivity of urban heat island to urban green space- A model-based study of City of Colombo, Sri Lanka. Atmosphere (Basel). 2019, 10, doi:10.3390/atmos10030151. 77. Sodoudi, S.; Zhang, H.; Chi, X.; Mü ller, F.; Li, H. The influence of spatial configuration of green areas on microclimate and thermal comfort. Urban For. Urban Green. 2018, 34, 85–96, doi:10.1016/j.ufug.2018.06.002. 78. Amani-Beni, M.; Zhang, B.; Xie, G. di; Xu, J. Impact of urban park’s tree, grass and waterbody on microclimate in hot summer days: A case study of Olympic Park in Beijing, China. Urban For. Urban Green. 2018, 32, 1–6, doi:10.1016/j.ufug.2018.03.016. 79. Steeneveld, G.J.; Koopmans, S.; Heusinkveld, B.G.; Theeuwes, N.E. Refreshing the role of open water surfaces on mitigating the maximum urban heat island effect. Landsc. Urban Plan. 2014, 121, 92–96, doi:10.1016/j.landurbplan.2013.09.001. 80. Theeuwes, N.E.; Solcerová, A.; Steeneveld, G.J. Modeling the influence of open water surfaces on the summertime temperature and thermal comfort in the city. J. Geophys. Res. Atmos. 2013, 118, 8881–8896, doi:10.1002/jgrd.50704. 81. Chen, Y.C.; Tan, C.H.; Wei, C.; Su, Z.W. Cooling effect of rivers on metropolitan Taipei using remote sensing. Int. J. Environ. Res. Public Health 2014, 11, 1195–1210, doi:10.3390/ijerph110201195. Land 2021, 10, 218 27 of 27 82. Gunawardena, K.R.; Wells, M.J.; Kershaw, T. Utilising green and bluespace to mitigate urban heat island intensity. Sci. Total Environ. 2017, 584–585, 1040–1055, doi:10.1016/j.scitotenv.2017.01.158. 83. Alves, A.; Vojinovic, Z.; Kapelan, Z.; Sanchez, A.; Gersonius, B. Exploring trade-offs among the multiple benefits of green- blue-grey infrastructure for urban flood mitigation. Sci. Total Environ. 2020, 703, 134980, doi:10.1016/j.scitotenv.2019.134980. 84. Yang, B.; Lee, D.K.; Heo, H.K.; Biging, G. The effects of tree characteristics on rainfall interception in urban areas. Landsc. Ecol. Eng. 2019, 15, 289–296, doi:10.1007/s11355-019-00383-w. 85. Yao, L.; Chen, L.; Wei, W.; Sun, R. Potential reduction in urban runoff by green spaces in Beijing: A scenario analysis. Urban For. Urban Green. 2015, 14, 300–308, doi:10.1016/j.ufug.2015.02.014. 86. ASCE Hydrology Handbook; Second Edi.; American Society of Civil Engineers: New York, 1996; ISBN 0-7844-013 8-1. 87. Alves, P.L.; Formiga, K.T.M.; Traldi, M.A.B. Rainfall interception capacity of tree species used in urban afforestation. Urban Ecosyst. 2018, 21, 697–706, doi:10.1007/s11252-018-0753-y. 88. Highfield, W.E. Section 404 permitting in coastal texas: A longitudinal analysis of the relationship between peak streamflow and wetland alteration. Environ. Manage. 2011, 49, 892–901, doi:10.1007/s00267-012-9832-7. 89. Te Chow, V.; Maidment, D.R.; Mays, L.W. Applied Hydrology; McGraw-Hill series in water resources and environmental engineering; McGraw-Hill: Singapore, 1988; ISBN 0070108102. 90. Yang, X.; You, X.Y.; Ji, M.; Nima, C. Influence factors and prediction of stormwater runoff of urban green space in Tianjin, China: Laboratory experiment and quantitative theory model. Water Sci. Technol. 2013, 67, 869–876, doi:10.2166/wst.2012.600. 91. Shang, H.; Zhang, K.; Wang, Z.; Yang, J.; He, M.; Pan, X.; Fang, C. Effect of varying wheatgrass density on resistance to overland flow. J. Hydrol. 2020, doi:10.1016/j.jhydrol.2020.125594. 92. Velasco, E.; Roth, M.; Norford, L.; Molina, L.T. Does urban vegetation enhance carbon sequestration? Landsc. Urban Plan. 2016, 148, 99–107, doi:10.1016/j.landurbplan.2015.12.003. 93. Tang, Y.; Chen, A.; Zhao, S. Carbon storage and sequestration of urban street trees in Beijing, China. Front. Ecol. Evol. 2016, 4, 1–8, doi:10.3389/fevo.2016.00053. 94. Pansit, N.R. Carbon storage and sequestration potential of urban trees in Cebu City, Philippines. Mindanao J. Sci. Technol. 2019, 17, 98–111. 95. Padawangi, R.; Douglass, M. Water, water everywhere: Toward participatory solutions to chronic urban flooding in Jakarta; 2015; Vol. 88; ISBN 2015883517. 96. Kim, H.; Lee, D.K.; Sung, S. Effect of urban green spaces and flooded area type on flooding probability. Sustain. 2016, 8, doi:10.3390/su8020134. 97. Kim, H.W.; Park, Y. Urban green infrastructure and local flooding: The impact of landscape patterns on peak runoff in four Texas MSAs. Appl. Geogr. 2016, 77, 72–81, doi:10.1016/j.apgeog.2016.10.008.
Land – Multidisciplinary Digital Publishing Institute
Published: Feb 21, 2021
Access the full text.
Sign up today, get DeepDyve free for 14 days.