Examining lag effects between industrial land development and regional economic changes: The Netherlands experience

In most empirical applications, forecasting models for the analysis of industrial land focus on the relationship between current values of economic parameters and industrial land use. This paper aims to test this assumption by focusing on the dynamic relationship between current and lagged values of the ‘economic fundamentals’ and industrial land development. Not much effort has yet been attributed to develop land forecasting models to predict the demand for industrial land except those applying static regressions or other statistical measures. In this research, we estimated a dynamic panel data model across 40 regions from 2000 to 2008 for the Netherlands to uncover the relationship between current and lagged values of economic parameters and industrial land development. Land-use regulations such as land zoning policies, and other land-use restrictions like natural protection areas, geographical limitations in the form of water bodies or sludge areas are expected to affect supply of land, which will in turn be reflected in industrial land market outcomes. Our results suggest that gross domestic product (GDP), industrial employment, gross value added (GVA), property price, and other parameters representing demand and supply conditions in the industrial market explain industrial land developments with high significance levels. It is also shown that contrary to the current values, lagged values of the economic parameters have more sound relationships with the industrial developments in the Netherlands. The findings suggest use of lags between selected economic parameters and industrial land use in land forecasting applications.


Introduction
Industrial land constitutes a significant part of the land value in urban areas and regions. Considering that land use reflects decisions of the private and public agents, the patterns of land development observed in an urban environment reflect the decision-making processes of those who manage scarce land resources. Industrial and commercial development policies such as the availability of subsidies and tax benefits applied to new development schemes, socio-economic growth and technological development, external trade policies, fixing of a1111111111 a1111111111 a1111111111 a1111111111 a1111111111 other markets. These are imperfect market information, cyclical nature of property developments resulting in oversupply or undersupply in the market, high transaction costs, a number of development constraints, irreversibility, and illiquidity [25][26][27]. Considering existence of such failures, there are lags between the decision to begin real property development and its completion [27,28].
From government's perspective, it is important to get the information on lags between the time when an irreversible decision is undertaken and when the development's first revenues are received. These lags may be quite long [29]. For instance, in the Netherlands, there is a gap of six to eight years between development of a land-use plan and provision of the industrial sites for development [30]. Following the development of an industrial site, the land use has almost a permanent character considering that urban re-development is very time consuming and expensive [30]. Therefore, it is important that any substantial industrial capital investment programme should be supported by various considerations such as societal needs, budgetary constraints and potential impacts of development lags on social welfare. The underdeveloped land has negative social impacts e.g. potential lack of business estates for investors, vacant lands, and soaring land and industrial property prices while overdeveloped land may result in vacancies in industrial property, overconsumption of land, and declining property prices. Without an estimation of the relationship between development lags and industrial property development, any capital investment policy could result in a non-reversible and unsustainable consequence in the long run.
From policy makers' and planners' perspective, it is necessary to consider how sensitive the industrial property market to economic fluctuations and development lags. However, this is a difficult task without developing a reliable forecast of industrial land demand. Therefore, it is essential to focus on statistical and econometric techniques to explain the interaction between industrial development and real estate development time lags. In this study, a dynamic panel regression approach is adopted to predict industrial land demand. The model can enable policymakers and planners to produce more effective policies and regulations to guide the real property market, tackle existing and potential problems and promote social welfare. The model can also facilitate construction stakeholders to forecast the industrial land through incorporating development lags into the model in the short to medium term to enable the industrial construction market to plan forward [31]. This would enable the construction industry and other construction stakeholders to optimise their returns from industrial investment and sustain economic development.
The research on forecasting demand for industrial or commercial land is scarce as the literature is mainly dominated by studies that focus on the drivers of urban land-use change [32,33] and the factors governing variations in industrial property rents (or prices) [34,35] or property developments [36]. In the scarce literature, there are studies of forecasting work for the industrial property market utilising supply and demand dynamics and prices (or rents) [37,38]. Some examples of research focusing on the dynamics of industrial demand in relation to economic and sectorial variables are Beckers and Schuur [18], Ramirez [39] and Sing [40]. The research by Sing [40] found significant coefficients for the first three lags of the GDP variable on industrial space demand in Singapore. By contrast, Beckers and Schuur [18] focused only on the relationship between employment and industrial land demand in the Netherlands. Though the study found significant non-linear relationship between sectorial employment and land use, it does not provide any evidence of time lag effects of the employment variables on industrial land use.
The studies by Hoymann [13], and Reginster and Rounsevell [41] utilised regression models to estimate demand for urban area. The former study compared results of the residential demand predictions with the observed data through application of different statistical approaches. The latter research applied regression techniques for the prediction of urban demand, which was fed into their studied land-use model to simulate future land-use changes. Despite estimations of significant coefficients for the selected predictors of urban land use, none of their regressions considered development lags concerning the urban land demand. Among few studies of forecasting work of the industrial land, Batista e Silva et al. [10] estimated demand for industrial and commercial land for the EU countries by using a number of statistical approaches and compared prediction results with the observed data following the research by Hoymann [13]. However, the study did not consider lag effects of the economic development on the growth of industrial and commercial land use in their studied countries.
Except few studies, time lag effects of economic changes on land use have rarely been empirically assessed particularly focusing on industrial land developments. In this respect, our study attempts to contribute to the existing research through providing a dynamic analysis example to uncover the relationship between lagged values of selected economic parameters and industrial land development across 40 regions spanning 8 years (2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008) in the Netherlands. We use a dynamic panel data approach which enables us to integrate the spatial characteristics of the regions into the time-series analysis, and can improve the robustness of the estimates.
The remainder of the paper is organised as follows: Section 2 reviews the literature on real property development lags and construction demand forecasting methods. The data and methodology are described in section 3, and estimation results are presented in section 4. Section 5 compares the actual and the model estimation results of industrial land developments. Section 6 concludes the paper.

