The Impact of Climate Change on Regional Government Expenditures: Evidence from Russia

This paper explores an almost untouched topic in the fast-growing climate econometrics literature—the implications of climate change for government expenditures. Using a rich sub-national dataset for Russia covering 1995–2009, we estimate the impacts of changes in climatic conditions through short-term variation and medium-term changes in average regional temperatures and precipitation. We show a strong and robust negative (but non-linear) relation between regional budget expenditures and population-weighted temperature. The results indicate that an increase in temperature results in a decrease in public expenditures and that the magnitude of this effect diminishes the warmer the region. Further, our results suggest that the benefits from warming accumulate and that adaptation measures could help leverage those benefits. The estimated decreases in regional government expenditure are, however, quite small. In mild warming scenarios, according to our estimates Russia saves roughly USD 3–4 billion in regional government expenditures between 2000 and the 2020s without undertaking adaptation measures, depending on the scenario. It should be noted that our results are estimated for a scenario of mild temperature increase (1–2 ∘\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{\circ }$$\end{document}C). Larger temperature increases are likely to have dramatic consequences e.g. from loss of permafrost and methane release that are impossible to predict with available historical data.


Introduction
There is a strong consensus that the earth is experiencing man-made climate change caused by such factors as accumulation of greenhouse gases in the atmosphere and black carbon deposition on land and sea surfaces. The recent report by the IPCC (2014) documents warming of the atmosphere and oceans, rising sea level, as well as the loss of snowpack and ice sheets. Climate change is likely to have major impacts on natural and human systems on all continents in coming years, eventually reflecting into countries' economies and consequently on their public finances. The academic discussion around fiscal impacts of these huge changes so far is largely policy-based and focused on two large topics: mitigation costs like emission taxes and direct budgetary effects of extreme weather events that are likely to become increasingly intense and frequent due to global warming (See Jones et al. 2012 for a survey).
Even modest and gradual climate disruption, however, may affect public finances through e.g. shifts in economic structures, weather-induced changes in public health, revised heating/cooling expenditures and adaptation costs related to public infrastructure. The research bias towards mitigation may arise from an assumption that adaptation is largely a private sector issue. For example, Tol (2005) argues that national governments and international organizations need not participate in climate change adaptation efforts due to the local nature of the problem.
On the other hand, governments, as providers of public goods and services and holders of broad powers, implicitly have a stake in adaptation as far as it impacts the ability of the government to function effectively and address market failures. Governmental adaptation costs may arise on many fronts, e.g. through public transport networks, developing public health responses or securing coastal areas. Notable uncertainties and imperfect information (e.g. in distinguishing weather variations from permanent climate shifts) could prevent efficient private-sector adaptation and response . Osberghaus and Reif (2010) discuss local externalities (e.g. overdrafts of groundwater from stepped up farm irrigation) and the production of local public goods (e.g. sea dikes) in response to climate change. Moreover, poor countries may lack private sector with adequate financial resources for adaptation due to financial market imperfections and other factors. Finally, governments may have to move ahead with adaptation measures even against some possible theoretically efficient allocation of responsibilities and costs, if political pressure from voters and interest groups becomes too compelling.
This study is an early step in filling the research gap in econometric analysis of public adaptation costs of climate change. By focusing on sub-national data for an enormous and climatically heterogeneous country, Russia, we highlight the effects caused by moderate changes in temperature and precipitation on regional government expenditure. We want to emphasize that our aim is not to achieve a complete picture of climate change impacts although temperature and precipitation are the key variables in climate change analysis and discussions. Instead, we use the two most widely used indicators (temperature and rainfall) to proxy for the magnitudes of climate change impacts as is often done in the climate econometrics literature (e.g. Schlenker and Roberts 2009). Thus, our results should be considered only to concern climate change with respect to changes in temperature and precipitation although we refer simply to climate change in our analysis for brevity.
Using annual data from 1995 to 2009 to analyse changes in regional government expenditure across 78 Russian regions, our estimations show a significant negative relation between temperature and expenditure per capita. Our main result is that temperature rise over the short term (a limited adaptation setting) reduces regional government expenditures in cold regions and that this effect attenuates in a non-linear manner the warmer the region. While the results become less trustworthy over a longer time horizon due to the scantiness of observations, we find evidence suggesting that weather benefits accumulate, i.e. warming tends to have a larger medium-term effect on decreasing regional government expenditures than in the short term. We also find some evidence that housing and communal expenditures are one mechanism through which climate and weather affect total expenditures. The benefits are fairly small although they might increase with proper adaptation. Under our mild warming scenario, Russia saves between just over USD 2 billion to USD 4 billion in regional government expenditures between 2000 and the 2020s without any adaptation measures. We find support to our results with the aggregate temperature variable from higher frequency degree days variables.
We believe Russia serves as a useful benchmark in quantifying potential fiscal effects of global changes in climate. It is the world's largest country in terms of land area and has a highly versatile climate that provides a good basis for empirical analysis. The public sector in Russia clearly also has a non-negligible effect on the economy with its regional government expenditures equalling almost 20 % of GDP. Further, the average temperature in Russia has increased considerably faster than global temperature in recent decades, suggesting that warming of climate is truly happening there. If moderate, but persistent, climate change has any fiscal effects, Russia would be a natural candidate for a place where they would materialize. One needs to note, however, that this discussion only reflects regional government expenditure impacts in the case of mild climate change and reactive adaptation. Radical changes (e.g. partial melting of the permafrost, which covers roughly two-thirds of Russia) would have large, hard-to-quantify effects.
The paper is structured as follows. In Sect. 2, we elaborate our choice of Russia as an ideal setting for examining expenditure impacts of climate change. Section 3 discusses the existing literature. Section 4 outlines the methodological issues on estimating climate change effects and presents our estimation strategy. Section 5 describes our dataset and the variables used. The estimation results are presented and discussed in Sect. 6. Section 7 concludes.

