Next Article in Journal
Improvements in Sample Selection Methods for Image Classification
Next Article in Special Issue
Characterizing Spatio-Temporal Dynamics of Urbanization in China Using Time Series of DMSP/OLS Night Light Data
Previous Article in Journal
On Recovering Missing Ground Penetrating Radar Traces by Statistical Interpolation Methods
Previous Article in Special Issue
Estimation of Gross Domestic Product Using Multi-Sensor Remote Sensing Data: A Case Study in Zhejiang Province, East China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Industrial Wastewater Discharge Retrieval Based on Stable Nighttime Light Imagery in China from 1992 to 2010

1
College of Resources and Environmental Science, Nanjing Agricultural University, Nanjing 210095, China
2
College of Geography, Nanjing Normal University, Nanjing 210023, China
3
Nanjing Institute of Geography & Limnology, Chinese Academy of Science, Nanjing 210008, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2014, 6(8), 7566-7579; https://0-doi-org.brum.beds.ac.uk/10.3390/rs6087566
Submission received: 14 May 2014 / Revised: 8 August 2014 / Accepted: 11 August 2014 / Published: 14 August 2014
(This article belongs to the Special Issue Remote Sensing with Nighttime Lights)

Abstract

: Industrial wastewater (IW) discharge, which is a known point source of pollution, is a major water pollution source. Increasing IW discharge has imposed considerable pressure on regional or global water environments. It is important to estimate the IW distribution in grid units to improve basin-scale hydrological processes and water quality modeling. For the first time, we use the nighttime light imagery produced by the Defense Meteorological Satellite Program’s Operational Linescan System (DMSP/OLS) to estimate the spatial and temporal variations in the IW distribution from 1992 to 2010 in China. The digital number values per unit area (DNP) of each stable light image were calculated using nighttime light imagery and were regressed against the IW per unit area (IWP) to estimate the total industrial wastewater (TIW) for each province. The results indicated strong positive correlations between the DNP and the IWP for each province during different years. The fitted linear regression models were used to estimate IW discharge in China with reliable accuracy. The IW estimation using the satellite data was consistent with the statistical results. The results also revealed that the IW discharge coverage expanded, whereas the IW discharge intensity decreased from 1992 to 2010 in China.

1. Introduction

Economic development and population growth have caused depletion of global water resources and serious problems for water environments [13]. A key challenge to sustaining future water environments is estimating both current and future water pollution. Population growth and increasing food requirements have caused increasingly serious domestic sewage and agricultural diffuse source pollution [4]. Moreover, rapid economic development, especially in industrial enterprises, creates more industrial wastewater [5]. With increasing point and non-point source pollution, the threat to the global water environment has become more severe [6]. Water environmental problems are particularly challenging in China, because China has the largest population, a fast-growing economy, a rising water demand, relatively scarce water resources, a dated infrastructure and inadequate governance [7]. Understanding the spatial characteristics and differences in wastewater discharge is important for the development of a water conservancy plan in China.

Industrial wastewater (IW) has been the most important source of water pollution since the Industrial Revolution [8]. A country’s IW is the total volume from each province, which is typically determined by adding the wastewater discharge from each administrative region. Currently, IW data are obtained from official statistical data that are collected by different administrative divisions, i.e., the province, county, township, the prefecture under the provinces and the village under the townships. However, water environmental problems generally occur at the basin or watershed scales, which correspond to the best natural unit for performing water pollution research and management [9,10]. Because hydrological processes and the establishment of water quality occur within individual basins, a distributed model that is based on a raster grid or hydrological response units (HRUs) is the most efficient method for examining problems related to water environments [11,12]. IW statistical data that are based on an administrative unit cannot be applied to hydrological processes or be used in a water quality model on the basin scale. Therefore, the use of basin-scale statistical IW data is an important issue.

Considering that population, gross domestic product (GDP) or other socio-economic factors can be spatialized to grid cells retrieved from remote sensing data [13,14], we intended to establish a relationship between the IW statistical data and the remote sensing data. The Operational Linescan System (OLS) is a sensor controlled by the Defense Meteorological Satellite Program (DMSP). This instrument can detect nocturnal artificial lighting under clear night sky conditions without moonlight due to its low-light imaging capability [15]. The nighttime light (NTL) data collected by the DMSP/OLS measure the light on the Earth’s surface, such as that generated by human settlements, gas flares, fires and illuminated marine vessels [16]. Previous studies have shown that the DMSP/OLS nighttime light imagery is a reliable tool for monitoring the spatial distribution of human settlements and for evaluating socioeconomic parameters [1720]. Moreover, several studies have shown a strong relationship between the NTL data and important socioeconomic variables, such as urbanization dynamics [16,2124], urban population estimates [25,26], population density [27,28], economic activity [29], GDP [14,30], energy consumption [20,3133] and the distribution of released anthropogenic heat [34]. Recently, progress has been achieved in environmental issues, such as CO2 emissions [33], the global distribution of anthropogenic nitrogen oxide emissions [35] and the global water footprint [3], using DMSP/OLS nighttime light imagery. However, no study has applied the NTL data to examine industrial wastewater, the chemical oxygen demand (COD) or discharged ammonia nitrogen (NH3N) at national or global scales.