Literature review Real property development lags
The overbuilding (and underbuilding) of real property development is well documented in the literature, particularly for commercial and residential property markets. Wheaton [42] provides evidence on office property over-buildings in the past decades in US (see also King and McCue [43]; Grebler and Burns [44]; Hendershott [45]). Barras and Ferguson [46] showed that new commercial and industrial property demand coincides with the business cycle, however there is a lag of three to four quarters between new orders and construction completions. Barras [47] examined the variables affecting demand for office space such as output prices, user cost of capital and profitability. His results suggested that lagged output is the most important variable. He also stated that the long production period between new building orders and completions are the crucial endogenous factors causing building cycles. DiPasquale and Wheaton [25] have demonstrated that stock adjustment through new construction is much slower than price adjustment in real property markets when there is an external shock in the economy. In the long run, the stock adjusts gradually considering the lags in the delivery of new stock. Capital investment decisions are based on the forecasts of property prices at the time of new construction completions [42]. This implies that any external shock to the economy is first reflected in the real property prices and rents but not in physical assets.
Grenadier [48] (cited in Ball et al. [49]: 205) argues that 'the lag explanations assume myopic behaviour of developers and therefore they should take into account the future state of the market'. Ganesan and Tse [50] examined the lead-lag relationship between construction activity and the aggregate economy in Hong Kong. The study concluded that economic output results in construction flow. However, the initial impact of a change in GDP would be on the change in real estate investment projects but not on the level of construction output considering that construction activity is very sensitive to credit conditions and land development restrictions. Wheaton [42] examined the future expectations of agents in a stock-flow model where development lags, the degree of durability and market elasticities can vary. He concluded that stock-flow models with myopic behaviour present high sensitivities to all of these parameters and that rational model can also generate cyclical behaviour based on exogenous structure on the market's operation that leads to historical dependence.
The existence of real property construction lags was examined by Spiegel [51] in the framework of a general equilibrium model for the residential housing market to explain the housing construction cycles in relation to expected housing returns. Barras [52] alternatively constructed a dynamic model of building cycles for the London office market. The adjustment process in his model is determined by three types of lags: the occupiers' response to change in rents, the developers' response to changes in demand through rent adjustment process, and the construction delay between starts and completions (see also Spiegel [51]). The existence of lags in the demand and supply relationship of the office market in Stockholm is considered through application of error correction models in the study of Englund et al. [53]. The lagged values of rents, vacancy rates and employment variables were found to be significant in their three-equation system model. A similar study by Hendershott et al. [37] analysed long run and short-run rent, vacancy and supply adjustments in the US retail market. In the model estimations, various lags were considered for the variables utilised in systems of demand and supply equations. Some other recent examples on the existence of lagged variables utilised in explaining the adjustment processes in demand and supply relationship in commercial and industrial property markets can be found in Fuerst and Grandy [54], White and Ke [55], Mesthrige [56], Bruneau and Cherfauh [57].

Construction demand forecasting methods
Various statistical models have been established and tested for forecasting the demand for construction activities. There are two types of models, namely, the univariate and the causal models [58]. The univariate models are based on the past values of time series of the selected variable. The common univariate models applied in a wide range of studies to forecast construction demand are Box-Jenkins [16] and exponential smoothing techniques (see Jiang and Liu [59] for detailed review). The causal modelling techniques, by contrast, consider a number of related variables to explain the construction demand while forecasting the future construction activity. The causal models cover the classical multiple regression and advance multivariate regression techniques. Multiple regression models were applied for various scenarios to forecast sectorial construction demand in Thailand [60] to investigate the relationship between space demand and residential property investment [61] and private housing demand [62] in Hong Kong.
The advanced multivariate models to predict construction demand are vector autoregressive and vector error correction models, which can estimate a selected regressand based on its own lags and the lags of other variables. There are few studies in the literature applying these techniques for the construction demand forecasting as these techniques are mainly utilised for the estimation of various economic indicators including unemployment, interest rate, exchange rate, and housing prices (see Ng et al. [63] for a review). There are recent studies applying vector error correction modelling techniques to forecast construction demand including private construction investment [63,64] and labour demand [65] in the Hong Kong construction market, and changes in construction demand following the economic crises in Australia [59].
Despite the fact that vector error correction models captures both the short-term dynamics as well as the long-term equilibrium adjustment process among the variables, these models cannot explain spatial and market-based differences across the regions. These techniques assume general equilibrium throughout the entire land market, implying homogeneity of the estimated coefficients across different submarkets and different regions. The existence of submarkets and locational differences of the regions can indicate structural variations of demand for industrial land considering that each submarket has specific demand and supply conditions, as well as different locational characteristics. The consideration of these factors is complex in time-series analysis and further improvements are necessary for the estimation of dynamic models. In our model, we have utilised a dynamic panel data approach by including a number of location and area specific characteristics as potential instruments to consider spatial heterogeneity in the study area. In this, our model provides a significant improvement over time-series models (i.e. univariate and multivariate models) which do not consider such regional and spatial variations.