Temperature Changes and Climate Policy in Russia
Russia is an interesting subject for studying potential fiscal effects of climate change. First, the country's geographical size and versatile climatic conditions make it an ideal candidate for an empirical study of climate effects. Russia encompasses Arctic, sub-Arctic, moderate and subtropical climatic zones. The average annual temperatures differ by as much as 15 • C between the southwest and northeast corners of Russia (Kotlyakov 2002). The permafrost, as mentioned, covers roughly two-thirds of Russia's land area.  (1995)(1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009) in Russian regions, ( • C). Temperatures coded at the regional level are based on administrative borders. The approximate location of the Ural Mountains is shown with the black dashed line. The two regions are excluded from our analysis: the Chukotka region (square) as an outlier and the Chechnya Republic (dot) due to the unreliability of available data. Source: meteo.ru database of weather stations, authors' calculations. Figure 1 shows population-weighted temperature patterns among Russian regions based on our data (description in Sect. 5.2.). The general pattern of warm southern and western regions contrasts nicely with cold northern and far-eastern regions. This rich climatic variation also improves the identification properties of our climate data compared to climatically homogenous countries such as the UK or Japan.
Temperatures in Russia have been rising faster than the global average for several decades. Temperatures in the Arctic have been rising at almost double the global average. According to Russia's Federal Service for Hydrometeorology and Environmental Monitoring (Roshydromet 2008), the average annual temperature has increased by 1.29 • C during the last roughly one hundred years in Russia, while global warming for the same period was 0.74 • C. Furthermore, average warming in Russia has apparently intensified in recent decades; it was +1.33 • C for the period 1976-2006 and the trend has continued. These trend differences are visible in Fig. 2, which shows the simple 10-year moving average of annual temperature deviations from the 1951-1980 average globally and in Russia from 1910 to 2014. Due to its extreme northern location, temperature increases have clearly been more drastic in Russia than globally since the 1980s.
The expected climate change impacts for Russia are versatile as evidenced in the impact assessment by Roshydromet (2014). For example, growing season for crops will become longer but at the same time some pests are expanding northward and eastward. Higher thermal comfort is expected in Russian north but the adverse health effects of heat waves will likely increase. Energy demand for heating will decrease due to shorter heating season but correspondingly the need for cooling in the warm season will increase. If climate change and global warming are likely to have had an effect on government finances somewhere, Russia is a natural candidate.
Second, Russia plays a key political role in climate change negotiations. The world's fourth largest CO 2 emitter after China, the US and India (Turkowski 2012), Russia is also the largest national terrestrial carbon sink (Lioubimtseva 2010). Despite the strong evidence of climate warming and active involvement in international climate negotiations, however, Russia has  (2014) been reluctant to commit to concrete greenhouse gas emission reductions. The reticence to get involved with climate change mitigation reflects several factors. Climate scenarios for Russia, a relatively cold and dry country, at least superficially imply beneficial trends like increasing average temperatures and precipitation. Russia also hopes to reap big economic benefits from climate change such as the "Northern Sea Route" (once the legendary Northeast Passage), as well as offshore drilling opportunities and high-latitude agriculture. As a result, many Russian scientists are convinced Russia is a net winner from climate change. Even the International Panel of Climate Change (2007) has noted the initially positive economic effects for the Russian economy.
With low enthusiasm for taking the lead in dealing with climate change, Russia, which is classified as a high-income country by the World Bank, prefers to align itself with developing economies. Further, the Russian public generally opposes the use of public funds on climate change mitigation (Korppoo 2008). Nevertheless, GAIN's (2012) climate change vulnerability ranking indicates Russia faces serious challenges from climate change, ranking it 39 th most susceptible out of some 180 countries surveyed. Moreover, any gains from climate change will be unevenly distributed across Russia's regions. Many regions face devastating consequences over the longer term. Like the rest of the world, climate change is likely to increase the frequency of extreme weather events. Forest and peat fires could increase, as well as flooding, land erosion and public health challenges (Korppoo 2008;Turkowski 2012). Near-surface thawing of the permafrost due to climate warming could be catastrophic for Russia's northern infrastructure, which is largely built on permafrost foundations (Pynnönniemi 2012). Furthermore, thawing of the permafrost could lead to huge increase in earth-based methane emissions that would exacerbate climate change dramatically (Lioubimtseva 2010).
At the time of writing Russia has opted out from international climate agreements, choosing instead to establish its own program. President Vladimir Putin signed a decree in September 2013 committing Russia to keep its greenhouse gas emissions at least 25 % below the 1990 level by 2020. There are also other signs that Russia has increased interest in dealing with climate change. A presidential advisor on climate change has been appointed, new legislation enhancing energy efficiency and renewable energy usage has been accepted and Russia is planning its own carbon credit trading system. At the moment, Russia's targets could hardly be described as ambitious. Its 2011 carbon emission levels were roughly a third below the year 1990 baseline level. Furthermore, its outdated industrial base is highly energy inefficient and Russia's energy efficiency will automatically improve with almost any fixed capital investment. Thus, the current emission reduction goals will probably be met without any actual mitigation efforts. If Russia's greenhouse gas emissions were to increase at their recent pace, they would return to the 1990 level around 2025 (Lioubimtseva 2010;Turkowski 2012). The current economic recession has made it even easier to hit stated emission targets.

