Advanced Modelling Tools to Support Planning for Sand/Gravel Quarries

: Sand and gravel quarry planning must guarantee the public interest in the procurement of raw materials while ensuring environmental sustainability. An Analyzing Planning Support System for sand and gravel quarry plan can assist decision-makers during the planning process. The proposed Analyzing PSS uses integrating geologic, economic, environmental, and geographic information to quantify raw materials and the size of quarries. This kind of tool is useful to support public authority decisions. The study provides the results of an experience conducted in the province of Brescia (NUT 3 in Northern Italy).


Introduction
Sand and gravel quarrying is needed to obtain natural raw materials mainly used by the construction sector, and it is a strategic activity for the economy. This activity is conditioned both by the presence of raw materials in the subsoil and by sites' accessibility and suffers market pressures caused by the demand for the material. Quarrying causes strong environmental impacts and must be limited, mitigated, and compensated as much as possible. Natural aggregate materials (i.e., sand and gravel) are considered "the most consumed bulk product in the world after water" [1] and, although the use of recycled aggregate is promising [2], the worldwide primary sand/gravel production is estimated to grow in the next decades [3]. Worldwide sand and gravel excavation is not monitored, and the use of aggregates is strictly correlated to cement production [4]. In the European Union there was a decline of raw material needs due to the economic crisis of the construction sector in 2008. The raw material needs increased again from 2013 to 2018 when the total production of aggregates was estimated at 2.894 million tons, of which 1.135 million tons was sand and gravel [5].
Administration authorities are called upon to govern the aggregate quarrying process, considering its multiple (and sometimes conflicting) dimensions, such as the economic and environmental ones. In civil society, stakeholders hold contradictory theses, as there is a public interest both in the procurement of raw materials and in environmental protection. Sand and gravel quarry planning tackles all the dimensions of the quarrying process, building an integrated system of information and (geo)data, developing specific strategies, and establishing rules. Although authorities have to make complex decisions during the quarry planning process, there is a lack of objective tools to turn the planning goals into possible scenarios related to the strategy.
Considering the growing pressure put on the mining industry facing environmental challenges [6], literature analysis shows how research on quarry planning is mainly dealing with environmental problems. Some studies concentrate on how to assess the impacts associated with the extraction and processing of raw materials so as to prevent, mitigate,

