602
Views
3
CrossRef citations to date
0
Altmetric
Original Articles

Spatial surrogate for domestic combustion’s air emissions: A case study from Silesian Metropolis, Poland

&
Pages 1012-1019 | Received 20 Jan 2017, Accepted 30 Mar 2017, Published online: 31 Jul 2017

ABSTRACT

The aim of this study is improving currently applied methodology for spatial disaggregation, as well as mapping air emission inventories by taking into account the auxiliary spatial data on population density. District heating infrastructure occurring in more populated areas changes distinctly the spatial distribution of estimated air emission; however, it does not change the initial estimate. That means the total, disaggregated value is constant. Considered sources of domestic combustion are located in the central part of the Silesian Metropolis, in the southern part of Poland. A large part of this area is strongly urbanized and supplied with heat (hot water) from the district heating system. Data on population density help to determine the area within which the dwellers use heat energy and hot water supplied by the heating infrastructure, apart from heating with small domestic boilers and stoves. This causes the domestic combustion’s emission impact within the distinguished area to be significantly lower in comparison to the official guidelines on air emission inventories. The important differences in spatial air emissions distributions calculated using a top-down approach are found for strongly urbanized areas supplied partly with heat and hot water from the district heating network. This fact should be taken into account when preparing detailed, high-resolution emission inventories for air regional and local quality modeling.

Implications: The spatial issues connected with elaboration of the high-resolution emission inventories are presented for the example of the populated area of the Silesian Metropolis (Poland). Spatial distribution of the population density is used to determine the area supplied with heat and hot water from the district heating system. It changes distinctly the spatial distribution of the air emission from small residential combustion sources.

Introduction

A top-down approach is widely applied for compilation of the emission inventories of air pollutants (European Environment Agency [EEA], 2016) and greenhouse gases (Intergovernmental Panel on Climate Change [IPCC], 2006). The general idea of this technique is to estimate air emission at the national level and disaggregate to a desired scale (size of grid cell; Horabik-Pyzel and Nahorski, Citation2016) using spatial surrogates.

Spatial surrogates, described widely by Maes et al. (Citation2009), are spatially referenced maps of selected statistics, in which spatial distributions are proportional to unknown, real spatial distribution of estimated air emission. This methodology is qualitatively cheap (Boychuk et al., Citation2012) and can significantly help to reduce the uncertainty in (national) total emission inventory (Bun et al., Citation2010).

A number of studies, including Fameli and Assimakopoulos (Citation2016), Kryza et al. (Citation2010), Kuenen et al. (Citation2014), Maes et al. (Citation2009), Markakis et al. (Citation2010), and Wang and Li (Citation2016), frequently apply various spatial surrogates to disaggregate global estimates; however, the dependency between the real unknown spatial emission distribution and particular surrogate is not always precise enough.

The quality of the applied spatial surrogate strongly influences the quality of the obtained spatial emission distribution. The result of applying inappropriate surrogates is frequently the misinterpretation of the emission map (result). Furthermore, this affects the air quality modeling results and provokes various misunderstandings among policymakers at all levels of national administration.

In this paper we present an emission surrogate for air emissions from residential combustion, including domestic heating and cooking primarily. The initial approach using population density (Gallego, Citation2010; Kryza et al., Citation2010) is improved with geostatistical methods to determine the area within which the residential emission’s impact is lower due to using the heat energy and hot water supplied from the district heating network.

The analysis performed concerns small domestic combustion sources, located in the Katowice Urban Zone (KUZ), which is a part of the Silesian Metropolis situated in the southern part of Poland ().

Figure 1. (a) Location of the Katowice Urban Zone in Poland and Europe. (b) Detailed map of the KUZ’s districts with the population density (inhabitants km−2).

Figure 1. (a) Location of the Katowice Urban Zone in Poland and Europe. (b) Detailed map of the KUZ’s districts with the population density (inhabitants km−2).

Materials and methods

Statistical data