Data and methodology Methodology
Modelling industrial land market dynamics. The growth of economic activities (i.e. economic output and employment) in cities and regions results in physical changes in industrial land in the medium to long-term following an increase in demand for new industrial property developments. A considerable correlation between industrial and residential property and the economy is confirmed by the study of Wheaton [42]. A similar study by Ball and Wood [66] found evidence between aggregate economy and construction activity in the real property sector (see also Barras [52]). Therefore, the demand for industrial space can be represented as a function of economic activity in urban areas and regions. The increase in demand for new industrial space increases rental prices of the industrial estates in the short-run. Ball et al. [49] argued that changes in the nature of production and output generally results in more or less intensive use of space for a given level of output. Therefore, in the short-run firms do not change their demand for physical space but they adjust the existing space according to their new production needs.
In the economic literature, the effect of transport infrastructure on industrial land development is incorporated as an additional input in the production function since accessibility to certain activities raises productivity or reduces cost of production [67,68]. Therefore, accessibility can also be considered as an additional factor determining the demand for industrial land through influencing the profitability of industrial production. Given this literature, we can suggest that changes in economic activities, accessibility and rental prices of real property influence demand for industrial space in the market where the demand function can be stated as: where E i,t represents economic parameters (i.e. GDP, industrial GVA and industrial employment), A i,t is the accessibility indicator which will be elaborated in future sections, and R i,t is the industrial rental price of the i th region at time t. If there is excess demand for the existing industrial estates, it will rise the prices and rental prices of the industrial properties. The increase in rents and price for industrial space provides an incentive for developers to initiate new investments in the market. Therefore, developers will generate new construction investments in the long run. Property developers generally make investment decisions based on the prices in the recent past period [69]. Following Wheaton [42], we can state that supply depends on historical prices, P i,t-s, and asset prices at the time of the delivery, P i,t , as industrial properties are occupied by multiple tenants. The accessibility to industrial activities (i.e. jobs) is an important determinant of industrial land price, which consequently influences supply of land considering that land developers take into account the issues of accessibility of land uses and their price for their investment decisions in the market. A zone which has higher accessibility to industrial jobs imply that the subject zone has higher industrial land prices compared to others having lower accessibility values. Therefore, an accessibility measure, A i,t can be included in the supply function as an independent variable.
In order to obtain an understanding of the impact of industrial land conversions to other uses (e.g. social facilities, residential uses) on the development of industrial land use, we include a variable, C i,t , in the supply equation indicating the areas of industrial land existing at time (t-s) and converted to other uses at time t. Considering that land supply is also dependent on the developed space in the previous periods (t-s), land-use restrictions such as land zoning policies imposed by local or national authorities, Z i,t , and natural limitations such as protected areas, water bodies or sludge areas, N i,t , land supply function can be provided as: It can be suggested that industrial property supply is also determined by other factors such as construction costs, lending rates (i.e. cost of finance) and vacancy rates in the industrial property market. The focus of our study is to examine the key economic variables as well as supplyside parameters including zoning, natural restrictions and industrial land converted to other uses. The commonly used economic parameters are: GDP, industrial GVA, employment and property price (or rental price). Other factors determining the supply of land are out of the focus of the study.
In the equilibrium, demand is equal to supply and there is no incentive in the property market to change the activities of the economic agents. Solving the equilibrium of demand and supply equations of Eqs (1) and (2), a reduced form of new equation of developed industrial land area, LAND i,t representing the conditions of demand and supply in the industrial property market, can be formed as follows: As prices and rental prices are linked through capitalisation rates, there is a direct relationship between prices and rents in the market. Assuming that capitalisation rate remains unchanged, property price can be written as a function of rental price. However, in this study rental prices of Dutch industrial properties are not included in the analysis as the time-series of rental price data is unavailable at the regional level. As a future focus, we recommend rental prices to be considered in the analysis based on the availability of the subject data. Estimation methodology. Based on construction demand forecasting literature, the process determining the variable of industrial land depends on the set of some related past variables in the same way as industrial land depend on the set of its past values. In our sample, data for all the related variables is specified at the regional level. In this, we can proceed to generate a dynamic panel data model considering the spatial variations in our pooled time-series data. Eq (3) can be adapted to a dynamic panel data model with the lagged effects of the dependent variable i.e. LAND i,t-s incorporated in the error term (see Maddala and Lahiri [70] for details) where y it is the industrial land where the industrial construction is taking place, k represents the time lag, γ i is the panel-level effect; x i,t-k is a 1×k 1 vector of exogenous variables for each individual region i and at time t; α i (k) is k 1 ×1 vector of parameters to be estimated; ω i are I.I.D.
error terms over the whole sample with a variance σ 2 ε . The γ i and ω i are assumed to be independent for each i over all t.
As the explanatory variables may be correlated with the unobserved panel-level effects, the estimation of Eq (4) by using standard OLS methods makes the estimators inconsistent. Following Arellano and Bond [71], Arellano and Bover [72], and Blundell and Bond [73], the above-mentioned issues have been covered through the application of Generalised Methods of Moments (GMM) approach. Consequently, the estimators are constructed by first differencing Eq (4) to remove panel-level effects and using instruments to form moment conditions. The first-differenced version of (4) can be written as: Dε it ¼ rDε i;tÀ 1 þ Do it ; jrj < 1 in which the panel-level effects are eliminated by the first differencing operation. The GMM approach utilises instruments for the estimation of Eq (5). The details on instrumental variables are provided in the Appendix.

