Milk and cheese production are important contributors of the Trentino A. A. (a mountain region in Nord-East Italy) economy. In 2015, 135.094 tons of milk were produced, mostly curled and strained for cheese production by a large number of small cheese plants scattered around the region, processing on average 14,000 liters of milk per day. The region is very sensitive to the circular economy paradigm and since a long time is trying to recycle the waste from agri-food activities. An important regional project is dedicated to recycle the whey, a polluting by-product of cheese production, representing the 85% - 90% of the milk transformed in cheese. Using appropriate technologies is possible to separate the whey components and sell in different market channels: animal feedstock, proteins, vitamins for human consumption, lactose for PHA and others. Simulation about milk quota removal in Italy [1] , estimate a loss of the producers’ incomes about 4 billion €, caused by the fresh milk prices decrease between 5 and a 10% (under the baseline scenario) and a consumers’ surplus gain of 3.7 billion €. [1] Figure 1 reports some historical changes in dairy sector in Italy: the number of dairy farms declined during the period 1995-2015 from 97,044 to 32,994 and the production increased from 10.4 million tons to 11.4 million tons. The 35% of largest dairy farms produce between 200 and more than 2000 tons per year, and cover the 88% of the total production [2] . Despite these events, the milk production still represents an interesting solution for many farmers; however consistent structural and organizational changes are required to increase the competitive advantage versus international competitors [3] .

Milk, cheese and whey productions are strictly correlated: 10 kg of fresh milk produce approximately 1 kg of cheese and 9 kg of whey; the total cheese production in Italy it is estimated a quantity of whey between 8 and 9 million tons; lactose is the most important component (40-45 gr. per Kg of whey), and is responsible of the high values of BOD (Biochemical Oxygen Demand) (BOD: 40,000 - 60,000 ppm) and COD (Chemical Oxygen Demand) (COD: 50,000 - 80,000 ppm) if released into water bodies [4] . According to Wissmann et al. [5] the pollution caused by 50,000 liters (13,209 gallons) of whey is equivalent to a city settlement of 25,000 inhabitants. Nevertheless the whey can be a potential economic resource for the quantity of valuable components, at this moment only a small fraction of them are separated with ultrafiltration, fermentation, inverse osmosis methods and sold through different market channels [6] . The 53% of the whey produced in Italy is dried, transformed in powder and exported in Germany and France excluding WPC-Whey Protein Concentrate and WPI-Whey Protein Isolate. In

Italy the largest whey quantity is used for animal feedstock (65% of the total consumption), another 20% is sold as infant formulas and the remaining 15% is used in chocolate, ice cream, bakery and confectionery industry [9] [10] .

In USA the whey powder used for animal feeding has a lower incidence (estimated 45%); most of the whey is sold to the dairy industry. A growth of the whey consumption is expected in the nutritional segment; used in nutritional formulations such as whey powder, demineralized whey, WPC and WPI, whose demand in the health, pharmaceutic and nutritional sectors is expected to grow in next years, an interesting development is expected also in the Biopolymer industry. The derived whey products are growing at a rate of 3% per year, mostly for whey powder and lactose. With the progress in whey processing technologies, new market opportunities are disclosed to operators and the logistic of transport, packaging, storage, conservation and the environmental impact are becoming growingly important for the competitiveness of the dairy chain.

The purpose is to afford the transport problem due to shipping growing volumes of liquid whey at minimum transport cost from CP (cheese production) to WPP (whey processing) by selecting the optimal route to reduce cost and environmental impact. Preliminary information is requested about the dimension of the whey collection basin, transport costs, distance from (CP) to processing plant, (WPP), type of road (high way, state, trigonal or provincial and communal roads), road conditions, traffic intensity during the day, number of city crossed, physical obstacles orographic nature¹. This information will be used to minimize both the transport costs, and the environmental impact caused by CO₂ emissions [9] [10] [11] [12] [13] .

This logistic problem requires to define:

1) The algebraic formulation of the objective function and constraints;

2) The balance condition that is the sum of the supplies of all the sources equal to the sum of the demands for all the destinations;

3) The selection of an iterative process to emulate the simplex algorithm, by starting with an initial basic feasible solution (IBFS),check for the demand- supply conditions and proceed iteratively to find the final optimal solution [10] [14] [15] .