Russia's Regional Public Finances
Russia is quite heterogeneous in economic terms, which makes it an ideal candidate for fiscal decentralization. Institutionally, Russia is a federation consisting of the federal government and 80 subjects of federation, i.e. regional governments. Over roughly the past twenty years, the Russian Federation has travelled a long path from a unitary state, through extreme decentralization, to becoming a well-defined, highly centralized, federation. Broadly speaking, the past twenty years can be divided into two periods; haphazard and spontaneous fiscal and political decentralization in the 1990s and recentralization in the 2000s. The evolution of Russia's fiscal federalist arrangements (i.e. how fiscal rights and responsibilities are divided between levels of government) has been actively researched and documented in numerous studies including De Silva et al. (2009), Desai et al. (2003, Solanko and Tekoniemi (2005) and Zhuravskaya (2010). Four broad issues characterize Russian fiscal federalism.
First, the Russian regions are heterogeneous in almost every aspect of economic development. The differences in e.g. per capita income, in public health, wages, industrial structure and tax revenues are huge. In 2009, the fiscal capacity gap difference between the richest and the poorest region was 48 times (World Bank 2011). Over the last ten years, the federal fiscal transfers have become better formalized and focused. Thus, federal equalization transfers today iron out the largest inequalities, although a seven-fold gap in regional fiscal capacity still remains.
Second, formal revenue autonomy is negligible. Even in the 1990s, regional authorities had only limited powers to decide on tax bases or tax rates. Apart from the corporate income tax in the 1990s and the property tax in the 2000s, all tax rates and bases are centrally determined. Further, possibilities to finance budgets with sub-national debt are limited (Kurlianskaya 2013). Third, autonomy on the expenditure side is quite limited. Most social spending is implemented through regional budgets. These include outlays for health, education and housing that jointly constitute almost 70 % of consolidated regional expenditures (Rosstat, Regioni Rossii). Most of all social expenditures are mandated in federal laws and regulations. Regions have somewhat more freedom in deciding on various benefits to regional (public and private) enterprises under the heading of "national economy." These expenditures constituted roughly 15 % of regional expenditure in 2006-the midpoint of our time-series.
Fourth, minimal formal autonomy notwithstanding, Russian regions have considerable powers in implementation of federal and regional legislation. Despite fairly uniform rules and regulations on social expenditure items, actual outcomes differ widely (World Bank 2011). Over the past twenty years, regions have proven successful in influencing both tax bases and effective rates e.g. by affecting regional tax collection, tolerating tax arrears and giving preferential treatment through regional regulations (Slinko et al. 2005).

Literature Review
While there are many mechanisms through which climate change adaptation costs might reflect on government expenditures 1 , the literature on fiscal costs and benefits of climate change is scarce. A typical limitation is that such studies must rely on assumptions about future adaptation behaviour. Among the few available studies, Osberghaus and Reif (2010) provide estimates for European countries with a methodology largely based on "guesstimates" of public sector involvement in adaptation investments. The World Bank (2010) estimates the cost of hard adaptation measures (i.e. those requiring engineering responses) in developing countries by analysing additional investment needs in various economic sectors. Both studies estimate the effects to be relatively low, although there are big differences, especially among the poorer countries. Bräuer et al. (2009) use a combination of qualitative-quantitative analysis of climate change impacts on several sectors to estimate the direct and indirect effects of climate change on government finances in Germany. They estimate that by 2100 the decrease in revenues and increase in government spending could result in a negative impact on the government expenditures equalling a GDP loss of 0.6-2.5 % compared to a reference scenario. For 2050, the results are less clear and might even have a slightly positive effect. These methods have utility, but the underlying behavioural assumptions are always open to criticism.
A relatively new, but rapidly growing, body of literature on fiscal effects of natural disasters complements adaptation literature through the use of historical data on past natural disasters. There is strong scientific consensus that climate change will increase the severity and occurrence of weather related disasters, and that such disasters have significant economic and fiscal repercussions. Ouattara and Strobl (2013) use data for Caribbean countries and find a positive government spending reaction persisting up to two years to large hurricanes. Schuknecht (1999) finds a significant negative effect of catastrophes on fiscal balances in a fixed effects setting with developing country data. Borensztein et al. (2008) find similar results. Based on EU and US case studies, Heipertz and Nickel (2008) conclude that the combined direct and indirect impact of extreme weather events on public finances was between 0.3-1.1 % of GDP. Lis and Nickel (2009) estimate a panel data model of 138 countries, concluding that extreme weather events have a fiscal impact of 0.23-1.1 % of GDP depending on the country group. Noy and Nualsri (2011) assert that developed countries exhibit counter-cyclical, and developing countries pro-cyclical, fiscal behaviours in the aftermath of natural disasters. Melecky and Raddatz (2011), using data of high-and middle-income countries for 1975-2008, find that government expenditure increases and budget deficits worsen after natural disasters.
These studies consistently indicate that extreme-weather-related phenomena affect public expenditures. What is missing is in the literature is evidence on how moderate, but persistent, changes in climate might affect public finances. This is precisely where we wish to make a contribution to the literature by offering an alternative way to estimate fiscal costs that is based on realized, historical data from a large, heterogeneous economy that has been experiencing warming.
The existing literature on determinants of public expenditure is deep and has its roots already in the political economy discussions of the late nineteenth century, but climaterelated variables are all but non-existent in the discussions. Broadly speaking, factors shaping public expenditure can be grouped into macroeconomic, demographic and institutional categories. The literature often argues that there is a positive relation between income levels and most expenditure categories and that government expenditures tend to increase in economic upturns, but adjust slowly in downturns (e.g. Shelton 2007). However, both results have been seriously challenged by e.g. Durevall and Henrekson (2011) who use long historical datasets from the UK and Sweden.
On the contrary, there is a broad consensus that the demographic structure of a country matters for the level and structure of its public expenditure. Population density tends to produce negative elasticities for defence, transport and communication and merit goods. Also, population density would seem to be positively linked with housing expenditures (Sanz and Velázques 2002). Age structure and dependency ratios affect public expenditures greatly (Sanz and Velázques 2002;Shelton 2007;Durevall and Henrekson 2011).
A country's political system and institutional arrangements bear on its public expenditures. Since most of the existing literature on public expenditure determinants uses cross-country data, the effects of the institutional differences are not easily analysed. However, a rich body of political economy literature considers institutional and political determinants of fiscal outcomes, mainly using regional datasets. In his seminal paper based on the US state government data, Merrifield (2000) argues that a whole range of political variables, including political fragmentation and turnover rates of state-level politicians, may greatly shape statelevel public finances. There is a rapidly growing body of literature taking advantage of the variation in Russia's regional institutions. To cite an example, a recent paper of Schultz and Libman (2015) finds that regional responses to the massive forest fires across Russia in the summer 2010 depended on characteristics of regional governors.