Data
The data used in the analysis consists of 40 regions at the NUTS3 (nomenclature of terrestrial units for statistics) level in the Netherlands, covering each year between 2000 and 2008. The NUTS3 coincide with the so-called COROP regions in the Netherlands, which has been originally designed to reduce cross-border commuting and provides a rough approximation of functional labour market regions in the Country [74]. Therefore, the NUTS3 level is selected as the spatial unit in the analysis considering that the NUTS3 level is the closest to the spatial structure of the Dutch labour and industrial land markets [18,75]. The industrial land-use data was obtained from land-use maps for Netherlands for the years 2000, 2003, 2006, and 2008 provided by the PBL Netherlands Environmental Assessment Agency (http://www.pbl. nl/en/) and were utilised to compute the areas of industrial land, natural protection and naturally restricted areas for land development, and large-scale urban development zones which incorporates industrial land uses [76] for each corresponding year. The areas of industrial land, protected and naturally restricted areas, and industrial land in development zones for those years where data is unavailable were calculated through applying a linear interpolation between known data points for all regions. This was done to provide continuity in the timeseries data regarding the areas of industrial land. The existing data points for the years 2000, 2003, 2006 and 2008 were individually checked for each NUTS 3 region, and a linear relationship was observed between the existing data points to a large extent. The areas of industrial land converted to other land uses were computed through application of a GIS-based analysis. In this approach, two spatial datasets (vector datasets) for the two subsequent years were compared regarding the industrial land uses existing in period t-s and converted to other uses in the subsequent period t. It is noted that industrial uses were mostly converted to social and recreational activities as well as residential uses. From this analysis, the areas of industrial land converted to other uses were obtained for the years 2000, 2003, 2006 and 2008. The unknown data points between these years were developed through linear interpolation.
The accessibility indicators were obtained from the work of Jacobs-Crisioni et al. [77] where potential accessibility measures were computed throughout Europe for each municipality within the modelled country and for the NUTS2 regions outside of the modelled countries (see also Jacobs-Crisioni and Koomens [78]). The potential accessibility measure is computed based on the formula: where A is the accessibility at origin municipality i; P is the population count in destination municipalities j in decade t; M represents the distance-decayed travel times in minutes; γ is assumed to be -1 as the real value of distance-decay is unknown for the study area. To compute the accessibility measures, domestic and foreign destinations (see Jacobs-Crisioni and Koomen [78] for the domestic and foreign accessibility component formulas) were taken into account; and the intra-zonal destinations were considered based on the Frost and Spence [109] approach in such a way that internal distance d j is assumed to be d j ¼ 0:5 The formula implies that intra-zonal distances are half of the radius of a hypothetically constructed circular zone. Travel times were obtained from the TRANS-TOOLS road network using the shortest path algorithm assuming free-flow travel times (for the details, we refer to Jacobs-Crisioni and Koomen [78]; Jacobs-Crisioni et al. [77]). The accessibility measures computed for the Netherlands are originally provided at the municipality level and were aggregated to the NUTS3 level to be included in the current study. As the accessibility measures were computed for each of the decade between 1961 and 2011 [78], the values for the years from 2000 to 2008 were developed through applying a linear interpolation between 1991, 2001 and 2011.
This study utilises gross domestic product (GDP), sectorial gross value added (GVA aggregated for each of the industrial sectors considered in this study), sectorial employment (EMP-in industry) and industrial property price (PRICE) data for each NUTS3 in the Netherlands from the online database provided by Eurostat [79] and the Netherlands Bureau of Statistics [80]. Some basic statistics on economic parameters is summarised in Table 1 for each NUTS3 region.
From the table, it can be seen that there are considerable differences across regions concerning annual average values of GDP, industrial property price, industrial GVA, and employment. The highest values are observed in the regions of Greater-Amsterdam, Rijnmond, Utrecht, Rest of Groningen, and South-East Brabant, which are housing majority of the economic activities in the Netherlands. An overview of the study area and its regional sub-divisions can be seen in Fig 1.

The Dutch industrial land market and land-use planning system
The lags between industrial demand and supply are sound in regions where there is scarcity of land for industrial development, strong regulations in real property markets, and strict institutional structures influencing the property transactions in the land markets. Netherlands is such an example of a region where there is scarcity of land ready-to-develop considering particularly the physical limitations in the Country (most of the land needs to be processed prior to development for the reasons of high water levels and low capacity of the soil to support the built structures in the country) strict land development restrictions imposed by various governmental layers, and a strict administrative system involving in various land-use planning activities. Land-use plans are highly restrictive in the Netherlands due to the fact that it is impossible to develop new urban areas without making changes in the land-use plans. Though land-use restrictions are important in the Country, it is also possible to have changes in the Examining lag effects between industrial land development and regional economic changes Examining lag effects between industrial land development and regional economic changes land-use plans. However, this requires a time consuming effort with the administrative procedures that may lead to uncertain outcomes. With respect to the institutional structure, there are three levels of government to lay down a strategic plan. These consist of: '(a) the national spatial planning decisions, (b) the provincial regional plans, (c) the municipal structure plans' [81]. While municipal structure plans are dictated by higher governmental levels, these structure plans prepared by the Dutch municipalities (complying with provincial and national plans) are the only legally binding and highly influential plans in determining the conditions in the land market (see Hajer and Zonneveld [81]; Needham and Louw [82] for a review of the institutional structure). The local land-use plans are also subject to a certain level of judicial flexibility and rarely overturned by higher planning authorities [30,83]. The growth of industrial development is facilitated by large-scale plans to develop new industrial land uses, and the municipalities are the main responsible bodies for the supply of industrial land. In Needham and Louw's ( [82]: 84) explanation, "in the Netherlands, the public provision of industrial land is seen as a means of implementing spatial, economic and employment policy" (see also Badcock [84]; Louw et al. [30]).
Considering this highly restrictive land-use planning system, an increase in land values would be expected as the supply of land diminishes as a result of land-use restrictions. The evidence shows that its final effect on urban development is mostly negative [85]. However, it is also evident that 'the spatial context and the competition between municipalities and jurisdictions can determine the outcomes of land-use regulation on urban development' [85]. Spatial competition for urban jurisdictions is sound in the Netherlands, where spatial planning is characterised by 'compact development' policies resulting high development densities and 'clustered de-concentration' aiming at concentrating suburbanisation at specific locations to increase urban densities in those particular locations [74,86].
While keeping high densities through successful management of urban development, the restrictive planning system has also led to successful open space preservation in the Country in the last five decades [87]. Nature conservation policies comprising, NATURA 2000-a European network of protected areas in the EU, and National Ecological Network (NEN)-a network of natural reserves areas put restriction on urban development to avoid conversion of protected land [88]. Koomen and Dekkers [89] suggest that there are more restrictive Buffer zone and Green Heart policies in the Netherlands, which cover important national landscapes as some of them were approved as world heritage sites by UNESCO. Regarding the preservation of Buffer zone and Green Heart, Koomen et al. [86] showed that between 1995 and 2003, the propensity of transformation from open space to urban use was significantly lower in these protected zones compared to other parts of the Randstad Area. Concerning the success in the application of open space preservation policies, a variable, namely natural limitations (see Eq 2 in the methodology section), is included in the supply equation to account for the impact of natural restrictions on urban development at the regional level. This variable represents the areas of natural protection zones such as NATURA2000 and NEN, and other natural restrictions for urban development such as water bodies and sludge areas in each NUTS3 region.
Besides open space preservation, large-scale urban development zones [76] introduced by the national spatial planning authorities in the Netherlands are highly successful in facilitating urban growth as these neighbourhoods also incorporate industrial and business areas aiming at reducing the commuting distance between new neighbourhoods and employment centres [85,90]. Though these large-scale plans were introduced by national government, they are mostly implemented at the regional and particularly municipal level with the involvement of private parties to invest in these new urban development schemes. Therefore, municipalities compete with each other to attract residents and developers of industrial and business areas to their territories [90]. There is not much difference among the regional authorities and the municipalities in the application of spatial plans. It is mentioned in Broitman and Koomen [87] that 'the previous national planning report has provided more freedom to regional and local authorities to meet their objectives' (they referred to: VROM et al. [91]). The industrial land in each of these large-scale development zones can be considered as a planning restriction that should be implemented by the corresponding local authorities. Therefore, we use industrial land area in each urban development zone for each corresponding NUTS3 region in the study area (see Eq 2) to represent the spatial planning policies of the government authorities implemented at the regional and municipal levels (see Levkovich and Rouwendal [85]; Broitman and Koomen [87]) for the inclusion of spatial planning policies in the models of the land markets in the Netherlands).
In the Netherlands, there is a large gap between development of a land-use plan and provision of the industrial land for development [30]. Given this strict institutional and regulatory structure, municipalities have often developed their own methods of forecasting the demand for industrial land for the provision of serviced land without delay. Considering the competition among the municipalities in attracting higher numbers of firms than their neighbours, the supply of industrial land is higher than the sale/lease of land on industrial estates in the Netherlands [83]. The over-supply of industrial land can also be explained by the existence of time lags between the industrial land use and economic parameters, and an asymmetric relationship between land development and economic change. This can be observed from our data on the percentage changes of industrial land, GDP, GVA, property price and employment in the industrial sectors (Fig 2). From Fig 2, it is clear that developments in the industrial land Examining lag effects between industrial land development and regional economic changes respond slowly to the fluctuations in the economic parameters. For instance, there are considerable changes in GDP, GVA, industrial property price and the number of industrial employment between 2000 and 2008, which are associated with modest changes in the industrial land developments.
Considering the discussed characteristics of the industrial land market, long-term forecasts of industrial land use have become an important issue in the Netherlands. Therefore, various methods are used in Netherlands as well as other countries for the estimation of future industrial land use (for detailed review, see Beckers and Schuur [18]). Knoben and Traa [92] (cited in: Beckers and Schuur [18]) stated that employment-based approach is the most commonly used methodology in Netherlands among others. A similar approach based on employment forecasts is also applied in European and non-European countries such as Ireland, UK, and Australia [93,94]. This methodology is based on future projection of employment by sector and multiplying these numbers with a land-use parameter to get the average size of land requirements per worker [18,94]. There are also other variables utilised in forecasting of industrial land such as economic output and income [10] and real property values [16]. The influence of a set of common economic fundamentals on the value of real property can be found in Hendershott et al. [95] and Wilson and Zurbruegg [96] (see Beckers and Schuur [18] for a detailed review of industrial land forecasting methods and the use of employment-based method in the Netherlands). Considering the importance of key economic parameters in forecasting industrial land in Netherlands and internationally, this study focuses on the relationship between industrial land development and these indicators to uncover the influence of lags of the subject indicators on the development of industrial land use.

