一、 研究背景
 
2004-2009年，中国陆续在东北地区、中部地区、内蒙古东部地区和汶川地震受灾严重地区开展了一系列增值税转型改革试点工作，将过去的生产型增值税彻底转型为消费型增值税。增值税转型改革是对税基作出调整，旨在消除固定资产投资重复征税、打通固定资产投资上下游增值税抵扣链条，降低企业投资税负成本，从而激发企业投资活力、鼓励企业进行设备更新和技术改造升级。因此，增值税转型改革将发挥典型税收激励作用；增值税转型改革激励了企业进行设备更新换代，完成技术改造升级，从而增强了企业的可持续发展能力。同时，增值税转型改革试点政策最先只在部分省份和地区实行，经过多年陆续在不同地区和省份的实践，于2009年推广到全国。因此，该政策为本文识别税收激励对企业污染排放的因果效应提供了外生冲击。
 
本文采用2000-2012年中国工业企业数据库和中国企业污染排放数据库的合并微观数据，借助2004-2009年逐步推开的增值税转型改革这一自然实验，采用多时点双重差分估计的方法实证研究税收激励如何影响企业的污染排放问题。相比已有文献，本文在如下方面做出贡献：①已有文献鲜有关注税收激励如何影响企业减排这一重要问题，因此，本文基于微观企业数据，精准识别了税收激励与企业减排之间的因果关系，并通过深入的理论分析，全面考查了税收激励如何影响企业污染排放的内在机理，这为更加精准的绿色税收工具使用，尤其是在协同环境保护与经济发展方面提供了切实可行的决策支持。②本文构建了一个税收激励下的企业污染排放决策模型，从而得以从企业投资行为角度，厘清税收激励通过何种机制作用于企业污染排放决策，即分解出研发创新效应、清洁生产效应和规模扩张效应三条传导路径，并通过扎实的实证研究进行了全面的检验和分析，在一定程度上丰富了现有企业污染决定因素的研究文献。
 
二、研究设计
 
(一)模型构建
 
本文参考Beck et al.(2010),将中国增值税转型改革作为自然实验，构建多时点双重差分模型，以识别税收激励对企业污染排放的因果效应：
 

 

其中，i代表企业，j代表行业，k代表城市，t代表年份，α为截距项，η_i为企业固定效应，θjt为行业-年份交互固定效应，δkt为城市-年份交互固定效应，μijkt为随机误差项。由于增值税转型改革在城市一行业层面实施，本文将标准误聚类到城市一行业维度。
 
被解释变量ln(pol)it代表取自然对数的企业污染排放规模，测度方法借鉴包群等(2013),将中国企业污染排放数据库中汇报的二氧化硫排放量、烟粉尘排放量、化学需氧排放量、氨氮排放量、工业废水排放量、工业废气排放量这六项污染物排放物指标降维，采用主成分分析法构建企业综合污染排放强度，该指标能够更为全面地评估企业的污染排放量。核心解释变量VATjkt代表该企业所在城市一行业是否施行增值税转型改革。控制变量