Estimation Through Annual Weather Variation and Mid-Term Climate
Change Dell et al. (2014) note that the word climate is reserved in the literature for the distribution of weather outcomes (e.g. annual average temperature, precipitation sum, etc.) and can be summarized by weather averages over several decades, while the word weather describes a particular realization from the climate distribution and can vary notably in the short term. In this paper, we will use the same terminology and conceptual approach. There is one notable econometric problem in estimating climate impacts. This follows from the very definition of climate, which essentially makes it a fixed regional variable in an econometric sense. Thus, we cannot include a general regional fixed effect in a climate model (as it would be perfectly collinear with the climate variable) and a model without a regional fixed effect might suffer from omitted variable bias. This is most likely the case in Russia, where there is considerable heterogeneity among regions. We approach the problem of estimating climate with two distinct estimation strategies. First, we will run a fixed-effects (FE) model: where Y is regional expenditures in year t in real terms, X is annual temperature realization, P similarly for precipitation, Z j the control variables, and μ i and θ t regional and time dummies. Regional fixed effects control for any unobserved region-specific time-invariant variables that may affect our dependent variable while the time dummy controls for country-wide trends (e.g. oil price changes).
As we cannot estimate the climate effect directly in an FE model, the strategy is to identify the effects of climate change through annual variation in weather. Deschênes and Greenstone (2007) use a similar approach with an agricultural application. However, as Seo (2013) argues, the FE model fails to grasp potential adaptation due to its short-term aspect. Thus, the FE approach basically gives us the effect of climate change with very limited adaptation (that would be applicable within a year) and might be interpreted as a "no-adaptation" benchmark result. The intuition behind this strategy is that what would happen if the annual weather variation from the average regional weather (i.e. climate) in a particular year would become permanent.
To tackle jointly the adaptation and omitted variable problems, we employ our second estimation strategy, which, Dell et al. (2014) argue, comes closest to estimating medium-term impacts of climate change. Consider the model (presented along lines of Dell et al. 2014, where our time subscript is now d indicating decade. Thus, C id being the climate in period d obtained through averaging annual weather (temperature and precipitation separately) figures and taking the first difference: Burke and Emerick (2013) 1970-1985 and 1985-2000. Due to the limited data span available for Russia, we will use two non-overlapping periods of 1995-2001 and 2002-2009. We acknowledge that weather averages of less than ten years probably cannot credibly depict local climate, which somewhat undermines the analysis. Despite of this shortcoming, we believe that the long-difference model can serve as a guideline for whether the short-term effects obtained from the FE model will increase or decrease in the longer term-especially taking into account that the Russian mean temperature has steadily increased during our time-series as evidenced in Fig. 2 above. 2 Dell et al. (2014) and Burke and Emerick (2013) discuss the various benefits of the long-difference model approach. In the mid-to long-term, climate change might induce adaptation in behaviour of economic agents (e.g. technology and policies), intensification of climate change effects (e.g. desertification) and general equilibrium effects (price and factor reallocation) that cannot be captured with an FE model. Long-difference models are not a panacea though. For example, if the two periods are too far apart, e.g. several decades, it may well be that the transmission mechanism between climate and economy has changed so much that it is not plausible to assume common regression coefficient for the two periods. Such setting could happen due to above-mentioned technological advancements used for adaptation. However, comparing LD model with other options such as hedonic models (e.g. Mendelsohn et al. 1994;Dell et al. 2014) Burke and Emerick (2013) argue for the superiority of the LD model as it removes omitted region fixed effects.
It is important to note that generally fixed effect estimations is preferred to first differences due to its better efficiency (Angrist and Pischke 2009). However, because FE tends to amplify the R 2 especially with small data sets and the deviations from the mean is algebraically the same as differencing when T = 2 (Angrist and Pischke 2009) as in our case, we choose to use the LD model here to follow the earlier literature. As can be seen in Table 7 in the Appendix where we have used fixed effects, the results follow the ones in Table 3 obtained through differencing.

Model Specification
In climate-related applications, the issue of non-linearity of weather (especially temperature) often arises. Several health-related studies (e.g. Deschênes and Greenstone 2011) confirm the intuitive observation that extreme heat and cold have negative effects on human health, while an optimum lies somewhere in between. Schlenker and Roberts (2009) demonstrate similar results for agriculture. Due to heterogeneity of Russia's regions, it is quite possible that non-linearities exist concerning the temperature effect on expenditures. Non-linearities in climate-related studies are often captured by examining temperature bins (as in the two studies mentioned above). For the bin approach to work properly, however, the data need to be highly disaggregated at the regional level (e.g. US county-level data that provide enough observations for each bin). This is not possible with the Russian data. Thus, we test for non-linearities using both linear and quadratic specifications for the climate variables. This is the reason why the form of the temperature and precipitation variables is not explicitly specified in Eq. (1).
Climate and weather are ideal independent variables for econometric analysis. As Dell et al. (2014) note, weather shocks can be considered as random draws from the climate distribution and thus have strong causal identification properties. Despite the causal exogeneity of climate and weather, endogeneity problems can arise especially in a cross-sectional setting when we omit time-invariant variables that correlate with climate. Burke and Emerick (2013) argue that, like the FE model, the LD model has good identification properties concerning the climate variable. This is because climate trends (i.e. average weather over a longer period at various points of time) can also be interpreted as randomly assigned treatments to different regions.
Biased estimations could arise also in a fixed effects setting where there are time-varying omitted variables. To make sure we specify the model correctly, we therefore include several control variables. Dell et al. (2012) and Burke and Emerick (2013) rely only on region and time fixed-effect controls to avoid the so called over-controlling problem, which arises when we include control variables that are, in effect, the mechanism through which our variable of interest (climate) affects the dependent variable. To be prudent in including time-varying control variables, we therefore limit ourselves to variables already mentioned in the literature as important determinants of government expenditures.
In a setting where our statistical units (Russian regions) are geographically correlated, the error terms could be correlated as well and thereby violating the assumption of non-correlated error terms. While spatial autocorrelation leaves estimated coefficients unbiased, it can lead to incorrect standard errors undermining inference. This problem might arise in our case where neighbouring regions are correlated climatically in Russia. Because of this, we use Driscoll and Kraay (1998) standard errors in the FE model. These standard errors are robust to heteroskedasticity and general forms of spatial and temporal dependence (Hoechle 2007). For the LD model, we use typical heteroskedasticity robust standard errors because Driscoll and Kraay errors work poorly when the cross-sectional dimension of the data is much larger than the time dimension (i = 78 and t = 2 in our long-difference setting).
Another issue to resolve is whether to include the lagged variables in the FE estimations. As Dell et al. (2014) note, inclusion of lags of the dependent variables as explanatory variables creates bias in panel models with short time periods. Thus, this is a potential problem with our relatively short panel of T = 15 (1995-2009). In addition to shortness of our time dimension, there is no certainty on the underlying data generation process for the expenditure variables and lag of y is not typically introduced into estimation in the econometric climate literature. For these reasons, we leave it out of our analysis as well. Concerning lags of the independent variables, the situation changes. There is no clear theory for the proper time dimension in which expenditure variables could react to weather. Thus, we test the best fit of weather variables and find that expenditures seem to have the most robust fit with the first lag of the temperature variable. This is reasonable taking into consideration that effects can easily be reflected into the expenditures with a lag due to budget planning and sequencing. In the LD model, lags are not relevant because we are looking at a change in the mid-term climate defined as annual averages over 1995-2001 compared to 2002-2009.

Dependent Variable
Our main dependent variable is total regional government expenditures per capita as reported by Russia's Ministry of Finance and Federal State Statistical Service (Rosstat). Additionally, to detect the potential mechanisms through which weather or climate can affect expenditures, we run the analysis separately for per capita housing and communal expenditures (shortened as housing expenditures in the following analysis). This is among the few regional budget classes that are comparable and available throughout the 1995-2009 period in Russia. Although data-wise we cannot separate housing and communal expenditures from each other, housing expenditures at least provide an intuitive transmission mechanism from climate to expenditures (i.e. in a cold climate, higher temperature means lower heating costs). Roshydromet (2014) estimates that energy demand for heating will reduce notably but at the same time energy demand for cooling will increase in the warm season. According to Solanko (2006), heating is perceived as a semi-public good in Russia due to the heavy historical reliance on district heating typically provided by municipal heat and power plants. Furthermore, Rosstat figures show that in 2009 still roughly one-sixth of the housing stock in Russia was under state or municipal ownership (in 2000 roughly a third), with ownership shares ranging regionally from 2.5 to 35 %. Public ownership of housing is much higher in colder regions, and indeed the correlation between average temperature and share of public ownership in our data is −0.44. This suggests housing as a potential channel for impacts of temperature changes on regional budget expenditures.
Our dependent variables are recorded in annual rubles terms, so Russia's high and volatile inflation rate has to be properly taken into account. To do this, we deflate regional expenditure figures with regional consumer price indices (cpi). As a robustness check, we also deflate expenditures using the annual price of the regional consumption basket (i.e. basically a purchasing power parity approach). To save space, we report here only the results from the cpi approach; they are easier to interpret and the results from the consumption basket specification are highly similar (available on request). Actual expenditure-per-capita variables are introduced in log form. Thus, the interpretation of the continuous variables' estimation coefficients is β times 100 %. Our dataset comprises annual regional level data across Russian regions for 1995-2009. We exclude Chechnya from the sample due to data unreliability and the Chukotka region due to its political peculiarities. 3 Chukotka, with a population of 50,000 people, clearly has the lowest population density of all Russian regions (average population density in Russia is over 2000 times that of Chukotka) and it is the second coldest region in Russia. This leaves us with 78 regions in our final sample. Descriptive statistics of all independent and dependent variables are presented in Table 1.

Regional Climate Data
Our climate data are obtained from the online database of the All-Russian Research Institute of Hydro-meteorological Information (meteo.ru), which is funded by Roshydromet. The database provides daily and monthly temperature and precipitation data over several decades for more than five hundred weather stations located across the Russian Federation. For our estimations, we calculate the annual average temperature and sum of precipitation for the years 1995-2009 for 78 Russian regions in our main specification. Climate variable specification is highly context dependent. In agricultural studies, it is common to use "growing degree days" within certain temperature thresholds during the growing season. This approach tries to capture the biological relation between weather and vegetation growth. In some applications, seasonality may also matter (e.g. tourism). Because there is no clear intuition on seasonality or potential opposite effects concerning budget expenditures, we use annual temperature figures in our main specification. This also makes it easier to interpret the results in the framework of climate predictions.
However, for robustness check we construct higher frequency climate variables reflecting potential impacts through agriculture and energy use from heating and cooling. These variables are growing degree days (GDD) as described above as well as heating and cooling degree days (HDD and CDD, respectively). GDD variable is calculated as the sum of daily average temperatures between 8 and 32 • C and truncated to 32 • C during the growing season of April to August (similarly to Deschênes and Greenstone 2007;Guiteras 2009). HDD and CDD on the other hand are calculated as annual sum of daily average outdoor temperature differences to indoor comfort temperature of 18.3 • C (65 • F) following EPA (2014). For example, if outdoor temperature is 10 • C we have a HDD of 8.3 • C for the particular day and in similar fashion for outdoor temperature of 23.3 • C we get a CDD of 5 • C. Our quadratic model take into account potential non-linearities concerning HDD and CDD, which according to Dell et al. (2014) is important as extreme temperatures provoke much stronger energy demand increases.
The weather data is calculated from the weather stations located within each region. Typically, a region has several weather stations, so the weather data are averaged from these stations to get a regional figure. For the large northern and eastern regions, we included only weather stations situated relatively close to regional capitals. We also exclude weather stations situated at very high altitudes. On average, we use data from 3.5 stations per region to calculate the annual regional temperature and precipitation variables of our estimations. Our data-generating approach creates an upward bias in temperature data with respect to pure geographic dimensions, as the weights of large and cold northern and eastern regions in the weather data are reduced. However, as discussed in Dell et al. (2014), there are various criteria for aggregating weather data. In economic applications, the main options are spatialand population-weighted aggregation and should be selected depending on context. Spatial weighting is more appropriate for e.g. agricultural studies, while topics related to human activity are likely best captured by population-weighted weather data. Our data generation approach can be considered as largely population-weighted (due to weather station selection), but still allows a bit more weight for geographical size than a pure population-weighted aggregation. For robustness, we also run the regressions with data where 5 % of observations are removed from both tails of the temperature and precipitation variables. This diminishes the role of potential weather outlier regions.

Control Variables
Our control variables were selected to represent the three categories of variables (macro, demographic and institutional), which, as discussed in Sect. 3, have been found important in the literature on determinants of public expenditure. We performed a robustness check on covariate selection in each regression by altering the included control variables (not reported here to save space). For the FE model, both the significance and coefficients of the weather variables were highly robust and quite strongly robust for the long-difference estimation.
To control for regional business cycles that could have notable effects on budget, we introduce regional industrial growth rate (Industrial_growth) as a covariate. Following Sanz and Velázques (2002), we include separate dependency ratios (dependents to working-age population) for younger and older population (Underage and Overage, respectively). We also control for regional population density (Pop_density). In the housing expenditure estimations we include the share of public houses in total regional housing stock (Public_housing). The source for these variables is Rosstat's Regioni Rossii publications. Variables reflecting regional political and institutional circumstances are always hard to come by. Such variables are particularly scarce for Russia at the regional level, especially those that would consistently cover the full period 1995-2009. We chose here the regional share of the Communist Party vote (Communist) in elections of the Duma (Russia's lower house of parliament) as a proxy for regional political activism. 4 This is a fairly straightforward measure of regional political dimensions and the Communist Party is the only party that has held a significant number of seats in the Duma throughout Russia's transition. This measure can also be considered as a proxy for regional antipathy towards the ruling party. The data was retrieved from the website of the Central Election Commission of the Russian Federation.