Results from panel unit root tests
Pooled time-series data tend to exhibit a time trend implying non-stationary data processes due to the existence of time dependency across means, variances and co-variances of the variables. Therefore, research on testing for unit roots in panels have been developed [97,98]. It is demonstrated in the literature that the application of regression analysis to non-stationary data results in misspecification errors and misleading statistics such as high R 2 's and t-statistics [99]. Therefore, it is important to test stationary of the data prior to regression applications. To this end, we applied three different tests in order to test the stationary feature of our dynamic panel data. These consist of: Levin, Lin and Chu (LLC) test, Fisher-Dickey-Fuller (FDF) test, and Im, Pesaran and Shin (IPS) test. The first test assumes that there is a common unit root process across the cross sections. This follows an Augmented Dickey Fuller (ADF) specification to test the hypothesis of a unit root against the alternative that each time series is stationary [100]: where y it refers to the pooled variable, x mt represents exogenous variables in the model (i.e. individual time trends, spatial fixed effects), δ mi are the corresponding coefficients, p i is the lag order, and ε it are the error terms which are assumed to be mutually independent. By contrast, Fisher-Dickey-Fuller, and Im, Pesaran and Shin tests estimate a separate ADF regression for each of the cross sections to allow for individual unit root processes across the cross sections. Table 2 summarises the results from unit root tests on the subject variables embedded in the model. The variables were tested by using both the first and the second lag specifications Examining lag effects between industrial land development and regional economic changes (Eq 6), and the first differences of the variables were also tested. The results from these tests suggest that except the variable LAND, most of the tests concerning other variables fails to reject the existence of a unit root following the use of the first lag specifications. The first differences of some variables become stationary as reported by the LLC, Fisher-DF and IPS tests; however, this is not the case for the LAND variable. Regarding LAND, both the first and the second lag specifications rejects the null hypothesis across all three tests implying a stationary process of the panel data. Therefore, we employ the LAND variable in the model without considering its first difference as it does not follow a stationary process. Since the results from the three test statistics suggest that GVA, GDP, PRICE and EMP should be taken by following the second lag (as it is stationary by contrast to the first lag specifications), we continue our analysis by using the second lag specifications of the subject variables in the model.