The primary objective of this study was to evaluate whether DMSP/OLS NTL data can be used as a potential remote sensing data source for estimating IW, COD and NH3N discharge spatial distributions in China. Statistical relationships between the NTL data digital number (DN) values and the industrial wastewater parameters were first determined for all provinces in China. Then, we interpreted the differences between the estimated IW discharge determined from the nighttime light imagery and the official statistical data. Finally, we analyzed the spatial and temporal variations in the industrial wastewater discharge from 1992 to 2010 in China.

2. Data and Methods

Nearly twenty years of DMSP/OLS stable light data (from 1992 to 2010) were used in this study. The data were obtained from version 4 of the NTL data and downloaded from the website of the National Geophysical Data Center (NGDC) at the National Oceanic and Atmospheric Administration (NOAA) [36]. The version 4 data contain 3000-km swath width NTL scan lines that are divided into an array of grids with a 0.55-km spatial resolution, which are aggregated and composited into 30 arc second grids. Each grid in the composite data contains a digital number (DN) that indicates the average nighttime light intensity observed within each year. Because the OLS passes over each location every day, multiple observations may be available for any single location over the course of a year. Several constraints were considered to identify the best quality nighttime light data for producing the composite image [37]. For example, ephemeral events, such as wildfires, were discarded. The final dataset contained lights from urban areas and other sites with persistent lighting, including gas flares. Background noise in the composite image was identified and replaced with a value of zero; the final DN values ranged from 0 to 63 [16]. Previous studies have shown that there was no substantial difference among the various satellite datasets in different years [34]. Therefore, it is feasible for the data from different nighttime light detection satellites to be applied in this study.

A series of datasets were utilized to analyze the relationship between the light DN values and the IW discharge in China. The amount of discharged industrial wastewater and the gross industrial production (GIP) data for each province from 1992 to 2010 were obtained from the China Statistical Yearbooks from 1993 to 2011 [38]. These yearbooks were compiled by the National Bureau of Statistics of China and contain detailed statistics on all societal aspects in China, including the economy, population and resources. The statistics were separately compiled for each province and district.

To obtain the correct area for analyzing the statistical relationship between the NTL and IW datasets, the ordinary spatial reference of the geographic coordinate system that was related to the WGS 84 datum was transformed to the Albers equal-area projection in China. Moreover, the spatial resolution of the NTL data was resampled to a 1 km × 1 km area. Vector data for the Chinese provincial boundaries were also transformed to an equal-area projection. The provincial boundary map was overlaid on the nighttime light images, and the DN values per unit area (DNP) for all provinces were calculated for each of the 19 years from 1992 to 2010. The total IW (TIW) discharge for each province divided by the area of the province represents the amount of discharged IW per unit area (IWP). Then, statistical relationships between the IWP and the DNP were analyzed. The relative difference (RD) between the predicted TIW (PTIW) and the statistical TIW (STIW) was calculated as follows to analyze the accuracy of the fitted models:

RD = ( PTIW STIW ) / STIW
where PTIW is the predicted total IW from the nighttime light image and STIW is the statistical total IW. For the COD and NH3N discharge calculation and analysis, similar methods were adopted. The acronyms of the variables involved in the calculations and analyses are shown in Table 1.

3. Results and Discussion

3.1. Relationship and Fitted Model between the DMSP/OLS Data and the Industrial Wastewater Discharge

The industrial wastewater discharge per unit area (IWP) from 1992 to 2010 was calculated by dividing the TIW by the area of each province in China. The NTL digital number values per unit area (DNP) for all of the provinces were calculated using the ArcGIS zonal statistics tool. The relationship between the IWP and the DNP for each year was analyzed. Moreover, a linear model was fitted to the data and used to estimate the industrial wastewater discharge from the NTL light imagery data.

The primary model criterion is that when the digital number (DN) of the nighttime image increases, the discharged IW should also increase. Therefore, the IW discharge was expected to increase with an increasing DN. Based on this assumption, different regression models, including the exponential, power function and quadratic polynomial models, were used. In addition, a statistical hypothesis test was conducted. A comparison among the different model results suggested that the linear model, i.e., y = kx + b, performed the best. Here, k is the slope and b is the intercept of the linear function.