The considered area (approximate extent: 100 km W–E and 60 km N–S) consists of 41 communities, which are the Polish local administrative units (LAUs). The most populated LAUs (population density more than 2,000 inhabitants km−2) are located in the central part of the KUZ (red areas, ). A significant number of inhabitants settled in urban areas are provided with heat and hot water from the district heating system, named the Silesian Heating System (SHS). The simplified map of the district heating infrastructure, derived from expertise of Plebankiewicz and Jankowski (Citation2007), and Nikodem (Citation2008), is presented in .

Figure 2. The simplified scheme of the SHS infrastructure: triangles, plants; lines, main heating arteries.

Figure 2. The simplified scheme of the SHS infrastructure: triangles, plants; lines, main heating arteries.

Comparing the map of population density (, and further: ) with the simplified scheme of the SHS (), the widest access to the SHS infrastructure occurs in the most populated areas. Moreover, the high-resolution population density distribution () is compliant with the spatial distribution of the district heating infrastructure (), due to higher density accumulated in the center of the KUZ.

Figure 3. Population density (inhabitants·km−2) in the KUZ aggregated in 2,000 m × 2,000 m grid (Gallego, Citation2010).

Figure 3. Population density (inhabitants·km−2) in the KUZ aggregated in 2,000 m × 2,000 m grid (Gallego, Citation2010).

Apart from using the heat and hot water from the district heating system, a significant part of the settled people use small domestic heating sources (boilers, stoves). This fact results in a considerable change of emission structure (share of particular emission sources) and its spatial distribution, because at the highly populated areas the air emission from small stoves is not proportional to the population density. Thus, to determine the area supplied from the SHS infrastructure, where the share of individual heating sources is smaller, it is necessary to take into account also auxiliary data, such as the percentage of population using the public heating system. This kind of data is usually scarce, affected by the methodological issues (simplifications), hard to obtain, or derived from experts’ judgments.

Spatial approach

The spatial approach applied in this paper uses population density to determine the area supplied with the heat and hot water from the SHS system. presents the spatial distribution of population density at the area of KUZ. The data, originally mapped in fine grid (100 m × 100 m) (Gallego Citation2010), is aggregated to coarse: 2,000 m × 2,000 m.

For observing spatial variability of the logarithmized population density, empirical semivariograms (Journel, Citation1989; Oh and Hwang, Citation2015; Pebesma, Citation2004; Zawadzki, Citation2011) were considered as follows:

(1)

where γ(h) is the semivariance for distance h; h the distance (m); N(h), the number of grid pairs within distance h; and Z(xi), Z(xi + h), population densities in points (grids) xi and xi + h (Bohling, Citation2005; Journel, Citation1989; Zawadzki, Citation2011).

Considering the population density as point data focused in centroids of particular pixels (regularly gridded data), the estimated at point is weighted using surrounding data points according to spatial covariance (Creesie, Citation1990; Deutsch, Citation2002; Zawadzki, Citation2011):

(2)

The ordinary kriging (OK) estimator is obtained by minimizing the variance over λ1, λ2, .., λN (Creesie, Citation1990; Stein, Citation2013; Zawadzki, Citation2011):

(3)

applying the method of Lagrange multipliers. That is,

(4)
(5)

The sequential Gaussian simulation (SGS) is the most straighforward method of simulating a multivariate Gaussian field (Webster and Oliver, Citation2007). Geostatistical simulation uses previously kriged or simulated data (eq 3) to reproduce the covariance between all simulated values (eqs 4 and 5). The algorithm is typically used to generate a set of realizations based on the same input distribution functions, heterogeneity, and the semivariogram (Deutsch, Citation2002; Soltani et al., Citation2013). Each realization defines particular spatial distribution of population density. The general scheme for the SGS (Manchuk and Deutsch, Citation2012) is listed next:

  1. Generate a random path through a set of points.

  2. For each point:

    • Search for conditioning data and previously simulated values.

    • Krige for the conditional distribution.

    • Simulate a Gaussian value from the distribution (as in the preceding).

    • Add the result to the pool of previously simulated values.

  3. Return to step 1 to generate more realizations.