Materials and Methods
The methodological approach of the paper (Figure 1) relies on a four-step linear analysis. The first and fundamental phase is represented by the input data selection to implement the PSS. The input data must be site-specific, and they must be chosen in accordance with the public authority and with the stakeholders involved in the quarry planning process.
poses a research study focus on performance of PSS and highlights research gap importance of PSS taxonomy (Informing, Communicating, and Analyzing PSS) thermore, in the 21st century, the complexity of PSS traditional technology has b come by GIS-based Planning Support technologies [33], and they are actually use tical planning to inform and improve policies and decisions [34].
This study proposes an Analyzing PSS (GIS-based Planning Support techno ibrated for sand and gravel quarry planning, providing an advanced modellin support public decision-makers during the planning process. The study will foc swering the following research questions: Is it possible to assist public authority ing the proper amount of sand and gravel to be extracted and the dimension quarry through a GIS-based Planning Support technology?

Materials and Methods
The methodological approach of the paper (Figure 1) relies on a four-step lin ysis. The first and fundamental phase is represented by the input data selection ment the PSS. The input data must be site-specific, and they must be chosen in ac with the public authority and with the stakeholders involved in the quarry plan cess.
Then, the Analyzing PSS can be set in order to state both the proper amoun and gravel to be extracted and the dimension of each quarry. The PPS is applied case study, and the results are presented. In the paper, finally, the discussion re potential and the limitation of the MAA proposed, suggesting recommendation authorities (decision-makers) involved in the redaction of a sand and gravel qua

Input Data
In order to implement an Analyzing PSS and calculate the volume of raw m be extracted from each quarry, the following data are needed: 1. Estimation of the overall sand and gravel demand. The value is expressed meters (m 3 ). 2. Monitoring database of the previous quarry plan. This includes data on the of raw material to be extracted from each quarry according to the old plan effectively extracted volumes (m 3 ), and the quarries' perimeters. 3. Territorial database in GIS. This includes geographical, geological, infras environmental, and landscape data. The data are acquired from geoportals authorities who are competent about the territory, or data can be acquire field using on-site inspections. The database also includes existing quarries ters. 4. Data about the mining companies. This includes data about environmenta certificates, the number of employees, and the availability of plants for the of aggregates in support of the circular economy. 5. Information and data on the mining companies' business development p data can be collected through a survey promoted by the public authority t prepare the new quarry plan. This information is useful in the decision-ma cess. Then, the Analyzing PSS can be set in order to state both the proper amount of sand and gravel to be extracted and the dimension of each quarry. The PPS is applied in a real case study, and the results are presented. In the paper, finally, the discussion reveals the potential and the limitation of the MAA proposed, suggesting recommendations for the authorities (decision-makers) involved in the redaction of a sand and gravel quarry plan.

Input Data
In order to implement an Analyzing PSS and calculate the volume of raw material to be extracted from each quarry, the following data are needed:

1.
Estimation of the overall sand and gravel demand. The value is expressed in cubic meters (m 3 ).

2.
Monitoring database of the previous quarry plan. This includes data on the volumes of raw material to be extracted from each quarry according to the old plan, data on effectively extracted volumes (m 3 ), and the quarries' perimeters.

3.
Territorial database in GIS. This includes geographical, geological, infrastructural, environmental, and landscape data. The data are acquired from geoportals of public authorities who are competent about the territory, or data can be acquired on the field using on-site inspections. The database also includes existing quarries' perimeters.

4.
Data about the mining companies. This includes data about environmental quality certificates, the number of employees, and the availability of plants for the recovery of aggregates in support of the circular economy.

5.
Information and data on the mining companies' business development plans. The data can be collected through a survey promoted by the public authority that must prepare the new quarry plan. This information is useful in the decision-making process.

Method
The goal of the Analyzing PSS is to define the volume of raw material to be extracted from each quarry so as to satisfy the overall sand and gravel demand established in an update of a quarry plan. Overall sand and gravel demand estimation is assumed as input data and is not a goal of the PSS. The PSS proposed and the methodology are applicable in a territory where an upgrade of the sand and gravel quarry plan is needed and is based on five main steps ( Figure 2):

1.
Exclusion of existing non-exploitable quarries

Method
The goal of the Analyzing PSS is to define the volume of raw material to be extracted from each quarry so as to satisfy the overall sand and gravel demand established in an update of a quarry plan. Overall sand and gravel demand estimation is assumed as input data and is not a goal of the PSS. The PSS proposed and the methodology are applicable in a territory where an upgrade of the sand and gravel quarry plan is needed and is based on five main steps ( Figure 2  Step 1. An assessment of the potential of existing quarries is carried out. If a quarry has been depleted or if for any reason it is no longer possible to go on extracting raw material, the quarry is excluded from the new quarry plan. If conversely, the deposit still is exploitable and there are no critical issues, the quarry is confirmed in the new quarry plan. The assessment is based on the data present in the monitoring database of the pre- Step 1. An assessment of the potential of existing quarries is carried out. If a quarry has been depleted or if for any reason it is no longer possible to go on extracting raw Sustainability 2021, 13, 6380 5 of 14 material, the quarry is excluded from the new quarry plan. If conversely, the deposit still is exploitable and there are no critical issues, the quarry is confirmed in the new quarry plan. The assessment is based on the data present in the monitoring database of the previous quarry plan and in the mining companies' survey results. Step 2. A volume quota called "Base Quota" (BQ) is assigned to each quarry as a fixed percentage of an "Intermediate Quota" (IQ). The IQ is proportional to the volumes assigned to the quarries by the previous quarry plan. Proportionality is calculated both at the single quarry level and also at the Local Quarries Network (LQN) level.
LQNs are basins, spatially defined, which aggregates quarries according to the relational fields and socio-economic, geographic, historical, and culturally homogeneous features. In quotas allocation, proportionality with the previous plan is maintained to respect former companies' investments, and LQNs are used to maintain the territorial balance. Therefore, the IQ is calculated by the equation below, Equation (1).
where D is the sand and gravel demand; Vpq is the previous quarry plan volume assigned to the single quarry; Vcq is the volume assigned by the previous quarry plan in the quarries confirmed by the new plan; and the summation i-th is the overall value for all the quarries or the overall value for the quarries included in a single LQN. Furthermore, the quotas are carried out using a Multi-criteria Analysis Approach (MAA), which requires the input data of the estimated sand and gravel demand data, the monitoring database of the previous quarry plan data, and the territorial database in GIS.
As already explained before, the "Base Quota" (BQ) is a fixed percentage of IQ. The remaining part of IQ constitutes "quota s" (Qs), which is the input data for the estimation of the "Premium Quota" (PQ) in the following step 3. Therefore, the relation among quotas can be synthesized as follows (Equation (2)): The relation between QB e Qs is strictly related to the importance of the "premium" given by the decision-maker, so an example of the assigned weight at "base" and "premium" quotas is presented in the case study.
Step 3. A volume quota, called "Premium Quota" (PQ), is assigned to each quarry to balance territorial/environmental constraints with mining business sustainability. The calculation of the PQ is carried out using a Multi-criterial Analysis Approach (MAA). As input data, it requires the estimation from the overall sand and gravel demand data, the monitoring database of the previous quarry plan, the territorial database in GIS, and the data about the mining companies, as well as "quota s", which is the part of the "Intermediate Quota" (IQ) not previously assigned. A score based on territorial/environmental aspects and the mining companies' features becomes the base for assigning PQ and given to each quarry.
Such criteria are described below: • Territorial/environmental features of where the quarry is. Aspects such as environmental impairment, interference with the ecological network, and interferences with floodable areas due to flooding institutions are assessed. The score is assigned according to a weighted linear function between 0 and 2, taking into account the score given to each feature assessed (high critical issues = 0, medium critical issues = 1, low critical issue = 2). • Internal quarry features. Aspects such as the perception of cultivation activity, the presence of landfills within the quarry area, the presence of degraded situations, the perception of the recoveries and prescribed mitigations, and the accessibility of the quarry are assessed. The score is assigned according to a weighted linear function between 0 and 2, taking into account the score given to each feature assessed (high critical issues = 0, medium critical issues = 1, low critical issue = 2). • Mining companies' possession of environmental (ISO 14001), quality (ISO 9001), and safety in the workplace (ISO 45001) certifications. The score is assigned according to Boolean logic (no certificate = 0, possession of certificate = 1).

•
Mining companies' employer number. The score is assigned according to a linear function between 0 (the company with the lowest employee number) and 1 (the company with the highest employee number).

•
Presence at the quarry of reusable inert waste plants (to support circular economy). The score is assigned according to Boolean logic (no reusable inert waste plant = 0, presence of plant = 1).
To assign the "Premium Quota" (PQ) to each quarry, each quarry's score (P) and the sum of the score assigned to all the quarries must be known. To assign a quota that is proportional to the size of the quarry, it is necessary to maintain a proper relation between the assigned base quota (BQ) and the premium one (PQ) to be assigned. The equation below has been used.
where PQ is the "Premium Quota"; P is the score obtained by the quarry for the territorial/environmental aspects and the mining quota features criteria, Qs is "quota s" coming from the previous step; and the summation i-th is the overall value for all the quarries. With equal "Base Quota", a hypothetical quarry to which a high score is attributed will obtain a higher PQ compared to another quarry to which a lower score is attributed. Adding the "Base Quota" (defined in the previous step) to the "Premium Quota", one obtains the volume quota to be extracted from each quarry, called "Maximum Quota" (MQ). The feasibility of the excavation of MQ for each quarry is then carried out. Any limits to the expansion or deepening of the quarries are also checked. The data for this assessment are selected from the territorial database in GIS and in the results of the survey on the mining companies. Finally, the maximum excavable volume is determined for each quarry. If there are any MQ that exceed the maximum excavable volume, it is necessary to proceed as foreseen in the next step.
In the fourth step, the way to reassign those MQ that cannot be extracted from the deposits of existing quarries is defined. The first possibility is to reassign those quotas (volumes) to the other quarries that belong to the same LQN using a criterion of proportionality, respecting the "maximum quota" already assigned. Those quarries that can receive those additional volumes (quotas) must have a deposit with a quantity of raw material higher than MQ. A second possibility is to reassign the volumes to other quarries belonging to the same mining company. This second chance may be evaluated by using the mining companies' survey results. In the event that some quarries still have excess volumes, it is necessary to go back to point 3.2 of the method and proceed with the reassignment of excess volumes again.
The reassignment of volumes is based on the data obtained from the previous steps, the data obtained from the mining companies survey, and the data from the territorial database in GIS where the LQN are spatially defined.
The fifth and last step definitively assigns the proper amount of sand and gravel to be extracted in each quarry, called "Final Quota" (FQ). The assignment of "Final Quota" is based on the data obtained from the previous steps of the method. Each step of the method needs the input of the data selected and collected in the preliminary step of the study. The input data needed in each step of the method are represented in Table 1. Step 2 Step 3 Step 4 Step 5 All steps, but the first, use the results from previous steps as input data.
Sustainability 2021, 13, x FOR PEER REVIEW 7 of 14 Step 2 Step 3 Step 4 Step 5 All steps, but the first, use the results from previous steps as input data.  Data is required;  Data is optional;  Data is not required; Input data 1: estimates of overall sand and gravel demand. Input data 2: database of the previous quarry plan monitoring. Input data 3: mining companies survey results. Input data 4: territorial database in GIS. Input data 5: mining companies' data.
Step 3: Assignment of sand/gravel premium volumes and overall check.

Presentation of the Case Study
The province of Brescia is an administrative district of the Lombardy Region, located in northern Italy, between the Alps and the Po river plain. According to the Statistical regions of the Eurostat Union [35], it is classified as "NUTS level 3", and, from an international point of view, it is considered as the eastern part of the metropolitan area of Milan, which is the most developed in Italy (Figure 3). Statistics show that Brescia is one of the most prosperous provinces in Europe and Italy [36]: it has more than 1.2 million inhabitants, and it is among the most populous provinces in Italy (5th out of 110, Eurostat 2016), with a Gross Domestic Product of 40,806.86 million euro (5th out of 111, Eurostat 2016), which however have a bad performance for air pollution [37] and soil consumption phenomena [38]. A vast deposit of sand and gravel is present in the plain. The raw material is first-rate. Lombardy is the Italian Region in which the highest amount of sand and gravel is extracted (8,350,000 cubic meters in 2018, 27% of the national total amount [39]), and Brescia is the Lombard province with the highest number of active quarries [40]. Extraction activities represent 3% of provincial GDP. In 2017, the local authority of the Province of Brescia decided to draw up the new Plan for sand and gravel quarries (new PPC), as well as a Strategic Environmental Assessment (SEA). According to Lombardy Regional Law no. 14 of 1998, the plan established the location, quality, and quantity of raw materials as well as the rules for quarrying activities. Step 2 Step 3 Step 4 Step 5 All steps, but the first, use the results from previous steps as input data.  Data is required;  Data is optional;  Data is not required; Input data 1: estimates of overall sand and gravel demand. Input data 2: database of the previous quarry plan monitoring. Input data 3: mining companies survey results. Input data 4: territorial database in GIS. Input data 5: mining companies' data.
Step 3: Assignment of sand/gravel premium volumes and overall check.

Presentation of the Case Study
The province of Brescia is an administrative district of the Lombardy Region, located in northern Italy, between the Alps and the Po river plain. According to the Statistical regions of the Eurostat Union [35], it is classified as "NUTS level 3", and, from an international point of view, it is considered as the eastern part of the metropolitan area of Milan, which is the most developed in Italy (Figure 3). Statistics show that Brescia is one of the most prosperous provinces in Europe and Italy [36]: it has more than 1.2 million inhabitants, and it is among the most populous provinces in Italy (5th out of 110, Eurostat 2016), with a Gross Domestic Product of 40,806.86 million euro (5th out of 111, Eurostat 2016), which however have a bad performance for air pollution [37] and soil consumption phenomena [38]. A vast deposit of sand and gravel is present in the plain. The raw material is first-rate. Lombardy is the Italian Region in which the highest amount of sand and gravel is extracted (8,350,000 cubic meters in 2018, 27% of the national total amount [39]), and Brescia is the Lombard province with the highest number of active quarries [40]. Extraction activities represent 3% of provincial GDP. In 2017, the local authority of the Province of Brescia decided to draw up the new Plan for sand and gravel quarries (new PPC), as well as a Strategic Environmental Assessment (SEA). According to Lombardy Regional Law no. 14 of 1998, the plan established the location, quality, and quantity of raw materials as well as the rules for quarrying activities. Step 2 Step 3 Step 4 Step 5 All steps, but the first, use the results from previous steps as input data.  Data is required;  Data is optional;  Data is not required; Input data 1: estimates of overall sand and gravel demand. Input data 2: database of the previous quarry plan monitoring. Input data 3: mining companies survey results. Input data 4: territorial database in GIS. Input data 5: mining companies' data.
Step 3: Assignment of sand/gravel premium volumes and overall check.

Presentation of the Case Study
The province of Brescia is an administrative district of the Lombardy Region, located in northern Italy, between the Alps and the Po river plain. According to the Statistical regions of the Eurostat Union [35], it is classified as "NUTS level 3", and, from an international point of view, it is considered as the eastern part of the metropolitan area of Milan, which is the most developed in Italy (Figure 3). Statistics show that Brescia is one of the most prosperous provinces in Europe and Italy [36]: it has more than 1.2 million inhabitants, and it is among the most populous provinces in Italy (5th out of 110, Eurostat 2016), with a Gross Domestic Product of 40,806.86 million euro (5th out of 111, Eurostat 2016), which however have a bad performance for air pollution [37] and soil consumption phenomena [38]. A vast deposit of sand and gravel is present in the plain. The raw material is first-rate. Lombardy is the Italian Region in which the highest amount of sand and gravel is extracted (8,350,000 cubic meters in 2018, 27% of the national total amount [39]), and Brescia is the Lombard province with the highest number of active quarries [40]. Extraction activities represent 3% of provincial GDP. In 2017, the local authority of the Province of Brescia decided to draw up the new Plan for sand and gravel quarries (new PPC), as well as a Strategic Environmental Assessment (SEA). According to Lombardy Regional Law no. 14 of 1998, the plan established the location, quality, and quantity of raw materials as well as the rules for quarrying activities.
Data is not required; Input data 1: estimates of overall sand and gravel demand. Input data 2: database of the previous quarry plan monitoring. Input data 3: mining companies survey results. Input data 4: territorial database in GIS. Input data 5: mining companies' data.
Step 3: Assignment of sand/gravel premium volumes and overall check.

Presentation of the Case Study
The province of Brescia is an administrative district of the Lombardy Region, located in northern Italy, between the Alps and the Po river plain. According to the Statistical regions of the Eurostat Union [35], it is classified as "NUTS level 3", and, from an international point of view, it is considered as the eastern part of the metropolitan area of Milan, which is the most developed in Italy (Figure 3). Statistics show that Brescia is one of the most prosperous provinces in Europe and Italy [36]: it has more than 1.2 million inhabitants, and it is among the most populous provinces in Italy (5th out of 110, Eurostat 2016), with a Gross Domestic Product of 40,806.86 million euro (5th out of 111, Eurostat 2016), which however have a bad performance for air pollution [37] and soil consumption phenomena [38]. A vast deposit of sand and gravel is present in the plain. The raw material is first-rate. Lombardy is the Italian Region in which the highest amount of sand and gravel is extracted (8,350,000 cubic meters in 2018, 27% of the national total amount [39]), and Brescia is the Lombard province with the highest number of active quarries [40]. Extraction activities represent 3% of provincial GDP.
Step 2 Step 3 Step 4 Step 5 All steps, but the first, use the results from previous steps as input data.  Data is required;  Data is optional;  Data is not required; Input data 1: estimates of overall sand and gravel demand. Input data 2: database of the previous quarry plan monitoring. Input data 3: mining companies survey results. Input data 4: territorial database in GIS. Input data 5: mining companies' data.
Step 3: Assignment of sand/gravel premium volumes and overall check.

Presentation of the Case Study
The province of Brescia is an administrative district of the Lombardy Region, located in northern Italy, between the Alps and the Po river plain. According to the Statistical regions of the Eurostat Union [35], it is classified as "NUTS level 3", and, from an international point of view, it is considered as the eastern part of the metropolitan area of Milan, which is the most developed in Italy (Figure 3). Statistics show that Brescia is one of the most prosperous provinces in Europe and Italy [36]: it has more than 1.2 million inhabitants, and it is among the most populous provinces in Italy (5th out of 110, Eurostat 2016), with a Gross Domestic Product of 40,806.86 million euro (5th out of 111, Eurostat 2016), which however have a bad performance for air pollution [37] and soil consumption phenomena [38]. A vast deposit of sand and gravel is present in the plain. The raw material is first-rate. Lombardy is the Italian Region in which the highest amount of sand and gravel is extracted (8,350,000 cubic meters in 2018, 27% of the national total amount [39]), and Brescia is the Lombard province with the highest number of active quarries [40]. Extraction activities represent 3% of provincial GDP. In 2017, the local authority of the Province of Brescia decided to draw up the new Plan for sand and gravel quarries (new PPC), as well as a Strategic Environmental Assessment (SEA). According to Lombardy Regional Law no. 14 of 1998, the plan established the  In 2017, the local authority of the Province of Brescia decided to draw up the new Plan for sand and gravel quarries (new PPC), as well as a Strategic Environmental Assessment (SEA). According to Lombardy Regional Law no. 14 of 1998, the plan established the location, quality, and quantity of raw materials as well as the rules for quarrying activities. The plan, valid for ten years, was proposed by provincial authorities and approved by the regional council.
The study case is emblematic in that it concerns one of the Italian provinces that has a very high number of active quarries, a territory with first-rate sand and gravel, and a local public administration that has decided to draw up a new sand and gravel quarry plan.

Towards the New Quarry Plan of the Province of Brescia
The Analyzing PSS has been used to define the volume of raw material to be excavated from each quarry so as to satisfy the overall sand and gravel demand calculated for drawing up the new quarry plan of the Province of Brescia. The input data for the case study are: Mining companies' data. The database contains data from specialized databases (AIDA Van Dijk, Accredia, consulted in 2020) and information about waste recovery plants (source Province of Brescia).

5.
Mining companies' survey results, promoted by the public authority within the Strategic Environmental Assessment of the plan. Seventy-seven responses to the survey were collected.
The overall demand for sand and gravel was estimated using a model provided by the regional deliberation no. 8/11347 of 2010. The application of the case study of the five steps of the Analyzing PSS will be reported in the subparagraphs below.

Exclusion of Existing No Longer Exploitable Quarries
In the case of Brescia, the potentiality of the 53 existing quarries, present in the PPC 2005 monitoring database, was assessed. The assessment was also carried out using the results of the mining companies survey. It emerges that 11 quarries out of 53 cannot go on providing sand and gravel. The reason that led to the extension of the quarries are (they may become one or more than one):

•
Depletion of the deposit and start of the quarry recovery (nine quarries).

•
The mining company could no longer go on extracting raw material from the quarry (four quarries).

•
The raw materials present in the deposit are of poor quality (one quarry). • Therefore, the confirmed quarries are 42.

Assignment of Sand/Gravel Base Volume (Pre-Dimensioning)
In the case of Brescia, the BQ is 80% of IQ (36,982,226 m 3 ). This percentage was fixed by the decision-maker to gradually introduce the "premium" concept (that consequently is 20% of IQ equal to 9,245,557 m 3 ).
The "Intermediate Quota" (IQ) proportional to the quotas assigned in PPC 2005 was calculated for each quarry. Six LQNs were identified. LQNs are expressions of the relational fields, socio-economic, geographic, historical, and culturally homogeneous features, and they coincide with the homogeneous territorial areas already defined by Lombardy Region in the Regional Territorial Plan. To maintain territorial balance and respect companies' investments, the sum of the "Base Quotas" assigned to each LQN by Sustainability 2021, 13, 6380 9 of 14 PPC 2005 is proportional to the sum of each "Base Quota" assigned to each LQN by the new quarry plan. The shutting down of quarries in a LQN is compensated by an increased production from the remaining quarries. In the new quarry plan, each LQN has a variable number of quarries and of assigned volumes (Table 2) based also on territorial features, the presence of sand and gravel, and the accessibility level.

Assignment of Sand/Gravel Premium Volumes and Overall Check
In the case of Brescia, "Premium Quotas" (PQ) were calculated both on a reward score and on "s quotas" (calculated in the previous step).
For what concerns the premium gain, a score from 0 to 7 based on five criteria was assigned. For the two territorial criteria, scores from 0 to 4 were assigned, while for the three company criteria, scores from 0 to 3 were assigned. In the case in which more than one company operates in a quarry, the score was calculated for each mining company and it was then weighed on the amount of volumes foreseen by the authorization issued to each company.
The choice was to give greater importance to territorial features than to companies' ones.
In the case of Brescia, a total of 137.69 scores were assigned. For what concern "s quotas", they were calculated to maintain a ratio of proportionality with the quarries dimension (indirectly through the proportionality with the "Base Quota"). The sum of "s quotas" was 9,245,557 m 3 , which is 20% of the sand and gravel demand estimate. To apply Equation (3) to the case study, it was also necessary to calculate the value of the sum of the results deriving from the multiplication between "s quotas" values and score values assigned to each quarry. This value, equal to 30,503,782, was used as the denominator to maintain a correct relationship between the BQ assigned and the PQ to be assigned. To assign the PQ to each quarry, the following equation was applied.

PQ =
Qs i × P i × 9, 245, 557 m 3 30, 503, 782 The "Premium Quotas" distributed in LQN are shown in Table 3. Adding the "Base Quota" (defined in the previous step) to the "Premium Quota", we have the volume quota to be extracted from each quarry, called "Maximum Quota" (MQ). The feasibility of MQ was then checked for each quarry.
In the case study, a series of criticalities regarding the constraints of enlargement and deepening of the existing quarries were detected following the surveys and supported by existing territorial database evidence in GIS. Such constraints are both to the nature of the deposit and to the coexistence of environmental pressures caused by other activities with an impact on the territory, particularly by landfills and airport activities (Montichiari airport development Plan). Therefore, after the overall check (point 3.4 of the method), n • 4 quarries were not confirmed, and 12 others were deemed not expandable in width and depth.

Reassignment of Non-Assignable Volumes
In the study case, a share equal to 4,664,865 m 3 of the volume (~10% of the total) could not be extracted in the assigned quarries. Such volumes were reassigned to other quarries belonging to the same LQN in the following cases:
In the case study, it was not possible to extract all the assigned quotas from LQN A quarries owing to territorial and morphological conditions (Piedmont area, river valley). The exceeding volumes, equal to 65,909 m 3 , were redistributed to quarries belonging to all the LQNs.
The results of the mining companies' surveys highlighted the opportunity for seven quarries to transfer 4,780,968 m 3 (~10% of the total) to quarries owned by the same company.

Final Volume Assignment
The FQ distributed in 38 quarries are shown in Table 4. Table 4. Sum of quarries' "Final Quota" (FQ) for each LQN.

LQN n • Quarries Volumes (m 3 )
A "Camonica Valley" 2 620,000 B "Trompia Valley and Sabbia Valley" 3 1,370,000 C "Sebino Franciacorta" 3 3,320,000 D "Gardesan Riviera" 3 1,739,000 E "Hills and high Brescia plain" 20 32,375,000 F "Low Brescia plain" 7 6,790,000 Figure 4 shows the territorial features of the LQNs (at the top) and a scheme of the Analyzing PSS applied to the case study of Brescia (at the bottom). Sustainability 2021, 13, x FOR PEER REVIEW 11 of 14

Discussion
The proposed Analyzing PSS is usable also without a previous plan; in fact, the previous plan monitoring database can be replaced by the simulation of the distribution of volumes in existing quarries or in not yet exploited deposits.

Discussion
The proposed Analyzing PSS is usable also without a previous plan; in fact, the previous plan monitoring database can be replaced by the simulation of the distribution of volumes in existing quarries or in not yet exploited deposits.
The assigning of a "Premium Quota" is the innovative part of the methodology, as it leads to reflection on the duality of the quarrying action concerning its public interest: on one hand the provision of raw materials, on the other the protection of the environment. The quantitative criterion of proportionality is functional not to excessively change the size of available deposits and to ensure a balanced supply of material throughout the territory: quarrying will go on less in a small deposit than in a very large one. Additionally, the proportionality criterion guarantees business continuity to those companies that have made investments by installing extraction and processing plants for raw materials.
Furthermore, it must be underlined that "the extractive industry of sand, gravel and rock differs in at least two important aspects from most other industrial sectors. First, the location of the industry and the quality of the material produced are determined by geology in relation to areas of high demand, i.e., large cities. Second, the success of the sector is dependent on the success and competitiveness of the downstream industries. Economic success, therefore, is controlled by two factors: quality and quantity of deposits of natural resources for aggregates and the political, legal, administrative, social, and economic environment in which extraction takes place (Wagner/Tiess, 2004: 27). For both reasons, aggregate companies are used to taking a long-term perspective" [41].
The premium boosts mining companies' proper behavior. LQN identification is a strength, because it allows the balancing of the response for territorial needs, according to the actual availability of the deposits. From the theoretical point of view, considering LQN, the distribution of volumes reduces the negatives deriving from the transport of raw materials. The presence of LQN when redistributing volumes favors a fair distribution on the territory of the quotas of excavable volumes.
Among the limits, the Analyzing PSS requires an ad hoc structured database. The preliminary step of selecting and acquiring the base information is particularly toilsome. Since it is an iterative method, to reach the definition of the "Final Quota" assigned to each quarry, it is necessary that each step is faced and resolved, as each variation affects the whole system.
In the proposal of the new plan, 38 quarries were confirmed to satisfy a demand of 46,000,000 m 3 of primary sand and gravel (input data 1 and 34% lower than PPC 2005), no new quarries were planned, and sixteen quarries were cancelled. Among the 38 confirmed quarries, nineteen were expanded and eight were reduced.
The construction of the new plan allowed an important shift from the traditional planning process to a GIS-based Analyzing PSS approach based on MAA.
In the case of Brescia, defining the LQN causes severe limitations due to the criteria with which the units were defined: they are only partially based on territorial and socioeconomic features (such as urban density, self-contained labor market, population size, and accessibility) in favor of corporate nature agreements. Therefore, political geography prevails over economic geography. Additionally, the political decision-maker did not demand setting limits to quarries' expansion where there is an exploitable quality deposit: it would have been advisable to identify maximum limits in the volume redistribution mechanism, so as to avoid excessive pressures on the territory and excessively large areas.

Conclusions
The objective of this study was to understand if it is possible to assist the public authority in defining both the proper amount of sand and gravel to be extracted and the size dimension of each quarry. This goal can be achieved through the application of the proposed GIS-based Planning Support technology. However, we must be aware that the final performance of the Analyzing PSS is affected by decision-maker choices and conditioned by economic interests at the regional level. Planning action is complex and subject to economic and environmental pressures. In the case of raw materials quarrying, economic interests (both mining companies' ones and in relation to the economy of the sector) often prevail over environmental protection. A non-trivial pressure is represented by the induced activity in the labor sector, where in the case study, the quarries represent 3% of the provincial GDP. Environmental pressure is determined by the concentration of quarry activities due to the raw materials in the subsoil. In Italy, such impacts are partially offset by the charge on aggregates (to compensate for the environmental costs caused by a quarry, and to protect the landscape) that "is only one element of a very complex planning, authorization, and regulation system related to quarrying activities" [41].
Finally, the proposed PSS does not automatically build alternative scenarios from which to choose the best one, but it implies making political decisions during the various steps. The Analyzing PSS objectifies volumes' assignment, but it implies a well-defined governing direction. In the study case, the decision-maker's choice to identify a percentual of the total premium quote equal to 20% of the sand and gravel demand did not allow the adequate exploitation of the premium potentiality.  Data Availability Statement: Restrictions apply to the availability of these data. The data are not publicly available due to the cooperation agreement submitted by authors with the Province of Brescia.

Conflicts of Interest:
The authors declare no conflict of interest.