Figure 1a–c depicts the linear relationships between the DNP and the IWP, the IWCODP and the IWNH3NP for 31 provinces in 2010, respectively. A high correlation was observed between the DNP and the IWP (Figure 1a) and between the DNP and the IWNH3NP (Figure 1c); the corresponding R2 values were 0.798 and 0.756, respectively. Although the correlation between the DNP and the IWCODP was the weakest among the three relationships (Figure 1b) with an R2 value of 0.512, all of the correlations were significant at the 0.01 level. These results of quantitative analysis between DNP and IWP are consistence with the current situation in China, because industry is mainly distributed in urban regions with higher NTL DN values.

Furthermore, the model coefficients and the fitted model parameters for the relationships between the DNP and the IWP, IWCODP and IWNH3NP for all the provinces were calculated for each year. Table 2 shows the Pearson correlation coefficient (r), the slope (k) and the determination coefficient of the linear regression models that were constructed between the DNP and the IWP, IWCODP and IWNH3NP for 31 provinces from 1992 to 2010.

Significant positive correlations between the DNP and the IWP, IWCODP and IWNH3NP are shown in Table 2. All of the Pearson correlation coefficients for the relationship between the DNP and the IWP exceeded 0.85, while the correlation coefficients range from 0.716 to 0.849 for the relationship between the DNP and the IWCODP in different years. Although the correlation coefficients between the DNP and the IWNH3NP were relatively low, i.e., ranging from 0.627 to 0.870, all of the values exhibited a significant correlation at the 0.01 level according to a two-tailed test. Generally, the linear relationships between the DNP and the IWP, the IWCODP and the IWNH3NP were significant in each year.

Based on the previous assumption and the analyses between the DNP and the IW discharge variables, the slope and determination coefficients of the linear regression functions were obtained for each year (Table 2). These parameters and functions can be used to estimate the IW spatial distribution from the NTL light data in China.

3.2. Differences between the TIW Predicted from the Nighttime Lights Data and the Statistical TIW

The IWP was estimated using the linear regression model (as Table 2 shown) for each grid in China. The total IW (TIW) for each province was calculated as the total of the values from all of the grids within each province using the ArcGIS zonal statistics tool. Then, the predicted TIW (PTIW) was compared with the official statistical TIW (STIW) for all of the provinces.

Table 3 shows the differences between the PTIW estimated from the nighttime lights data and the STIW for 31 provinces in 2010. As shown in Table 3, the linear model results were reasonable in most provinces, whereas significant errors occurred in specific areas. The relative difference ranged from 0 to 0.5 in 16 of the 31 provinces, whereas the model exhibited significant errors in five provinces (i.e., the difference rate exceeded one). Despite these errors, the DMSP/OLS data provided useful estimates of the IW, IWCOD and IWNH3N spatial distributions in China, exhibiting a total error of −8.2%. Theoretically, industrial wastewater discharge is closely correlated with industrial development and the DMSP/OLS NTL data. The large relative difference may primarily be related to the significant differences in economic development in different areas.

Previous studies have shown that the DMSP/OLS data are closely correlated with the social and economic activities, such as the GDP and the population distribution [14,19]. For regions in which economic development was accompanied by industrial wastewater discharge, Table 3 demonstrates that the linear regression model results are generally credible. However, for regions in which economic development did not significantly rely on industrial activities, a significant estimation error occurred in the model. For example, the estimates were more accurate for the provinces with a higher percentage of STIW in China, such as Shandong, Henan, Guangdong and Jiangsu. The areas in which the model results overestimated the actual statistical data were typically the provinces that exhibited primary and tertiary industrial growth; these activities are typically related to less industrial wastewater being discharged in the process of economic development. These areas included Xizang, Beijing, Hainan, Qinghai, Heilongjiang, Yunnan and Jilin. Although the IW discharge was overestimated, the percentages of the TIW in China are very low in these provinces. Moreover, the error was not obvious in the national IW discharge amount. The model underestimated the results for the developing areas in which industry is important and the production efficiency is relatively low, such as in Guangxi, Hunan and Fujian provinces. In these regions, the IW discharge per unit gross industrial production (GIP) was higher.

Consequently, regional differences in industrial productivity for economic development may be a major reason for the model errors. Industrial wastewater discharge is closely correlated with industry development and the DN values obtained from the DMSP/OLS NTL data. However, the relationship between the industrial wastewater discharge efficiency and the NTL DN value is complex.

3.3. Spatial and Temporal Variations in the Chinese IW Discharge Distribution

The spatial distribution of the grid-scale IW discharge from 1992 to 2010 in China is shown in Figure 2. The annual evolution of the IW discharge was comparable to economic development and urban expansion in China.

The IW discharge was spatially concentrated in eastern, northeastern and southeastern China, especially in the Jing-Jin-Ji urban circle, the Yangtze River Delta, the Pearl River Delta and other coastal cities of China. The results show that the inland cities appear as spots on the IW distribution map, whereas the eastern portion of China and the coastal cities exhibit an intermittent distribution pattern (Figure 2).