Short-Term Analysis Results
We next present the results for our fixed-effects models for both total and housing expenditures. As will be shown, it appears that a typical linear model fails to sufficiently grasp the relation between weather and expenditures although the nationwide non-linear estimate is relatively similar. We will present the linear results as a reference point to the methodologies used in other climate econometrics literature. The results of the fixed effect estimates for both total and housing expenditures are shown in Table 2. The table presents results for both linear and quadratic models for total and housing expenditures with full data, as well as for robustness check with data where 5 % from both tails of the temperature and precipitation variables are removed.
For total expenditures, our basic linear model (1) suggests that a rise in the annual (lagged) temperature decreases regional expenditures, i.e. a 1 • C increase leads to a roughly 3 % decrease in regional expenditures per capita. No such significant relation seems to exist for precipitation. Underage and Overage are positively associated with expenditures, which is intuitive in the sense that education and healthcare are mostly sub-national responsibilities in Russia. Industrial growth also shows a significant positive relation and a larger role of the Communist Party would seem to be reflected as smaller expenditures. The causality of the Communist variable is unclear, however. It is possible that the Communist Party enjoys greater support in poor regions that lack money for public services. Indeed, there is a strong negative correlation (−0.58) between regional income per capita and support of the Communist Party in the data. In the basic model for housing expenditures (5), the control variable results seem similar with the difference that Population density is now highly significant while Underage is not. However, there would seem to be no significant temperature or precipitation effect.
In the linear estimations, where the tails are cut for climate variables (models 3 and 7), our control variables seem fairly robust, while climate variables are insignificant. However, the linear model appears to lack sufficient flexibility to capture the effect of the climate variables. A quadratic relation seems to fit for temperature effect on total expenditures and is even more robust when tails are removed. There is also evidence in the cut tails regression (4) that precipitation has a similar (although less robust) non-linear effect.
To get a better sense of the non-linear climate relation, we build a graph where we show the temperature effects on total expenditures at different regional temperature levels based on the results in Table 2. In Fig. 3, the y-axis shows how many per cent expenditures would drop for a temperature increase of 1 • C in regions with particular average temperature. The quadratic specification suggests that increase in temperature has a steep negative impact on expenditures in colder regions, while the effect approaches zero in warmer regions. In our full data results, the impact remains negative for even the hottest regions (but reverses at an annual average temperature of around 12 • C), while the cut data show additional warming starts to raise expenditures in those regions where average temperature is around 8 • C or higher. This supports the intuitive hypothesis that cold regions benefit from warmer weather, while the benefits diminish in warmer regions (and could even reverse for the hottest regions). We find some support for such U-shape effect for housing expenditures in the cut data (see Table 2 model (8), not graphed) suggesting that savings in housing expenditures due to warming become net costs already at around an average regional temperature of 4 • C, which is roughly the Russian average temperature in our data.
To obtain a rough nationwide effect, we calculate a population-weighted average effect of regional effects. For total expenditures, we get a nationwide coefficient of −0.026 (i.e. 2.6 % decline per one degree increase) with full data and −0.015 with cut data. Thus, the nonlinear analysis confirms the decreasing impact of warming on expenditures, while the effect is clearly smaller. This suggests that non-linearities are helpful in understanding regional differences in climate impacts. As these might be important for policymaking, they should not be sidestepped.

