Fully Dynamic Input-Output / System Dynamics Modeling for Ecological-Economic System Analysis

The complexity of ecological-economic systems significantly reduces our ability to investigate their behavior and propose policies aimed at various environmental and/or economic objectives. Following recent suggestions for integrating nonlinear dynamic modeling with input-output (IO) modeling, we develop a fully dynamic ecological-economic model by integrating IO with system dynamics (SD) for better capturing critical attributes of ecological-economic systems. We also develop and evaluate various scenarios using policy impact and policy sensitivity analyses. The model and analysis are applied to the degradation of fish nursery habitats by industrial harbors in the Seine estuary (Haute-Normandie region, France). The modeling technique, dynamization, and scenarios allow us to show trade-offs between economic and ecological outcomes and evaluate the impacts of restoration scenarios and water quality improvement on the fish population.

This complexity significantly reduces our ability to understand the way ecological-economic systems behave and provide useful policy insights [3,8].Furthermore, because ecological-economic systems are often viewed as non-separable [1], splitting them in two would be "a poor choice of boundary" [2] and could lead to severe misperceptions, policy failures, and undesirable or counterproductive outcomes [2,7].
Tackling such complexities is only possible via modeling and computer simulations.The ecological-economic model used in this paper captures the complexity of both ecological and economic systems.Practically speaking, dynamic modeling is used to capture the nonlinear feedback dynamics [9].Recently, integration of such nonlinear dynamic modeling (e.g., SD) with an IO analysis [10][11][12] has been suggested.While dynamic modeling captures the nonlinear dynamics of ecological-economic systems, IO enables the study of sectoral impacts [13,14].However, there have been few attempts to integrate these two modeling approaches [15][16][17].
Obviously, depending on the type of analysis undertaken, other modeling techniques can be used.For example, if the focus is on individual heterogeneity, an agent-based (AB) model [18] may be preferred to a more aggregate-level dynamic modeling approach such as SD [6,[19][20][21].A hybrid AB-SD model can also be considered [22,23].
Our study develops a fully dynamic ecological-economic model by extending an IO/SD model initially developed by [16] for the ecological-economic system surrounding the Seine estuary.The primary contribution of the previous paper lies in the synchronization of IO and SD, whereas ours is to fully dynamize the model.More specifically, technical coefficients within the IO model will become dynamic along with the relationships between the ecological and economic systems and within the ecological system.Constant technical coefficients have been subject to important criticisms against the use of IO modeling [13].A fully dynamic ecological-economic model better capturing important nonlinear dynamic behaviors and feedback is in line with a suggestion by [12].In addition, to capture a qualitative aspect of the ecological system and reflect a critical uncertainty outside the system boundary (i.e., not endogenously determined but externally given in the model), our model adds water quality parameters for the Seine River.
In addition to setting up a model, we also propose scenario developments comprising policy impact and policy sensitivity analyses as methods to examine the model in a more systematic way rather than provide an ad hoc analysis as in [16].It is another primary contribution of our study.
The remainder of the paper is organized as follows.Section 2 develops the model (IO and SD), details its dynamization, and considers various policy scenarios.Section 3 outlines the results, assesses the impacts of various policies, and undertakes sensitivity analyses.The last section concludes and proposes topics for further study.

Study Area
We apply the IO/SD model to the case of degradation of fish nursery habitats by industrial harbors in the Seine estuary (Haute-Normandie region, France).It is well documented that fish nurseries are at risk when harbors develop [24,25].From 2002From -2004, the harbor of Le Havre (Grand Port Maritime du Havre) added 10 km of dykes to extend new infrastructures on the sea to adapt to current worldwide development of huge capacity container ships.The Grand Port Maritime du Havre is the largest harbor in France in terms of container ship traffic, which would have probably been difficult to achieve without the 2002-2004 extension project; and it is France's second largest harbor for crude oil imports.Upstream in the Seine estuary, the harbor of Rouen is also very important (Europe's first for cereal exports and France's second for refined petroleum product transportation).Those two harbors therefore provide a large amount of direct employment and hence, have a critical economic impact in France [26].Several studies have analyzed ecological impacts of harbors and other human activities on nursery areas in the Seine estuary [16,27,28].

Model Development
As this is an extension of an earlier model [16], readers interested in the technical details can refer to the previous paper, which also includes the full model description of the SD component in Powersim language in the supplementary information (S4).The model in digital format is also available from the authors upon request.Here we restrict our explanation to the main features of the model and concentrate on its dynamic extension.Following the convention of the SD approach [29], we conduct various model tests as part of the model development (i.e., boundary adequacy, structure assessment, dimensional consistency, extreme condition, integration error, and sensitivity analysis).

Model Overview
Figure 1 is a simplified representation of how the model captures the ecological-economic system of the study site.It displays the key variables only to highlight the main relationships within the ecological-economic model.
The economic system is embedded in the ecological system as a sub-system.The ecological system is modeled with Powersim (Powersim Studio 10, http://www.powersim.com/),an SD software, and the economic sub-system (IO) is modeled with Microsoft Excel.We integrate the economic sub-system into the ecological system by utilizing Powersim's function to connect to various datasets, including Excel.The integration is not only capable of transferring data but also allows both Powersim and Excel to run computations internally at each time step of the simulation.Although most economic variables are modeled in Excel (inside the dashed box), some are modeled in SD (Powersim) for technical efficiency.
"Restoration rate" and "Soles caught originating from the internal part of the Seine estuary" are the key variables that connect the ecological system and economic sub-system.Although a higher "Restoration rate" quantitatively improves "Nursery areas (in the internal part of the Seine)", it incurs a "Cost of environmental measures"."Soles caught originating from the internal part of the Seine estuary" involves various feedback loops that are a source of nonlinear dynamics of the model.B1 and B2 stand for negative or balancing feedback loops [29].For instance, a larger "Sole stock from the internal part of the Seine" increases "Catchable stock", which results in larger "Intermediate domestic consumptions" that drives more "Soles caught originating from the internal part of the Seine estuary".This increase, however, dampens "Sole stock from the internal part of the Seine", and so on."Soles caught originating from the internal part of the Seine estuary" involves other feedback loops such as that of "Dynamic IO".It is, however, not clear whether the feedback is positive or negative because the impacts of "Final demand for Sole (Foreign and Domestic)" on IO are mixed.1 is a simplified representation of how the model captures the ecological-economic system of the study site.It displays the key variables only to highlight the main relationships within the ecological-economic model.
The economic system is embedded in the ecological system as a sub-system.The ecological system is modeled with Powersim (Powersim Studio 10, http://www.powersim.com/),an SD software, and the economic sub-system (IO) is modeled with Microsoft Excel.We integrate the economic sub-system into the ecological system by utilizing Powersim's function to connect to various datasets, including Excel.The integration is not only capable of transferring data but also allows both Powersim and Excel to run computations internally at each time step of the simulation.Although most economic variables are modeled in Excel (inside the dashed box), some are modeled in SD (Powersim) for technical efficiency.
"Restoration rate" and "Soles caught originating from the internal part of the Seine estuary" are the key variables that connect the ecological system and economic sub-system.Although a higher "Restoration rate" quantitatively improves "Nursery areas (in the internal part of the Seine)", it incurs a "Cost of environmental measures"."Soles caught originating from the internal part of the Seine estuary" involves various feedback loops that are a source of nonlinear dynamics of the model.B1 and B2 stand for negative or balancing feedback loops [29].For instance, a larger "Sole stock from the internal part of the Seine" increases "Catchable stock", which results in larger "Intermediate domestic consumptions" that drives more "Soles caught originating from the internal part of the Seine estuary".This increase, however, dampens "Sole stock from the internal part of the Seine", and so on."Soles caught originating from the internal part of the Seine estuary" involves other feedback loops such as that of "Dynamic IO".It is, however, not clear whether the feedback is positive or negative because the impacts of "Final demand for Sole (Foreign and Domestic)" on IO are mixed.
We added "Water quality" as part of the system because of its significance, which is later discussed.It is, however, exogenously given without a feedback relationship in the model because the dynamics and effort of water quality improvement are outside the boundary of the ecologicaleconomic system in the study site (a large part of the contamination in the Seine estuary comes from the city of Paris located at 350 km upstream).Interconnections between the economic sub-system and ecological system.The figure is adapted from [16] with primary model extensions underlined in red.Note: "+" and "−" indicate variable changes in the "same" and "opposite" direction, respectively.Interconnections between the economic sub-system and ecological system.The figure is adapted from [16] with primary model extensions underlined in red.Note: "+" and "−" indicate variable changes in the "same" and "opposite" direction, respectively.We added "Water quality" as part of the system because of its significance, which is later discussed.It is, however, exogenously given without a feedback relationship in the model because the dynamics and effort of water quality improvement are outside the boundary of the ecological-economic system in the study site (a large part of the contamination in the Seine estuary comes from the city of Paris located at 350 km upstream).[16] used the commodity-by-industry IO table for the study area (Haute-Normandie region) comprising 12 commodities and 12 industries for 2007, whereas this paper uses the industry-by-industry IO table for the same area comprising 37 sectors for 2012.Since the regional table does not exist, we construct it for the area studied using the 2010 French national table [30], update data to 2012, and regionalize the table with techniques developed by Jackson [31], Lahr [32], and McDonald [33].
There are several reasons for using industry-by-industry IO tables and 37 sectors.First, detailed sectors allow more precise (sectoral) policy analyses.Second, regionalization of the industry-by-industry IO table is more robust and less time consuming than that of the commodity-by-industry IO table.Finally, the calculation of technical coefficients in regionalized industry-by-industry IO tables can help avoid inconsistencies (i.e., negative values of coefficients) that may occur with commodity-by-industry IO tables.
The IO table (Table 1) comprises three matrices-X, the intermediate sales matrix; F, the final demand matrix; and V, the value added payments matrix-five vectors-x, x , v, mi , mf , representing total industry output, its transpose, total value added payments, a row vector of imports consumed by industries, and a row vector of imports consumed as final demand-and one scalar m for total imports.The time notation (t) is suppressed when it is not necessary.Selling Sector (i = 1, . . ., n; n = 37) We can derive the following relationship from the industry-by-industry IO table (Table 1), which calculates sectoral output (x) based on static technical coefficients: I and i are respectively an identity matrix and column vector of 1's known as a summation vector; f = Fi; (I − A) −1 is renamed L for conciseness and known as the Leontief inverse or total requirement matrix; and A is the matrix of technical (or IO, or direct input) coefficients made of elements a ij = x ij /x j = x ij /x i [14].Equation (1) calculates the direct and indirect impacts of changes in the final demand (F) on the industry outputs (x).The calculations of the final demand for each sector follow [16] with three adaptations.
First, the final domestic demand for sole is calculated separately from the demand for the 36 other sectors i to relate sole consumption to environmental conditions and measures.The final domestic demand for sole is calculated in tons in the economic sub-system as follows: where f t i=sole, 1 tons enters the SD model and depends on the income elasticity (The income elasticities are adopted from Gohin [34].)(e i ) as well as on changes in household disposable income (Y t ) from t − 1 to t, which in turn is a function of the cost of environmental measures (ψ t i ) paid by industries.Second, the cost of restoration is allocated among economic sectors as follows.The investment (k = 4) is defined as where is the total amount of the gross operating surplus used for investment (other than investments in nursery restoration), µ j is the fixed coefficient of the share of gross operating surplus of sector j used for investment in the output of sector j at the reference year, 1 2 ψ t−1 α j inv is the part of the total annual restoration cost paid by sector j through a reduction of the part of its gross operating surplus that would have been used for non-restoration investments otherwise, ψ t is the total annual restoration cost at time t, ψ t i is the restoration cost paid at time t by other sectors to sector i = 18 (the construction sector) to build infrastructures required for nursery restoration; ψ t i = 0 for i = 18, α j is the share of the total annual restoration cost paid by sector j, and ĉap i is the fixed capital formation coefficient calculated in the IO table at the reference year as the ratio of investment in sector i on total gross operating surplus used for investment.
The third point worth mentioning here is that the sectoral output x t−1 j is not calculated via Equation (1) as in most IO modeling publications.Next, Section 2.2.2.2 explains how we dynamize x t−1 j using Equation (7).This allows us to increase the dynamic property of the equations.In the scenario development detailed next, we set two cost allocation rules by changing the combination of α j .

Dynamization of the IO Table
Changes in output ∆x caused by changes in final demand ∆f are expressed as: If constant technical coefficients could be a reasonable assumption for a short-term period (~5-10 years), it is much less so in the longer term.To overcome this issue, several approaches have been proposed [35]) and we use a classical method adopted in various IO-EC models [13,[35][36][37].
The key idea is to capture the difference between expected and actual output.To do that, we need to differentiate sectoral output x from Table 1 in two types: (1) the expected (or predicted) output, a vector of sectoral expected output z i (vector z) that is conditional on the base year IO table and, hence, contains the deterministic structure of this table (i.e., x ij values in Table 1, which do not change over time); and (2) the actual output, a vector of sectoral actual output x i (x) that is the historical or forecasted value and expressed as a function of constant price expected output (z).Please note that in dynamic IO-EC models, x is a vector of historical sectoral outputs at the stage the model is being built, that is, when modelers use observed time series data to build up statistical regressions that compute x in the current year (time t) as a function of the difference in the previous year (time t − 1) between observed (also named actual) values of x and values estimated (also named expected) by a static IO model.However, once the building stage is completed, the dynamic IO model can be run for forecasting to simulate future sectoral outputs x.This is why x is considered a historical (observed) value at the building stage and a forecasted value for forecasting purposes.The term "forecasted" is voluntarily used here to make it clear that it includes dynamic changes in the structure of the economy, which is not the case when we use the terms "estimated" or "expected." Both sectoral outputs z and x are computed at the constant price.Unless otherwise mentioned, all prices in this paper are expressed in 2012 Euros.Following Kim et al. [13] and Israilevich et al. [37,38], z and x are expressed as follows: A is a matrix of constant technical coefficients, which is taken from the base year; f is the final demand vector; and z is identical to x in the base year for which technical coefficients (i.e., A) are known but then generally differs from x over time as coefficients change.
The Equation ( 6) relationship can be rewritten as (see [37,38] for detailed description; f in Equation ( 6) represents the general symbol for mathematical functions, not to be confused with the final demand vector f.): where βt is a diagonal matrix with the elements βi,t , computed as: This transforms the static IO equation into a dynamically determined relationship.The coefficients for βi,t (α 0 , α z and α g ) are estimated using generalized least squares (GLS) to correct for first-order autocorrelation in the residuals in Equation ( 8) based on historical values (i.e., time series of observed data for x t , f t , and A) and on z i,t−1 computed in Equation (5).
g i,t is the set of exogenous explanatory variables selected by the modeler, which vary across i (e.g., sectoral value added, gross domestic product (GDP), or employment) and z i,t−1 is a lagged input-output-generated predicted output.Equation ( 9) explains the difference between the actual output (x i,t ) and the expected output (z i,t ) (ratio on the equation's left side) as a function of the ratio of actual and expected outputs of the previous period ( ), exogenous variables (g i,t ), and a stochastic component, ε i,t .In other words, Equations ( 6) and ( 7) capture overall changes in technical coefficients over time.
Figure 2 shows overall changes in technical coefficients (i.e., the share of an intermediate input, x ij , consumed by a sector j in its total input consumption, x j ) from 2012 to 2032.The changes are calculated by subtracting the value of technical coefficients in the reference year 2012 (also called the base year) from their value in 2032.Values for the reference year 2012 are obtained through the regionalization technique of the national IO table as mentioned previously.Technical coefficients for 2032 are estimated using dynamization.
employment assistance; travel agencies and reservation services; building services, etc.).This also suggests that the dynamization equation of technical coefficients captures the increasing importance of the functionality economy with respect to the conventional economy.That is, industries today prefer renting their equipment to buying them (e.g., car leasing, photocopy renting, computer leasing, etc.) or subcontracting services to other companies rather than supplying them on their own.The model has two stocks: "Nursery areas" and "Sole stock from the internal part of the Seine."The nursery areas include 21 categories with different sole abundance to represent the spatial heterogeneities; the categorization is based on the sediment type-gravel, sand, or silt-and the depth.These areas are assumed to be independent, as we do not know how they interact with each other.Another simplification of the ecological model stands with the fact that we do not develop the physical, chemical, and biological conditions required by soles to reproduce, e.g., the spawn is performed in the bottoms where soles inhabit and therefore depends on environmental conditions such as water temperature and acidity, water quality and food availability [42].The sole stock uses a cohort structure, ages 1 through 10.In line with the findings from a similar study on the Chicago metropolitan region by [38], changes are modest (mean: −0.0001; standard deviation: 0.0121).Dark-and light-colored cells respectively show positive and negative changes from 2012-2032.To highlight the changes, we set cells with very little change (i.e., absolute value of changes smaller than 10 −17 ) in white, along with cells in which no change occurs.
Negative changes are dominant in the top left industries of Figure 2 (cells from columns and rows 1 to 18).This suggests that our dynamization equation of technical coefficients succeeds in capturing the current trend of progressive greater degrees of replacement of French industries by those located in developing and emergent countries in which labor costs are much lower (China, central and Eastern Europe, Brazil, etc.).This is particularly the case on the diagonal showing self-consumption within the same sector.For instance, regarding the cell in row 5 and column 5 describing self-consumption within the "wood, paper, printing" sector, the decline is likely because the paper industry (in column 5) is importing ever more wood and paper pulp (from row 5) from Brazil and other emergent and developing countries, as suggested in time series statistics from FBCA [39], Copacel [40], and INSEE [41].
Positive changes are dominant in the bottom right industries of Figure 2 (cells from columns and rows 22 to 37).This suggests that the dynamic technical coefficients reflect the current trend of the increasing rate of services consumption by other service sectors.For example, most service sectors have increased their share of intermediate inputs supplied by sector 30, "rental and leasing activities, etc." (which also includes many other services to companies and businesses such as interim employment assistance; travel agencies and reservation services; building services, etc.).This also suggests that the dynamization equation of technical coefficients captures the increasing importance of the functionality economy with respect to the conventional economy.That is, industries today prefer renting their equipment to buying them (e.g., car leasing, photocopy renting, computer leasing, etc.) or subcontracting services to other companies rather than supplying them on their own.3 depicts the SD part of the IO/SD model in a stock and flow diagram.The model has two stocks: "Nursery areas" and "Sole stock from the internal part of the Seine."The nursery areas include 21 categories with different sole abundance to represent the spatial heterogeneities; the categorization is based on the sediment type-gravel, sand, or silt-and the depth.These areas are assumed to be independent, as we do not know how they interact with each other.Another simplification of the ecological model stands with the fact that we do not develop the physical, chemical, and biological conditions required by soles to reproduce, e.g., the spawn is performed in the bottoms where soles inhabit and therefore depends on environmental conditions such as water temperature and acidity, water quality and food availability [42].The sole stock uses a cohort structure, ages 1 through 10.The connection between SD and the economic sub-system is represented by bold arrows.Such an arrow leaving a circle (hence, leaving an auxiliary variable) means that the value of that variable is transferred to the economic sub-system.In addition, by correspondence, a bold arrow entering the circle means the value is transferred from the economic sub-system to the ecological system.
Based on Cordier et al.'s [16] work, we update parameter values and add two new components: varying restoration schedule and water quality.
First, concerning water quality, its improvement reflects the qualitative aspect of the nursery.The degradation of water quality is a major anthropogenic disturbance to soles in the Seine estuary [25].Using the backward predictive approach, Rochette et al. [25] estimated that the juvenile density in the Seine estuary is approximately 23% lower today compared to 1850, primarily because of water quality degradation.Consequently, we assume that water quality improvement can recover the juvenile abundance by 23% maximum as: * where Water Quality Index: [0.50, 1.00].
Here, we simply assume that the juvenile abundance is proportional to water quality by conducting a sensitivity analysis.The current (or business as usual) Water Quality Index is set at 0.77 (=100-23%).Further improvement could be expected, as the water quality of the Seine River has shown significant improvements regarding phosphate and ammonium pollution since 2000, when the European Water Framework Directive came into force [43].However, those authors also point  The connection between SD and the economic sub-system is represented by bold arrows.Such an arrow leaving a circle (hence, leaving an auxiliary variable) means that the value of that variable is transferred to the economic sub-system.In addition, by correspondence, a bold arrow entering the circle means the value is transferred from the economic sub-system to the ecological system.
Based on Cordier et al.'s [16] work, we update parameter values and add two new components: varying restoration schedule and water quality.
First, concerning water quality, its improvement reflects the qualitative aspect of the nursery.The degradation of water quality is a major anthropogenic disturbance to soles in the Seine estuary [25].Using the backward predictive approach, Rochette et al. [25] estimated that the juvenile density in the Seine estuary is approximately 23% lower today compared to 1850, primarily because of water quality degradation.Consequently, we assume that water quality improvement can recover the juvenile abundance by 23% maximum as: Abundance Multiplier i = Water Quality Index * Abundance Multiplier max i (10) where Water Quality Index: [0.50, 1.00].
Here, we simply assume that the juvenile abundance is proportional to water quality by conducting a sensitivity analysis.The current (or business as usual) Water Quality Index is set at 0.77 (=100-23%).Further improvement could be expected, as the water quality of the Seine River has shown significant improvements regarding phosphate and ammonium pollution since 2000, when the European Water Framework Directive came into force [43].However, those authors also point out that water quality is decreasing regarding nitrate pollution.Therefore, we set the plausible range of Water Quality Index to [0.50, 1.00], following Sterman's [29] warning not to be overconfident about uncertainty and therefore to include a relatively wide range.Our estimation of water quality seems sufficient for this type of analysis but hydrodynamics of water, flow, velocity, . . .have not been considered here and may potentially impact the analysis).
Second, concerning restoration schedules, we try to explore the impacts of the timing of restoration.Assuming restoration evolves at a constant rate for 10 years from 2013, to meet a predetermined target at the end date, restoration at time t is obtained from: where t −1 , t 0 , and b are a year before the beginning of the restoration, the year the restoration begins, and a parameter determining the speed of restoration, respectively (see S3 for the derivation process).For example, when b = 0 (Figure 4a), restoration per year increases as time passes.When b = 5 (Figure 4c), restoration decreases linearly over time.In all cases, however, the same amount of nursery (i.e., Target restoration) is restored at the end of the 10-year period.
Sustainability 2018, 10, x FOR PEER REVIEW 9 of 21 sufficient for this type of analysis but hydrodynamics of water, flow, velocity, … have not been considered here and may potentially impact the analysis).
Second, concerning restoration schedules, we try to explore the impacts of the timing of restoration.Assuming restoration evolves at a constant rate for 10 years from 2013, to meet a predetermined target at the end date, restoration at time is obtained from: where , , and are a year before the beginning of the restoration, the year the restoration begins, and a parameter determining the speed of restoration, respectively (see S3 for the derivation process).For example, when 0 (Figure 4a), restoration per year increases as time passes.When 5 (Figure 4c), restoration decreases linearly over time.In all cases, however, the same amount of nursery (i.e., Target restoration) is restored at the end of the 10-year period.

Scenario Development
As suggested by the Millennium Ecosystem Assessment [44], we develop scenarios to explore the complex ecological-economic system, inform planning and decision-making processes, and help bridge the gap between scientific understanding and policy needs.
Our scenarios combine two approaches: policy impact assessment and policy sensitivity analysis.Policy impact assessments examine the economic, social, and environmental impacts of public policy (e.g., OECD [45] and the European Commission [46] and are frequently used in IO [14].Policy sensitivity analysis is one of the three sensitivity analyses proposed in SD [29].It tests if policy implications change when assumptions about the model (e.g., parameters with uncertainty) vary over a plausible range [29].
As illustrated in Table 2, we develop scenarios using two policy options, one uncertainty parameter, and five outcome indicators.

Scenario Development
As suggested by the Millennium Ecosystem Assessment [44], we develop scenarios to explore the complex ecological-economic system, inform planning and decision-making processes, and help bridge the gap between scientific understanding and policy needs.
Our scenarios combine two approaches: policy impact assessment and policy sensitivity analysis.Policy impact assessments examine the economic, social, and environmental impacts of public policy (e.g., OECD [45] and the European Commission [46] and are frequently used in IO [14].Policy sensitivity analysis is one of the three sensitivity analyses proposed in SD [29].It tests if policy implications change when assumptions about the model (e.g., parameters with uncertainty) vary over a plausible range [29].
As illustrated in Table 2, we develop scenarios using two policy options, one uncertainty parameter, and five outcome indicators.Since the timing of restoration matters, we test the three restoration timings proposed in Figure 4; that is, increasing, decreasing, and constant restoration over time.All three options restore the same surface (23.71513 km 2 ) over the same time frame (10 years), which allows the total surface area of nurseries to recover to its 1979-1980 level, as mentioned previously.In fact, it is possible to test infinite variations of timing by changing b ∈ [0, 5] in Equation ( 11), but for the purpose of this paper, we test the three timings "Increase (b = 0)," "Constant (b = 2.371513)," and "Decrease (b = 5)".
We consider two cost allocation rules across economic sectors and final demand categories.Technical details with alternative cost allocation rules are explained in [47].Technically, the two rules are reflected in the model using different α j values: the share of the total annual restoration cost paid by sector j in Equation ( 3); these are displayed in Figure 5.The first cost allocation rule-"No sharing"-follows the "polluter pays" principle developed by the OECD [48,49].According to this principle, the economic sector-harbors in our case study, which is included within the more general transport sector-directly responsible for environmental degradations pays to fix the situation.The second cost allocation rule-"Sharing"-follows the shared environmental responsibility principle developed by Gallego and Lenzen [50], Lenzen et al. [51], and Lenzen and Murray [52].According to this principle, a sector using products from a direct or indirect polluter in the supply chain should bear a share of its environmental responsibility.We calculate this share as a function of the number of commodities purchased by the polluter and of its ability to change production processes toward more ecological ones.Applying a shared environmental responsibility principle allows us to propose an alternative in which harbors do not bear the restoration cost alone.The aim is to alleviate the negative impacts on harbors' GOS and, hence, retain investment capacity and competitiveness.It is quite important given that harbors generate benefits for society, contribute to the public interest, and provide positive externalities relative to climate change mitigation.
We then conduct a policy sensitivity analysis to test if different water quality levels lead to different policy recommendations.We use the Latin hypercube method available in Powersim and did 40 runs with varying Water Quality Indexes (Equation ( 10)).Because we do not know the probability distribution of water quality, we simply adopt the uniform distribution with the range [0.50, 1.00], as described in Section 2.2.3.
Finally, we select five outcome indicators that represent the system's performance [9].The three economic indicators are GDP and disposable income-both aggregate economic values-as well as GOSs that represent sectoral impacts of restoration policies.The remaining two ecological indicators are soles caught, expressed in physical units, and nursery areas, expressed in physical and monetary units.TEV of the nursery area is estimated using the benefit transfer method [53].We adopt the value of the estuary in the TEEB valuation database [54].To avoid double counting of soles, the TEV includes flood prevention, material, recreation, spiritual, and cognitive services.The unit value is estimated at 85,749 Euro/km 2 /year.
The outcomes of the various scenarios have been computed with a zero discount rate, as varying the discount rate from 0 to 5% did not show meaningful insights.

Policy Impact Assessments
Figure 6a-d show the impacts of restoration schedules with cost allocation rule 1 on the outcome indicators over time except for GOSs, which will be shown later, separately.Results with cost allocation rule 2 are not presented because they are very similar, although not identical to those detailed here.In other words, the choice of the cost allocation rule does not have significant impact on these four outcome indicators.Economic outcome indicators (GDP and disposable income) are

Policy Impact Assessments
Figure 6a-d show the impacts of restoration schedules with cost allocation rule 1 on the outcome indicators over time except for GOSs, which will be shown later, separately.Results with cost allocation rule 2 are not presented because they are very similar, although not identical to those detailed here.In other words, the choice of the cost allocation rule does not have significant impact on these four outcome indicators.Economic outcome indicators (GDP and disposable income) are lower with restoration than at the BAU level-because of restoration costs-whereas the opposite holds for ecological outcome indicators (total surface of nursery areas and soles caught).Restoration schedule also plays a role, as increasing restoration rate over time (the "Increasing" scenario) is preferred for economic outcomes (Figure 6a,b), whereas faster restoration early in time (the "Decreasing" scenario) is preferred for ecological outcomes (Figure 6c,d).
Sustainability 2018, 10, x FOR PEER REVIEW 13 of 21 lower with restoration than at the BAU level-because of restoration costs-whereas the opposite holds for ecological outcome indicators (total surface of nursery areas and soles caught).Restoration schedule also plays a role, as increasing restoration rate over time (the "Increasing" scenario) is preferred for economic outcomes (Figure 6a,b), whereas faster restoration early in time (the "Decreasing" scenario) is preferred for ecological outcomes (Figures 6c,d).S1 for the cumulative values of the impacts of restoration over the simulation period displayed in Figure 6.
Figure 7 shows the impacts of restoration schedules on selected GOSs from 37 sectors to highlight the sectoral differences.Figure 7a-c suggest that, irrespective of the cost allocation rule, the total GOS earned by all companies in the region is lower with restoration than without restoration (BAU).S1 for the cumulative values of the impacts of restoration over the simulation period displayed in Figure 6.
Figure 7 shows the impacts of restoration schedules on selected GOSs from 37 sectors to highlight the sectoral differences.Figure 7a-c suggest that, irrespective of the cost allocation rule, the total GOS earned by all companies in the region is lower with restoration than without restoration (BAU).

Restoration Schedule Decreasing
Constant Increasing GOS: total  S2 for the cumulative values of the impacts displayed in Figure 7.   S2 for the cumulative values of the impacts displayed in Figure 7.
The cost allocation rule matters for the transport sector (Figure 7d-f) because it reduces the GOS, as harbors (which are included in the transport sector category) are bearing 100% and 49.75% of the cost in rules 1 and 2, respectively.Rule 2 ("Sharing") reduces the impact, as other sectors acting later in the supply chain share the restoration cost.For the other sectors and total GOS, the choice of the cost allocation rule does not make any visible difference so that lines for rules 1 and 2 overlap.
The construction sector benefits from the restoration because it implements it and thereby increases its GOS.However, negative impacts outweigh positive ones (Figure 7g-i).Indeed, during the restoration period, economic growth slows in the region, thereby reducing household consumption and inducing a slowdown for the construction sector through indirect linkages with the other sectors of the regional economy.
Figure 7j-l show a reduction in GOS for the textile sector across all three scenarios, in opposition with the aggregate figures (Figure 7a-c).This reflects the downward technical coefficients of the dynamized IO (Figure 2).This captures the past and current trend in the textile sector.French production has increasingly been substituted by imports from developing and emergent countries in which labor costs are much lower.Between 1994 and 2015, the output of the French textile industry decreased by 60% (own calculation in IO data from INSEE [55]).

Policy Sensitivity Analysis
We conducted a policy sensitivity analysis to explore how water quality level (WQ) influences policy recommendations.We only present the impact on soles caught because other outcome indicators were barely affected (see Supplementary Information for summary statistics) and there is no impact on nursery areas because their dynamics are exogenously determined by the restoration scenarios.
There are two important results.First, Figure 8 shows the impacts of water quality on soles caught when restoration is implemented.The results for High WQ and Low WQ correspond respectively to the largest and smallest values of soles caught computed by the sensitivity analysis.Since water quality improvement positively contributes to the abundance of juvenile soles (Equation (10) and Figure 3), it is reasonable to assume that better WQ leads to larger quantities of soles caught (Cons.Rest.and High WQ) than BAU WQ (Cons.Rest.and BAU WQ) throughout the simulation period.In total, the cumulated amount of soles caught is 8084 tons for high water quality, which is 23.9% higher than the BAU level (6527 tons).The cost allocation rule matters for the transport sector (Figure 7d-f) because it reduces the GOS, as harbors (which are included in the transport sector category) are bearing 100% and 49.75% of the cost in rules 1 and 2, respectively.Rule 2 ("Sharing") reduces the impact, as other sectors acting later in the supply chain share the restoration cost.For the other sectors and total GOS, the choice of the cost allocation rule does not make any visible difference so that lines for rules 1 and 2 overlap.
The construction sector benefits from the restoration because it implements it and thereby increases its GOS.However, negative impacts outweigh positive ones (Figure 7g-i).Indeed, during the restoration period, economic growth slows in the region, thereby reducing household consumption and inducing a slowdown for the construction sector through indirect linkages with the other sectors of the regional economy.
Figure 7j-l show a reduction in GOS for the textile sector across all three scenarios, in opposition with the aggregate figures (Figure 7a-c).This reflects the downward technical coefficients of the dynamized IO (Figure 2).This captures the past and current trend in the textile sector.French production has increasingly been substituted by imports from developing and emergent countries in which labor costs are much lower.Between 1994 and 2015, the output of the French textile industry decreased by 60% (own calculation in IO data from INSEE [55]).

Policy Sensitivity Analysis
We conducted a policy sensitivity analysis to explore how water quality level (WQ) influences policy recommendations.We only present the impact on soles caught because other outcome indicators were barely affected (see Supplementary Information for summary statistics) and there is no impact on nursery areas because their dynamics are exogenously determined by the restoration scenarios.
There are two important results.First, Figure 8 shows the impacts of water quality on soles caught when restoration is implemented.The results for High WQ and Low WQ correspond respectively to the largest and smallest values of soles caught computed by the sensitivity analysis.Since water quality improvement positively contributes to the abundance of juvenile soles (Equation (10) and Figure 3), it is reasonable to assume that better WQ leads to larger quantities of soles caught (Cons.Rest.and High WQ) than BAU WQ (Cons.Rest.and BAU WQ) throughout the simulation period.In total, the cumulated amount of soles caught is 8,084 tons for high water quality, which is 23.9% higher than the BAU level (6527 tons).While Figure 8 showed the combined impact of restoration and various WQ levels on soles caught, Figure 9 now compares the impact of restoration on soles caught with three scenarios ("Constant", "Decreasing", and "Increasing": Figure 4) and that of water quality improvement without restoration.The purpose is to explore how much water quality improvement contributes to the quantity of soles caught.As Figure 9 shows, if water quality greatly improves, without restoration, soles caught could be higher than with any restoration scenarios (e.g., High WQ).In cumulative values, although better water quality (High WQ) can improve soles caught by 23.3% (7332 tons) compared to BAU levels (5945 tons), the "Decreasing" restoration scenario (Dec.Rest.)improves the catch by 11.6%.There is, of course, a possibility that low water quality (e.g., Low WQ) leads to less catch (27.5% lower).
The best policy recommendation in terms of soles caught would therefore seem to hope for the best WQ level possible.However, as this seems unrealistic or at least difficult to achieve, an interesting policy option candidate (whose end result is close to high WQ level) would be restoration at a decreasing rate; that is, undertake more effort today and less in the future to arrive at the predetermined restoration target.While Figure 8 showed the combined impact of restoration and various WQ levels on soles caught, Figure 9 now compares the impact of restoration on soles caught with three scenarios ("Constant", "Decreasing", and "Increasing": Figure 4) and that of water quality improvement without restoration.The purpose is to explore how much water quality improvement contributes to the quantity of soles caught.As Figure 9 shows, if water quality greatly improves, without restoration, soles caught could be higher than with any restoration scenarios (e.g., High WQ).In cumulative values, although better water quality (High WQ) can improve soles caught by 23.3% (7332 tons) compared to BAU levels (5945 tons), the "Decreasing" restoration scenario (Dec.Rest.)improves the catch by 11.6%.There is, of course, a possibility that low water quality (e.g., Low WQ) leads to less catch (27.5% lower).
The best policy recommendation in terms of soles caught would therefore seem to hope for the best WQ level possible.However, as this seems unrealistic or at least difficult to achieve, an interesting policy option candidate (whose end result is close to high WQ level) would be restoration at a decreasing rate; that is, undertake more effort today and less in the future to arrive at the predetermined restoration target.

Model Development and Analysis
We developed a fully dynamic IO/SD ecological-economic model by extending [16] with dynamization of the technical coefficients in the IO component so the model captures dynamic relationships within and between the ecological and economic systems.
Our fully dynamic model is an answer to criticisms about constant technical coefficients, showing that changes in technical coefficients differ by sector and that their magnitude [37] and directions of change (positive or negative) seem to be reasonable.The dynamization of IO shows a sharp contrast within various sectors when restoration policies are implemented (Figure 7).
We used econometrics, more precisely GLS, to evolve from a static to dynamic IO (see Equation ( 8)).Other techniques could have been investigated for capturing changes that are not an extension of the past, including using experts' opinions on future technological changes [56] and the RAS method [57], for example.

Model Development and Analysis
We developed a fully dynamic IO/SD ecological-economic model by extending [16] with dynamization of the technical coefficients in the IO component so the model captures dynamic relationships within and between the ecological and economic systems.
Our fully dynamic model is an answer to criticisms about constant technical coefficients, showing that changes in technical coefficients differ by sector and that their magnitude [37] and directions of change (positive or negative) seem to be reasonable.The dynamization of IO shows a sharp contrast within various sectors when restoration policies are implemented (Figure 7).
We used econometrics, more precisely GLS, to evolve from a static to dynamic IO (see Equation ( 8)).Other techniques could have been investigated for capturing changes that are not an extension of the past, including using experts' opinions on future technological changes [56] and the RAS method [57], for example.
We then analyzed the IO/SD model by developing various scenarios that explore the complex behavior of the ecological-economic system as applied in the Millennium Ecosystem Assessment [44].By taking advantage of the SD software, we combined policy impact and policy sensitivity analyses to develop various scenarios.Whether assessment indicators are barely or significantly impacted by restoration and water quality improvement, all cases reflect dynamic interactions between and within the ecological-economic model.Hence, both findings seem relevant.
Varying restoration schedules revealed trade-offs between economic impacts and ecological benefits.Also, the sectoral analyses undertaken in this paper showed that the overall negative impacts of the economic slowdown due to restoration costs outweigh the benefits for the regional economy, mostly materialized by advantages for the construction sector.
The policy sensitivity analysis using water quality improvement identified the potential importance of examining exogenous factors outside the system boundary.Water quality is exogenously given (it is beyond the regional authorities' control) and our analyses show that its improvement could increase soles caught by a greater amount than any restoration policy.Hence, water quality should be a priority for research and public policy aimed at restoring fish nursery habitats.Obviously, improving water quality will result in economic costs for the region.Also, because the Seine River's water quality depends on neighboring river bodies [43], coordination between regions and countries seems essential and economic costs may be shared.To quantify the economic and ecological consequences of water quality improvement resulting from changes in industrial processes, agricultural practices and urban waste water treatment plants, a full hydrological and biogeochemical model of the Seine-Normandie water basin should be developed for many contaminants (HAPs, Nitrates, Phosphates, organic matter, heavy metals, pesticides, residual medical drugs, plastic chemicals such as Bisphenol-A, etc.).However, this is a real scientific challenge given the huge surface area of the Seine-Normandie (95,000 km 2 ) extending from the Seine estuary, downstream, to Paris city, upstream, and including almost 200,000 industries and 18 million inhabitants.However, developing such a model is beyond current scientific capacities and should be developed in future research [58].

Future Research Topics
Despite the potential importance of IO/SD models [10][11][12], the development of such models has been rare.Our study paves the way for a new generation of ecological-economic models based on IO/SD, leaving room for future research on modeling and environmental policy analysis.
First, research on sustainability indicators and ecological-economic modeling should be undertaken together [10,59].Our model captured the correlations and trade-offs between sustainability indicators, but research on sustainability indicators [60] often does not reflect the relationships between them, ignoring their interdependencies.In turn, indicator selection can guide ecological-economic modeling [61], as a model should be built for a specific purpose [29], such as capturing the dynamics and relationships between sustainability indicators.Our model was not systematically guided by indicator selection.
Second, improving the degree of system closure in economic components by adding an econometric analysis could further improve modeling.Indeed, keeping final demands exogenous, such as in most IO models, may undervalue socio-economic changes [36], but regional IO-econometric models could improve inter-industry impact studies [13,[36][37][38].In our model, two final demand categories are endogenous: household consumption (computed as a function of incomes) and investments (computed as a function of profit).However, both are partly based on an exogenous parameter: household consumption includes exogenous income elasticities (e i ) and investments are based on a fixed capital formation coefficient ( ĉap i ) calculated in the IO table at the reference year.Future research should replace exogenous parameters with statistical regressions computed from time series data related to the studied area.

Figure 1 .
Figure 1.Interconnections between the economic sub-system and ecological system.The figure is adapted from[16] with primary model extensions underlined in red.Note: "+" and "−" indicate variable changes in the "same" and "opposite" direction, respectively.

Figure 1 .
Figure 1.Interconnections between the economic sub-system and ecological system.The figure is adapted from[16] with primary model extensions underlined in red.Note: "+" and "−" indicate variable changes in the "same" and "opposite" direction, respectively.

Figure 3
Figure3depicts the SD part of the IO/SD model in a stock and flow diagram.The model has two stocks: "Nursery areas" and "Sole stock from the internal part of the Seine."The nursery areas include 21 categories with different sole abundance to represent the spatial heterogeneities; the categorization is based on the sediment type-gravel, sand, or silt-and the depth.These areas are assumed to be independent, as we do not know how they interact with each other.Another simplification of the ecological model stands with the fact that we do not develop the physical, chemical, and biological conditions required by soles to reproduce, e.g., the spawn is performed in the bottoms where soles inhabit and therefore depends on environmental conditions such as water temperature and acidity, water quality and food availability[42].The sole stock uses a cohort structure, ages 1 through 10.

Figure 2 .
Figure 2. Changes in technical coefficients for the 37 sectors (positive changes are in dark color, negative ones in lighter color, and white cells mean no change).See S1 for regression models used for creating the figure.See S2 for the full descriptions of sector names.

Sustainability 2018 , 21 Figure 3 .
Figure 3.The IO/SD model of the ecological system and its economic sub-system (Stocks and flows are represented by boxes and double arrows respectively.Circles and diamonds denote auxiliary variables and constants).

Figure 3 .
Figure 3.The IO/SD model of the ecological system and its economic sub-system (Stocks and flows are represented by boxes and double arrows respectively.Circles and diamonds denote auxiliary variables and constants).

Figure 6 .
Figure 6.Impacts of restoration schedules with cost allocation rule 1. See TableS1for the cumulative values of the impacts of restoration over the simulation period displayed in Figure6.

Figure 6 .
Figure 6.Impacts of restoration schedules with cost allocation rule 1. See TableS1for the cumulative values of the impacts of restoration over the simulation period displayed in Figure6.

Figure 7 .
Figure 7. Impact of restoration schedule and cost allocation rule on selected GOSs.See TableS2for the cumulative values of the impacts displayed in Figure7.

Figure 7 .
Figure 7. Impact of restoration schedule and cost allocation rule on selected GOSs.See TableS2for the cumulative values of the impacts displayed in Figure7.

Figure 9 .
Figure 9. Soles caught to compare the contributions of restoration vs. water quality.High WQ and Low WQ are computed by the sensitivity analysis with varying WQ index.See Tables S6-S8 for more detailed results of the sensitivity analysis.

Figure 8 .
Soles caught.Cons.Rest.stands for restoration at constant rate.WQ stands for water quality.Cons.Rest.and High WQ and Cons.Rest.and Low WQ are computed by the sensitivity analysis with varying WQ index.Cost allocation rule 1 is applied.See Tables S3-S5 for more detailed results of the sensitivity analysis.
Figure 8. Soles caught.Cons.Rest.stands for restoration at constant rate.WQ stands for water quality.Cons.Rest.and High WQ and Cons.Rest.and Low WQ are computed by the sensitivity analysis with varying WQ index.Cost allocation rule 1 is applied.See Tables S3-S5 for more detailed results of the sensitivity analysis.
Soles caught to compare the contributions of restoration vs. water quality.High WQ and Low WQ are computed by the sensitivity analysis with varying WQ index.See Tables S6-S8 for more detailed results of the sensitivity analysis.