The increasing IW discharge area and decreasing IW discharge intensity (IW discharge quantity per km2) trends are clearly depicted in Figure 2. The IW distribution in 1992 (Figure 2a) suggests that the IW discharge area was not extensive, except in economically developed areas, in which the discharge intensity was high. Several years later, the IW discharge region expanded beyond the developed areas. A significant enhancement in the IW discharge in eastern China is shown in the IW discharge distribution map in 1998 (Figure 2b). The trend in the IW distribution in 2004 (Figure 2c) and 2010 (Figure 2d) in China suggests that the IW discharge region continued to expand. For the IWCOD and the IWNH3N distributions, the spatial patterns and the temporal trends were consistent with the IW distribution. The extent of the distribution clearly broadened from 2003 (Figure 2e,g) to 2010 (Figure 2f,h).

The temporal changes in the IW discharge intensity exhibit a good trend. Figure 2 shows that the trend in the IW distribution in China suggests that the IW discharge intensity gradually decreased, although the IW discharge region continued to expand. This trend is also shown in Table 2 in which the slope of the linear regression model decreases from 1992 to 2010. A decreasing slope indicates that the IW discharge intensity decreased even though the light DN values remained constant each year. The highest predicted IW discharge intensity decreased from 50.071 × 104 t∙km−2 in 1992 to 6.332 × 104 t∙km−2 in 2010. The highest IWCOD discharge intensity decreased from 12.815 t∙km−2 in 2003 to 4.429 t∙km−2 in 2010. Moreover, the highest IWNH3N discharge intensity also decreased from 1057.05 kg∙km−2 in 2003 to 562.79 kg∙km−2 in 2010. According to these findings, the total IW discharge increased in conjunction with economic growth. However, the IW discharge intensity decreased due to economic development.

Figure 3 shows the ratio of the TIW to the total DN (TDN) and the TIW to the GIP from 1992 to 2010 in China. An obvious decreasing trend was observed. The IW discharge quantity per light DN value decreased gradually over the 20-year period. The trend was similar to the variations in the GIP to the TIW ratio. Similarly, the IW discharge per unit GIP decreased in conjunction with increased production efficiency in China.

4. Conclusions

This study highlighted the possibility of using the Defense Meteorological Satellite Program’s Operational Linescan System (DMSP/OLS) nighttime light (NTL) imagery to retrieve both spatial and temporal variations in the industrial wastewater (IW) distribution from 1992 to 2010 in China. We found strong positive correlations between the digital number values per unit area (DNP) and the IW per unit area (IWP) for each province in different years. The Pearson correlation coefficients (r) between the DNP and the IWP are greater than 0.85, while the r ranges from 0.716 to 0.849 for the relationship between DNP and the chemical oxygen demand discharge amount per unit area (IWCODP). Moreover, the r between DNP and the ammonia nitrogen discharge amount per unit area (IWNH3NP) ranges from 0.627 to 0.870, and all of the values show a significant correlation at the 0.01 level.

Based on the fitted linear model between DNP and IWP, the DMSP/OLS NTL data provided a dependable estimation of the IW discharge with a total error of −8.2% in 2010. The estimation was more accurate for the provinces with a higher percentage of IW discharge in China, such as Shandong, Henan, Guangdong, Jiangsu, etc. Fitted linear regression models could be used to estimate the IWP, the IWCODP and the IWNH3NP from 1992 to 2010 in China with reliable accuracy. The amount of IW estimated using the satellite data was consistent with the statistical results.

The industrial wastewater discharge distribution maps for individual grid units, which were produced using the satellite data, were significantly more credible and detailed compared to the statistical maps. The IW discharge was spatially concentrated in eastern, northeastern and southeastern China, especially in the Jing-Jin-Ji urban circle, the Yangtze River Delta, the Pearl River Delta and other coastal developed cities, while the IW discharge is lower in undeveloped areas, such as northwestern and southwestern provinces in China. The grid-scale industrial wastewater discharge data retrieved from the NTL data can be used in large-scale basin water quality models. The results also revealed that the IW discharge coverage expanded while the IW discharge intensity decreased from 1992 to 2010 in China. The maximum of IWP changed from 50.071 × 104 t∙km−2 in 1992 to 6.332 × 104 t∙km−2 in 2010. The maximum of IWCODP decreased from 12.815 t∙km−2 in 2003 to 4.429 t∙km−2 in 2010, and the maximum of IWNH3NP also decreased from 1,057.05 kg∙km−2 in 2003 to 562.79 kg∙km−2 in 2010.