Fig. 3
Linear and quadratic fit for weather variation at different temperatures (x-axis is annual average temperature of a region in • C, y-axis is the temperature effect at different temperatures, 100 * β = %)

Longer-Run Effects
We next address the question of adaptation through the LD model (Eq. 2), as well as the results with respect to climate scenarios. As noted in Sect. 4, the FE model only captures short-term adaptation (within a year). Thus, its ability to estimate longer-run climate impacts including adaptation is weak. In Table 3, we present our estimates for the LD model, which should better account for adaptation and/or potential intensification effects than the FE model. However, it should be emphasized that the long-difference results should be treated with great caution and taken mainly as indicating direction as to whether the benchmark effects of the FE model are likely to increase or decrease over the longer term. This is due to the small amount of available observations and relatively short time span for weather averages in our data. We do not find any significant non-linear impacts with the LD model and thus present only the linear model.
The results of Table 3 suggest that temperature has a cumulative effect on total expenditures; the longer-run coefficient is ten times larger in the absolute sense than the short-term FE estimate. The estimation suggests that a +1 • C increase in regional average temperature leads to a decrease of roughly a third in real expenditures per capita. This magnitude seems implausibly high and we have to take into account that both expenditures and temperature differences in our data consists of solely positive values-due to economic development and general warming of the whole country-with the exception of two regions with temperature change of −0.01 • C . Thus, the LD model effectively captures the fact that in regions of greater temperature increase, expenditures have increased less. Also, there might be some trend or non-fixed variable that we are unable to control for in the LD model. Nevertheless, it would seem that the temperature effect suggested by the FE model is probably larger when adaptation and intensification come into play. While temperature is significant and precipitation is not in the case of total expenditures, the situation is reversed for housing expenditures.