Estimation of parameters is performed using the gstat package from the statistical software R (Bivand et al., Citation2013; Pebesma, Citation2004; Pebesma and Graeler, Citation2015).

Results and discussion

The population density hot spot is located in the central part of considered region (), where the SHS system is strongly developed (). Studies by Meeker (Citation1985), Persson (Citation2015), Persson and Werner (Citation2010), and Volkova et al. (Citation2012) confirm that the population density is one of the main parameters affecting the linear heat density, which is the main parameter characterizing district heating systems (Nussbaumer and Thalmann, Citation2016).

The suggestion that domestic (residential) air emissions at the most populated areas are proportional to the population density (Fameli and Assimakopoulos, Citation2016; Kryza et al., Citation2010; Maes et al., Citation2009; Markakis et al., Citation2010, Citation2013) can be misleading. Due to the occurring heating network infrastructure, a substantial part of the settled population does not use individual sources of heat (e.g., small domestic boilers and stoves). The emissions estimated using top-down approach should be reduced within the area supplied with the heat energy from the district heating system.

Earlier investigations (Nikodem, Citation2008; Plebankiewicz and Jankowski, Citation2007) emphasized that the SHS system supplies heat to approximately 50% inhabitants. That means the outlying, hot-spot values should be corrected. The exact area of correction is determined using the SGS.

Variography

To determine the area covered by the SHS infrastructure, spatial continuity of population density was analyzed. To mitigate significantly outlying values, the initial population density was logarithmized (log10); however, the variogram cloud () did not suggest global outliers (Haslett et al., Citation1991; Ploner, Citation1999). The most “trend-free” direction was observed for azimuth 90º (E) (, top left).

Figure 4. (a) Variogram cloud for logarithmized (log10) population density in the KUZ: semivariance = (log10(inhabitants km−2))2. (b) Anisotropy of the population density in the KUZ, from left to right, bottom: 0°, 45°, top: 90°, 135°.

Figure 4. (a) Variogram cloud for logarithmized (log10) population density in the KUZ: semivariance = (log10(inhabitants km−2))2. (b) Anisotropy of the population density in the KUZ, from left to right, bottom: 0°, 45°, top: 90°, 135°.

The theoretical variogram models considered for the data are exponential (eq 6), spherical (eq 7), and Gaussian (eq 8). The h parameter represents the lag distance, a practical range, and c sill:

(6)
(7)
(8)

Variogram models applied in this paper were taken into account earlier in work by Oh and Hwang (Citation2015), for purposes of investigation of population density spatial distribution.

To choose the best fitted semivariogram model and its parameters (range, sill, and nugget), we compared the sums of squared errors (SSE) for each model. The least SSE value indicates the best fit. Results are presented in and .

Table 1. Semivariogram models for logarithmized population density in the KUZ.

Figure 5. The theoretical semivariogram models for population density in the KUZ: exponential (solid line), spherical (dashed), and Gaussian (dotted). Semivariance = (log10(inhabitants · km−2))2.

Figure 5. The theoretical semivariogram models for population density in the KUZ: exponential (solid line), spherical (dashed), and Gaussian (dotted). Semivariance = (log10(inhabitants · km−2))2.

The value of the nugget effect for Spherical and Gaussian variograms was significantly bigger than for the exponential, which suggests they should not be used as the best fits (Zhou et al., Citation2013). On the other hand, the SSE value for the exponential model was significantly smaller. Comparing this result with the previous studies (Oh and Hwang, Citation2015; Zhang et al., Citation2016), the exponential variogram was accepted for investigation of spatial continuity of the population density.

Ordinary kriging estimation and the simulation (SGS)