Future endeavors should improve the prediction accuracy by considering other factors that affect the IW discharge in the fitted regression models. The use of satellite imagery by integrating industrial wastewater data, domestic sewage and agricultural non-point source pollution is another challenge for water environment pollution.

Acknowledgments

We sincerely appreciate the editor and the five anonymous reviewers for their critical and valuable comments to help improve this manuscript. This study was jointly supported by the National Natural Science Foundation of China (No. 41171071), the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD), the Key Project of Chinese Academy of Sciences (KZZD-EW-10-4) and the Scientific Research Foundation for the Returned Overseas Chinese Scholars, State Education Ministry. Special thanks are extended to the Earth Observation Group, National Geophysical Data Center, and National Oceanic and Atmospheric Administration of the United States for providing the nighttime light imagery.

Author Contributions

Zhaofu Li and Hongyu Liu designed and performed the data analysis and authored the manuscript. Zhengqin Xiong and Hengpeng Li revised the manuscript and provided several important suggestions. Chuan Luo and Panpan Li collected the statistical data. All of the authors participated in editing and revising the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Vörösmarty, C.J.; Green, P.; Salisbury, J.; Lammers, R.B. Global water resources: Vulnerability from climate change and population growth. Science 2000, 289. [Google Scholar] [CrossRef]
  2. Rockström, J.; Lannerstad, M.; Falkenmark, M. Assessing the water challenge of a new green revolution in developing countries. Proc. Natl. Acad. Sci. USA 2007, 104, 6253–6260. [Google Scholar]
  3. Zhao, N.; Ghosh, T.; Currit, N.A.; Elvidge, C.D. Relationships between satellite observed lit area and water footprints. Water Resour. Manag 2011, 25, 2241–2250. [Google Scholar]
  4. Novotny, V. Diffuse pollution from agriculture—A worldwide outlook. Water Sci. Technol 1999, 39, 1–13. [Google Scholar]
  5. Liu, J.; Diamond, J. China’s environment in a globalizing world. Nature 2005, 435, 1179–1186. [Google Scholar]
  6. Loucks, D.P.; van Beek, E. Water Resources Systems Planning and Management—Facts about Water; UNESCO: Paris, France, 2005. [Google Scholar]
  7. Liu, J.; Yang, W. Water sustainability for china and beyond. Science 2012, 337, 649–650. [Google Scholar]
  8. Rogers, P.; Silva, R.; Bhatia, R. Water is an economic good: How to use prices to promote equity, efficiency, and sustainability. Water Policy 2002, 4, 1–17. [Google Scholar]
  9. Mitsch, W.J.; Day, J.W.; Gilliam, J.W.; Groffman, P.M.; Hey, D.L.; Randall, G.W.; Wang, N. Reducing nitrogen loading to the gulf of mexico from the mississippi river basin: Strategies to counter a persistent ecological problem ecotechnology—The use of natural ecosystems to solve environmental problems—Should be a part of efforts to shrink the zone of hypoxia in the gulf of Mexico. BioScience 2001, 51, 373–388. [Google Scholar]
  10. Withers, P.J.; Lord, E.I. Agricultural nutrient inputs to rivers and groundwaters in the UK: Policy, environmental management and research needs. Sci. Total Environ 2002, 282, 9–24. [Google Scholar]
  11. Krysanova, V.; Müller-Wohlfeil, D.-I.; Becker, A. Development and test of a spatially distributed hydrological/water quality model for mesoscale watersheds. Ecol. Model 1998, 106, 261–289. [Google Scholar]
  12. Di Luzio, M.; Srinivasan, R.; Arnold, J.G. A GIS-coupled hydrological model system for the watershed assessment of agricultural nonpoint and point sources of pollution. Trans. GIS 2004, 8, 113–136. [Google Scholar]
  13. Imran, M.; Stein, A.; Zurita-Milla, R. Investigating rural poverty and marginality in Burkina Faso using remote sensing-based products. Int. J. Appl. Earth Obs. Geoinf 2014, 26, 322–334. [Google Scholar]
  14. Wu, J.; Wang, Z.; Li, W.; Peng, J. Exploring factors affecting the relationship between light consumption and gdp based on DMSP/OLS nighttime satellite imagery. Remote Sens. Environ 2013, 134, 111–119. [Google Scholar]
  15. Elvidge, C.D.; Baugh, K.E.; Hobson, V.R.; Kihn, E.A.; Kroehl, H.W.; Davis, E.R.; Cocero, D. Satellite inventory of human settlements using nocturnal radiation emissions: A contribution for the global toolchest. Glob. Chang. Biol 1997, 3, 387–395. [Google Scholar]
  16. Zhang, Q.; Seto, K.C. Mapping urbanization dynamics at regional and global scales using multi-temporal DMSP/OLS nighttime light data. Remote Sens. Environ 2011, 115, 2320–2329. [Google Scholar]
  17. Small, C.; Elvidge, C.D. Night on earth: Mapping decadal changes of anthropogenic night light in Asia. Int. J. Appl. Earth Obs. Geoinf 2013, 22, 40–52. [Google Scholar]
  18. Small, C.; Pozzi, F.; Elvidge, C.D. Spatial analysis of global urban extent from DMSP-OLS night lights. Remote Sens. Environ 2005, 96, 277–291. [Google Scholar]
  19. Elvidge, C.D.; Baugh, K.E.; Kihn, E.A.; Kroehl, H.W.; Davis, E.R.; Davis, C.W. Relation between satellite observed visible-near infrared emissions, population, economic activity and electric power consumption. Int. J. Remote Sens 1997, 18, 1373–1379. [Google Scholar]
  20. Shi, K.; Yu, B.; Huang, Y.; Hu, Y.; Yin, B.; Chen, Z.; Chen, L.; Wu, J. Evaluating the ability of NPP-VIIRS nighttime light data to estimate the gross domestic product and the electric power consumption of China at multiple scales: A comparison with DMSP-OLS data. Remote Sens 2014, 6, 1705–1724. [Google Scholar]
  21. Yu, B.; Shu, S.; Liu, H.; Song, W.; Wu, J.; Wang, L.; Chen, Z. Object-based spatial cluster analysis of urban landscape pattern using nighttime light satellite images: A case study of China. Int. J. Geogr. Inf. Sci 2014, 5, 1–28. [Google Scholar]
  22. Wei, Y.; Liu, H.; Song, W.; Yu, B.; Xiu, C. Normalization of time series DMSP-OLS nighttime light images for urban growth analysis with Pseudo invariant features. Landsc. Urban Plan 2014, 128, 1–13. [Google Scholar]
  23. Sutton, P.C. A scale-adjusted measure of “urban sprawl” using nighttime satellite imagery. Remote Sens. Environ 2003, 86, 353–369. [Google Scholar]
  24. Shi, K.; Huang, C.; Yu, B.; Yin, B.; Huang, Y.; Wu, J. Evaluation of NPP-VIIRS night-time light composite data for extracting built-up urban areas. Remote Sens. Lett 2014, 5, 358–366. [Google Scholar]
  25. Yang, X.; Yue, W.; Gao, D. Spatial improvement of human population distribution based on multi-sensor remote-sensing data: An input for exposure assessment. Int. J. Remote Sens 2013, 34, 5569–5583. [Google Scholar]
  26. Amaral, S.; Monteiro, A.M.V.; Camara, G.; Quintanilha, J.A. DMSP/OLS night-time light imagery for urban population estimates in the Brazilian Amazon. Int. J. Remote Sens 2006, 27, 855–870. [Google Scholar]
  27. Zhuo, L.; Ichinose, T.; Zheng, J.; Chen, J.; Shi, P.J.; Li, X. Modelling the population density of china at the pixel level based on DMSP/OLS non-radiance-calibrated night-time light images. Int. J. Remote Sens 2009, 30, 1003–1018. [Google Scholar]
  28. Sutton, P.C.; Elvidge, C.; Obremski, T. Building and evaluating models to estimate ambient population density. Photogramm. Eng. Remote Sens 2003, 69, 545–553. [Google Scholar]
  29. Elvidge, C.D.; Imhoff, M.L.; Baugh, K.E.; Hobson, V.R.; Nelson, I.; Safran, J.; Dietz, J.B.; Tuttle, B.T. Night-time lights of the world: 1994–1995. ISPRS J. Photogramm. Remote Sens 2001, 56, 81–99. [Google Scholar]
  30. Triantakonstantis, D.; Stathakis, D. Cokriging areal interpolation for estimating economic activity using night-time light satellite data. Lect. Notes Comput. Sci 2014, 8582, 243–252. [Google Scholar]
  31. Chand, T.R.K.; Badarinath, K.V.S.; Elvidge, C.D.; Tuttle, B.T. Spatial characterization of electrical power consumption patterns over india using temporal DMSP-OLS night-time satellite data. Int. J. Remote Sens 2009, 30, 647–661. [Google Scholar]
  32. Letu, H.; Hara, M.; Yagi, H.; Naoki, K.; Tana, G.; Nishio, F.; Shuhei, O. Estimating energy consumption from night-time DMPS/OLS imagery after correcting for saturation effects. Int. J. Remote Sens 2010, 31, 4443–4458. [Google Scholar]
  33. He, C.; Ma, Q.; Li, T.; Yang, Y.; Liu, Z. Spatiotemporal dynamics of electric power consumption in Chinese mainland from 1995 to 2008 modeled using DMSP/OLS stable nighttime lights data. J. Geogr. Sci 2012, 22, 125–136. [Google Scholar]
  34. Chen, B.; Shi, G.; Wang, B.; Zhao, J.; Tan, S. Estimation of the anthropogenic heat release distribution in China from 1992 to 2009. Acta Meteorol. Sin 2012, 26, 507–515. [Google Scholar]
  35. Toenges-Schuller, N.; Stein, O.; Rohrer, F.; Wahner, A.; Richter, A.; Burrows, J.P.; Beirle, S.; Wagner, T.; Platt, U.; Elvidge, C.D. Global distribution pattern of anthropogenic nitrogen oxide emissions: Correlation analysis of satellite measurements and model calculations. J. Geophys. Res. Atmos 2006, 111. [Google Scholar] [CrossRef]
  36. NOAA National Geophysical Data Center. Version 4 DMSP-OLS Nighttime Lights Time Series. Available online: http://ngdc.noaa.gov/eog/dmsp/downloadV4composites.html (accessed on 17 March, 2014).
  37. Elvidge, C.D.; Erwin, E.H.; Baugh, K.E.; Ziskin, D.; Tuttle, B.T.; Ghosh, T.; Sutton, P.C. Overview of DMSP nightime lights and future possibilities. Proceedings of 2009 Joint Urban Remote Sensing Event, Shanghai, China, 20–22 May 2009; pp. 1665–1669.
  38. National Bureau of Statistics of China. China Statistical Yearbook (1993–2011); China Statistical Press: Beijing, China. (In Chinese)
