Next Article in Journal
Interpreting the Fuzzy Semantics of Natural-Language Spatial Relation Terms with the Fuzzy Random Forest Algorithm
Previous Article in Journal
Examining the Association of Economic Development with Intercity Multimodal Transport Demand in China: A Focus on Spatial Autoregressive Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Changes in Gross Primary Production (GPP) over the Past Two Decades Due to Land Use Conversion in a Tourism City

by
I Wayan Nuarsa
1,2,
Abd. Rahman As-syakur
2,3,*,
I Gusti Alit Gunadi
1,2 and
I Made Sukewijaya
1
1
Faculty of Agriculture, Udayana University, Kampus Bukit Jimbaran, Bali 80361, Indonesia
2
Center for Environmental Research (PPLH), Udayana University, PB Sudirman Street, Denpasar, Bali 80232, Indonesia
3
Center for Remote Sensing and Ocean Science (CReSOS), Udayana University, Postgraduate Building, Floor 3, PB Sudirman Street, Denpasar, Bali 80232, Indonesia
*
Author to whom correspondence should be addressed.
ISPRS Int. J. Geo-Inf. 2018, 7(2), 57; https://doi.org/10.3390/ijgi7020057
Submission received: 8 November 2017 / Revised: 2 February 2018 / Accepted: 5 February 2018 / Published: 7 February 2018

Abstract

:
Understanding the changes in gross primary production (GPP), which is the total carbon fixation by terrestrial ecosystems through vegetation photosynthesis, due to land use conversion in a tourism city is important for carbon cycle studies. Satellite data from Landsat 5, Landsat 7 and Landsat 8 and meteorological data are used to calculate annual GPP for 1995, 2003 and 2014, respectively, using the vegetation production model (VPM) in the tourism city Denpasar, Bali, Indonesia. Five land use types generated from topographic maps in three different years over the past two decades are used to quantify the impacts of land use changes on GPP estimation values. Analysis was performed for two periods to determine changes in land use and GPP value as well as their speed. The results demonstrated that urban land development, namely, the increase of settlement areas due to tourism activity, had overall negative effects on terrestrial GPP. The total GPP of the whole area decreased by 7793.96 tC year−1 (12.65%) during the study period. The decline is due to the conversion of agriculture and grassland area into settlements, which caused the city to lose half of its ability to uptake carbon through vegetation. However, although forest area is declining, forest maintenance and restoration by making them protection areas has been helpful in preventing a drastic decline in GPP value over the past two decades. This study provides information that is useful for carbon resource management, tourism, policy making and scholars concerned about carbon reduction in a tourism city.

1. Introduction

Around the globe, urbanization is expected to increase significantly in the coming decades and is projected to continue [1,2]. Rapid urbanization and urban sprawl have a significant impact on conditions of urban land use changes due to the increasing demands for the construction of new residential, commercial, utility, and transport infrastructures [3,4]. Urbanization is a process related to population increase in urban regions that is associated with settlement increases and automatically reduces vegetation areas due to limited land [5]. Population growth and economic changes are two important factors influencing land use changes and the availability of agricultural land, especially in urban areas [6]. One of the drivers of economic and population growth is tourism activities [7,8]. Tourism has a strong impact on land use, both at particular sites and in adjoining zones and travel corridors [9]. In Bali, Indonesia, an increase in the tourism industry and the population is creating huge social and environmental problems, and the main threats to irrigated agriculture and the related livelihoods are land conversion and the transfer of water from agriculture to non-agricultural uses [10].
Most land use changes in Indonesia at the local and regional scales are conversions from agricultural land into residential land and from forest to agriculture [11,12,13,14]. Land use change can have a large impact on ecosystems, altering their composition and structure as well as their function, including matter and energy cycles [15]. Several studies have explained that land use and land cover changes in urban areas produce global consequences, including effects on the climate, e.g., [16], hydrology, e.g., [17], and biodiversity, e.g., [18], and that these changes alter ecosystem services and affect the biological carbon uptake of urban vegetation, e.g., [19]. Land use changes in urban areas play a significant role in determining local, regional, and global carbon emissions [20]. Meanwhile, urban ecosystems can also account for a significant portion of terrestrial carbon storage at both local and regional scales [21].
Human activities are highly dependent on ecosystem goods and services, especially net ecosystem CO2 exchange (NEE), the balance between gross primary production (GPP) and ecosystem respiration [5]. GPP is a key measure of carbon mass flux in carbon cycle studies, which is the total carbon fixation by terrestrial ecosystems through vegetation photosynthesis [22]. In addition, GPP can be used to quantify vegetation productivity and to understand temporal variability in productivity [23]. The quantification of the magnitude of carbon uptake by terrestrial vegetation and how it varies due to land use changes in urban areas is an important question, as future potential sequestration is influenced by both increased atmospheric CO2 and the changing climate [24], especially for tourism cities. Therefore, the estimation of GPP in an urban area is essential for the quantification of net terrestrial carbon and thus is of particular importance for global carbon cycle research. Previous studies on GPP impacts from land cover/land use changes associated with urban growth have been documented in Southeastern Michigan [19,25]. However, changes in GPP due to land use conversion as a result of urbanization in the tourism area have not been reported in previous work.
The remote sensing of vegetation GPP is an important step in analysing terrestrial carbon cycles in response to changing climate [26]. Repetitive and systematic satellite remote-sensing observations of vegetation dynamics and ecosystems allow us to characterize vegetation structure, including leaf area index (LAI) and the absorptivity of photosynthetically active radiation (APAR) [27,28]. Remote-sensing techniques are advantageous because they allow monitoring of the GPP on a regional and local scale through its relationship with the fraction of absorbed photosynthetically active radiation (fAPAR) driven by vegetation indices [29,30]. Moreover, remote sensing is the best choice for urban areas to provide a better understanding of the effects of land use change on vegetation dynamics [31,32,33]. Moderate-resolution Imaging Spectroradiometer (MODIS) satellite remote sensing data are commonly used to calculate GPP, e.g., [34]. However, for urban regions that have heterogeneous landscapes, Landsat data are more frequently used [19,35].
GPP is generally estimated using models because it is difficult to measure directly [36]. Satellite-based, light-use efficiency (LUE) models have been widely used because they rely on simple algorithms to estimate GPP [37]. The LUE model states that carbon exchange is a function of the amount of light energy absorbed by vegetation and the efficiency with which that light energy is used to fix carbon [38]. One of the LUE models that is commonly used for estimated GPP is the vegetation production model (VPM) [39]. The VPM is a light-use efficiency model that utilizes remote sensing imagery to estimate GPP based on the impacts of temperature, water stress, and phenology [40]. GPP from Landsat data based on VPM estimates has been validated by several studies that indicate good agreement with CO2 eddy flux towers [41,42]. The potential of using the VPM model to scale up GPP estimation in Indonesia has also been evaluated for CO2 flux tower sites that use MODIS data [43].
Based on these conditions, in this work, we attempt to use Landsat data and land use information from topographic maps to estimate vegetation carbon uptake by urban tourism areas and their changes over the past two decades due to land use conversion. Three types of Landsat data, including Landsat 5, Landsat 7 and Landsat 8, are used to calculate annual GPP for 1995, 2003 and 2014, respectively. Annual changes in GPP values estimated using VPM in Denpasar, an urban tourism city on Bali Island, Indonesia, were analysed and extracted for each land use type. It is expected that the results can provide information related to the impact of land use change on ecosystem services and particularly to the biological carbon uptake capability of urban vegetation.

2. Data and Methods

2.1. Research Location

Denpasar (Bali, Indonesia) is located at 08°35’31”–08°44’49” south latitude and 115°10′23″–115°16′27″ east longitude. Denpasar city has four districts: North Denpasar, South Denpasar, West Denpasar, and East Denpasar. In 1995, Denpasar had an area of 120.47 km2, which increased to 124.68 km2 because of land reclamation at Serangan Island and Benoa port in 1996, which is located in South Denpasar. A map of Denpasar city, including boundaries before and after reclamation, is presented in Figure 1. Denpasar is located between two neighbourhood regencies, namely, the Badung Regency to the north and west, the Gianyar Regency to the east, and the Badung Strait area to the south. The three districts, including Denpasar, have the highest tourism activity in Bali. In Denpasar, the centre of tourism activity is in the southeastern part, namely, the Sanur region.
The population of Denpasar reached 335,196 people in 1992 and increased to 880,600 people in 2015, an increase of more than 250%. Topographically, the study area has a flat relief with zero elevation in the south and southeast, a highest elevation of 89 m a.s.l. in the north, and an average elevation of 23 m a.s.l. In general, the rainfall patterns of the Denpasar area are influenced by monsoons, with the maximum amount of precipitation occurring during the peak of the wet season from December to February and decreasing to a minimum during the valley of the dry season from June to August [44,45]. Moreover, the annual average temperature of Denpasar is 27.6 °C, the average annual rainfall is 1803.62 mm year−1 [44], and settlement was the dominant land use in Denpasar city in 2006 followed by rice field areas [46].

2.2. Data Use