Ordinary kriging for logarithmized population density in the KUZ is given in . Clearly, the observable population density hot spot coincides with the locations supplied with heat and hot water from the SHS network (see and ). The OK estimation determines the most probable spatial distribution of population density (in this case, logarithmized values). The result can help to determine the area where the population density should be reduced by nearly 50% (Plebankiewicz and Jankowski, Citation2007). Similarly to the study by Leopold et al. (Citation2012), the Gaussian simulation was used to determine hot-spot values.

Figure 6. (a) Ordinary kriging (OK) for the logarithmized population density in the KUZ: estimation. (b) OK variance.

Figure 6. (a) Ordinary kriging (OK) for the logarithmized population density in the KUZ: estimation. (b) OK variance.

For investigation of population density spatial distribution, the sequential Gaussian simulations were generated. presents 25 realizations of the logarithmized (log10) population density’s spatial distribution as well as the average of the realizations ().

Figure 7. (a) Twenty-five realizations of the SGS of logarithmized population density in the KUZ. (b) The mean value of the realizations.

Figure 7. (a) Twenty-five realizations of the SGS of logarithmized population density in the KUZ. (b) The mean value of the realizations.

Basing on the data derived from the earlier analysis (Nikodem, Citation2008; Plebankiewicz and Jankowski, Citation2007), the logarithmized values of population density were reduced by log102 (the initial data on population density are reduced by half). Using information derived from experts’ analysis also from the national statistics (CSO, Citation2015) the hot-spot values were assessed as bigger than 2.75 (102.75 ≈ 562.3 [inhabitants]). Analyzing spatial distribution of the population density in Poland (Zasina and Zawadzki, Citation2016), this is the approximated value, connected with supplying 50% of dwellers with heat and hot water. The final results (direct values of population density) are shown in .

Figure 8. Population density (inhabitants km−2) in the KUZ as the emission surrogate, corrected using the SGS realizations.

Figure 8. Population density (inhabitants km−2) in the KUZ as the emission surrogate, corrected using the SGS realizations.

Using the available statistics, including the official database (CSO, Citation2015), and Polish data on air emissions submitted for international purposes (NCEM, Citation2017), it was found that the top-down estimated annual value of mercury (Hg) released into the air (annual emission flux) from residential combustion sources located in the analyzed area (KUZ) is nearly 83 kg. Taking advantage of the estimated emission data, we prepared two spatial distributions of the Hg air emission, that is, using a “standard approach” based only on the population density (Gallego, Citation2010) (), and an improved approach using the geostatistical techniques (SGS simulations, ). The data distributed using the standard approach indicate that the emission’s hot spot is mostly concentrated in the center of the area, underestimating values located at the peripheries, where district heating is not available (see ).

Figure 9. (a) Annual Hg air emission from residential combustion (kg), distribution based on population density, Gauss–Krüger projection. (b) Annual Hg air emission from residential combustion (kg), distribution based on geostatistical simulations, Gauss–Krüger projection.

Figure 9. (a) Annual Hg air emission from residential combustion (kg), distribution based on population density, Gauss–Krüger projection. (b) Annual Hg air emission from residential combustion (kg), distribution based on geostatistical simulations, Gauss–Krüger projection.

The geostatistically improved spatial distribution () increases the importance and emission from less populated areas, situated in the rim of the region, primarily. The comparison of values in particular grids is given in the histograms in . The proposed methodology decreases the mercury emission associated with the most populated areas.

Figure 10. (a) Histogram of annual Hg air emission from residential combustion (kg), distribution based on population density. (b) Histogram of annual Hg air emission from residential combustion (kg), distribution based on geostatistical simulations.

Figure 10. (a) Histogram of annual Hg air emission from residential combustion (kg), distribution based on population density. (b) Histogram of annual Hg air emission from residential combustion (kg), distribution based on geostatistical simulations.

Conclusions

The methodology presented in this paper uses spatial distribution of the population density, as well as indirect information about heating network development (Nikodem, Citation2008; Plebankiewicz and Jankowski, Citation2007), to improve spatial disaggregation of air pollutant emission estimated using a top-down approach (Maes et al., Citation2009).