Figure 1. (a) The relationship between the DNP and the IWP (for 31 provinces); (b) the relationship between the DNP and the IWCODP (for 31 provinces); (c) the relationship between the DNP and the IWNH3NP (for 31 provinces) in 2010.
Figure 1. (a) The relationship between the DNP and the IWP (for 31 provinces); (b) the relationship between the DNP and the IWCODP (for 31 provinces); (c) the relationship between the DNP and the IWNH3NP (for 31 provinces) in 2010.
Remotesensing 06 07566f1 1024
Figure 2. Estimated industrial wastewater discharge distribution in China. (ad) The predicted IWP (104 t/km2) distribution in 1992, 1998, 2004 and 2010, respectively; (ef) the predicted IWCODP (t/km2) distribution in 2003 and 2010, respectively; and (gh) the predicted IWNH3NP (kg/km) distribution in 2003 and 2010, respectively.
Figure 2. Estimated industrial wastewater discharge distribution in China. (ad) The predicted IWP (104 t/km2) distribution in 1992, 1998, 2004 and 2010, respectively; (ef) the predicted IWCODP (t/km2) distribution in 2003 and 2010, respectively; and (gh) the predicted IWNH3NP (kg/km) distribution in 2003 and 2010, respectively.
Remotesensing 06 07566f2a 1024Remotesensing 06 07566f2b 1024
Figure 3. Temporal trends in the TIW to the TDN ratio and the TIW to the GIP ratio from 1992 to 2010 in China.
Figure 3. Temporal trends in the TIW to the TDN ratio and the TIW to the GIP ratio from 1992 to 2010 in China.
Remotesensing 06 07566f3 1024
Table 1. Acronyms of the variables involved in the calculations and analyses. COD, chemical oxygen demand; NTL, nighttime light.
Table 1. Acronyms of the variables involved in the calculations and analyses. COD, chemical oxygen demand; NTL, nighttime light.
AcronymTIWPTIWSTIWIWPDNPTDN