Robustness Check with Alternative Temperature Variables
Our results for the three degree day variables, working as a robustness check, are presented in Tables 5 and 6 in the Appendix. The tables replicate the models of Tables 2 and 3 but we have left out the control variables in the table to save space as they were robust to the new temperature variables. As can be seen, we do not find evidence on temperature effects through agriculture as the GDD variable is insignificant in both linear and quadratic models and in short and longer-run analysis. This could be due to relatively small role of regional governments in supporting Russia's agricultural sector. For example, in 2004 regional budget expenses to the agricultural sector in Russia were roughly 1.6 bln. USD (less than 2 % of all budget expenses). On the other hand, HDD would seem to support the results of our main specification. As presented graphically in Fig. 4 of Appendix regions that have low HDD (i.e. warm regions) face decreasing expenditures if the HDD increases (i.e. temperature goes lower). Intuition behind this result is that in the warmest regions the expenditure drop due to colder weather in general is larger than possible increase for required additional heating. Conversely, for cold regions that have already high heating needs, further increase in HDD accrues increasingly high expenditures. The HDD estimations are very consistent with our main specification also in the sense that LD estimate is roughly ten times larger the FE estimate. We can find similarly intuitive results for CDD and housing expenses as in the case of HDD but these results are clearly less robust. This could be a result from Russia being generally a cold country and thus not needing much of air conditioning and also Russia is still relatively poor, especially its warmest regions, and thus might not have a high coverage of air conditioning.