The remainder of this paper is structured as follows: Section 2 introduces the network theory, to find the optimum transport solution; Section 3 describes the whey supply in the basin with distances among CP and road network condition; Section 4 describes the case study based on the optimization approach; Section 5 reports the comments about results of simulation, policy implications and suggestions to improve the whey collection.

The contribution of the present paper to the transport problem is twofold. The first one is to introduce a novel formulation that extends a globally inclusive facility hierarchy problem [16] applied to whey collection and processing. The second is the contribution to the circular economy by reducing the impact caused by whey dispersion, and CO₂ emission caused by transport distance².

²The WPP will use different technologies as concentration, demineralization, ultrafiltration, crystallization, PHA to separate the whey ingredients.

Definition: k is the homogeneous commodity produced by a given CP and delivered through a route to the WPP that can be located in one of the CP of the network, assuming that any CP can be also a potential location of the WPP. The solution consists in finding the optimal route to minimize the transport costs and environmental impact caused by CO₂ emissions [8] . The environmental impact is measured by assuming that the emission is a fixed quantity per km, then the total emission is correlated to the route distance.

The following data are required to define this problem:

k i quantity of whey produced by a given CP and delivered to the next CP to the end WPP, x_i,j, variable that indicate if the link between nodes i and j is open; c_ij is the unit cost of shipping the whey along the a given branch of the network.

Figure 2 describes the whey route from node i (departure) to node j (destination); with the assumption that CP = WPP, i = 1 , ⋯ , n and j = 1 , ⋯ , n .

The logistic of transport is an allocation problem illustrated in Figure 2: m sources (CP) and n destinations (WPP) in this case n = m because every destination can

be located in one of the sources; the whey route from node i (departure) to node j (destination). The m sources can ship the product to any of the n destinations at per unit carrying cost c i j (unit transportation cost from source i to destination j). The transport problem requires to minimize the shipping cost of the whey commodity from a source CP defined the supply node i to a destination node j for j = 1..n defined the demand node. Along the transitory nodes it is completed the route from source i to destination node j. As the CP (cheese plant) are interchangeable with WPP (whey processing plant), the possible WPP locations are equal to CP plants [17] so that the number of sources corresponds to the number of destination. The cost depends on the number of truckloads (x_ij) shipping a given quantity of whey from source i to destination j and c_ij is the cost per unit of commodity shipped from supply i to destination j; this means that the transport cost depends on the distance. The accumulation of whey delivered by the all CP_i represents the total supply then s i = ∑ i = 1 m C P i and the cumulated demand d_i of the whey from all different CP_j is d j = ∑ j = 1 n C P j . In Figure 3 is reported an example the oriented graph network; node 1 is the origin (departure) and node 5 is the destination; two different routes are hypothesized in the example reported in Figure 3.

The O.F is targeted to find the minimum transport cost of the product that passes through the branches from the initial node i to the final node j in a route constrained by the quantity of product shipped from origin to the destination. The Figure 3 describes the network formed by 5 nodes and 8 branches: node 1, the origin is connected to the other nodes with 3 branches, node 2 has 4 optional branches, node 3 has 3 alternative branches, node 4 has 4 possible branches, node 5 the destination has 2 branches. The node 1 is the origin and node 5 is the destination; the example allows two routes: the route A is crossing the nodes 1, 2, 5 and gives the solution 1; the route B is crossing the nodes 1, 3, 4, 5 and gives the Solution 2. The branches in red represent the route pattern, the branches in bleu are alternative choices and proceed toward the final node or

in opposite directions meaning that the flow direction is inverted .

With the route A, the branch connecting nodes 1 and 2 measures 6 and the flow is equal 5; the branch connecting nodes 2 and 5 measures 4 and the flow is 7. The final solution is given by a total flow equal to 12, the total distance is 10 and the total transport cost, assuming the unit transport cost c = 1 the total cost is 58 that is the value of the route A.

With route B the length of the branch 1 - 3 is 1 and allows a flow equal 5, the branch 3 - 4 measures 2 and allows a flow equal 4; the branch 4 - 5 measures 3 and allows a flow equal 3. The total flow is 12 equal to the solution 1 but the total distance is 22 and assuming a cost c = 1, the solution is a transport cost = 22 inferior to solution 1. Other routes can be hypothesized as the 1-3-4-2-5 or 1-3-4-1-2-5 but they are less efficient in normal route conditions to Solution 2.