Full NameTotal industrial water dischargePredicted TIW from nighttime light imageStatistical TIW by National Bureau of Statistics of ChinaIndustrial wastewater discharge per unit areaDigital number per unit areaTotal nighttime light digital number values

AcronymGIPIWCODIWNH3NIWCODPIWNH3NP

Full NameGross Industrial ProductionQuantity of COD in industrial wastewaterQuantity of ammonia nitrogen in industrial wastewaterCOD discharge amount per unit area or unit NTL light valuesAmmonia nitrogen discharge amount per unit area
Table 2. The Pearson correlation coefficient (r), linear regression model slope (k) and determination coefficient (R2) between the DNP and the IWP (unit: 104 t/km2), the IWCODP (unit: t/km2) and the IWNH3NP (unit: kg/km2) for 31 provinces from 1992 to 2010.
Table 2. The Pearson correlation coefficient (r), linear regression model slope (k) and determination coefficient (R2) between the DNP and the IWP (unit: 104 t/km2), the IWCODP (unit: t/km2) and the IWNH3NP (unit: kg/km2) for 31 provinces from 1992 to 2010.
Yearr between the DNP and the IWPk of the IWP Regression ModelR2 of the IWP Regression Modelr between the DNP and the IWCODPk of the IWCODP Regression ModelR2 of the IWCODP Regression Modelr between the DNP and the IWNH3NPk of the IWNH3NP Regression ModelR2 of the IWNH3NP Regression Model
19920.865 **0.8070.747------
19930.856 **0.7320.732------
19940.883 **0.5290.779------
19950.873 **0.5100.762------
19960.852 **0.5350.726------
19970.866 **0.4150.750------
19980.869 **0.3800.754------
19990.868 **0.3540.753------
20000.887 **0.2970.786------
20010.879 **0.2770.772------
20020.891 **0.2420.793------
20030.914 **0.2600.8350.849 **0.1950.7210.774 **15.990.599
20040.914 **0.2170.8350.819 **0.1530.6700.675 **11.100.455
20050.926 **0.1990.8570.793 **0.1570.6290.627 **11.760.393
20060.913 **0.1800.8330.790 **0.1340.6230.772 **11.980.595
20070.909 **0.1770.8260.789 **0.1260.6220.815 **11.450.664
20080.908 **0.1580.8190.763 **0.1030.5650.835 **9.9810.679
20090.873 **0.1420.7610.753 **0.0980.5670.809 **7.6930.654
20100.894 **0.1010.7980.716 **0.0620.5120.870 **8.7500.756