Results from dynamic panel regressions
Regression analysis was conducted to estimate Eq (5) applying GMM. Prior to the regression analysis, partial correlation coefficients across explanatory variables were computed for the pooled data and those having significant correlation coefficients are presented in Table 3. Considering the collinearity between 'GDP' and 'GVA', GVA to GDP ratio was calculated and included in the analysis as an additional variable. There is significant collinearity between pairs of some other variables, which are also presented in the Table. These are various combinations of the variables including: 'ACCESS', 'EMP', 'PRICE' and 'CONV_LAND'. Therefore, these variables were included in the regressions either separately or they were treated with caution.
We have specified industrial land development as a log-linear function of exogenous variables as defined in Eq (3). In order to compare the impact of lagged values of predictors with the contemporaneous values, we estimated two sets of model specifications given as Where 'X' stands for economic parameters, 'Z' are the other parameters and 'k' represents the time lags as described previously. Tables 4 and 5 present the results for GMM system equations and the associated diagnostic test statistics for various combinations of variables selected as instruments. The tables specifically demonstrate the results concerning the regressions which do not show any misspecification errors following the diagnostic test statistics. From Tables 4  and 5, it can be noted that Eq (3) representing both supply and demand conditions is estimated in regressions B1, C1, B2, C2 and reduced forms of the equation are also estimated as presented  (2). We first discuss the selection of instrumental variables utilised in the estimations of Models 1 and 2 based on various combinations of explanatory variables. Basically, an instrumental variable is a variable that is uncorrelated with the error term but correlated with the explanatory variables in the estimated model [101]. A set of instrumental variables, which is correlated with the explanatory variables in Models 1 and 2, but uncorrelated with the error term in our dynamic panel data model was developed. These include the variables representing 'areas of mining sites, public facilities, airport and rail infrastructure', 'areas of land converted from industry to other uses', 'GDP', 'accessibility', 'price', and 'population' which were computed for each NUTS3 region within the period 1996-2008. The areas of mining sites, public facilities, and airport and rail infrastructure were derived from the land-use maps for Netherlands for the years 1996, 2000, 2003, 2006, and 2008. The data for missing years were linearly interpolated by using the known data points for all the NUTS3 regions.
The above outlined instrumental variables can represent area, industry and location specific characteristics in our study area. This can be explained, for instance, by the existence of a metropolitan or an urban area in a region implying that the region accommodates large number Wald Chi (2)   Examining lag effects between industrial land development and regional economic changes of population, highly skilled labour force and intensive physical infrastructure, directly linked with the growth of industrial properties. In the first difference equations regarding Models 1 and 2, we consider either two-period lagged values of the variables including mining sites, airport, rail, public facilities, price, population; one-period lagged values of accessibility, areas of land converted from industry to other uses, GDP, and/or the dependent variable i.e. LAND were lagged two periods, as potential Wald Chi (2)   Examining lag effects between industrial land development and regional economic changes instruments. This points to the assumption of weak exogeneity (Eq A2 in Appendix) concerning all exogenous variables. From the Sargan statistics given at the bottom of Table 5, the weak  exogeneity of the different variable sets (see the Table notes) is not rejected for all the regressions estimated for various combinations of the explanatory variables. These test results indicate that there are no misspecification errors in the selection and the use of the specified instruments in the estimations of Models 2. Table 4 reports the Arellano-Bond statistics to detect zero first-order and second-order correlations in the level residuals. The moment conditions (Eq A1 in Appendix) set for the estimation of dynamic panel model utilise the orthogonality conditions between the differenced errors and the lagged values of the dependent variable [100]. This is based on the assumption that the error terms in our dynamic panel model (Eq 5) are serially uncorrelated and the differenced error is with unit root [100]. These conditions are satisfied if the null hypothesis of the absence of first-order serial correlation is rejected while the null of the absence of secondorder serial correlation is not rejected. In Table 4, we present the regressions which satisfy the stated conditions. From Table 4, almost all the explanatory variables have expected signs except CONV_-LAND, and most of them are significant at the 1% level. The sign of the CONV_LAND variable is ambiguous as it has a positive sign in the regressions B1 and D1 and a negative sign in C1. It may stem from possible correlations between CONV_LAND and other explanatory variables included in the regressions. Land zoning policies and accessibility have a positive impact while natural restrictions for development have an expected negative impact on industrial land developments. Regarding estimations of Model 1, the coefficients for economic output, output-to-income ratio and price are larger than other predictors. For instance, the positive sign of GVA/GDP implies that high share of industrial GVA in national income results in an increase in industrial space. The B1 outcomes show that variations in industrial demand in Netherlands are mainly explained by price that is followed by manufacturing output to income ratio and then employment. When we exclude economic income and output from the regressions as in D1, industrial property price is considered as highly influential in explaining industrial land developments. The coefficients indicate the elasticity of estimated variables on the industrial land development, that is, a 1% increase in industrial employment corresponds with 0.05% increase in land use in Model A1, and with 0.34% in Model B1.
The second model specifications in Table 5 explain the impact of lagged economic predictors on industrial land use. We keep all the contemporaneous predictors the same as specified in the first group of models (Table 4) and added different lags of economic parameters to the corresponding models (see Eq 3) as shown in Table 5. From A2, B2 and C2 it is followed that the contemporaneous value of EMP and GVA/GDP becomes insignificant in contrast to Models A1, B1 and C1 (Table 4) when the lagged values of these predictors are added to the regressions. This points to the existence of development lags in the industrial property market when there is an upturn (or down turn) in the economy. The increase in economic activity results in an increase in user demand. However, firms tend to adjust the existing space according to the requirements of new production activities as suggested by Ball et al. [49]. Additional demand is absorbed by the vacant industrial properties in the short run and this will result in increasing land values and property prices. This may encourage developers to initiate new industrial property investments. Considering development constraints and long time period required for the delivery of the new construction projects, there is a lag between industrial land development and the growth in economic activities [25,49]. The insignificant coefficients of EMP and GVA/GDP confirm this fact. On the other hand, from A2 and B2, we observe an increasing influence of the lagged GVA on industrial land. The coefficients of GVA and GVA/GDP with two-period to eighth-period lags lie in the range of 0.008 to 0.13 in Model A2 and 0.027 to 0.30 in B2, respectively; the sixth-period lag indicating a strong reinforcing effect on the dependent variable in both models.
To analyse the possible impact of employment on industrial land, a separate model with current and lagged values of employment variable was also estimated. The coefficients of this model were found insignificant and did not satisfy the model specification test statistics addressing the bias in the estimated coefficients. Therefore, we estimated employment variable with the other variables to have robust estimates for the subject variable. The outcomes of these models indicate that the coefficients of EMP are quite significant for the two, four and six-period lags where the four-period lags in A2 and B2 demonstrate the highest elasticity on industrial land developments.
As can be followed from Table 5, most of the explanatory variables in the regressions have the expected signs with significant coefficients. Accessibility and land zoning policies have positive impacts on industrial land developments as can be followed from the Models presented in the Table. By contrast, A2 and C2 indicate insignificant coefficients for either contemporaneous value of EMP and GVA/GDP or second-period lag of GVA/GDP, which has an unexpected negative coefficient in the latter case. Natural restrictions have an expected negative impact on industrial land development in Models B2 and D2. However, it has an unexpected positive impact in Model C2. From C2 and D2, we note that the amount of industrial land converted to other uses has a negative impact on industrial land. The models incorporating PRICE in the regressions show that its current value has the largest impact on industrial land which is followed by four-period and six-period lags as highlighted by Models B2 and D2. This confirms the fact that property developers generally make investment decisions based on the prices at the time of delivery as well as the prices in the recent past periods [42,69].
To demonstrate the individual effect of GVA, EMP and PRICE on industrial space, we estimated various dynamic panel regressions for individual lags of the subject variables on industrial land. The results from these regressions are provided in Tables 6, 7 and 8. From Table 6, we observe a decreasing influence of GVA from two-period lag to the eight-period lag; twoperiod and four-period lags having the highest impact on industrial land developments  In parenthesis are the standard errors. GMM-Robust: Diff. Eq. Industrial Land (-2); Mining(-2), Airport(-2), Public Facilities(-2), Population(-2), all differences are forward-orthogonal deviations. https://doi.org/10.1371/journal.pone.0183285.t006 Examining lag effects between industrial land development and regional economic changes indicating the existence of lags between economic growth and land development in Netherlands. A similar pattern is observed with PRICE (Table 8) where the coefficients range from 0.15 to 0.34; the contemporaneous value of the variable having the strongest impact on industrial land. This coincides with the results presented in Table 5. As suggested by the theory, developers take into account the property price information in the very recent periods for their new investment decisions. Therefore, we suggest use of the contemporaneous value of PRICE as an indicator for the developers in the market regarding new industrial property investments. By contrast to previous findings, EMP (Table 7) has an increasing effect on industrial land, the eight-period lag indicating a strong reinforcing effect on the dependant variable.
GMM-Robust: Diff. Eq. Industrial Land (-2); Mining(-2), Airport(-2), Public Facilities(-2), Population(-2), all differences are forward-orthogonal deviations. https://doi.org/10.1371/journal.pone.0183285.t008 Examining lag effects between industrial land development and regional economic changes This implies that firms initially do not demand for new industrial land following an increase in number of employment but they prefer adjusting their existing industrial space to accommodate new workers. If there is excess demand for new space following the initial adjustments, firms demand for additional industrial land in the long run. The results in Table 7 demonstrate existence of lags between growth in employment and land development in the long run. These issues should be taken into account in the industrial land forecasting work which utilises employment-based methods as in the case of Netherlands.