The presented algorithm changes the spatial distribution of population density, treated as an air emission surrogate, using local reclassification, which decreases the value (weight) of emission surrogate at more populated areas and increases them in less populated areas.

The applied methodology is based on geostatistical determination (variography), ordinary kriging (OK) estimation, and geostatistical simulation (SGS). The exponential variogram model is chosen as the most effective for population density. It confirms the results of previous studies in this matter (Oh and Hwang, Citation2015; Zhang et al., Citation2016; Zhou et al., Citation2013).

Carrying out this analysis, it is found that population density is an appropriate emission surrogate (Maes et al., Citation2009); however, it can be improved, considering spatial variability of this value. This kind of correction can be introduced easily and can help to analyze in detail emission impacts, air emissions spatial distributions, and further to improve emission inventories for air quality modeling.

The presented analysis is suitable for all kind of air pollutants or greenhouse gases emissions estimated using a top-down approach. Our results should be useful for policymakers and scientists working with analysis on various levels.

Funding

The authors thank the China Section of the Air & Waste Management Association for the generous scholarship they received to cover the cost of page charges and make the publication of this paper possible.

Additional information

Funding

The authors thank the China Section of the Air & Waste Management Association for the generous scholarship they received to cover the cost of page charges and make the publication of this paper possible.

Notes on contributors

Damian Zasina

Damian Zasina works in the Institute of Environmental Protection–National Research Institute, Warsaw, Poland. He is also a PhD student and teaching assistant in the Faculty of Building Services, Hydro and Environmental Engineering, at the Warsaw University of Technology. His research interests focus on using of geospatial methods for disaggregation and uncertainty analysis in air emission inventories.

Jaroslaw Zawadzki

Jaroslaw Zawadzki is a professor of statistics and geostatistics in the Faculty of Building Services, Hydro and Environmental Engineering, at the Warsaw University of Technology. He specializes in environmental studies and production engineering. His research interests focus on environmental risk assessment, environmetrics, geostatistics, remote sensing, and environmental

physics.