Warming Effect in the Framework of Climate Scenarios
Finally, the use of simplistic proxies (temperature and precipitation) to describe climate change allows us to discuss our results in the context of climate change scenarios. Utilizing the GAEZ database (http://gaez.fao.org), which takes climate data from the Climate Research Unit of East Anglia (CRU-A, -B, -C and -D), we retrieve two climate-change predictions from the database to obtain the change in average temperature from 1995-1999 to 2020s. These predictions are simulated with the commonly used HadCM3 climate prediction model for the scenarios A2 and B1. 5,6 Naturally, there is no objective reason to rely on these particular scenarios. However, our intention is to merely illustrate our results in a framework of actual climate predictions and these scenarios are regularly used in relevant discussions (e.g. the World Bank's climate change portal) and full presentation of prediction uncertainties would not serve the purpose of the paper. Although the scenarios initially provide quite similar climate predictions for Russia as a whole (before starting to deviate in decades further out), certain regional differences are apparent, so we present both for robustness. Due to relatively weak robustness of our precipitation estimates and higher uncertainty of the precipitation predictions we will focus here on temperature change scenarios. We provide estimates only up to the 2020s (i.e. the average of 2011-2040 reference point given by the climate scenarios), because temperatures are expected to continue to further increase (>2 • C), i.e. the benefits from warming become increasingly uncertain. Moreover, it is problematic to derive predictions outside the variation available in our data. From the scenarios, we calculated regional temperature predictions and multiplied these with the regional temperature gradients (i.e. derivate of the obtained quadratic function) to get regional impact estimates. We calculated a population-weighted average of the regional impacts to get an estimate for whole Russia.
To get some insight into the monetary impacts, we take Russia's regional government expenditures in year 2000 (roughly 1 trillion rubles). Inflating this figure using the consumer price index gives roughly 4 trillion rubles in 2013 prices, or USD 125 billion using the average RUB/USD exchange rate for 2013. The results are summarized in Table 4. In the B1 scenario, Russian regional governments save between just over USD 2 billion to USD 4 billion (in real terms, non-discounted), depending on the data set, between 2000 and the 2020s due to warming. In the A2 scenario, the savings from milder temperatures ranges between USD 2 billion to roughly USD 3.5 billion over the period. Compared to, say, Russia's consolidated regional budget expenditures of USD 275 billion in 2013, the benefits are quite small. Barring any major uncertainties in the mid-scenario climate comparison, our weather data suggests that roughly half of the predicted 1.1-1.2 • C increase in annual mean temperature had taken place by the end of the 2000s i.e. half of the expenditure benefits had already accrued.
An intriguing question is whether these results are applicable in a global perspective. Russia offers a wide spectrum of temperature zones. This improves generality of the results compared to climatically homogenous countries like Japan or the UK. On the other hand, Russia is a cold country compared to most. The population-weighted mean annual temperature for Russia based on our data is around 4-5 • C, leaving it with a small number of peers such as the Nordic countries and Canada (as presented in Dell et al. 2012). The fact that our results suggest a diminishing convex-shaped effect of warming on expenditures in such a cold country as Russia would indicate that similar benefits are probably not available in warmer countries. However, our longer-term results suggest that, with proper adaptation and potential intensification effects, the benefits could well be larger for Russia than the short-term analysis indicates.

Conclusions
In this paper, we examined the effect of climate change (concerning temperature and precipitation) on regional government expenditures using a rich regional-level dataset. This effect was estimated through short-term variation and mid-term change in average temperature and precipitation. The former approach conceptually gives us a no-adaptation benchmark effect while the latter approach captures potential adaptation better, but is less reliable due to the scantiness of observations and data variation.
Our short-term approach suggests that a rise in temperature reduces regional public expenditures in cold regions, becomes negligible in warmer regions and reverses in the hottest regions. Thus, the effect is non-linear. We find similar relation for precipitation but the result is not very robust. We also find some evidence that housing and communal expenditures provide a mechanism through which temperature affects total expenditures. We found support for the temperature effect result from estimations with higher frequency data using degree day variables.
Our non-linear and population-weighted estimate for the Russia-wide impact of a 1 • C increase in annual temperature is a decrease of 2.6 % in regional budget expenditures in real terms. To estimate the monetary impact of warming we limited the duration of our climate scenarios up to 2020s, when the expected warming still only slightly exceeds 1 • C. Our estimates indicate that warming effects save the Russian regional governments somewhere over USD 2 billion to USD 4 billion, depending on the climate scenario, in expenditures in non-discounted US dollars between 2000 and the 2020s. Roughly half of these savings had already accrued by the end of the 2000s. Our mid-term model suggests that these benefits could well be larger when adaptation and intensification effects come into play. However, the mid-term result is less robust and our results cannot account for a large increase in temperature, which would likely have a dramatic effect through e.g. thawing of the permafrost that covers roughly two-thirds of Russia's land area.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Table 5
Estimations on total expenditures and housing expenditures with additional temperature variables (growing, heating and cooling degree days); linear and quadratic fixed effects, full data models 1-2 and 5-6, models 3-4 and 7-8 with data excluding 5 % of both tails from temperature and precipitation variables Log of total expenditures,   Driscoll-Kraay standard errors in parentheses. Control variables same as in Table 2, but not reported. *, **, *** indicate significance at ≈10, 10, 5 and 1 % levels All coefficients multiplied by 1000 to improve presentation  Median HDD Fig. 4 Graphical result of the heating degree days estimations (a1) and (a2), x-axis is HDD, y-axis is the HDD effect at different HDD (100 * β = %)