**Represents that the correlation is significant at the 0.01 level (two-tailed); “-” indicates no data; n = 31.

Table 3. Difference rate between the PTIW and the STIW for 31 provinces in China in 2010.
Table 3. Difference rate between the PTIW and the STIW for 31 provinces in China in 2010.
RegionsThe PTIW from Nighttime Light Images (104 t/year)STIW (104 t/year)Difference RatePercentage of the STIW in China (%)Gross Industrial Production (GIP) (108 ¥)The STIW per 108 Yuan of GIP (104 t/108 ¥)
Jiangsu193,209263,760−0.2711.121,753.9312.12
Zhejiang122,004217,426−0.449.214,297.9315.21
Shandong209,951208,2570.018.821,238.499.81
Guangdong169,951187,031−0.097.923,014.538.13
Guangxi57,405165,211−0.657.04511.6836.62
Henan137,566150,406−0.096.313,226.3811.37
Fujian61,161124,168−0.515.27522.8316.51
Hebei161,895114,2320.424.810,707.6810.67
Hunan45,01395,605−0.534.07343.1913.02
Hubei59,48794,593−0.374.07767.2412.18
Sichuan59,20793,444−0.373.98672.1810.78
Jiangxi40,40872,526−0.443.15122.8814.16
Liaoning101,17771,5210.413.09976.827.17
Anhui83,14570,9710.173.06436.6211.03
Shanxi92,87749,8810.862.15234.009.53
Shaanxi77,14745,4870.701.95446.108.35
Chongqing23,66845,180−0.481.94359.1210.36
Inner Mongolia51,15239,5360.291.76367.696.21
Heilongjiang128,01238,9212.291.65204.117.48
Jilin73,12238,6560.891.64506.318.58
Shanghai29,96636,696−0.181.57218.325.08
Yunnan59,21830,9260.911.33223.499.59
Xinjiang30,44025,4130.201.12592.159.80
Ningxia15,76421,977−0.280.9827.9126.55
Tianjin33,44419,6800.700.84840.234.07
Gansu24,66615,3520.610.61984.977.73
Guizhou24,71114,1300.750.61800.067.85
Qinghai−9,6169031−2.060.4744.6312.13
Beijing36,14181983.410.33388.382.42
Hainan19,49457822.370.2571.0010.13
Xizang−32,445736−45.080.0163.924.49
China2,179,3382,374,733−0.082100.0220,064.7710.79

Share and Cite

MDPI and ACS Style

Li, Z.; Liu, H.; Luo, C.; Li, P.; Li, H.; Xiong, Z. Industrial Wastewater Discharge Retrieval Based on Stable Nighttime Light Imagery in China from 1992 to 2010. Remote Sens. 2014, 6, 7566-7579. https://0-doi-org.brum.beds.ac.uk/10.3390/rs6087566

AMA Style

Li Z, Liu H, Luo C, Li P, Li H, Xiong Z. Industrial Wastewater Discharge Retrieval Based on Stable Nighttime Light Imagery in China from 1992 to 2010. Remote Sensing. 2014; 6(8):7566-7579. https://0-doi-org.brum.beds.ac.uk/10.3390/rs6087566

Chicago/Turabian Style

Li, Zhaofu, Hongyu Liu, Chuan Luo, Panpan Li, Hengpeng Li, and Zhengqin Xiong. 2014. "Industrial Wastewater Discharge Retrieval Based on Stable Nighttime Light Imagery in China from 1992 to 2010" Remote Sensing 6, no. 8: 7566-7579. https://0-doi-org.brum.beds.ac.uk/10.3390/rs6087566

Article Metrics

Back to TopTop