The main remote sensing data used in this study were level-1 geometrically corrected (L1G) Landsat 5, Landsat 7, and Landsat 8 data acquired on 2 February 1995, 21 March 2003, and 27 March 2014, respectively, which were obtained from the United States Geological Survey (USGS) via the EarthExplorer database (https://earthexplorer.usgs.gov). However, clouds and cloud shadows obstructed the view of the land surface during dataset acquisition in the study area. Clouds and their shadows are one of the main problems when using Landsat data in the tropics [47,48]. The brightening effect of the clouds and the darkening effect of cloud shadows create significant sources of noise in the Landsat data and cause problems in cloud detection, which is an initial step in most analyses, including error estimation for vegetation indexes [49]. Clouds and their shadows covered 1.77%, 2.87% and 0.13% of the entire area of Denpasar for Landsat 5 on 2 February 1995, Landsat 7 on 21 March 2003, and Landsat 8 on 27 March 2014. Most clouds and cloud shadows cover the agriculture and grass, settlement, and forest land use types. The area of clouds and cloud shadows for the three main Landsat datasets was determined using the automated cloud-cover assessment (ACCA) algorithm [50]. The application of the ACCA algorithm requires green, red, near-infrared (NIR), shortwave infrared (SWIR), and thermal infrared (TIR) spectral bands in surface reflectance form. To replace the data lost due to clouds and cloud shadows, Landsat data acquired adjacent to the main data are employed to eliminate the impact of pixels lost in GPP estimation. As a consideration before using these data to replace the lost data, first, the Landsat data must have the date of acquisition closest to the primary data to avoid land surface changes due to agricultural processes. Second, the land surface in the main data and the secondary data used as a patch for lost pixels must have similar spectral characteristics, especially for unchanged land surfaces such as built-up areas. Third, the Landsat data used as the missing pixel cover must be cloud- and cloud shadow-free in the area to be patched. Based on these three considerations, a Landsat 5 acquisition on 26 May 1995, a Landsat 7 acquisition on 28 March 2003, and a Landsat 8 acquisition on 11 March 2014 were used as patches for pixels lost in 1995, 2003 and 2014, respectively. However, only 1.11% of the 2.87% pixels lost in 2003 are covered by secondary data, and the Landsat 7 acquisition on 22 April 2003 was utilized to cover the nearly 1.76% missing pixels that remained.
The spectral similarity between the main data and secondary data is derived by correlating the remote sensing indices that are used to estimate GPP from VPM. The two remote sensing indices are the enhanced vegetation index (EVI) and the land surface water index (LSWI). EVI was proposed based on a feedback-based approach that incorporates both background adjustment and atmospheric resistance concepts into the normalized difference vegetation index (NDVI) [51]. The EVI employs an additional blue band for atmospheric correction and variable soil and canopy background reflectance [52] and is calculated from blue, red, and NIR bands. The LSWI, which is proposed by Xiao et al. [53], is the normalized difference between NIR and SWIR spectral bands. Because the SWIR spectral band is sensitive to vegetation water content and soil moisture, a combination of the SWIR and NIR spectral bands has been utilized to derive water-sensitive vegetation indices. For instance, when the leaf water content or soil moisture increases, SWIR absorption increases and SWIR reflectance decreases, resulting in an increase in the LSWI value [40]. The EVI and LSWI are calculated as:
EVI   = 2.5 × ρ n i r ρ r e d ρ n i r + ( 6 × ρ r e d 7.5 × ρ b l u e ) + 1
LSWI = ρ n i r ρ s w i r ρ n i r + ρ s w i r
where ρblue, ρred, ρnir and ρSWIR are the Landsat (5, 7 and 8) surface reflectance for the blue, red, NIR and SWIR bands, respectively. The relationship between the two indices from the primary and secondary Landsat datasets for the study site represented very good agreement, with correlation (r) values ranging from 0.90 to 0.96 for EVI and from 0.90 to 0.95 for LSWI, indicating that secondary data can be used as a patch. The scatter plots between the indices from the primary and secondary Landsat datasets are presented in Figure 2.
Topographic maps for 2003 and 2014 with a scale of 1:25,000, interpreted by the Japan International Cooperation Agency (JICA) [54] and Regional Planning and Development Board (Bappeda) of Denpasar city, were used to generate land use types for 2003 and 2014, respectively. However, due to the absence of land use data for 1995, the topographic maps in 1992 from the Indonesian Geospatial Information Agency (BIG, previously known as the Indonesian National Coordination for Survey and Mapping, BAKOSURTANAL) were used to produce land use types for 1995. Land use data from 1992 and 2014 were derived from a basic map that is used by stakeholders and policy makers published by the Indonesian government through central and local governments, respectively. Land use data for 2003 were interpreted by JICA based on Landsat data from 2003 that were utilized by the Public Works Department of Bali Province. Through the Indonesian Geospatial Information Agency, the Indonesian government has a regulation that when a map is used as a base map for the government, it should have an accuracy ≥85%. Therefore, all land use data used in this study have an accuracy level greater than or equal to that allowed. All land use data used in this study are in the form of vector data and are subsequently converted into raster data with the same spatial resolution as the Landsat data. To address the complexity of land use types, land use was labelled, and the clusters were combined into five land use types, similar to the classification scheme used by Zhao et al. [25], i.e., settlement, forest, agriculture and grass, others, and water (Table 1). Water and other land use types were not included for analysing changes in GPP values due to land use conversion in the current study as these cover types were assumed to be relatively constant with respect to GPP [25].

2.3. Methods

GPP is estimated by the VPM. The VPM is an LUE model based on the conceptual partitioning of light absorption by chlorophyll pigments (FPARchl) and non-photosynthetic vegetation (NPV such as branches, trunks, stems or senescent leaves), where only FPARchl is used for photosynthesis [39]. A brief description of VPM is given below:
GPP   =   ε g × FPAR chl   ×   PAR
where εg is light-use efficiency (μmol m−2 s−1, photosynthetic photon flux density, PPFD); FPARchl is the fraction of photosynthetically active radiation (PAR) absorbed by chlorophyll; and PAR is the photosynthetically active radiation (μmol photosynthetic photon flux density, PPFD) calculated as 0.45 × R (R, incoming global solar radiation). In the current study, daily PAR is summed in a year to obtain an annual PAR value; thus, annual GPP can be obtained.
Light-use efficiency (εg) is a function of the maximum light-use efficiency (ε0), temperature (Tscalar), water condition (Wscalar), and leaf phenology (Pscalar) and is defined as:
ε g = ε 0 × T s c a l a r × W s c a l a r × P s c a l a r
A maximum light-use efficiency, ε0, value of 0.40 gC mol−1 PPFD [55] was used in this study.
The parameter Tscalar represents the temperature sensitivity of photosynthesis and is calculated at each time step using the equation developed for the terrestrial ecosystem model [56]:
T s c a l a r = ( T T min ) × ( T T max ) [ ( T T min ) × ( T T max ) ] ( T T o p t ) 2
where T is the monthly average air temperature (when PAR > 1 μmol m−2 s−1) obtained from the Indonesian Meteorology, Climatology and Geophysics Agency (BMKG) and the parameters for Tmin, Tmax, and Topt for photosynthesis were 0, 48, and 28 °C, respectively, based on Vetrita et al. [43].
The Wscalar takes into account the complex impact of water stress on photosynthesis (i.e., changes in stomatal conductance, leaf water potential, etc.) caused by soil moisture and/or atmospheric water deficits [57]. Wscalar is calculated as follows [39]:
W s c a l a r = 1 + LSWI 1 + LSWI max
where LSWImax is the maximum LSWI during the growing season for each pixel.
In the VPM, Pscalar is included to account for the effect of leaf age on photosynthesis at the canopy level [58]. The calculation of Pscalar is dependent upon the longevity of leaves (deciduous, versus evergreen). For a canopy dominated by leaves with a life expectancy of one year, Pscalar is calculated in two different phases as a linear function [39]. Because urban areas have very complex and heterogeneous plant structures and always have new leaves emerging during the plant-growing season, Pscalar is set to be 1.0 in this study.
The FPARchl within the photosynthetically active period in the VPM is assumed to be a linear function of EVI, and the coefficient a is set to be 1.0 [39]:
FPAR chl   =   EVI   ×   a
PAR irradiance measurements were not routinely carried out at radiometric sites, such as standard actinometric stations, where broadband solar irradiance components are routinely measured [59]. PAR is the region of incoming global solar radiation reaching the earth surface ranging from wavelengths of 400 to 700 nanometers [60]. Therefore, PAR is assumed to be approximately half of the global solar radiation [61]. Based on the literature compiled by Escobedo et al. [62], the ratio between PAR and global solar radiation ranges from 0.42 to 0.52. For this study, a ratio of 0.45 is used to calculate PAR from global solar radiation because this factor is employed in most studies related to GPP estimation, e.g., [63]. Global solar radiation was calculated using the modified Sayigh universal formula for the Indonesian region [64] by utilizing daily data on bright sunshine hours received during the day, maximum air temperature, and relative humidity obtained from the BMKG in a year to obtain annual PAR for 1995, 2003 and 2014. Details on the global solar radiation computation method using the modified Sayigh universal formula can be found in Halawa and Sugiyatno [64].
In this study, the remote sensing indices EVI and LSWI are assumed to be the same throughout the year. Consequently, the monthly and seasonal variability of vegetation cover and water stress cannot be known throughout the year. Only PAR and Tscalar have monthly and seasonal variability, because they are calculated from daily meteorological data. However, the primary remote sensing Landsat data used in this study were taken during the end of March, which is the beginning of the transition season (March–April–May) between the northwest monsoon (rainy season) and the southeast monsoon (dry season) for the study area [45,65], except for the data for 1995, which were acquired in February. During the transition season, the water content is not too high because rainfall begins to decrease and is not too low because it is not the dry season. In Denpasar, the vegetation in settlement areas is dominated by annual and perennial plants, such as roadside trees, shade trees, and annual flowering plants. The forest land use type is also dominated by perennial plants and mangroves. Nuarsa et al. [66] showed that annual, perennial and mangrove plants in the tropics have relatively stable vegetation indices throughout the year and are unaffected by the monsoon seasons. The agriculture land use type exhibits a seasonal pattern that is influenced by monsoons due to the planting period in which the dominant plant are paddy crops [67]. However, Jayanti et al. [68] found that paddy crops in Bali, Indonesia exhibit peaked vegetation index values in late June and early December (cropping period (reproductive stage)) with the lowest values in mid-September (fallow and flooding/transplanting period) and are at an average value in March and late October (cropping period (vegetative growth stage)). Moreover, As-syakur et al. [69] also showed that the vegetation health index (VHI), which can describe variations in vegetation indexes and water stress other than drought conditions, does not have high fluctuations throughout the year in the tropics except during El Niño events. Therefore, it can be assumed that the vegetation cover and water content in the transition season can represent the average annual condition for the study area.
To understand the impact of land use changes on GPP estimation in an urban tourism area of Denpasar, GPP values are extracted for each type of land use during the study year. The definition of land use types is based on Table 1. Previously, the area of each land use type and the speed of change were calculated. The total value of GPP is obtained from the sum of the GPP values in each pixel multiplied by the pixel area in all research areas or in each land use type, while the average GPP value is obtained by dividing the total values of GPP by the number of image pixels. In this study, the speed of GPP changes was also calculated in all research areas and in each land use type.
The analysis of land use change was performed for two periods (1992–2003 and 2003–2014) with the topographic data sources previously described. The land use change can be measured using Equation (8) in hectares and Equation (9) in percent, and the equations are defined as follows:
L U C a ( h a ) = ( L U a , t L U a , t 1 )
L U C a ( % ) = L U C a ( h a ) L U a , t 1 × 100
where LUCa (ha) is changes in land use type a in hectares from time t (current year) to time t−1 (former year) and LUCa (%) is land use change in percent. LUa,t and LUa,t−1 are the total land area of land use type a in hectares at time t (current year) and time t−1 (former year). In the current study, land use expansion in each type was also calculated and can be measured using a simple equation of the land use expansion index (LUEI). The LUEI can reflect the temporal patterns of expansion or reduction in each type of land use and is defined as follows:
L U E I a = L U C a ( h a ) / L U a , t N t N t 1
where LUEIa (ha year−1) is the annual expansion index for land use type a in hectares and can be converted to percent per year (% year−1) by multiplying by 100 and can also be converted to hectares per day (ha day−1) after dividing by the number of days in the year (356). N is the total number of years from time t (current year) to time t−1 (former year). A positive LUEI value indicates expansion, and a negative value indicates reduction.

3. Results and Discussion

3.1. Land Use Changes

Table 2 lists the land use area and proportion of the total area in each land cover type in 1992, 2003 and 2014. In 1992, agriculture and grassland comprised the largest proportion of the total area, approximately 65.31% of the total area, while the water and other land use types covered less than 5% of the study area by 1992. Figure 3a shows that agriculture and grassland were located mainly in the areas surrounding settlements, except in the southeastern and southern parts of Denpasar. Interestingly, during 2003, the proportion of agriculture and grassland began to equalize with that of settlements, with settlement areas and agriculture and grass areas accounting for approximately 43.49% and 42.03% of the total area, respectively. During this period, the settlement land use type began to spread in all directions from the city centre and clustered in the southeast region. It is noteworthy that compared to the spatial pattern in 1992, agriculture and grassland became more fragmented and scattered between settlements in 2003 (Figure 3b), and a large area of forest began to establish in the southern part of Denpasar. The forest in the southern part of Denpasar is dominated by mangroves, which began to be replanted in the pond area in 1995, although the project itself began in 1992 [70]. Settlement became the most dominant land cover type in Denpasar city during 2014, accounting for 58.57% of the study area by 2014, followed by agriculture and grassland, accounting for 28.21% and 9.68% of the study area, respectively. The two other types of land use covered less than 5% of the total area. The area of settlement expanded broadly and resulted in clear visibility of the agriculture and grass land use type, especially in the southern part of the study area, which is adjacent to the centre of tourism activity (Figure 3c). Tourism centres are not only located in the southeastern part of Denpasar but also in another district neighbouring the southwestern side of Denpasar, namely, the Kuta tourism centre.
Land use has changed in the tourism area of Denpasar from 1992 to 2014. Based on the area in 1992, settlement increased by approximately 204.75% in 2014, and the forest and agriculture and grass land cover types decreased to 12.68% and 55.29% ha in 2014, respectively (Table 2 and Table 3). The fastest change in land use in Denpasar was found from 1992 to 2003, especially for settlement and agriculture and grassland. Based on the LUEI value, the settlement area increased by approximately 0.69 ha day−1, and the agriculture and grass land cover type lost an area of approximately 0.60 ha day−1. Meanwhile, forests decreased by approximately 0.04 ha day−1 during the first nine years. During the period 2003 to 2014, settlement still showed an expanding pattern with an expansion rate of 0.43 ha day−1, and agriculture and grass declined with an LUEI rate of 0.39 ha day−1. In general, the rate of change in the settlement areas in Denpasar over the past two decades increased by approximately 0.67 ha day−1, and the forest and the agriculture and grass areas decreased by approximately 0.02 ha day−1 and 0.60 ha day−1, respectively.
The characteristics of the settlement land use type in Denpasar became increasingly complex from 1992 to 2014, with the geographical variations characterized mainly by the transformation from greenspaces into buildings for residence and government activity as well as into multi-functional buildings that provide general tourism services, including accommodation, shopping, restaurants, and entertainment. The changes in land use in the Denpasar region not only show expanding settlement areas but also expanding forestland due to the conversion of water and other land use types in the southern part of Denpasar and reclamation areas in Serangan Island. The clear land use change in this tourism city is consistent with the studies of Kytzia et al. [71], Xi et al. [72], and Mao et al. [73]. However, the process of greenspaces changing into settlement areas in Denpasar is very fast compared to that in other cities in Indonesia and the world. In Jakarta, the capital city of Indonesia, the settlement area has grown to approximately 0.02 ha day−1 over the past four decades (1970 to 2008) [13], and in Beijing, China, the urban land cover from 1986 to 2001 expanded at a rate of 0.06 ha day−1 [74].
Denpasar has been the capital of Bali Province since 1958. Therefore, all activities of the provincial government and activities related to education and business are centred in Denpasar. Although Denpasar City is not a major tourist destination in Bali, being the centre of government and the presence of tourism destinations in the Sanur region have led to growth in tourism, which can have a direct impact on land use change. As shown in Figure 3a,b, most of the forest in southeastern Denpasar, which is the Sanur region, has was converted to settlement from 1992 to 2003, and most of the change was from forest to multi-functional buildings for tourism activities. Bali has been a tourism destination since before Indonesia became independent [75]. Tourism in Bali has continued to grow, especially after the construction of the Bali Beach Hotel (now known as the Inna Grand Bali Beach Hotel), the first five-star accommodation for tourism located in the Sanur region of Denpasar. This hotel has been in operation since 1966, has a height of 10 floors, which was the tallest building at the time of its construction and is still the tallest in Bali today. For example, over 11,278 foreign tourists per year visited the island in 1969, and by 2014, the number had grown to 3,766,638. From 1990 to 2000, foreign tourist visits to Bali grew an average of 17.08% per year, and from 2000 to 2014, the average increase was only 11.11% per year. The rapid growth of foreign tourists visiting the island before the new millennium may have led to rapid changes in land use in Denpasar from 1992 to 2003 by increasing the number of multi-functional buildings for tourism services and the incomes of employees in the tourism industry. In addition, the change in land use resulting from the indirect impact of the development of the tourism industry on population growth through migration in the southern part of Bali, including Denpasar, has increased the need for residential buildings. Denpasar City is the most populous area in Bali Province, and in 1992, the population density reached 2704 people per km2 and increased to 6622 people per km2 in 2014. Bali has experienced the most rapid economic development of all the islands of Indonesia [76]. One of the main reasons for this growth is the many income potential created directly by tourism and indirectly by industries that produce goods for the tourism industry. Consequently, tourism activities can ensure income increases and improvements to the standard of living. In total, 51 percent of people’s income and 38 percent of Bali’s employment opportunities are directly linked with expenditures by tourists and tourism investment [77]. On average, the tourism sector has accounted for 30% of the gross domestic product (GDP) of Bali Province over the past two decades, where economic growth is largely contributed by the southern regions as centres of tourism activity, including Denpasar city [14]. The growth in economic income, living standards, and population are common causes of land use change in urban areas [78,79]. The relationship between population growth and regional GDP and land use change was also described by Wu et al. [80] in the Yangtze River Basin, China. Moreover, an increase in the number of tourists can affect the land needs associated with their activities [72], and similar results were presented by Allen et al. [81] in South Carolina.
As previously described, the Inna Grand Bali Beach Hotel is the tallest building in Bali today with a height of 10 floors. In Bali, the Decree of the Governor of Bali Province No. 13/Perbang.1614/II/a/1971, which was made in 1971, sets the maximum building height to be 15 metres, typically only four floors. This regulation stems from local wisdom of the Hindu religion and Balinese culture on the urban and tourism development processes. This rule ensures that the construction of buildings such as residences, government offices, businesses, multi-functional buildings for tourism, and others do not spread skyward. Consequently, sideways expansion is a logical choice for obtaining large buildings with the largest volumes but directly and indirectly sacrifices greenspaces for settlement areas. This regulation may be one of the causes of the rapid change in land use in Denpasar, although government regulation can control and block the acceleration of land use change [82,83,84]. As mentioned earlier, Jakarta and Beijing did not experience any rapid land use change compared with Denpasar. Although both cities are capital cities, they allow the existence of taller buildings. Therefore, further studies related to the impacts of this regulation need to be undertaken to identify their linkage with declining greenspaces and the expansion of settlement.

3.2. Annual Changes in GPP

The VPM was used to estimate GPP in 1995, 2003 and 2014 in the tourism city of Denpasar. Table 4 reports the annual total GPP estimates for each of the land cover types and their percentage of the total. Meanwhile, the spatial distributions of annual GPP in Denpasar as analysed from the Landsat satellite data by the VPM method for 1995, 2003 and 2014 are presented in Figure 4. Because the current study only studied changes in GPP, the study only shows GPP values for land use types related to greenspaces and settlement and does not include GPP values for non-greenspaces and settlement land types. The annual GPP for the entire region in the years 1995, 2003 and 2014 was estimated to be 61,609.15, 55,621.81 and 53,815.19 tC year−1, respectively. In 1995, the estimated GPP in agriculture and grassland accounted for the largest portion, approximately 75.37% of the total estimated GPP, followed by forest with a proportion of 13.10%, and settlement accounted for 11.53%. It can be clearly seen in Figure 4a that the GPP was estimated to be higher than 550 gC m−2 year−1 in the agriculture and grass land use type. In 2003, the allocation of agriculture and grassland still accounted for the largest area, although it is lower than that in 1995, 55.33% of the total area, while the forest and settlement land use types exhibited increasing allocations of estimated GPP with proportions of 16.61% and 28.06% of the study area by the year 2003, respectively. There is an interesting spatial pattern of estimated GPP in 2003. Inside the area around the settlement, there is still a fairly high GPP value (>550 gC m−2 year−1), and the south part of Denpasar has been dominated by high GPP values as a result of the replanting of mangroves. However, the high estimated GPP regions in the southeastern part of Denpasar in 1995 became more fragmented and scattered with low values of GPP in 2003 (Figure 4b), while the same location in 2003 was dominated by multi-functional buildings that serve tourism activities. Interestingly, the GPP estimation for settled areas has overtaken that of agriculture and grass areas, with a value of 41.39% of the total GPP estimation in 2014. The GPP estimation from forest increased by 19.08%, while that of agriculture and grass areas also declined, although the proportion was still quite large with a share of 39.53% of the total area in 2014. It is noteworthy that compared to their spatial pattern in 2003, the medium to high value GPP values were more fragmented and scattered around the regions with low GPP estimates in 2014. High estimated GPP was clearly observed in the southern part of Denpasar and in small areas of the agriculture and grass land use type (Figure 4c).
Figure 5 shows the mean estimates of annual GPP as analysed from Landsat satellite data for each land use type for 1995, 2003 and 2014. The mean GPP estimation for settlements decreased slightly (339.23 to 304.99 gC m−2 year−1) during the 20-year period and increased marginally for the agriculture and grass land use type (593.53 to 605.65 gC m−2 year−1). Meanwhile, the average GPP estimation in forests over the past two decades increased dramatically compared to that for the other two types of land use from 1995 to 2014 (544.97 to 846.40 gC m−2 year−1). However, in general, the mean GPP estimated for the settlement and in agriculture and grassland land cover types was almost similar, and that for the forest increased over the past two decades. This study found that the estimated GPP for the settlement region is half of that for agriculture and grassland and that a change of 1 ha of agriculture and grassland into settlement will cause Denpasar city to lose 3 tC ha−1 year−1 of carbon uptake carbon by green vegetation (GPP).
The region-wide net decrease in GPP was approximately 7793.96 tC year−1 during the past two decades, a decrease of 12.65% from the GPP value in 1995 (Table 5). The estimated total GPP decreased during the 20-year period with an average decline of 0.63% per year. However, the fastest decline occurred in the first nine years compared to 2003 to 2014 (Table 6), which is similar to the decrease in the agriculture and grassland use type. Interestingly, based on land use data, GPP estimates increased for settlements and forests, although the forest area decreased, while a decline in the value of estimated GPP occurs only in the agriculture and grassland use type. The GPP estimate for settlements increased at a rate of 13.31% per year, while during 2003 to 2014, the increase in the GPP estimate was less than that for the first nine years, and over the past 20 years, the increase has been large, reaching 10.68% per year. The estimated GPP in forests over the past two decades has continued to grow, even though the forest area has declined, with the fastest growth occurring from 1995 to 2003. Decreasing agricultural and grass land use has also led to a decline in the value of GPP estimates, and the decline in the past 20 years reached more than 50%. The fastest change, a value of 3.75% per year, was found from 1995 to 2003, and from 2003 to 2014, the value of estimated GPP loss was approximately 2.57% per year.
The current study found that the total value of estimated GPP for the settlement land use type increased dramatically. From the lowest proportion in 1995, this land use type achieved the highest ability to uptake total carbon through vegetation of all land use types in 2014. The average value of estimated GPP is similar in the three observation periods, although the average value is lower than the result obtained by Zhao et al. [19] (390.24 gC m−2 year−1), which may be due to the different methods used in GPP calculation. The increase in the GPP estimation value is very high but is still less than the increase in the settlement area, where in general, there is a linear relationship between the increase in area and the total GPP in the settlement area (see Table 3 and Table 5). The linear relationship between changes in settlement area and primary production is consistent with the study of Lu et al. [5]. A linear relationship was also found in the agriculture and grass land use type, with decreases in area of approximately 55.29%, followed by a loss of ability to uptake total carbon of 54.19%. Although only approximately half of the agriculture and grass land use type has been converted, the area reached up to 4350.51 ha. The settlement areas increased to 4906.44 ha during the study period, which indicates that agriculture and grass have been converted to settlement areas. Consequently, the ability of cities to uptake carbon from the atmosphere is reduced, as displayed in Table 5. Settlement areas can uptake carbon due to the presence of vegetation. As-syakur et al. [85] have shown that not all settlement areas in Denpasar are covered by structures, but have various built-up land uses that include vegetated and non-vegetated areas such as backyards and front yards. However, this city is unique in that the majority of the population is Hindu, and the city is a “flat city” due to the absence of tall buildings. These characteristics are not found in other regions in Indonesia or in other regions in Southeast Asia. The city of Denpasar contains a “holy area” for Hindu religious activities at each residence that is characterized by the presence of vegetation with varying levels of cover and heights. However, the increase in settlement area has not impacted the average GPP estimates, unlike the results obtained by Zhao et al. [19], which show an increase in the mean annual GPP at the time of an increase in the built-up area in the eastern United States. Given that the development of settlement areas in the city is difficult to control due to high demand for tourism activities, an increased settlement area accompanied by good management of “holy areas”, backyards, and front yards should be able to increase the ability of settlement areas to uptake atmosphere carbon through vegetation.
Interestingly, although the forest areas are degraded, forest restoration has been helpful in preventing a drastic decline in GPP in Denpasar city over the past two decades. A forest area of approximately 175.41 ha was lost during the study period, but the annual estimated GPP increased by approximately 27.23% from 1995 to 2014. This increase is due to the restoration of mangrove forests in the southern part of Denpasar from 1995 to the present. Reforestation is conducted by replanting mangroves in water (pond) land use areas and maintaining these areas by making them protection areas. The mangrove forest in southern Denpasar is part of a protected area established by the Indonesian government, the Ngurah Rai Grand Forest Park, through decision No. 544/Kpts-II/1993 in 1993.
The current study indicated the importance of forest (mangrove) restoration to sustain the carbon cycle in urban areas that have experienced rapid land use change due to increased tourism activity and urbanization. Ecological restoration is important to increasing the total carbon uptake by vegetation. Successful ecological restoration associated with primary production has been reported previously by Yang et al. [86], but this result was an increase in the amount of carbon sequestration due to an increased terrestrial vegetation area, not under the decline of forest area. However, Briber et al. [87] found an increase in tree productivity during forest conversion to urban areas. Forest productivity is a function of many factors such as species, age, resource availability, growing season length, and competition [88,89]. Mangroves are among the most carbon-rich forests in the tropics [90]. Duarte et al. [91] also suggest that mangrove carbon production is more rapid than other estuarine and marine primary producers. On the other hand, Landsat data are some of the primary sources for identifying areas of forests in space and time with good accuracy [92], including in urban areas [93] and mangroves [94].
However, this study still has some uncertainties. The land use data used in this study do not match the year, especially for GPP analysis in 1995, which may cause a mismatch in extracting the estimated GPP. This mismatch will lead to an underestimation of GPP values in settlements and an overestimation in agriculture and grass areas because most of the settlement areas come from agriculture and grass areas. From 1992 to 1995, Denpasar had an average population growth rate of 2.85% per year, and foreign tourist visits to Bali grew an average of 13.47% per year in the same period. One consequence is the increasing number of multi-functional buildings for tourism services and the number of residences needed to accommodate the population growth. However, the four-year difference may not have a significant impact on a study that had a long span period. In addition, the use of 30-m Landsat data can cause errors in detecting vegetation in heterogeneous settlement areas. As shown by As-syakur et al. [46], different spatial resolutions in the satellite data used for analyses of GPP in urban areas will lead to different averages for settlement areas, and more detailed spatial resolutions will give larger average values. Increased pixel sizes (or decreased spatial resolution) result in the loss of image detail [95]. Satellite data with high spatial resolutions may be able to reduce the difficulties with remotely sensed data from coarse-resolution satellites [96]. Finally, this is the first study that used the VPM method with medium-spatial-resolution satellite imagery to calculate GPP in urban areas; therefore, the results are questionable. Previously, Ciu et al. [97] tried to compare GPP-based VPM estimates (using MODIS data) and MODIS GPP data products (MOD17A2) with solar-induced chlorophyll fluorescence (SIF) over the most populous megacity area with better results compared to those produced with the MOD17A2 product. On the other hand, the total annual GPP estimated using the VPM in the current study was not very different from the GPP estimate produced by other LUE models or from the MODIS GPP product, as presented by As-syakur et al. [46] in the same location. This study can be used as an initial source of information related to land use changes and their impact on terrestrial carbon uptake by vegetation. However, urban areas are highly heterogeneous landscapes that have different or even opposing effects on overall urban vegetation productivity [98]. Urban areas also have complex climate systems that are affected by complex socio-ecological systems [99] and the urban–rural proportion as related to urban heat islands [100], which directly and indirectly affect the calculation of light-use efficiency. Therefore, a validation of the estimated GPP results of the VPM method with eddy flux towers needs to be performed to advance our quantitative understanding of the capability of this method to contribute to analyses of vegetation carbon uptake in urban areas.

4. Conclusions

A recent study indicates that land use has changed in the tourism area of Denpasar over the past two decades and that this change has had overall negative effects on terrestrial GPP. The total GPP of the entire study area decreased by 7793.96 tC year−1 during the 20-year period, with an average decline of 0.63% per year. The fastest decline in GPP values was found in the first nine years, especially for the agriculture and grass land use type, and a rapid increase in the settlement area was found at the same time. Although the growth in the settlement area can increase the total GPP value, the amount is not sufficient to cover the total amount of carbon that can be taken up by other types of greenspaces. Based on the mean GPP value for each land use type, a change of 1 ha from agriculture and grass area into a settlement will cause a loss of 3 tC ha−1 year−1 in the ability of the city to uptake carbon through vegetation. Moreover, despite losses in forest area of approximately 12.68%, forest areas will still be able to increase the ability to uptake carbon from the atmosphere as long as the forest is maintained and restored by making them protection areas.
Further studies related to the direct impacts of tourism activity on terrestrial carbon uptake by vegetation need to be undertaken in the future. Moreover, seasonal analysis and validation of the GPP calculated using VPM also requires future research using eddy flux towers to determine the capability of the method to estimate GPP in urban areas. Finally, this study provides information that is useful for carbon resource management, the tourism industry, policy makers and scholars. In this way, the adverse impacts of land use change due to tourism activities in urban areas on ecosystem services can be minimized.

Acknowledgments

Support for this work provided was by the Udayana University through research grants number 641-11/UN14.2/PNL.01.03.00/2016 and 673-53/UN14.4.A/LT/2017. We gratefully acknowledge the data received from the following organizations: Landsat data from the U.S. Geological Survey; Meteorological data from the Indonesian Meteorology, Climatology, and Geophysics Agency (BMKG); and topographic maps for 1992, 2003 and 2014 from the Indonesian Geospatial Information Agency (BIG), the Japan International Cooperation Agency (JICA), and the Regional Planning and Development Board (Bappeda) of Denpasar City, respectively.

Author Contributions

I.W.N., A.R.A., I.G.A.G. and I.M.S. designed research, wrote the paper and discussed it together.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Foley, J.A.; DeFries, R.; Asner, G.P.; Barford, C.; Bonan, G.; Carpenter, S.R.; Chapin, F.S.; Coe, M.T.; Daily, G.C.; Gibbs, H.K.; et al. Global consequences of land use. Science 2005, 309, 570–574. [Google Scholar] [CrossRef] [PubMed]
  2. Hoymann, J.; Goetzke, R. Simulation and Evaluation of Urban Growth for Germany Including Climate Change Mitigation and Adaptation Measures. ISPRS Int. J. Geo-Inf. 2016, 5, 101. [Google Scholar] [CrossRef]
  3. Yang, L.; Xian, G.; Klaver, J.M.; Deal, B. Urban land-cover change detection through sub-pixel imperviousness mapping using remotely sensed data. Photogramm. Eng. Remote Sens. 2003, 69, 1003–1010. [Google Scholar] [CrossRef]
  4. Rahman, M.T. Detection of Land Use/Land Cover Changes and Urban Sprawl in Al-Khobar, Saudi Arabia: An Analysis of Multi-Temporal Remote Sensing Data. ISPRS Int. J. Geo-Inf. 2016, 5, 15. [Google Scholar] [CrossRef]
  5. Lu, D.; Xu, X.; Tian, H.; Moran, E.; Zhao, M.; Running, S. The effects of urbanization on net primary productivity in southeastern China. Environ. Manag. 2010, 46, 404–410. [Google Scholar] [CrossRef] [PubMed]
  6. Bilsborrow, R.E.; Ogendo, H.W. Population-driven changes in land use in developing countries. Ambio 1992, 21, 37–45. [Google Scholar]
  7. Ivanov, S.; Webster, C. Measuring the impact of tourism on economic growth. Tour. Econ. 2007, 13, 379–388. [Google Scholar] [CrossRef]
  8. Surugiu, C.; Surugiu, M.R. Is the tourism sector supportive of economic growth? Empirical evidence on Romanian tourism. Tour. Econ. 2013, 19, 115–132. [Google Scholar] [CrossRef]
  9. Williams, A.M.; Shaw, G. Future play: Tourism, recreation and land use. Land Use Policy 2009, 26, S326–S335. [Google Scholar] [CrossRef] [Green Version]
  10. Roth, D. Environmental sustainability and legal plurality in irrigation: The Balinese subak. Curr. Opin. Environ. Sustain. 2014, 11, 1–9. [Google Scholar] [CrossRef]
  11. Wijaya, A.; Budiharto, R.S.; Tosiani, A.; Murdiyarso, D.; Verchot, L.V. Assessment of large scale land cover change classifications and drivers of deforestation in Indonesia. ISPRS Arch. 2015, 40, 557. [Google Scholar] [CrossRef]
  12. Ramdani, F.; Putra, A.P.; Utomo, B.N. Historical Urban Land Use Transformation in Virtual Geo-Library. ISPRS Int. J. Geo-Inf. 2015, 4, 1500–1511. [Google Scholar] [CrossRef]
  13. Kunu, P.J.; Lelolterry, H. Penggunaan Lahan dan Evolusi Penggunaan Lahan di Provinsi DKI Jakarta. J. Agroforestr. 2010, 4, 203–2007. [Google Scholar]
  14. As-syakur, A.R. Perubahan penggunaan lahan di Provinsi Bali. Ecotrophic J. Environ. Sci. 2011, 6, 1–7. [Google Scholar]
  15. Chen, T.; Huang, Q.; Liu, M.; Li, M.; Qu, L.; Deng, S.; Chen, D. Decreasing Net Primary Productivity in Response to Urbanization in Liaoning Province, China. Sustainability 2017, 9, 162. [Google Scholar] [CrossRef]
  16. Kalnay, E.; Cai, M. Impact of urbanization and land-use change on climate. Nature 2003, 423, 528–531. [Google Scholar] [CrossRef] [PubMed]
  17. Nie, W.; Yuan, Y.; Kepner, W.; Nash, M.S.; Jackson, M.; Erickson, C. Assessing impacts of landuse and landcover changes on hydrology for the upper San Pedro watershed. J. Hydrol. 2011, 407, 105–114. [Google Scholar] [CrossRef]
  18. Poschlod, P.; Bakker, J.P.; Kahmen, S. Changing land use and its impact on biodiversity. Basic Appl. Ecol. 2005, 6, 93–98. [Google Scholar] [CrossRef]
  19. Zhao, T.; Brown, D.G.; Fang, H.; Theobald, D.M.; Liu, T.; Zhang, T. Vegetation productivity consequences of human settlement growth in the eastern united states. Landsc. Ecol. 2012, 27, 1149–1165. [Google Scholar] [CrossRef]
  20. Hutyra, L.R.; Yoon, B.; Alberti, M. Terrestrial carbon stocks across a gradient of urbanization: A study of the Seattle, WA region. Glob. Chang. Biol. 2011, 17, 783–797. [Google Scholar] [CrossRef]
  21. Yan, Y.; Zhang, C.; Hu, Y.; Kuang, W. Urban Land-Cover Change and Its Impact on the Ecosystem Carbon Storage in a Dryland City. Remote Sens. 2016, 8, 6. [Google Scholar] [CrossRef]
  22. Wang, X.; Ma, M.; Huang, G.; Veroustraete, F.; Zhang, Z.; Song, Y.; Tan, J. Vegetation primary production estimation at maize and alpine meadow over the Heihe River Basin, China. Int. J. Appl. Earth Obs. Geoinf. 2012, 17, 94–101. [Google Scholar] [CrossRef]
  23. Wagle, P.; Xiao, X.; Torn, M.S.; Cook, D.R.; Matamala, R.; Fischer, M.L.; Jin, C.; Dong, J.; Biradar, C. Sensitivity of vegetation indices and gross primary production of tallgrass prairie to severe drought. Remote Sens. Environ. 2014, 152, 1–14. [Google Scholar] [CrossRef]
  24. Nemani, R.R.; Keeling, C.D.; Hashimoto, H.; Jolly, W.M.; Piper, S.C.; Tucker, C.J.; Myneni, R.B.; Running, S.W. Climate-driven increases in global terrestrial net primary production from 1982 to 1999. Science 2003, 300, 1560–1563. [Google Scholar] [CrossRef] [PubMed]
  25. Zhao, T.; Brown, D.G.; Bergen, K.M. Increasing gross primary production (GPP) in the urbanizing landscapes of southeastern Michigan. Photogramm. Eng. Remote Sens. 2007, 73, 1159–1167. [Google Scholar] [CrossRef]
  26. Wu, C.Y.; Gonsamo, A.; Zhang, F.M.; Chen, J.M. The potential of the greenness and radiation (GR) model to interpret 8-day gross primary production of vegetation. ISPRS J. Photogramm. Remote Sens. 2014, 88, 69–79. [Google Scholar] [CrossRef]
  27. Inoue, Y.; Peñuelas, J.; Miyata, A.; Mano, M. Normalized difference spectral indices for estimating photosynthetic hyperspectral and CO2 flux measurements in rice. Remote Sens. Environ. 2008, 112, 156–172. [Google Scholar] [CrossRef]
  28. Turner, D.P.; Gowerb, S.T.; Cohenc, W.B.; Gregorya, M.; Maiersperger, T.K. Effects of spatial variability in light use efficiency on Satellite-Based NPP monitoring. Remote Sens. Environ. 2002, 80, 397–405. [Google Scholar] [CrossRef]
  29. Myneni, R.B.; Williams, D.L. On the relationship between FAPAR and NDVI. Remote Sens. Environ. 1994, 49, 200–211. [Google Scholar] [CrossRef]
  30. Iwata, H.; Ueyama, M.; Iwama, C.; Harazono, Y. Variations in fraction of absorbed photosynthetically active radiation and comparisons with MODIS data in burned black spruce forests of Interior Alaska. Pol. Sci. 2013, 7, 113–124. [Google Scholar] [CrossRef]
  31. Green, K.; Kempka, D.; Lackey, L. Using remote sensing to detect and monitor land-cover and land-use change. Photogramm. Eng. Remote Sens. 1994, 60, 331–337. [Google Scholar]
  32. Rogan, J.; Chen, D. Remote sensing technology for mapping and monitoring land-cover and land-use change. Prog. Plan. 2004, 61, 301–325. [Google Scholar] [CrossRef]
  33. Gómez, C.; White, J.C.; Wulder, M.A. Optical remotely sensed time series data for land cover classification: A review. ISPRS J. Photogramm. Remote Sens. 2016, 116, 55–72. [Google Scholar] [CrossRef]
  34. Turner, D.P.; Ritts, W.D.; Cohen, W.B.; Gower, S.T.; Zhao, M.; Running, S.W.; Wofsy, S.C.; Urbanski, S.; Dunn, A.L.; Munger, J.W. Scaling gross primary production (GPP) over boreal and deciduous forest landscapes in support of MODIS GPP product validation. Remote Sens. Environ. 2003, 88, 256–270. [Google Scholar] [CrossRef]
  35. Gitelson, A.A.; Peng, Y.; Masek, J.G.; Rundquist, D.C.; Verma, S.; Suyker, A.; Baker, J.M.; Hatfield, J.L.; Meyers, T. Remote estimation of crop gross primary production with Landsat data. Remote Sens. Environ. 2012, 121, 404–414. [Google Scholar] [CrossRef]
  36. Chapin, F.S., III; Matson, P.A.; Vitousek, P. Principles of Terrestrial Ecosystem Ecology, 2nd ed.; Springer: New York, NY, USA, 2011. [Google Scholar]
  37. Yuan, W.; Cai, W.; Liu, S.; Dong, W.; Chen, J.; Arain, M.A.; Blanken, P.D.; Cescatti, A.; Wohlfahrt, G.; Georgiadis, T.; et al. Vegetation-specific model parameters are not required for estimating gross primary production. Ecol. Model. 2014, 292, 1–10. [Google Scholar] [CrossRef]
  38. Sims, D.A.; Luo, H.; Hastings, S.; Oechel, W.C.; Rahman, A.F.; Gamon, J.A. Parallel adjustments in vegetation greenness and ecosystem CO2 exchange in response to drought in a southern California chaparral ecosystem. Remote Sens. Environ. 2006, 103, 289–303. [Google Scholar] [CrossRef]
  39. Xiao, X.; Hollinger, D.; Aber, J.; Goltz, M.; Zhang, Q. Satellite-based modeling of gross primary production in an evergreen needle leaf forest. Remote Sens. Environ. 2004, 89, 519–534. [Google Scholar] [CrossRef]
  40. Xiao, X.; Zhang, Q.; Saleska, S.; Hutyra, L.; de Camargo, P.; Wofsy, S.; Frolking, S.; Boles, S.; Keller, M.; Moore, B. Satellite-based modeling of Gross Primary Production in a seasonally moist tropical evergreen forest. Remote Sens. Environ. 2005, 94, 105–122. [Google Scholar] [CrossRef]
  41. Madugundu, R.; Al-Gaadi, K.A.; Tola, E.; Kayad, A.G.; Jha, C.S. Estimation of gross primary production of irrigated maize using Landsat-8 imagery and Eddy Covariance data. Saudi J. Biol. Sci. 2017, 24, 410–420. [Google Scholar] [CrossRef] [PubMed]
  42. Danelichen, V.H.; Biudes, M.S.; Velasque, M.; Machado, N.G.; Gomes, R.S.; Vourlitis, G.L.; Nogueira, J.S. Estimating of gross primary production in an Amazon-Cerrado transitional forest using MODIS and Landsat imagery. An. Acad. Bras. Cienc. 2015, 87, 1545–1564. [Google Scholar] [CrossRef] [PubMed]
  43. Vetrita, Y.; Chaoyang, W.; Zheng, N.; Hirano, T. Evaluation of light use efficiency model using modis in tropical peat swamp forest, Indonesia. In Proceedings of the Second CReSOS International Symposium on Southeast Asia Environmental Problems and Satellite Remote Sensing, Bali, Indonesia, 21–22 February 2011; pp. 127–134. [Google Scholar]
  44. As-syakur, A.R.; Tanaka, T.; Osawa, T.; Mahendra, M.S. Indonesian rainfall variability observation using TRMM multi-satellite data. Int. J. Remote Sens. 2013, 34, 7723–7738. [Google Scholar] [CrossRef]
  45. Aldrian, E.; Susanto, D.R. Identification of three dominant rainfall regions within Indonesia and their relationship to sea surface temperature. Int. J. Climatol. 2003, 23, 1435–1452. [Google Scholar] [CrossRef]
  46. As-syakur, A.R.; Osawa, T.; Adnyana, I.W.S. Medium Spatial Resolution Satellite Imagery to Estimate Gross Primary Production in an Urban Area. Remote Sens. 2010, 2, 1496–1507. [Google Scholar] [CrossRef]
  47. Davis, A.; Marshak, A.; Cahalan, R.; Wiscombe, W. The Landsat scale break in stratocumulus as a three-dimensional radiative transfer effect: Implications for cloud remote sensing. J. Atmos. Sci. 1997, 54, 241–260. [Google Scholar] [CrossRef]
  48. Goodwin, N.R.; Collett, L.J.; Denham, R.J.; Flood, N.; Tindall, D. Cloud and cloud shadow screening across Queensland, Australia: An automated method for Landsat TM/ETM+ time series. Remote Sens. Environ. 2013, 134, 50–65. [Google Scholar] [CrossRef]
  49. Zhu, Z.; Woodcock, C.E. Object-based cloud and cloud shadow detection in Landsat imagery. Remote Sens. Environ. 2012, 118, 83–94. [Google Scholar] [CrossRef]
  50. Irish, R.R.; Barker, J.L.; Goward, S.N.; Arvidson, T. Characterization of the Landsat-7 ETM+ automated cloud-cover assessment (ACCA) algorithm. Photogramm. Eng. Remote Sens. 2006, 72, 1179–1188. [Google Scholar] [CrossRef]
  51. Liu, H.Q.; Huete, A.R. A feedback based modification of the NDV I to minimize canopy background and atmospheric noise. IEEE Trans. Geosci. Remote Sens. 1995, 33, 457–465. [Google Scholar]
  52. Huete, A.R.; Liu, H.Q.; Batchily, K.; Van Leeuwen, W.J.D.A. A comparison of vegetation indices over a global set of TM images for EOS-MODIS. Remote Sens. Environ. 1997, 59, 440–451. [Google Scholar] [CrossRef]
  53. Xiao, X.; Boles, S.; Frolking, S.; Salas, W.; Moore, B., III; Li, C.; He, L.; Zhao, R. Observation of flooding and rice transplanting of paddy rice fields at the site to landscape scales in China using VEGETATION sensor data. Int. J. Remote Sens. 2002, 23, 3009–3022. [Google Scholar] [CrossRef]
  54. Japan International Cooperation Agency (JICA). The Comprehensive Study on Water Resources Development and Management in Bali Province, in the Republic of Indonesia: Final Report; Japan International Cooperation Agency and Directorate General of Water Resources, Indonesian Ministry of Public Works: Jakarta, Indonesia, 2006.
  55. Hirano, T.; Segah, H.; Harada, T.; Limin, S.; June, T.; Hirata, R.; Osaki, M. Carbon dioxide balance of a tropical peat swamp forest in Kalimantan, Indonesia. Glob. Chang. Biol. 2007, 13, 412–425. [Google Scholar] [CrossRef]
  56. Raich, J.W.; Rastetter, E.B.; Melillo, J.M.; Kicklighter, D.W.; Steudler, P.A.; Peterson, B.J.; Grace, A.L.; Moore, B.; Vorosmarty, C.J. Potential net primary productivity in South America: Application of a global model. Ecol. Appl. 1991, 1, 399–429. [Google Scholar] [CrossRef] [PubMed]
  57. Nguy-Robertson, A.; Suyker, A.; Xiao, X. Modeling gross primary production of maize and soybean croplands using light quality, temperature, water stress, and phenology. Agric. For. Meteorol. 2015, 213, 160–172. [Google Scholar] [CrossRef]
  58. Kalfas, J.L.; Xiao, X.; Vanegas, D.X.; Verma, S.B.; Suyker, A.E. Modeling gross primary production of irrigated and rain-fed maize using MODIS imagery and CO2 flux tower data. Agric. For. Meteorol. 2011, 151, 1514–1528. [Google Scholar] [CrossRef]
  59. Jacovides, C.; Tymvios, F.; Asimakopoulos, D.; Theofilou, K.; Pashiardes, S. Global photosynthetically active radiation and its relationship with global solar radiation in the Eastern Mediterranean basin. Theor. Appl. Climatol. 2003, 74, 227–233. [Google Scholar] [CrossRef]
  60. Rao, C.N. Photosynthetically active components of global solar radiation: Measurements and model computations. Arch. Meteorol. Geophys. B 1984, 34, 353–364. [Google Scholar]
  61. González, J.A.; Calbó, J. Modelled and measured ratio of PAR to global radiation under cloudless skies. Agric. For. Meteorol. 2002, 110, 319–325. [Google Scholar] [CrossRef]
  62. Escobedo, J.F.; Gomes, E.N.; Oliveira, A.P.; Soares, J. Ratios of UV, PAR and NIR components to global solar radiation measured at Botucatu site in Brazil. Renew. Energy 2011, 36, 169–178. [Google Scholar] [CrossRef]
  63. Jin, C.; Xiao, X.; Wagle, P.; Griffis, T.; Dong, J.; Wu, C.; Qin, Y.; Cook, D.R. Effects of in-situ and reanalysis climate data on estimation of cropland gross primary production using the Vegetation Photosynthesis Model. Agric. For. Meteorol. 2015, 213, 240–250. [Google Scholar] [CrossRef]
  64. Halawa, E.E.H.; Sugiyatno, S. Estimation of global solar radiation in the Indonesian climatic region. Renew. Energy 2001, 24, 197–206. [Google Scholar] [CrossRef]
  65. As-syakur, A.R.; Osawa, T.; Miura, F.; Nuarsa, I.W.; Ekayanti, N.W.; Dharma, I.G.B.S.; Adnyana, I.W.S.; Arthana, I.W.; Tanaka, T. Maritime Continent rainfall variability during the TRMM era: The role of monsoon, topography and El Niño Modoki. Dyn. Atmos. Oceans. 2016, 75, 58–77. [Google Scholar] [CrossRef]
  66. Nuarsa, I.W.; Nishio, F.; Hongo, C.; Mahardika, I.G. Using variance analysis of multitemporal MODIS images for rice field mapping in Bali Province, Indonesia. Int. J. Remote Sens. 2012, 33, 5402–5417. [Google Scholar] [CrossRef]
  67. Takama, T.; Setyani, P.; Aldrian, E. Climate change vulnerability to rice paddy production in Bali, Indonesia. In Handbook of Climate Change Adaptation; Filho, W.L., Ed.; Springer: Berlin/Heidelberg, Germany, 2015; pp. 1731–1757. ISBN 978-3-642-38669-5. [Google Scholar]
  68. Jayanti, I.A.G.K.; Osawa, T.; Adnyana, I.W.S.; Tanaka, T.; Nuarsa, I.W.; As-syakur, A.R. Multitemporal MODIS Data to Mapping Rice Field Distribution in Bali Province of Indonesia Based on the Temporal Dynamic Characteristics of the Rice Plant. Earth Sci. Res. 2012, 1, 64–74. [Google Scholar]
  69. As-syakur, A.R.; Nuarsa, I.W.; Osawa, T. Impacts of El Niño on agricultural drought in Bali, Indonesia. In Proceedings of the 19th Symposium on Remote Sensing for Environment, Chiba, Japan, 16 February 2017; p. 54. [Google Scholar]
  70. Stevenson, N.J.; Lewis, R.R.; Burbridge, P.R. Disused shrimp ponds and mangrove rehabilitation. In An International Perspective on Wetland Rehabilitation; Streever, W.J., Ed.; Springer: Dordrecht, The Netherlands, 1999; pp. 277–297. ISBN 978-94-011-4683-8. [Google Scholar]
  71. Kytzia, S.; Walz, A.; Wegmann, M. How can tourism use land more efficiently? A model–based approach to land–use efficiency for tourist destinations. Tour. Manag. 2011, 32, 629–640. [Google Scholar] [CrossRef]
  72. Xi, J.; Zhao, M.; Ge, Q.; Kong, Q. Changes in land use of a village driven by over 25 years of tourism: The case of Gougezhuang village, China. Land Use Policy 2014, 40, 119–130. [Google Scholar] [CrossRef]
  73. Mao, X.; Meng, J.; Wang, Q. Modeling the effects of tourism and land regulation on land-use change in tourist regions: A case study of the Lijiang River Basin in Guilin, China. Land Use Policy 2014, 41, 368–377. [Google Scholar] [CrossRef]
  74. Wu, Q.; Li, H.Q.; Wang, R.S.; Paulussen, J.; He, Y.; Wang, M.; Wang, B.H.; Wang, Z. Monitoring and predicting land use change in Beijing using remote sensing and GIS. Landsc. Urban Plan. 2006, 78, 322–333. [Google Scholar] [CrossRef]
  75. Howe, L. The Changing World of Bali: Religion, Society and Tourism; Routledge: New York, NY, USA, 2006. [Google Scholar]
  76. Lorenzen, S.; Lorenzen, R.P. A Case Study of Balinese Irrigation Management: Institutional Dynamics and Challenges. In Proceedings of the Second Southeast Asia Water Forum, Bali, Indonesia, 29 August–3 September 2005. [Google Scholar]
  77. Pitana, I.G. Tri Hita Karana–The Local Wisdom of the Balinese in Managing Development. In Trends and Issues in Global Tourism 2010; Conrady, R., Buck, M., Eds.; Springer: Berlin/Heidelberg, Germany, 2010; pp. 139–150. ISBN 978-3-642-10829-7. [Google Scholar]
  78. Lambin, E.F.; Turner, B.L.; Geist, H.J.; Agbola, S.B.; Angelsen, A.; Bruce, J.W.; Coomes, O.T.; Dirzo, R.; Fischer, G.; Folke, C.; et al. The causes of land-use and land-cover change: Moving beyond the myths. Glob. Environ. Chang. Part A Hum. Policy Dimens. 2001, 11, 261–269. [Google Scholar] [CrossRef]
  79. Parker, D.C.; Manson, S.M.; Janssen, M.A.; Hoffmann, M.J.; Deadman, P. Multi-Agent Systems for the simulation of land-use and land-cover change: A review. Ann. Assoc. Am. Geogr. 2003, 93, 314–337. [Google Scholar] [CrossRef]
  80. Wu, X.; Shen, Z.; Liu, R.; Ding, X. Land Use/Cover Dynamics in Response to Changes in Environmental and Socio-Political Forces in the Upper Reaches of Yangtze River, China. Sensors 2008, 8, 8104–8122. [Google Scholar] [CrossRef] [PubMed]
  81. Allen, J.; Lu, K.S.; Potts, T.D. A GIS-based analysis and prediction of land-use change in a coastal tourism destination area. In Proceedings of the 1999 International Symposium on Coastal and Marine Tourism: Balancing Tourism and Conservation, Vancouver, BC, Canada, 26–29 April 1999; pp. 287–297. [Google Scholar]
  82. Conway, T.M.; Lathrop, R.G. Alternative land use regulations and environmental impacts: Assessing future land use in an urbanizing watershed. Landsc. Urban Plan. 2005, 71, 1–15. [Google Scholar] [CrossRef]
  83. Quigley, J.M.; Rosenthal, L.A. The effects of land use regulation on the price of housing: What do we know? What can we learn? Cityscape: J. Pol. Devel. Res. 2005, 8, 69–137. [Google Scholar]
  84. Mandelker, D.R. The role of the local comprehensive plan in land use regulation. Mich. Law Rev. 1976, 74, 899–973. [Google Scholar] [CrossRef]
  85. As-syakur, A.R.; Adnyana, I.W.S.; Arthana, I.W.; Nuarsa, I.W. Enhanced Built-Up and Bareness Index (EBBI) for Mapping Built-Up and Bare Land in an Urban Area. Remote Sens. 2012, 4, 2957–2970. [Google Scholar] [CrossRef]
  86. Yang, H.F.; Mu, S.J.; Li, J.L. Effects of ecological restoration projects on land use and land cover change and its influences on territorial NPP in Xinjiang, China. Catena 2014, 115, 85–95. [Google Scholar] [CrossRef]
  87. Briber, B.M.; Hutyra, L.R.; Reinmann, A.B.; Raciti, S.M.; Dearborn, V.K.; Holden, C.E.; Dunn, A.L. Tree productivity enhanced with conversion from forest to urban land covers. PLoS ONE 2015, 10, e0136237. [Google Scholar] [CrossRef] [PubMed]
  88. Iakovoglou, V.; Thompson, J.; Burras, L.; Kipper, R. Factors related to tree growth across urban-rural gradients in the Midwest, USA. Urban Ecosyst. 2001, 5, 71–85. [Google Scholar] [CrossRef]
  89. MCDonald, R.I.; Urban, D.L. Forest edges and tree growth rates in the North Carolina Piedmont. Ecology 2004, 85, 2258–2266. [Google Scholar] [CrossRef]
  90. Donato, D.C.; Kauffman, J.B.; Murdiyarso, D.; Kurnianto, S.; Stidham, M.; Kanninen, M. Mangroves among the most carbon-rich forests in the tropics. Nat. Geosci. 2011, 4, 293–297. [Google Scholar] [CrossRef]
  91. Duarte, C.M.; Middelburg, J.J.; Caraco, N. Major role of marine vegetation on the oceanic carbon cycle. Biogeosciences 2005, 2, 1–8. [Google Scholar] [CrossRef]
  92. Cohen, W.B.; Yang, Z.; Kennedy, R. Detecting trends in forest disturbance and recovery using yearly Landsat time series: 2. TimeSync—Tools for calibration and validation. Remote Sens. Environ. 2010, 114, 2911–2924. [Google Scholar] [CrossRef]
  93. Schneider, A. Monitoring land cover change in urban and peri-urban areas using dense time stacks of Landsat satellite data and a data mining approach. Remote Sens. Environ. 2012, 124, 689–704. [Google Scholar] [CrossRef]
  94. Lee, T.; Yeh, H. Applying remote sensing techniques to monitor shifting wetland vegetation: A case study of Danshui River estuary mangrove communities, Taiwan. Ecol. Eng. 2009, 35, 487–496. [Google Scholar] [CrossRef]
  95. Chen, D.; Stow, D.A.; Gong, P. Examining the effect of spatial resolution and texture window size on classification accuracy: An urban environment case. Int. J. Remote Sens. 2004, 25, 2177–2192. [Google Scholar] [CrossRef]
  96. Rocchini, D. Effects of spatial and spectral resolution in estimating ecosystem A-diversity by satellite imagery. Remote Sens. Environ. 2007, 111, 423–434. [Google Scholar] [CrossRef]
  97. Cui, Y.; Xiao, X.; Zhang, Y.; Dong, J.; Qin, Y.; Doughty, R.B.; Zhang, G.; Wang, J.; Wu, X.; Qin, Y.; et al. Temporal consistency between gross primary production and solar-induced chlorophyll fluorescence in the ten most populous megacity areas over years. Sci. Rep. 2017, 7, 14963. [Google Scholar] [CrossRef] [PubMed]
  98. Churkina, G. Modeling the carbon cycle of urban systems. Ecol. Modell. 2008, 216, 107–113. [Google Scholar] [CrossRef]
  99. Chowdhury, R.R.; Larson, K.; Grove, M.; Polsky, C.; Cook, E.; Onsted, J.; Ogden, L. A multi-scalar approach to theorizing socio-ecological dynamics of urban residential landscapes. Cities Environ. 2011, 4, 6. [Google Scholar] [CrossRef]
  100. Stewart, I.D.; Oke, T.R. Local climate zones for urban temperature studies. Bull. Am. Meteorol. Soc. 2012, 93, 1879–1900. [Google Scholar] [CrossRef]
Figure 1. The research location: Denpasar, Bali Province, Indonesia. Dark grey regions in southern Denpasar are Serangan Island and Benoa port after reclamation, and these areas are used for Denpasar boundaries in 2003 and 2014. The dotted line inside Denpasar is the district boundaries.
Figure 1. The research location: Denpasar, Bali Province, Indonesia. Dark grey regions in southern Denpasar are Serangan Island and Benoa port after reclamation, and these areas are used for Denpasar boundaries in 2003 and 2014. The dotted line inside Denpasar is the district boundaries.
Ijgi 07 00057 g001
Figure 2. Scatter plots of the Landsat base index from the main data with the data used for patching the no-data regions caused by clouds and their shadows. Top panels present enhanced vegetation index (EVI) and bottom panels present land surface water index (LWSI). (a,b) for Landsat 5 on 3 February 1995 with 26 May 1995; (c,d) for Landsat 7 on 21 March 2003 with 28 March 2003 (black dots) and for 21 March 2003 with 22 April 2003 (red dots); and (e,f) for Landsat 8 on 27 March 2014 with 11 March 2014.
Figure 2. Scatter plots of the Landsat base index from the main data with the data used for patching the no-data regions caused by clouds and their shadows. Top panels present enhanced vegetation index (EVI) and bottom panels present land surface water index (LWSI). (a,b) for Landsat 5 on 3 February 1995 with 26 May 1995; (c,d) for Landsat 7 on 21 March 2003 with 28 March 2003 (black dots) and for 21 March 2003 with 22 April 2003 (red dots); and (e,f) for Landsat 8 on 27 March 2014 with 11 March 2014.
Ijgi 07 00057 g002
Figure 3. Distribution of land use types in Denpasar, Indonesia from topographic maps. (a) 1992; (b) 2003; and (c) 2014.
Figure 3. Distribution of land use types in Denpasar, Indonesia from topographic maps. (a) 1992; (b) 2003; and (c) 2014.
Ijgi 07 00057 g003
Figure 4. Annual gross primary production (GPP) distribution in Denpasar from Landsat satellite data analysis. The grey colour indicates that the area is not counted due to land use types such as water and others. (a) 1995; (b) 2003; and (c) 2014.
Figure 4. Annual gross primary production (GPP) distribution in Denpasar from Landsat satellite data analysis. The grey colour indicates that the area is not counted due to land use types such as water and others. (a) 1995; (b) 2003; and (c) 2014.
Ijgi 07 00057 g004
Figure 5. Mean of annual GPP in Denpasar from Landsat satellite data analysis of each land use type.
Figure 5. Mean of annual GPP in Denpasar from Landsat satellite data analysis of each land use type.
Ijgi 07 00057 g005
Table 1. Definition of land use types for the current study.
Table 1. Definition of land use types for the current study.
Land UseDescription
SettlementCombines built-up areas, roads, backyards, and parking
ForestCombines mangroves, perennial plants, and shrubland
Agriculture and grassCombines agricultural farmlands (irrigated paddy fields), dry land, grassy fields, golf courses, and parks
OthersCombines bare land and sand
WaterCombines rivers, fishponds, and dams
Table 2. Area of each type of land use in Denpasar from 1992 to 2014 and their proportion of total area.
Table 2. Area of each type of land use in Denpasar from 1992 to 2014 and their proportion of total area.
Land use199220032014
Area (ha)%Area (ha)%Area (ha)%
Settlement2396.3419.89%5421.8743.49%7302.7858.57%
Forest1382.9411.48%1226.439.84%1207.539.68%
Agriculture and grass7868.0765.31%5239.9842.03%3517.5628.21%
Others62.550.52%359.552.88%210.061.68%
Water337.322.80%220.231.77%230.131.85%
Total12,047.22 12,468.06 12,468.06
Table 3. Changes in land use types in Denpasar from 1992 to 2014 in percent and the land use expansion index (LUEI) value (% year−1).
Table 3. Changes in land use types in Denpasar from 1992 to 2014 in percent and the land use expansion index (LUEI) value (% year−1).
Land Use1992 to 20032003 to 20141992 to 2014
ChangesLUEI (% year−1)ChangesLUEI (% year−1)ChangesLUEI (% year−1)
Settlement126.26%10.52%34.69%2.89%204.75%10.24%
Forest−11.32%−0.94%−1.54%−0.13%−12.68%−0.63%
Agriculture and grass−33.40%−2.78%−32.87%−2.74%−55.29%−2.76%
Others474.82%39.57%−41.58%−3.46%235.83%11.79%
Water−34.71%−2.89%4.50%0.37%−31.78%−1.59%
Table 4. Annual gross primary production (GPP) in Denpasar from 1995 to 2014 for each different land use as a total (tC year−1) and a percentage of the total.
Table 4. Annual gross primary production (GPP) in Denpasar from 1995 to 2014 for each different land use as a total (tC year−1) and a percentage of the total.
Land use199520032014
Total%Total%Total%
Settlement7102.0411.53%15,607.4928.06%22,274.9341.39%
Forest8071.1813.10%9238.8516.61%10,268.5819.08%
Agriculture and grass46,435.9375.37%30,775.4855.33%21,271.6839.53%
Total61,609.15 55,621.81 53,815.19
Table 5. GPP estimation changes in each type of land use in Denpasar from 1995 to 2014.
Table 5. GPP estimation changes in each type of land use in Denpasar from 1995 to 2014.
Land Use1995 to 20032003 to 20141995 to 2014
Settlement119.76%42.72%213.64%
Forest14.47%11.15%27.23%
Agriculture and grass−33.72%−30.88%−54.19%
Total−9.72%−3.25%−12.65%
Table 6. Changes in GPP estimation per year (% year−1) for each type of land use from 1995 to 2014.
Table 6. Changes in GPP estimation per year (% year−1) for each type of land use from 1995 to 2014.
Land Use1995 to 20032003 to 20141995 to 2014
Settlement13.31%3.56%10.68%
Forest1.61%0.93%1.36%
Agriculture and grass−3.75%−2.57%−2.71%
Total−1.08%−0.27%−0.63%

Share and Cite

MDPI and ACS Style

Nuarsa, I.W.; As-syakur, A.R.; Gunadi, I.G.A.; Sukewijaya, I.M. Changes in Gross Primary Production (GPP) over the Past Two Decades Due to Land Use Conversion in a Tourism City. ISPRS Int. J. Geo-Inf. 2018, 7, 57. https://doi.org/10.3390/ijgi7020057

AMA Style

Nuarsa IW, As-syakur AR, Gunadi IGA, Sukewijaya IM. Changes in Gross Primary Production (GPP) over the Past Two Decades Due to Land Use Conversion in a Tourism City. ISPRS International Journal of Geo-Information. 2018; 7(2):57. https://doi.org/10.3390/ijgi7020057

Chicago/Turabian Style

Nuarsa, I Wayan, Abd. Rahman As-syakur, I Gusti Alit Gunadi, and I Made Sukewijaya. 2018. "Changes in Gross Primary Production (GPP) over the Past Two Decades Due to Land Use Conversion in a Tourism City" ISPRS International Journal of Geo-Information 7, no. 2: 57. https://doi.org/10.3390/ijgi7020057

APA Style

Nuarsa, I. W., As-syakur, A. R., Gunadi, I. G. A., & Sukewijaya, I. M. (2018). Changes in Gross Primary Production (GPP) over the Past Two Decades Due to Land Use Conversion in a Tourism City. ISPRS International Journal of Geo-Information, 7(2), 57. https://doi.org/10.3390/ijgi7020057

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

Article Metrics

Back to TopTop