header-bear

Proudly serving Wesley Chapel, New Tampa and surrounding areas.

Increase in domestic electricity consumption from particulate air pollution

Last updated: 12-24-2020

Read original article here

Increase in domestic electricity consumption from particulate air pollution

Increase in domestic electricity consumption from particulate air pollution
Environmental impact
Abstract
Accurate assessment of environmental externalities of particulate air pollution is crucial to the design and evaluation of environmental policies. Current evaluations mainly focus on direct damages resulting from exposure, missing indirect co-damages that occur through interactions among the externalities, human behaviours and technologies. Our study provides an empirical assessment of such co-damages using customer-level daily and hourly electricity data of a large sample of residential and commercial consumers in Arizona, United States. We use an instrumental variable panel regression approach and find that particulate matter air pollution increases electricity consumption in residential buildings as well as in retail and recreation service industries. Air pollution also reduces the actual electricity generated by distributed-solar panels. Lower-income and minority ethnic groups are disproportionally impacted by air pollution and pay higher electricity bills associated with pollution avoidance, stressing the importance of incorporating the consideration of environmental justice in energy policy-making.
Download PDF
Main
Air pollution has resulted in many types of negative externalities, a situation that calls for policy intervention to address the associated damages. Policymakers and research are widely concerned with increases in mortality risk, which are direct damages induced by pollution, as well as co-damages in terms of other welfare losses. These damages are generated via different channels including the physical and mental health impact on human beings, decreases in labour productivity 1 , 2 , the decline in subjective well-being 3 , the harm of cognitive competence 4 , 5 , disturbances in ecosystem health 6 , a diminished value of local environmental amenities and properties 7 , increases in household medical expenditure and so on. Accurate assessment of such externalities is crucial to estimating the social cost of pollution for the design and evaluation of policies such as a Pigouvian tax imposed on polluters for such external costs or for pollution control, or a cap-and-trade programme that establishes a market issuing allowances to internalize such cost 8 . While direct pollution damages are often measured in existing studies, there are not many discussions in the literature about the magnitude of the co-damages. A key challenge to quantify these co-damages, however, is to understand the interactions among pollution, human behaviours 9 and technologies. People can mitigate exposure to environmental risks by taking various avoidance behaviours, such as adjusting outdoor activities 10 , 11 and purchasing face masks and air purification systems in the short term 9 , 12 , and migrating to new living locations in the longer term 13 . Avoidance behaviours alleviate the negative health impact of pollution 14 but come at a cost, for example, spending less time doing outdoor activities 10 , 15 , and may lead to further impacts such as increased energy consumption due to a shift from natural to mechanical ventilation 16 , and increased need for heating or air-conditioning or for other activities such as watching television 17 , 18 in residential buildings. Commercial buildings may also be affected via further complexities if individuals choose to work remotely due to air pollution to avoid exposure during commuting 19 . However, commercial buildings might have better indoor air quality due to better ventilation 20 so that people can stay in commercial buildings for longer period of time. These two effects can cancel out, and thus we hypothesize that air pollution does not have a statistically significant impact on commercial buildings as a whole. Such effects and the consequential extra environmental damage are, however, hardly addressed explicitly and quantitatively in current studies, which leads to biases in damage evaluation. Our paper fills in this gap in the literature.
While electricity demand is driven up by pollution-averting behaviours, air pollution can also reduce electricity supply. High concentrations of particulate matter reduce solar electricity generation due to the changed solar irradiance. The emission of aerosols can attenuate solar radiation by scattering and absorbing sunlight before it reaches the solar panel 21 , and thus reduces photovoltaic (PV) performance 22 , 23 . Large particles in particulate matter can also generate dust on top of solar panels. In areas with severe air pollution such as China, the potential of solar PV generation decreased on average by 11–15% between 1960 and 2015 (ref. 22 ); the decrease of point-of-array irradiance can even reach 35% in the most polluted areas 23 . Such interaction adds another dimension to the complexity of assessing pollution externalities. Existing studies take a predominantly engineering perspective that relies on computer simulations to calculate the change in solar irradiance due to air pollution or on field experiments to measure the changes in electricity generation of a few solar panels in response to air pollution. While providing critical estimation on the relationship between particulate pollution and solar electricity generation in certain refined meteorological and geographical conditions, these studies fall short in evaluating how much actual solar generation is affected at a large scale. Our paper contributes on empirical grounds and serves as a crucial reference for policy-making.
As pollution co-damages are closely related to both demand-side human behaviours and supply-side solar power generation, the distribution of these co-damages raises environmental justice concerns. Lower-income households or minority ethnic groups can be more vulnerable to the impact of air pollution. Individuals from these groups usually reside in locations with higher air pollution levels 24 . Moreover, they may live in affordable housing and buildings that are old, not well insulated and equipped with fewer energy-efficient appliances, all of which leads to higher energy-related expenditures 25 , 26 . Increased electricity bills due to more indoor hours, therefore, translates into a larger proportion of the household income for these groups, compared to their higher-income or non-minority counterparts. This constrains other essential expenditures such as on medical services by these lower-income and minority households, thus leading to further adverse health impacts 27 . Our analyses incorporate the equity aspects of pollution co-damages to provide necessary implications for policy design towards environmental justice.
This article demonstrates how the interactions among air pollution, human defensive behaviour and the energy supply system can influence the estimates of negative externalities caused by air pollution. Using consumer-level daily and hourly electricity consumption data and solar panel generation records in the city of Phoenix, Arizona, during the period 2013–2018, we show how particulate air pollution, indicated by concentrations of both PM10 and PM2.5 (that is, particulate matter 10 micrometres or less in diameter, and 2.5 micrometres or less, respectively), triggers consumer avoidance behaviours as well as lowers the generation of solar energy. Our sample covers 4,313 residential buildings and 17,422 commercial buildings. A variety of demographic and socio-economic characteristics are associated with the consumer dataset, based on which we further explore the heterogeneity of the co-damages associated with income and ethnicity. Estimates can be biased by endogeneity issues due to reverse causality (that is, air pollution induces changes in energy consumption as well as solar electricity generation, which in turn also affects the air quality) and missing variables (for example, unobservable characteristics of the local economy and physical environment can affect the air quality and energy consumption simultaneously). To address the endogenous biases, we use wind direction as an instrumental variable (IV) for pollution concentration. This IV has a direct impact on concentrations of pollutants but not on energy consumption, which creates variation in air quality that is exogenous to consumption, thus leading to an unbiased estimation of the pollutant coefficient. Our main results are based on daily average data. We also analyse hourly data to examine the intra-day heterogeneity in the impact on electricity usage. The study area of our analysis is the fifth-most populated city in the United States 28 and ranks among the top five most-polluted cities in the country 29 . This suggests that, even though they are based on a region of a developed country, our results can provide valuable insights and benchmark statistics when compared with studies in developing countries, with dense populations and low-ranked air quality.
Effect of air pollution on the demand sectors
Through an IV fixed-effects panel regression, we regress the individual household’s daily electricity consumption on air pollution level, while controlling for other confounding variables. Detailed modelling can be found in the Methods section. The validity of the IV estimation is also supported by the first-stage regression, which shows a significant positive correlation between the daily average cosine of the prevailing hourly wind direction angle and the concentration of air pollution, meaning that wind in the upwind direction of pollution sources would bring higher particulate concentration (columns 1 and 3 in Table 1 ). The considerable F statistics of far more than ten for testing the statistical significance of the excluded instrument indicate a strong IV in both the regressions for PM10 and PM2.5. We find that a higher concentration of particulate pollutant results in a statistically significant increase in residential electricity consumption. An increase of 1 µg m–3 in PM10 concentration raises the daily residential electricity consumption by 0.020 kW h (column 2 in Table 1 ). Residents are more sensitive to a change in PM2.5 concentration, as a 1 µg m–3 rise in PM2.5 concentration causes an 0.145 kW h (column 4 in Table 1 ) increase in daily electricity consumption. In this way, one more standard deviation of PM10 and PM2.5 would increase the daily residential electricity consumption by 0.85% and 1.74%, respectively, from the mean, based on the descriptive statistics in Supplementary Table 1 . Such effects are also seasonally heterogeneous (Supplementary Table 12 ) as the increased electricity consumption has a larger magnitude in the peak of summer (July and August), while the significance diminishes during the winter (November to April).
Table 1 Effect of air pollution on electricity consumption in residential buildings
Full size image
Effect of air pollution on the supply sector
We then used a similar panel IV regression to regress individual consumers’ daily solar electricity generation on air pollution level, while controlling for confounding variables ( Methods ). The IV of wind direction again proves powerful in explaining the variation of PM10 and PM2.5 with its positive significance in columns 1 and 3 in both Tables 3 and 4 , while the F statistics continue to verify it as strong. We find that particulate pollution also reduces the electricity generation of distributed-solar panels in both residential and commercial buildings. IV estimation shows that a 1 µg m–3 increase in PM10 concentration significantly reduces the electricity generated by solar panels by 0.435 kW h in residential buildings (column 2 in Table 3 ) and by 0.022 kW h in commercial buildings (column 2 in Table 4 ). PM2.5 has an even larger effect—a 1.888 kW h reduction per µg m–3 increase for residential buildings (column 4 in Table 3 ) and 0.093 kW h reduction per µg m–3 increase for commercial buildings (column 4 in Table 4 ). In terms of percentage change, one standard deviation increase of PM10 and PM2.5 would result in a 25.01% and 30.64% reduction, respectively, of solar electricity generation for residential buildings with solar panels from the mean solar electricity generation, and 0.13% and 0.15% reduction, respectively, for commercial buildings. The comparison also indicates that commercial buildings are much less affected if considering that the power of solar panels is on average larger in commercial buildings, referring to the descriptive statistics in Supplementary Tables 1 and 3 . A possible reason is that the solar panels in commercial buildings are better maintained, with dust cleaned off them in a timely manner.
Table 3 Effect of air pollution on solar energy generation in residential buildings
Full size table
Discussion
This study explores the co-damage of air quality degradation via human defensive behaviour on the demand side and the performance of clean energy techniques on the supply side, respectively. Our results show that particulate pollution, while exposing individuals to health risks with direct emissions, can further add to their loss with regenerative feedback, which boosts energy consumption due to longer times spent indoors and the downgraded performance of solar panels. While previous studies predominantly focus on the positive consequences of the defensive behaviours in alleviating health impacts 10 , 12 , this research shows the possible pathways in which air pollution generates extra damage by interacting with such defensive behaviours 9 . Our analysis also shows that residents from low-income or Hispanic groups are more heavily affected, highlighting the vulnerability of those of specific socio-economic status in responding to environmental change and the potential environmental justice issues that should be addressed by policy design 24 , 27 .
Several limitations should be noted. First, our analysis addresses the situation in the city of Phoenix, Arizona. In spite of its top rank for air pollution levels in US cities, the concentration of PM is still far less than that in many developing countries such as Mexico or China 34 , 35 . Meanwhile, response levels can also differ due to cultural differences. Therefore, our results should be extrapolated with caution. In addition, our dataset lacks information on specific household end-use activities (for example, heating and cooling, or air purification). Thus, we are not able to pinpoint exactly what appliance or appliances are more intensively used against higher particulate concentrations, for further details on the mechanisms that we discuss. We leave these for future research that draws on high-resolution data in various geographical areas.
Several critical policy implications stem from the findings of this research. First, when calculating the marginal damage factors from air pollution, policymakers need to explicitly consider the co-damages generated from the feedback of consumer behaviours and clean technology performance, which is insufficiently discussed in the current literature, as well as policy analysis and evaluation. Lack of consideration of these pollution co-damages will lead to an under-estimation of the welfare gains from pollution control policies. Our results also stress the necessity to investigate comprehensively the consequences of air quality alerting systems, for example, alleviated health risks 36 , changed automobile traffic flows as individuals endeavour to escape from pollution as a response 37 , decreased outdoor recreation 10 , 38 and so on. Second, the fact that air pollution disproportionally affects those of low socio-economic status threatens energy and environmental justice, and again stresses that air pollution control can not only result in health benefits as a whole, but also contribute to an equitable distribution of such benefit. The disproportional impact also highlights the importance of energy policies that can improve the home energy efficiency of lower-income and ethnic minority groups to accelerate the achievement of fairness and equity. Third, our findings provide one more justification for the need to clean the electricity grid and improve the efficiency of renewable energy generation techniques. In addition, the expansion of solar power should consider the effect of air pollution when setting reasonable development targets. The results comparing the impacts on commercial PV units and residential PV units suggest that there should be clear messages or incentives to communicate the importance of the cleaning and maintenance of PV units to residential consumers.
Methods
Data
The data were provided by Salt River Project, one of the two largest utility companies in Arizona. Hourly electricity consumption in kilowatt-hours was available for 4,313 residential units (spanning from May 2013 to April 2017) and 17,422 commercial units (spanning from May 2013 to April 2018). For the residential units in the sample, a Residential Equipment and Technology Survey was also conducted in 2014, which asked about detailed sociodemographic information, building characteristics, appliance and other energy technology attributes, and energy consumption behaviours. For the commercial units, a six-digit code in the North American Industry Classification System is available to identify the sector type of the building. We aggregate the electricity consumption to the daily level for analysis. The daily electricity price is constructed by taking the average of the hourly prices. For commercial consumers, both the electricity charge and demand charge are included as price variables. The zip code zone of each building is also available in the dataset, which enables a spatial match with the air quality and meteorological variables.
Salt River Project also has distributed-solar consumers in its service territory. These solar panels can be installed on the rooftop of buildings or can be ground-mounted. For each distributed-solar consumer, our dataset has information on the hourly electricity generated by the consumer’s solar panels, along with the installation dates of the solar panels. There are 260 residential distributed-solar consumers (6.03% of the residential sample) and 330 commercial distributed-solar consumers (1.89% of the commercial sample) in our dataset.
We combine meteorological observations from multiple sources. Records of air quality, including daily average concentrations of PM2.5 and PM10, are retrieved from pre-generated data files of the United States Environmental Protection Agency 39 . Climate factors including the daily average temperature, total precipitation and average wind speed are obtained from Global Surface Summary of the Day 40 . The hourly wind direction data come from the Environmental Protection Agency’s pre-generated data files. We obtain the solar irradiance data from the National Renewable Energy Laboratory’s National Solar Radiation Database 41 . For missing solar irradiance data for a given location in a given time period, we use the simulated solar irradiance by the National Renewable Energy Laboratory for a given day in that location in a typical meteorological year.
We adopt an inverse distance weighting interpolation that has commonly been used in previous literature 42 , 43 to match the air quality and meteorological records with the zip code zone of each building. First, the distance between each air quality monitoring station and the geometric centre of the zip code zone is calculated. Next, the daily records of all the stations less than 50 km away from the geometric centre are averaged together and weighted by their inversed distance to the centre. This weighted average is used as the matched air quality record for all the buildings within the zip code zone. The climate records are matched in a similar way. The inverse distance weighting is conducted in Stata 14.0 using the wtmean command with 34 meteorological stations and 67 air pollution monitoring stations. To test whether our analysis is sensitive to the radius of the inverse distance weighting procedure, we change the caliper to 10 km and 20 km and rerun the analysis, as shown in Supplementary Tables 16 – 19 (for 10 km) and Supplementary Tables 20 – 23 (for 20 km). The coefficients change only slightly in magnitude but their signs and statistical significance remain, indicating the robustness of our results.
Since datasets addressing individual travelling behaviour are rarely publicly available at the localized level for the study area, we resort to the COVID-19 Impact Analysis Platform by the University of Maryland 44 , 45 for a national-level exploration. Established for studies on COVID-19’s impact, this dataset includes the daily number of trips per person at the county level starting from 1 January 2020, which is further broken down into work and non-work trips. The information on trips comes from mobile device location data. Since the massive outbreak of COVID-19 in the United States took place no earlier than March, we adopt the records in January and February and match them with the air pollution and climate data from the above sources using a similar method.
Empirical strategies
We first estimate a generalized linear squared model on the panel dataset of residential and commercial units separately with the equation
$${\mathrm{Elec}}\_{\mathrm{Con}}_{it} = \beta _1{\mathrm{Pollution}}_{it} + {\mathbf{X}}_{{\boldsymbol{it}}} + \alpha _i + \tau _y + \delta _m + {\mathrm{Weekend}}_t + {\mathrm{Holiday}}_t + \varepsilon _{it}$$
(1)
where i indexes an individual residential or commercial consumer and t indexes the day of the sample. Elec_Conit refers to the daily electricity consumption of consumer i on day t. Pollutionit is the daily average concentration of either PM10 or PM2.5. Xit is a vector of control variables, including cooling degree days and heating degree days (estimated using daily average temperature), daily total precipitation, wind speed and electricity price (average daily electricity price for the residential consumers, and demand charge and energy charge for the commercial units). We also control for the concentration of ozone as another major pollutant that affects air quality and thus the outdoor activities of consumers. The variable αi is customer fixed effect, and it controls for the time-invariant attributes of the consumer such as square footage and the number of stories as well as environmental awareness of building occupants. The time fixed effects τy and δm include the year fixed effect and the month-of-year fixed effect. The time fixed effects capture the time-varying factors across years and seasons, such as economic development and change in local energy policies. Weekend and Holiday are dummy variables for holidays and weekends, respectively. The Holiday dummy is equal to one if the day belongs to the following US federal holidays: New Year’s Day, Martin Luther King Jr Day, Presidents’ Day, Memorial Day, Independence Day, Labor Day, Columbus Day, Veterans Day, Thanksgiving and Christmas. The εi,t is the error term. Standard errors are clustered at the building level. We are interested in β1, which indicates the electricity use increase per µg m–3 increase of particulate concentration, ceteris paribus.
We analyse how the impact of air pollution differs by different income groups. Using the available data on household characteristics, the sample is divided into three levels of per capita income: low, medium and high. The division, provided by the Pew Research Center, is based on the minimum household income level of different household sizes varying from one to five (US$24,042/34,000/41,641/48,083/53,759 for middle income, and US$72,126/102,001/124,925/144,251/161,277 for upper income in 2014 (ref. 46 )). Since the household size is recorded as 1.5, 3.5 and 5 persons, we take an average of the two adjacent minimum household income levels for the 1.5- and 3.5-person households.
We test whether and by how much the particulate pollution affects solar energy generation with,
$${\mathrm{Elec}}\_{\mathrm{Solar}}_{it} = \beta _1{\mathrm{Pollution}}_{it} + {\mathbf{X}}_{{\boldsymbol{it}}} + \alpha _i + \tau _y + \delta _m + {\mathrm{Weekend}}_t + {\mathrm{Holiday}}_t + \varepsilon _{it}$$
(2)
where Elec_Solarit refers to the daily electricity generated by solar for consumer i on day t, and the other terms are the same as in equation ( 1 ). The Xit is modified to adapt to factors that can affect the power generation of solar panels, including climate factors that can affect the performance of solar power (temperature, precipitation, wind speed and surface albedo) and electricity prices, which can affect the motivation of consumers to actively maintain solar panels in good condition (consumers are encouraged to do so if the price is higher). The distributed-solar consumers in our sample were on a net-metering plan under which they could sell excessive solar electricity at retail electricity prices.
The naïve general least squares (GLS) estimation (results shown in Supplementary Tables 6 – 9 ) suffers from endogeneity issues due to reverse causality and missing variables 47 , 48 . As air pollution changes the behaviour patterns and increases the energy consumption of consumers, the latter can result in more electricity generation and thus pollution emissions. Meanwhile, if consumers spend more time indoors, the demand for vehicle travel may also decrease and lead to reduced emissions from transportation 48 . Omitting such pathways would lead to a biased estimation of the effect of air pollution. Besides, air quality and individual socio-economic activities can be jointly affected by the same factors, such as the local economy and physical environment 47 . Since all such factors cannot be observed in our datasets, these missing variables could bias the estimation.
To address these issues, we resort to using wind direction for an IV estimation. Its validity has been verified by multiple existing air pollution studies 47 , 49 , 50 . The idea is that wind direction affects regional air quality as it transports pollutants in specific directions. As the wind direction fluctuates on a daily or even hourly basis, it can convert the study area upwind or downwind of the pollution. Other than this pathway, wind direction (while controlling for wind speed) can hardly affect electricity consumption or solar electricity generation, and thus can meet the exclusive restriction for a valid IV.
We use the daily average cosine of the angle between the prevailing wind direction and the hourly wind direction as our IV following the previous studies 47 , 51 with modifications to adapt to our daily-level data. We first plot the distribution of the hourly wind direction of all the climate stations to obtain the prevailing wind direction, which turns out to be 180°. We then calculate the cosine of the angle between each hourly wind direction observation and this prevailing direction, and finally obtain the daily average for each climate station, which matches with different zip code zones. In this way, we can conduct the first-stage regression before running equation ( 1 ) or ( 2 ) as
$${\mathrm{Pollution}}_{it} = \gamma _1{\mathrm{Wind}}\_{\mathrm{dir}}_{it} + {\mathbf{X}}_{{\boldsymbol{it}}} + \alpha _i + \tau _y + \delta _m + {\mathrm{Weekend}}_t + {\mathrm{Holiday}}_t + e_{it}$$
(3)
where Wind_dirit indicates the daily wind direction variable, eit is the error term and other terms are the same as in equations ( 1 ) and ( 2 ). The coefficient γ1, after we run the first-stage model, is statistically significant with an F value larger than ten, implying that the IV is relevant and strong. We then use the predicted values of pollution from equation ( 3 ) in the second-stage model when we run equation ( 1 ) or ( 2 ).
It should be noted that the maximum value of the electricity consumption of the commercial buildings in our sample is extraordinarily large (Supplementary Table 3 ). However, there is no way for us to rule out the possibility that this value is reasonable given the decent variation of daily electricity consumption in the commercial buildings that this value belongs to. Therefore, we keep these potential outliers for the main analysis but also rerun the regressions, dropping commercial buildings with a maximum daily electricity consumption over 500 kW h and 1,000 kW h. The results provided in Supplementary Tables 10 and 11 show that our key results remain robust after the change. Also, about 10% of the buildings have a constant daily electricity consumption of zero in the raw data. We regard them as shut-down buildings and remove them from our sample.
We further test how air pollution affects residential and commercial electricity consumption at the hourly level. The identification is similar to equation ( 1 ) but using the matched hourly data of electricity use and air quality (lagged for one hour). The electricity consumption and solar electricity generation of one particular hour will not influence the air quality of the previous hour, and thus there is no reverse causality issue. In addition, such an immediate hourly reaction of building energy use will not lead to an immediate change (within the same hour) in local PM pollution levels for the following reason. The hourly change in building electricity consumption leads to an hourly change in electricity generated at power plants. The coal-fired power plants surrounding the Phoenix metropolitan area are all located at least 100 miles away. This implies that the transmission of the PM pollution from these power plants to Phoenix will take time (considering that the average wind speed in Arizona cities is less than 23 miles per hour and the average wind speed in our sample is 2.66 m s–1 or 6 miles per hour), and thus will not influence the local PM pollution within an hour. The notable hourly variation in local PM pollution (such as in morning hours and late afternoon hours) in Arizona mostly comes from other sources such as motor vehicles and road dust, instead of from power plants, based on the study by Clements et al. (ref. 52 ). As a result, the hourly change in building energy consumption will not alter local PM pollution in the Phoenix metropolitan area immediately.
To examine whether individuals stay at home instead of commuting to work on polluted days, we conduct a regression analysis on personal trips with,
$${\mathrm{Trip}}_{jt} = \beta _1{\mathrm{Pollution}}_{jt} + {\mathbf{X}}_{{\boldsymbol{jt}}} + \pi _j + \delta _m + {\mathrm{dow}}_t + \varepsilon _{jt}$$
(4)
where Tripjt indicates the trips per person in county j on day t; πj and dowt denote the county and day-of-week fixed effects; and the other terms are similar to those in equations ( 1 ) and ( 2 ) but at the county level. On the basis of regressions using the total trips, we further test the effect of pollution concentration on the work and non-work trips. Due to a similar source of endogeneity, we are instrumenting the pollution using the wind direction with,
$${\mathrm{Pollution}}_{jt} = \gamma _1{\mathrm{Wind}}\_{\mathrm{dir}}_{jt} + {\mathbf{X}}_{{\boldsymbol{jt}}} + \pi _j + \delta _m + {\mathrm{dow}}_t + e_{jt}$$
(5)
where Wind_dirjt indicates the daily wind direction variable for county j on day t, and the other terms are the same as in equation ( 4 ). We calculate the daily average cosine of the angle between the prevailing wind direction and the hourly wind direction as our IV in a similar way to that described above. The prevailing wind direction is retrieved from the median of the wind angle of each county during the study period.
Data availability
Records of air quality and hourly wind direction were retrieved from pre-generated data files of the United States Environmental Protection Agency at https://aqs.epa.gov/aqsweb/airdata/download_files.html . Climate factors were obtained from Global Surface Summary of the Day at ftp://ftp.ncdc.noaa.gov/pub/data/gsod/ . The solar irradiance data from the National Renewable Energy Laboratory’s National Solar Radiation Database is at https://maps.nrel.gov/nsrdb-viewer . The high-frequency electricity data are from the Salt River Project. As they are restricted by a non-disclosure agreement, they are available from the authors upon reasonable request and with permission from the SRP. The county-level trip data are available upon request from the COVID-19 Impact Analysis Platform of the University of Maryland at https://data.covid.umd.edu/about/index.html . Source data are provided with this paper.
Code availability
All data and models are processed in Stata 14.0. The figures are produced in R studio (based on R 3.6.1). All custom code is available on GitHub at https://github.com/hepannju/Increase-in-domestic-electricity-consumption-from-particulate-air-pollution .
References
1.
Graff Zivin, J. & Neidell, M. The impact of pollution on worker productivity. Amer. Econ. Rev. 102, 3652–3673 (2012).


Read the rest of this article here