Evaluation of forecasting accuracy
The forecasts generated by Models 1 and 2 were evaluated on three measures of accuracy, namely, the average absolute error (AAE), the total absolute error (TAE), and the relative difference (RD). These measures can be explained as follows. Average absolute error Where E i is the estimated industrial land, A i is the actual industrial land at region i; and n is the total number of NUTS3 regions. Using these measures, we can compare the forecasting accuracy of the models for the years 2000, 2003, 2006 and 2008 for the first group of models, and 2006 to 2008 for the second group. The results obtained for the models are shown in Tables 9 and 10.
Results in Table 9 show that Model A1 which utilises contemporaneous values of EMP, GVA and ACCESS, and Model B1 which incorporates supply conditions in addition to parameters representing demand, performed better than other models. Regarding Models A1 and B1, the TAE is within the limit of around 10%, which is generally considered as acceptable [65] and the AAE is the smallest among all others. From Table 10, it is observed that Model A2 is the best performer considering the smallest values of TAE and AAE among others and also showing one of the lowest relative differences. Models B2 and C2 perform better than Model D2 considering the TAE measures for B2 and C2 are consistently within the acceptable limit of 10% in contrast to the TAE for D2 which is more than 15%. This indicates that the model incorporating supply-side parameters as the only driver for land-use change is not able to perform better than other models which consider other economic variables in the predictions. The inclusion of all economic parameters to the predictions as in Model B2 (that is embedding both supply and demand conditions) and the inclusion of demand conditions as in Model A2 improves the prediction accuracy as the Models perform better than others.
Finally, when we compare the results in Tables 9 and 10 (Table 10) compared to their counterparts provided for Models 1 (Table 9). For the same years, RD is considerably smaller in Models 2 ranging from -0.07 to 0.026 than those in Models 1which range from -4.31 to 2.21.
The models can also be evaluated based on the prediction interval forecasts. Basically, a prediction interval forecast is an interval of values that is calculated according to a specified level   Table 1 for the details of regions). Regarding all NUTS3 regions, Models A2 and B2 generated closer estimates of industrial land to the actual land values compared to Models C2 and D2 predictions.
Finally, we forecasted future land-use demand for the Netherlands based on two alternative scenarios to illustrate the differences of the use of different lagged values of economic parameters on future demand for industrial land. The scenarios consist of a baseline scenario where economic growth trends are derived from OECD baseline scenario [102] linked with the study 'The Netherlands in a Sustainable World', by the Netherlands Environmental Assessment Agency; and an alternative High Development Pressure scenario, which is based on the Global Economy scenario in the study 'Welfare, Prosperity and Quality of the Living Environment' [103]. The latter scenario assumes large numbers of economic migrants from new member states into the Netherlands. The details of the scenarios can be seen in 'The Netherlands in the Future Second Sustainability Outlook Report' [103]. From the Report, the forecasted annual growth rate of GDP is 1.9 percent for the baseline scenario and 2.6 percent for the High Development Pressure scenario between 2002 and 2040. These growth rates were applied to obtain the future values of GVA for the Netherlands for the 2009-2020 period, and utilised to forecast industrial land demand from the regressions that we estimated to identify the impact of individual lags of the GVA variable on industrial land use ( Table 6). The forecasted industrial land demand from this analysis is presented in Fig 5. From Fig 5(b), High Development Pressure scenario results in higher industrial land particularly when contemporaneous value, and second and fourth lags of GVA are used in the land demand forecasting. Both Fig 5(a) and 5(b) demonstrate that industrial land demand is the lowest when eight-period lag of GVA is utilised in the analysis. It is also confirmed in Table 6 that the estimated coefficient of eight-period lag of GVA is the smallest compared to others. The trends observed in Fig 5 indicate uncertainties in the projected industrial land use depending on the use of different lagged values of GVA in industrial land forecasting work. This implies that land-forecasting studies focusing on a single economic parameter for the projection of future land demand should consider the uncertainties stemming from the use of different lagged parameters in the forecasting analysis.
In future projection of industrial land demand, we have considered socio-economic forces as they were embedded in the alternative scenarios [103] that were developed to forecast the values of economic parameters (e.g. GDP). However, individual characteristics of industrial firms in the post-2016 period, which are different than existing firms, are not captured in the current analysis. During recent decades, advanced technologies accompanied by advanced knowledge and innovation systems have been emerging with an increasing pace. An example for an advanced production systems is 'additive manufacturing (or 3-D printing)' which is characterised by a process where new locations of economic activity and innovation emerge following an unprecedented change in technology [104]. This process is seen to be sensitive to provision of services, firm scale, clustering and networks [105]. Gress and Kalafsky [106] noted that for an effective production, additive manufacturing may be the catalyst for the Examining lag effects between industrial land development and regional economic changes movement of such manufacturing closer to the research sources. The existing advanced knowledge synergies and networking may reinforce the existing clusters, and together with new innovations, existing clusters could be transformed to house the activities of advanced manufacturing [107]. 'This is happening in geographically centralised manufacturing centres Examining lag effects between industrial land development and regional economic changes including Belgium, Germany and the Netherlands' [106]. Innovation and advanced manufacturing may take place in other locations as well through the construction of new knowledge synergies, innovation and networking.
EC ([110]: 6) identifies 'environmental' and 'organisational' factors that lead high-tech firms' location decisions with respect to production, R&D, and innovation. Environmental factors comprise: a) location characteristics (input cost factors, social factors, scientific and technological strengths, political stability); b) industry, sector, value chain characteristics (the value of being near an ecosystem with suppliers or other firms in the same sector); c) market considerations (proximity with important markets and key customers). Likewise, organisational factors consist of: a) strategic considerations (firms' strategy with respect to market entry and growth, product/market combinations, production organisation, distribution and services); b) technology and innovativeness (intensive collaboration between production, R&D and innovation and a more pro-active strategy towards knowledge and technology sourcing); c) product and production complexity (has a strong influence on the need to colocate production with development and research activities).
Concerning the representation of location decisions of firms using advanced manufacturing techniques in the analysis of industrial land development, it is useful to analyse the existing clusters of high-tech firms, research and innovation centres, their knowledge synergies and networking among them. The size and geographical location of these clusters and highlyskilled employment working in these high-tech industries are crucial parameters that could be included in the analysis to identify the location and land-use requirements of high-tech industries, research and innovation centres. Future projections on demand for advanced manufacturing outputs, geographies of consumption, required services to support advanced manufacturing industries, geographies of labour, innovation and R&D facilities are required to be incorporated in the projection of industrial land regarding the land requirements of high-tech production activities. We leave this as a future research topic depending on the accessibility of required data on the geographies of advanced manufacturing and innovation systems to be included in the analysis of industrial land use in the Netherlands.