The first step is to draw the graph of the road network for the collection basin reported in Figure 4, the 16 CP (cheese plants) are the 16 nodes of our network in the Trentino A. A., a northern region of Italy, delivering the whey every day. The optimal position of the WPP will minimizes both the transport costs and

environmental impact from CP to the WPP. Due to the orographic nature of this region, the road network is developed along vertical lines following the mountain compluvium lines: one is a highway and four state or province roads. The horizontal road development is represented by three state and province roads. This network configuration allows different route choices then the problem is to find the optimal route that will optimize the transport cost and environmental impact. The network indicates the concentration of ten CP in a restricted area of 45 × 31 square km that deliver the 55% of the total daily whey production; other six plants are distributed in an area of 93 × 61 square km., that is four times larger compared to the first one and offers only the 45% of the total daily whey supply. This non homogeneous distribution of the CP affects the transport costs and will be taken into account to select the optimal transport solution.

The map reported in Figure 4 indicates the CP distribution in the collection basin. The graph is oriented to Trento where is located the processing plant, then a hypothesis of solution is to find the optimal route connecting all CP to the WPP located in Trento. Other WPP locations will be simulated for the cost minimization due to the different CP concentration. Observing the CP distribution and the road map it is possible to select three alternative routes to optimize the transport costs; these are indicated in Figure 5. The first route is a circuit connecting ten CP concentrated in the 1^st quadrant (see Figure 4) starting with Mezzana, the whey is shipped to the final destination Trento. A second route is a circuit that includes the CP located in Val di Fassa, Predazzo and Cavalese (2nd quadrant of Figure 4) and final whey delivery to Trento; the third route is a circuit including CP located in Fiera Primiero and Lavarone (3rd quadrant) and delivery to Trento.

In Table 1 is reported the daily quantity (hundred Kilos) of whey collected from the 16 CP and collection cost from CP to WPP.

From the literature two empirical procedures are selected to solve the transport

problem: The first one is a numerical solution that emulate the simplex method and proceeds with successive iterations starting with an initial minimum cost value, proceeding to the next minimum cost and finally it is obtained the final minimum cost value. At the beginning the minimum cost cell is selected the corresponding row (supply) and column (demand) is selected at the crossing of row (supply) and column (demand) and the residual is calculated by finding the positive difference between demand and supply (or vice versa) and proceeds by finding the next minimum cost. [12] [13] .

The second procedure is more empirically oriented, and uses the observation of the network to choose the preferred route.

The advantage is that if the network features and road conditions are known it is simpler to solve the problem and allow to perform easily many simulations. Alternatively one can use some algorithm (i.e. the Dijkstra algorithm) to find the shortest route.

The network is composed by i = 1, ・・・, m sources (CP suppliers) and j = 1 , ⋯ , n destinations (customers): each sources CP_i supplies x i j quantity of whey to the next j destinations at per unit shipping cost C i j that is the transport cost per unit of whey shipped from source i to destination j. Each source i for 1 ≤ i ≤ m delivers si quantity to n destinations (customers) and the destination demands d_j for 1 ≤ j ≤ n . The supply constraints i is the flow limit of whey from the origin supply node i to the n destination (consumption) nodes j; the demand constraints (see column) indicate the quantity of product from the all origin nodes i allocated to one destination node j.

The total supply is the sum of the product delivered by a given CP to destination, then:

s i = ∑ i = 1 m s i j ; the demand is the sum of d_ij units, offered by the all CP_j, then d j = ∑ i = 1 n d i j this is the demand at every location. If origins and destinations coincide the problem can be represented by a square table with source i equal to destination j so that m rows equal to n columns.

The transport problem consists in the formulation of the objective function Z that is the transport cost minimization bounded with the supply and demand constraints as indicated below:

Min Z = ∑ i = 1 m ∑ j = 1 n c i j ∗ x i j (for all branches)

subject to

∑ j = 1 n x i j ≤ s i for i = 1 , 2 , ⋯ , m (supply constraint; row m)

∑ i = 1 m x i j ≤ d j for j = 1 , 2 , ⋯ , n for j = 1 , 2 , ⋯ , n (demand constraint; column n)

∑ i = 1 m s i = ∑ j = 1 n d j (equality constraint for balance problem)

X i j ≥ 0 for all i, j (positive or non negatives)

The OF with supply and demand constraints are reported in extended notation below:

Min Z = c 11 x 11 + c 12 x 12 + ⋯ + c 1 n x 1 n + c 21 x 21 + c 22 x 2 , 2 + ⋯ + c 2 n x 2 n + ⋯ + c m 1 x m 1 + c m 2 x m 2 + ⋯ + c m n x m n

subject to:

1) Supply constraints (row 1…16)

x 11 + x 12 + ⋯ + x 1 n ≤ s 1

x 21 + x 22 + ⋯ + x 2 n ≤ s 2

x m 1 + x m 2 + ⋯ + x m n ≤ s m

2) Demand constraints (column 1…16)

x 11 + x 21 + ⋯ + x m 1 ≤ d 1

x 12 + x 22 + ⋯ + x m 2 ≤ d 2

x 1 n + x 2 n + ⋯ + x m n ≤ d n

The Table 2 reports the information required to solve the problem:

1) For each branch i-, j is reported: The distance d i j in km from node i to node j, adjusted with time varying with local network conditions affecting the transport difficulty. The adjusted distance in km X_ij is included in the OF.

2) c i j is the cost of shipping one unit of whey from i to j then c i j = 0.67 ∗ d i j expressed in €/100Kg.

The data of Table 2 are used to elaborate the transport problem reported here:

F . O 10.05 x 12 + 26.13 x 13 + 22.78 x 14 + 30.82 x 15 + 32.83 x 16 + 33.5 x 17 + 30.15 x 18 + 26.80 x 19 + 22.78 x 1 , 10 + 43.55 x 1 , 11 + 64.32 x 1 , 12 + 75.04 x 1 , 13 + 88.44 x 1 , 14 + 93.8 x 11 , 5 + 76.38 x 11 , 6 + 10.05 x 2 , 1 + 18.76 x 2 , 3 + 18.09 x 2 , 4 + ⋯ + 68.34 x 16 , 15

s. to

supply constraint (nr = 16)

x 1 , 2 X 1 , 2 + X 1 , 3 + X 1 , 4 + X 1 , 5 + X 1 , 6 + X 1 , 7 + X 1 , 8 + X 1 , 9 + X 1 , 10 + X 1 , 11 + X 1 , 12 + X 1 , 13 + X 1 , 14 + X 1 , 15 + X 1 , 16 ≤ 310.63

x 16 , 1 + X 16 , 2 + X 16 , 3 + X 16 , 4 + X 16 , 5 + X 16 , 6 + ⋯ + X 16 , 55 ≤ 143.13

demand constraint (nr = 16)

x 1 , 1 + X 2 , 1 + X 3 , 1 + X 4 , 1 + X 5 , 1 + X 6 , 1 + X 7 , 1 + X 8 , 1 + X 9 , 1 + X 10 , 1 + X 11 , 1 + X 12 , 1 + ⋯ + X 13 , 1 + X 14 , 1 + ⋯ + X 16 , 1 ≤ 329.01

x 1 , 16 + X 2 , 16 + X 3 , 16 + X 4 , 16 + X 5 , 16 + ⋯ X 16 , 16 + ⋯ + ≤ 485.20

Procedure 1―To find the optimal solution is required the following steps:

Step 1―Define the OF: the transport cost minimization by selecting a route composed by nodes (CP) and branches, supply (row) and demand (column) constraints;

Step 2―Check for the balance condition: sum of row values equal to sum of column values;

Step 3―Find the minimum transport cost value c_ij in the transport Table 2 and select the corresponding cell;