References

  • Bivand, R.S., E. Pebesma, and V. Gómez-Rubio. 2013. Applied spatial data analysis with R. Springer–Verlag, New York, NY. ISBN 978-1-4614-7617-7. doi:10.1007/978-1-4614-7618-4.
  • Bohling, G. 2005. Introduction to geostatistics and variogram analysis. http://people.ku.edu/~gbohling/cpe940/Variograms.pdf (accessed September 21, 2015).
  • Boychuk, P., K. Boychuk, Z. Nahorski, and J. Horabik. 2012. Spatial inventory of greenhouse gas emissions from the road transport in Poland. ECONTECHMOD 1(4):9–15.
  • Bun, R., K. Hamal, M. Gusti, and A. Bun. 2010. Spatial GHG inventory at the regional level: Accounting for uncertainty. Climatic Change 103:227–44. doi:10.1007/s10584-010-9907-5
  • Creesie, N. 1990. The origins of kriging. Math. Geol. 22(3):239–52. doi:10.1007/BF00889887.
  • CSO 2015. Local data bank of the Polish Central Statistical Office. http://www.stat.gov.pl/bdl (accessed September 21, 2015).
  • Deutsch, C.V. 2002. Geostatistical reservoir modeling. Oxford, UK: Oxford University Press.
  • European Environment Agency. 2016. EMEP/EEA air pollutant emission inventory guidebook. EEA. Luxembourg: Publications Office of the European Union. doi:10.2800/247535.
  • Fameli, K.-M., and V.-D. Assimakopoulos. 2016. The new open Flexible Emission Inventory for Greece and the Greater Athens Area (FEI-GREGAA): Account of pollutant sources and their importance from 2006 to 2012. Atmos. Environ. 137:17–37. doi:10.1016/j.atmosenv.2016.04.004.
  • Gallego, F.J. 2010. A population density grid of the European Union. Population Environ. 31:460–73. doi:10.1007/s11111-010-0108-y.
  • Haslett, J., R. Bradley, P. Craig, A. Unwin, and G. Wills. 1991. Dynamic graphics for exploring spatial data with application to locating global and local anomalies. Am. Stat. 45(3):234–42. doi:10.2307/2684298
  • Horabik-Pyzel, J., and Z. Nahorski. 2016. Uncertainty of spatial disaggregation procedures: Conditional autoregressive versus geostatistical models. Proc. Fed. Conf. Comput. Sci. 8:449–57.
  • Intergovernmental Panel on Climate Change. 2006. 2006 IPCC guidelines for national greenhouse gas inventories. IPCC, Hayama, Kanagawa, Japan. http://www.ipcc-nggip.iges.or.jp/public/2006gl (accessed September 21, 2016).
  • Journel, A.G. 1989. Fundamentals of geostatistics in five lessons, vol. 8. Washington, DC: American Geophysical Union. http://pbadupws.nrc.gov/docs/ML0227/ML022770097.pdf (accessed July 5, 2016).
  • Kryza, M., M. Werner, M. Blas, A.J. Dore, and M. Sobik. 2010. The effect of emission from coal combustion in nonindustrial sources on deposition of sulfur and oxidized nitrogen in Poland. J. Air Waste Manage. Assoc. 60:856–66. doi:10.3155/1047-3289.60.7.856.
  • Kuenen, J.J.P., A.J.H. Visschedijk, M. Jozwicka, and H.A.C. Denier van der Gon. 2014. TNO-MACC_II emission inventory; A multi-year (2003–2009) consistent high-resolution European emission inventory for air quality modelling. Atmos. Chem. Phys. 14:10963–76. doi:10.5194/acp-14-10963-2014.
  • Leopold, U., G.B.M. Huevelink, L. Drouet, and D.S. Zachary. 2012. Modelling spatial uncertainties associated with emission disaggregation in an integrated energy air quality assessment model. In International Congress on Environmental Modelling and Software. Managing resources of a limited planet: Pathways and visions under uncertainty, Sixth biennial meeting, Leipzig, Germany, ed. S. Lange R. Seppelt, A.A. Voinov, and D. Bandkamp. pp. 966–72. ISBN 978-88-9035-742-8.
  • Maes, J., J. Vliegen, K. Van de Vel, S. Janssen, F. Deutsch, K. De Ridder, and C. Mensink. 2009. Spatial surrogates for the disaggregation of CORINAIR emission inventories. Atmos. Environ. 43(6):1246–54. doi:10.1016/j.atmosenv.2008.11.040.
  • Manchuk, J.-G., and C.V. Deutsch. 2012. Implementation aspects of sequential Gaussian simulation on irregular points. Comput. Geosci. 16:625–37. doi:10.1007/s10596-012-9274-6.
  • Markakis, K., A. Poupkou, D. Melas, and C. Zerefos. 2010. A GIS based anthropogenic PM10 emission inventory for Greece. Atmos. Pollut. Res. 1:71–81. doi:10.5094/APR.2010.010.
  • Markakis, K., E. Katragkou, A. Poupkou, and D. Melas. 2013. MOSESS: A new emission model for the compilation of model-ready emission inventories? Application in a coal mining area in northern Greece. Environ. Model. Assess. 18:509–21.
  • Meeker, D.O. 1985. District heating and cooling in the United States. Prospects and issues. Washington, DC: National Academy Press. ISBN 0-309-03537-6.
  • NCEM. 2017. Convention on long-range transboundary air pollution. Submissions 2017. http://ceip.at/ms/ceip_home1/ceip_home/status_reporting/2017_submissions (accessed March 1, 2017).
  • Nikodem, W. 2008. Rozwazania nad kompleksowa modernizacja systemu zaopatrzenia w cieplo Konurbacji Slaskiej. Energetyka 5:381–88[in Polish].
  • Nussbaumer, T., and S. Thalmann. 2016. Influence of system design on heat distribution costs in district heating. Energy 101:496–505. doi:10.1016/j.energy.2016.02.062
  • Oh, D.-R., and C.S. Hwang. 2015. A comparative study on the spatial statistical models for the estimation of population distribution. J. Korean Soc. Surveying Geodesy Photogram. Cartogr. 3(3):145–53. doi:10.7848/ksgpc.2015.33.3.145.
  • Pebesma, E. 2004. Multivariable geostatistics in S: The gstat package. Comput. Geosci. 30:683–91. doi:10.1016/j.cageo.2004.03.012.
  • Pebesma, E., and B. Graeler. 2015. Package ‘gstat’. https://cran.r-project.org/web/packages/gstat/gstat.pdf (accessed July 5, 2016).
  • Persson, U. 2015. District heating in future Europe: Modelling expansion potentials and mapping heat synergy regions. PhD thesis, Chalmers Institute of Technology, Gothenburg, Sweden.
  • Persson, U., and S. Werner. 2010. Effective width: The relative demand for district heating pipe lengths in City Areas. Presented to 12th International Symposium on District Heating and Cooling. Tallinn, Estonia. September 2010.
  • Plebankiewicz, M., and A. Jankowski. 2007. Cieplo dla aglomeracji miast slaskich do wsparcia z funduszy unijnych. Wokol Energet. 3:1–4. http://www.cire.pl/pliki/2/AGLOM.pdf ( accessed September 21, 2015).
  • Ploner, A. 1999. The use of the variogram cloud in geostatistical modelling. Environmetrics 10:413–37. doi:10.1002/(ISSN)1099-095X
  • Soltani, F., P. Afzal, and O. Asghari. 2013. Sequential Gaussian simulation in the Sungun Cu porphyry deposit and comparing the stationary reproduction with ordinary kriging. Universal J. Geosci. 1(3):106–13. doi:10.13189/ujg.2013.010210.
  • Stein, M.L. 2013. Interpolation of spatial data. Some theory for kriging. New York, NY: Springer-Verlag.
  • Volkova, A., V. Mashatin, A. Hlebnikov, and A. Siirde. 2012. Methodology for the improvement of large district heating networks. Environ. Climate Technol. 10:3945.
  • Wang, Y., and G. Li. 2016. Mapping urban CO2 emissions using DMSP/OLS ‘city lights’ satellite data in China. Environ. Plan. A 49:1–4.
  • Webster, R., and M.-A. Oliver. 2007. Geostatistics for environmental scientists. Chichester, UK John Wiley & Sons, Ltd. ISBN 978-0-470-02858-2.
  • Zasina, D., and J. Zawadzki. 2016. Spatial surrogate for air emissions from small residential combustion—Analysis using scarce top-down estimates. Presented at Systemy Wspomagania w Inzynierii Produkcji/Review of Problems and Solutions, September, Wisła, Poland, pp. 135–143.
  • Zawadzki, J. 2011. Metody geostatystyczne dla kierunków przyrodniczych i technicznych. Warsaw, Poland [in Polish].
  • Zhang, H., W. Qu, S. Niu, J. Qi, L. Ye, and G. Zhang. 2016. Geological statistics analysis of population distribution at township level in Henan Province, China. Int. Proc. Chem. Biol. Environ. Eng. 91:63–9.
  • Zhou, G., Y. Wu, X. Xie, D. Xing, F. Fan, and R. Yu. 2013. Research on characteristics in the spatial distribution of the whole county population of Sichuan province, China via ARCGIS analysis. In International conference on advances in social science, humanities, and management, pp. 1000–1006. Amsterdam, The Netherlands: Atlantis Press.

Reprints and Corporate Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

To request a reprint or corporate permissions for this article, please click on the relevant link below:

Academic Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

Obtain permissions instantly via Rightslink by clicking on the button below:

If you are unable to obtain permissions via Rightslink, please complete and submit this Permissions form. For more information, please visit our Permissions help page.