Discussion and conclusions
The findings of this study are policy relevant, as many forecasting models for real property demand analysis are using socio-economic indicators to predict industrial land development. The existence of time lag effects concerning particularly the socio-economic parameters is important to be considered in the models that focus on prediction of demand for industrial property. It is of particular significance for the land-use models where a reliable estimation of future land demand is required to be integrated into spatial land-use models. It is also important for policy makers and planners to formulate timely and forward-looking strategies to distribute and allocate scarce land resources efficiently. While previous research on real estate demand forecasting focused on cross-section data of the socio-economic parameters-that are limited to reflect cyclical fluctuations, it is important to develop reliable forecasting models which can provide the impact of economic fluctuations on land-use development. The results of this study should provide a basis for researchers, planners and policy makers to explore the applicability of our dynamic panel regression technique in modelling industrial land demand during the periods of economic boom or recession. The use of time lag effects by the model could serve as a useful source for construction and industry stakeholders to implement appropriate strategies during periods of cyclical fluctuations.
There are two key findings of the study: First, our models highlight the importance of economic output, national income, industrial property price and employment in influencing industrial land developments. In this, it is shown that fluctuations in the economic fundamentals are closely interrelated to similar fluctuations in the development of industrial land. This is in line with the findings of the economic literature outlined previously. Second, our findings suggest that industrial developments are mainly explained by the existence of time lags between the changes in industrial property price, industrial GVA (or GVA to GDP ratio) and employment, and the changes in industrial land developments. Considering time lags between demand and supply which can be related to inefficiencies of the land markets in the Netherlands, fluctuations in demand for industrial property are not immediately reflected in industrial developments.
Our dynamic panel regression models highlight the varying contribution of economic fundamentals to the industrial land development at different times. The models also successfully capture the regional and market-based heterogeneity at different times throughout the use of NUTS3 regions as the spatial unit of the analysis and through the utilisation of instrumental variables in the dynamic panel regressions that can represent socio-economic and location specific characteristics of the regions. Therefore, our modelling approach provides an improvement over the time-series models discussed in section 2.2. However, there are still issues to be covered in our modelling approach such as the inclusion of agglomeration, spillover effects and networking representing the knowledge synergies of high-tech industrial developments on land use. Future research can make a contribution to the model by focusing on such impacts to explain regional variations in industrial land developments.
Concerning industrial property price indicator, in contrast to the lagged values, the current value of property price is the key driver of industrial developments as it has larger coefficients compared to those with two to six-period lags. This implies that property price can be considered as a short-term indicator in the Dutch industrial property market in such a way that it gives signals to property developers whether to invest in a new industrial development scheme requiring expansion of industrial land. A further implication of the findings suggests that industrial property prices can be used to forecast short to medium-term demand for industrial land. However, the indicator is insufficiently robust to forecast the long-term future demand of industrial land use.
Although the first group of models find significant coefficients for the current values of output and income indicators, our results suggest caution with the consideration of industrial GVA and GVA/GDP in current values (as in Table 5) to forecast the demand for industrial space. A further issue is that contemporaneous values of these predictors in the second group models can explain only a small share of the variation in industrial land demand considering the existence of lagged effects of the predictors in industrial developments. Therefore, consideration of lagged values as well as the current values of economic indicators in the industrial land forecasting work will improve the robustness of the analysis. Furthermore, our models suggest that output and income indicators can appropriately explain short to long-term (up to eight years) relationship between the lagged values of subject indicators and industrial land use. This is crucial to be considered in the studies which focus on economic indicators (i.e. GVA, GDP) to forecast future land demand. The examples include Reginster and Rounsevell [41], Rounsevell et al. [14], Hoymann [13], and Batista e Silva et al. [10] among others.
A significant issue to be considered in the forecasting models is the influence of economic shocks on the development structure of industrial land. The post-2006 crisis has changed the structural relationship between economic parameters and land use. Therefore, the possibility of an economic shock in relation to the persistence of its effects to the time series plays a role in the demand-supply relationship of the industrial property. A further issue to be considered is the unprecedented reduction of population and employment in some regions in Netherlands which will come in the succeeding decades following a long period of uninterrupted growth [108]. This may lead to structural changes in labour and real property markets considering that business sectors may adopt their demand for industrial property and their location choice to the new economic conditions [18]. This would result in unprecedented effects on the longterm demand for industrial real estate. It is therefore important to update the models at various times in such a way that unforeseeable changes are fully reflected and the model forecasts are generated with high accuracy. For this reason, land forecasting applications should consider the influences of shocks in the economy and structural changes in labour and real property markets through the consideration of full time series data particularly for the post-2006 period.
The forecasting models have been increasingly used by national and local planning authorities to predict and allocate land uses on various sites in different regions. In particular, this is of great significance for the Netherlands due to the physical limitations of land development in the country, and due to the existence of strict institutional and regulatory structures influencing industrial land markets. We should note that the research presented here is not confined to the Netherlands and should be extended to cover other regions and countries in terms of further improvements in the methodology and applications with regard to industrial land analysis. However, it should be also noted that the application of current analysis requires longterm time-series spatial data regarding industrial land-use. Though long-term socio-economic data at the regional level is widely available in most countries, there are data accessibility issues for the time-series spatial data. This may limit application of the current study for other countries at the regional or local levels.
Though our model is able to incorporate spatial heterogeneities observed across the regions in the study area, the model cannot capture the spatial dependencies observed at neighbouring locations. As an alternative methodology for introducing spatial lag effects -that is based on the notion that the value of a variable at a given location is related to the values of the same variables measured at nearby locations, we refer to spatial error modelling approach, also known as spatial autoregressive model. This methodology can be used as an alternative approach in the future work for the prediction of demand for industrial land. However, it is important to mention that the model is unable to incorporate time lag effects of the dependant and explanatory variables. Therefore, a methodological development is necessary regarding the spatial error models to account for the impacts of time lags of the dependant and selected independent variables. This will be a challenging work as sophisticated techniques are required for methodological improvement.
This leads to the matrix Z i i.e. the matrix of instrumental variables: The GMM approach discussed so far utilises the moment conditions set by the Eqs (A1) and (A2) given the first-differenced Eq (5). The consistency of the GMM estimators relies on the correct specification of the model and the moment conditions, which in turn, are influenced by the first-differenced disturbances being serially uncorrelated and the explanatory variables being exogenous to the system of equations.
Supporting information S1 File. Panel data for the analysis of Netherlands. (XLS)