Step 4―First allocation: find the smaller value by comparing s_i and d_j referred to cell c_ij; (in our case min c_ij = 2.01 in cell 8, 7 and the smaller value between s i = 362.53 and d j = 239.37 is s i . Allocate X i j = supply value (in this case is 23,937 in the corresponding cell (i, j) and compute the difference: 362.53 - 239.37 = 123.16 that is the residual supply value.

Step 5―Second allocation: search for the next minimum cost C i j corresponding to value 3.35 in cell 7, 6 and s 7 = 185.73 and d 6 = 247.21 . The residual supply s 7 = 123.16 is is compared with d₆ that is greater so Residual d j = 247.21 − 123.16 = 124.05 and supply go to 0; this value is the new reduced demand allocated in the new cell c i j . Proceeding with these operations, the supply and demand requirements are progressively allocated, allowing to compute the partial costs of allocation at each step.

The min C_ij from the 16 columns reported in Table 2 are listed here:

c 8 , 7 = 2.01 ; c 7 , 6 = 3.35 ; c 6 , 8 = 5.36 ; c 6 , 5 = 6.03 ; c 2 , 1 = 10 , 05 c 10 , 9 = 10.72 ; c 4 , 10 = 11.39 ; c 4 , 3 = 14.74 ; c 10 , 2 = 18.09 ; c 9 , 11 = 24.79 ; c 13 , 12 = 24.80 ; c 11 , 16 = 36.85 ; c 13 , 15 = 42.88

The cost minimizations obtained with successive iterations are reported in Table 3; 15 iterations were required to find the optimal allocation for the most general network configuration that included all branches of the network

Simulation 2: For this simulation, the all CP of the first route are excluded from route 2 and 3 (value = 0) a priori because inefficient in term of distance and time as suggested by the Google map and roads previously observed in Figure 4 and Figure 5 then the new Table 4 is used for finding the optimal solution.

In Table 5 are reported the intermediate transport values with elimination of some road alternatives.

The optimal transport values from the two roads are the following:

Simulation 1 = 51295;

Simulation 2 = 46319

This procedure can be adopted preferably in case some information are available ex ante as graph map, distance, road condition, or preference about the route that could facilitate the search for the optimal transport solution.

The problem is the same, minimization of the transport cost: c_ij is the unit cost of transport; the 16 nodes and branches are forming the network that shows possible alternative routes to ship the whey from nodes i (origin) to nodes j (destination), X_ij is the quantity of whey shipped from i to j. We start by solving the problem previously illustrated in Figure 3 with 5 nodes and 8 branches. The problem is to find the minimum transport cost of the 8 branches subject to the 5 node constraints.

O F : 6 X 1 , 2 + X 1 , 3 + 4 X 2 , 5 − 3 X 3 , 2 + 2 X 3 , 4 + 4 X 4 , 1 + X 4 , 2 + 3 X 4 , 5

s. t. (maximum flow through the node-archs; the values are selected with the branch capacity) node-arch 1: X 1 , 2 + X 1 , 3 − X 1 , 4 = (1 branch with negative versus, max. capacity = 5)

Node arch 2: X 2 , 5 − X 2 , 1 − X 2 , 3 − X 2 , 4 = − 4 (2 branches with negative versus, max. capacity = −4)

Node arch 3: − X 3 , 2 + X 3 , 4 − X 3 , 1 = 0 (1 branch with negative versus, max. capacity = −4

Node arch 4: X 4 , 1 + X 4 , 2 + X 4 , 5 − − X 4 , 3 = 1 (1 branch with negative versus, max. capacity = 1

Node arch 5: − X 5 , 2 − X 5 , 4 = − 3 (2 branches with negative versus, max. capacity = −3)

Route options: The Figure 3 gives the information to create the incidence matrix node-arch A: for the 5 nodes v (1…5) and arch e (1…8) the corresponding entry A_ve is:

+1 if e exit from v (v is the tail of e = positive branch direction);

−1 if e entry in v ((v is the head of e = negative branch direction);

0, otherwise (not allowed flow).

The incidence matrix described in Table 6 reports the constraints: the value 1 means that the transport has the positive versus from the origin node i to the destination nodes j; the negative sign indicates the opposite direction required when one of the positive branches is not allowed.

Parameters

0 ⊆ V ; origin nodes

D ⊆ V ; destination nodes

T ⊆ V ; passage nodes

s_i, for iÎ0, supply product vertix i

d_i, for i ÎD, demand product vertix i

X_ij for i, jÎA flow of product on the arch i, j

u_ij for i, jÎA, capacity of the arch i, j (maximum admittable flow on the arch)

c i j ∗ x i j , (i, j)ÎA is the transport cost of the flow X_ij on the arch i, j

The second version of the transport problem is represented in Table 7 in LP notation:

The values at the right and bottom sides of the transportation Table 7 indicate respectively: The row reports for a given source of whey s_i for i = 1, ⋯ , 16 the number of possible destination (allocation) j for j = 1 , ⋯ , n ; for n = 16

A column indicates the whey supply from the sources s_i for i = 1, ⋯ , 16 to the m demand d_j for j = 1 , ⋯ , m that are the quantity supplied to each demand. In Table 8 is reported the numeric solution of the problem.

C i j for 1 ≤ i ≤ m and 1 ≤ j ≤ n is the unit cost of shipping the whey from the i-th source to the j-th destination;

X i j is the quantity of whey shipped from i-th source to j-th destination;

C i j ∗ X i j is the total transport cost from source (node i) to destination (node j) for all i, j. pair combinations (from the first to the last allocation).The cost cells are distributed in continuous³:

O F : C 1 , 1 X 1 , 1 = 464 X 11 ; C 1 , 2 X 1 , 2 = 513 X 1 , 2 ; C 3 , 4 = 685 X 3 , 4

For the equilibrium condition it is required:

∑ i = 1 m S i = ∑ i = 1 n d j then (see Table 8):

750 + 1250 + 1000 = 800 + 650 + 700 + 850 = 3000 i.e. demand = supply then this problem is balanced.

This procedure is used to solve the problem of the minimum transport cost from the CP origins to WPP destination with cost simulations of some predefined routes to fulfill specific objectives of the operators. The Table 5 previously described is adapted to this problem by using the incidence matrix to define the preferred routes. The optimal value is obtained from the following operation:

MATR. SOMMA. PRODUCT (B98:S115; B123:S140)

The first matrix B98:S115 reports the whey quantity multiplied by the adjusted distances (see Table 5), and transport cost from origin to destination; the 2nd matrix (B123:S140) reports the numeric coefficients of the incidence matrix: to generate the routes and calculate the corresponding costs.

The route is specified with coefficient value = 1 for positive direction and −1 for the opposite direction of the branch while the non activated branch is indicated with nul value 0. The objective function minimizes the transport cost between a given CP (origin) and the WPP (destination); the total cost will be calculated by

summing the costs of various branches composing the route: C T = ∑ ( C i j ∗ X i j ) (see Figure 3.1). Each node can only have one path to it and one path from it. The horizontal constraints are the supply (from) constraint and the vertical constraints are the demand (to) constraints. Forcing the sum of the “from” origin constraints to be equal to 1 will force the solver to choose only one path from each node. Forcing the sum of the “to” destination constraints to be equal to 1 will force the solver to choose only one path for a given node. The first group of simulations is computed by shipping the whey from the all CP to Trento destination that is also the main whey producer and where it is located the WPP location (Table 9).

Table 10 reports the incidence matrix with transportation cost and delivery to Trento following the routes designed by the operator. The first route is the circuit encompassing ten CP of the first quadrant in Figure 4 with starts to CP: Mezzana, Terzuolas, Rumo, Revo, Castelfondo, Fondo, Cavareno, Rumeno, Coredo Tuenno with end to Trento; the second route is composed by three CP: Cavalese, Predazzo and Val di Fassa and end to Trento; the third route is the circuit Fiera di Primiero and Lavarone and end to Trento. Assuming the whey delivery to Trento, the OF value is 62,211 €. Other three simulations are tried by shifting the final destination (F.D). The results are: F.D-Tuenno: OF = 58,257; F.D-Coredo: OF = 76,706; F.D-Cavalese: OF = 213,961. The results suggest these considerations: the optimal solution is found in the area where the higher supply of whey is compatible with minimum transport cost; Tuenno is the optimal location for whey delivery in absence of WPP contraints. The plants are presently located in Trento because is the highest whey supplier and the whey collection is easier. The 2nd procedure is quite interesting as it reveals the optimal position of WPP in Tuenno. Compared to Trento, the WPP located in Tuenno reduces the present transport cost of 6.3%; the WPP located in Coreno increases the transport cost of 23% and location in Cavalese increases the cost of 244%.

Compared the second with the first procedure, the transport cost with destination Trento are quite similar, the difference is 15.5% and is explained by the differences in route 2 and 3 followed to ship the product to Trento. Sustainability of the collection strategy is an important collateral effect of the transport: the available technologies to process the whey contribute to limit the dispersion of whey pollutants in the environment and demonstrate that economic and environmental targets of the whey processing can be simultaneously obtained. The production, processing and transportation of milk products, contributes with 2.7 percent to the global anthropogenic greenhouse gas emissions in Italy. According to the methodology of the IPCC [4] , the total annual CO₂ emissions from the transport sector calculated as follows:

The whey transport emission is computed by assuming that CO₂ emission of a normal diesel truck with capacity of 20 ton and emission of 20 g CO₂ eq per ton/ km. This coefficient has been used to compute the whey pollution due transport. The CO₂ emissions are: route 1 = 15.20; route 2 = 12.05; route 3 = 245.29; total