# The relationship of dairy farm eco-efficiency with intensification and self-sufficiency. Evidence from the french dairy sector using life cycle analysis, data envelopment analysis and partial least squares structural equation modelling

Andreas Diomedes Soteriades, Alistair William Stott, Sindy Moreau, Thierry Charroin, Mélanie Blanchard, Jiayi Liu, Philippe Faverdin

## ► To cite this version:

Andreas Diomedes Soteriades, Alistair William Stott, Sindy Moreau, Thierry Charroin, Mélanie Blanchard, et al.. The relationship of dairy farm eco-efficiency with intensification and self-sufficiency. Evidence from the french dairy sector using life cycle analysis, data envelopment analysis and partial least squares structural equation modelling. PLoS ONE, 2016, 11 (11), pp.e0166445. 10.1371/journal.pone.0166445. hal-02631798

HAL Id: hal-02631798

https://hal.inrae.fr/hal-02631798

Submitted on 27 May 2020

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

CrossMark
click for updates

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Soteriades AD, Stott AW, Moreau S, Charroin T, Blanchard M, Liu J, et al. (2016) The Relationship of Dairy Farm Eco-Efficiency with Intensification and Self-Sufficiency. Evidence from the French Dairy Sector Using Life Cycle Analysis, Data Envelopment Analysis and Partial Least Squares Structural Equation Modelling. PLoS ONE 11(11): e0166445. doi:10.1371/journal. pone.0166445

**Editor:** Nikolaos Georgantzis, University of Reading, UNITED KINGDOM

**Received:** March 23, 2016

**Accepted:** October 28, 2016

**Published:** November 10, 2016

**Copyright:** © 2016 Soteriades et al. This is an open access article distributed under the terms of the [Creative Commons Attribution License](http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** Funding was obtained from the following sources: 1. Scotland's Rural College (http://www.sruc.ac.uk/) (ADS); 2. Scottish Government Strategic Research Programme 2011-16 (http://www.gov.scot/Topics/Research/About/EBAR/StrategicResearch/future-research-strategy)

RESEARCH ARTICLE

# The Relationship of Dairy Farm Eco-Efficiency with Intensification and Self-Sufficiency. Evidence from the French Dairy Sector Using Life Cycle Analysis, Data Envelopment Analysis and Partial Least Squares Structural Equation Modelling

Andreas Diomedes Soteriades¹*, Alistair William Stott¹, Sindy Moreau², Thierry Charroin², Melanie Blanchard³, Jiayi Liu⁴, Philippe Faverdin⁵,⁶

1 Scotland's Rural College, Future Farming Systems Group, Edinburgh, United Kingdom, 2 Institut de l'Elevage, F-75000 Paris, France, 3 CIRAD - INRA - Montpellier SupAgro, UMR 0868 SEMLET Systèmes d'Elevage Méditerranéens et Tropicaux, F-34000, Montpellier, France, 4 Biomathematics and Statistics Scotland, Edinburgh, United Kingdom, 5 INRA, UMR 1348 PEGASE, F-35590 St-Gilles, France, 6 Agrocampus-Ouest, UMR 1348 PEGASE, F-35000 Rennes, France

* andreas.soteriades@sruc.ac.uk

## Abstract

We aimed at quantifying the extent to which agricultural management practices linked to animal production and land use affect environmental outcomes at a larger scale. Two practices closely linked to farm environmental performance at a larger scale are farming intensity, often resulting in greater off-farm environmental impacts (land, non-renewable energy use etc.) associated with the production of imported inputs (e.g. concentrates, fertilizer); and the degree of self-sufficiency, i.e. the farm's capacity to produce goods from its own resources, with higher control over nutrient recycling and thus minimization of losses to the environment, often resulting in greater on-farm impacts (eutrophication, acidification etc.). We explored the relationship of these practices with farm environmental performance for 185 French specialized dairy farms. We used Partial Least Squares Structural Equation Modeling to build, and relate, latent variables of environmental performance, intensification and self-sufficiency. Proxy indicators reflected the latent variables for intensification (milk yield/cow, use of maize silage etc.) and self-sufficiency (home-grown feed/total feed use, on-farm energy/total energy use etc.). Environmental performance was represented by an aggregate 'eco-efficiency' score per farm derived from a Data Envelopment Analysis model fed with LCA and farm output data. The dataset was split into two spatially heterogeneous (biophysical conditions, production patterns) regions. For both regions, eco-efficiency was significantly negatively related with milk yield/cow and the use of maize silage and imported concentrates. However, these results might not necessarily hold for intensive yet more self-sufficient farms. This requires further investigation with latent variables for intensification and self-sufficiency that do not largely overlap- a modelling challenge that occurred here. We conclude that the environmental 'sustainability' of intensive dairy farming depends on
(ADS); 3. Physiologie Animale et Système d'Elevage division of the Institut National de la Recherche Agronomique (French National Institute for Agricultural Research; http://www.inra.fr/en)
(ADS); 4. Project 'Efficience Environnementale et productions animales pour le développement durable', ANR-09-STRA-01-01, Agence Nationale de la Recherche (French National Research Agency; http://www.agence-nationale-recherche.fr/en/project-based-funding-to-advance-french-research/) (MB); 5. Scottish Government Strategic Research Programme 2011-16 (http://www.gov.scot/Topics/Research/About/EBAR/StrategicResearch/future-research-strategy)
Consultancy for Theme 4 WP4.1 Adaptation to change in land-based and other rural industries (JL). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

particular farming systems and circumstances, although we note that more self-sufficient farms may be preferable when they may benefit from relatively low land prices and agri-environment schemes aimed at maintaining grasslands.

## Introduction

Meeting the world's rapidly growing food demands in perpetuity while preserving the environment and the planet's natural resources is an enormous challenge for agriculture [1]. The European Union (EU) is attempting to address this challenge through the Common Agricultural Policy (CAP), whose earlier focus on market price support and recent shift to direct payments and abolition of milk quotas have resulted in an intensification of dairy farming in the past few decades [2,3]. However, several dairy studies find that a trade-off might exist between dairy farm intensification and environmental performance [4–7] because intensification often has externalities associated with the production of imported inputs.

Assessing the relationship between intensification and environmental performance requires 'global' metrics able to capture both on and off-farm environmental impacts of dairy farming. Global metrics can be calculated with Life Cycle Analysis (LCA), an internationally standardized method for estimating the environmental impacts of agricultural products from a whole-system perspective [8]. Numerous studies use LCA to assess dairy farm environmental performance by calculating 'eco-efficiency' ratios, that is, environmental impacts expressed per unit of milk or land area [5,6,8,9]. Yet eco-efficiency ratios have several drawbacks [10], for example the allocation of environmental impacts to several dairy farm products (milk, meat, crops) is challenging. Dairy studies are therefore increasingly coupling LCA indicators with the multiple-input, multiple-output method Data Envelopment Analysis (DEA [11]) to calculate single aggregated eco-efficiency scores per farm, by accounting for all LCA impacts (or carbon footprinting indicators), inputs and outputs simultaneously [10,12–19].

With global eco-efficiency indicators in hand, the next challenge is to identify farm management strategies that can improve eco-efficiency. The global efficiency indicator is based on a ratio between productions of goods and inputs (or environmental impacts). A wealth of dairy studies [5,6,8,10,13,19–22] show that eco-efficiency is generally influenced by two management strategies improving this ratio with an increase in outputs or a decrease in inputs: (i) the level of intensification at animal and farm-levels, i.e. higher production per unit of input so that increases in outputs outweigh potential increases in inputs and environmental impacts; and (ii) the farm's degree of self-sufficiency, i.e. its capacity to produce goods from its own resources [21], so that decreases in inputs and environmental impacts outweigh likely decreases in outputs. In fact, self-sufficiency, with higher control over nutrient recycling and thus minimization of losses to the environment, can be considered as a key agro-ecological principle, central to improving the sustainability of livestock systems [23].

Nevertheless, so far a few dairy studies have attempted to holistically measure self-sufficiency, production intensity and environmental performance and to study their relationships from a whole-system perspective [21]. Certainly, this reflects the challenge of obtaining sufficient data on all three aspects [21]. Moreover, specific relationships between these aspects cannot be predetermined; it is a managerial choice to intensify production while being less or more self-sufficient, with implications for the environment in either case. On the one hand, intensive and less self-sufficient systems may increase production in a very resource use-efficient manner. For example, increasing cow productivity and the use of concentrated feeds can
increase feed efficiency and milk yield per cow and per ha ($10^4 m^2$) on-farm land, and can decrease manure volume and nutrients [5,19,24]. Nonetheless, high efficiency may not be enough to reach environmental targets [25]. On the other hand, intensive and more self-sufficient systems recycle more elements on-farm, with lower environmental impacts per ha but not necessarily per unit of milk [5].

By coupling LCA data with DEA, Soteriades et al. [10] found that the less the proportion of maize silage in the total forage area of French specialized dairy farms, the higher the farms' eco-efficiency. This confirmed earlier findings that increasing maize silage in dairy farms generally implies higher animal and farm-level farming intensity and lower farm self-sufficiency (higher supplementation with protein-rich feeds, usually imported soybean meal) [5,8,26,27]. They also showed that eco-efficiency performance can be region-dependent. Their findings suggested a possible relationship between dairy farm eco-efficiency, intensification and self-sufficiency, which could differ between regions, although, as commented earlier, these relationships cannot be predetermined. Additionally, a more holistic analysis of this relationship should involve more indicators representing these three aspects. From a modelling perspective, such cases can be dealt with the multivariate method Partial Least Squares Structural Equation Modelling (PLS-SEM [28,29]). PLS-SEM can account for multiple indicators to search for latent patterns in the data when there is no or only little prior knowledge on how the variables are related, while accounting for (regional) heterogeneity [30].

The objective of this paper was to explore the relationships between dairy farm eco-efficiency, self-sufficiency, and farm and animal-level intensification, by extending the LCA-based dairy farm exercise of Soteriades et al. [10] that accounted for impacts both on and off-farm. That way it was possible to determine whether or not on-farm management practices promoting self-sufficiency and/or increasing input per ha and output per unit input- especially per ha- could explain dairy farm environmental performance at a larger scale. In other words, this study attempts to quantify the extent to which agricultural management practices linked to animal production and land use affect output and environmental outcomes (see [31]). Another contribution of this study is that it accounted for spatial heterogeneity between farms, which could result in different relationships between eco-efficiency, self-sufficiency and intensification in each region. The relationships were explored with PLS-SEM. From a policy perspective, this study provides information for the identification of dairy farming systems that can mitigate dairy farm impacts while ensuring food security.

## Materials and Methods

This study builds on the dairy farm exercise by Soteriades et al. [10] for French specialized dairy farms (see above) which combined LCA with DEA to calculate dairy farm eco-efficiency. A brief description of their exercise (data, model and results) is therefore given below before turning to the PLS-SEM framework employed here to study the relationships of eco-efficiency with intensification and self-sufficiency.

### The study of Soteriades et al. [10]

**LCA data and eco-efficiency.** Soteriades et al. [10] used main dairy farming impacts quantified by LCA that is, non-renewable energy use, land use, eutrophication, acidification and global warming potential [6,32]. Note that in LCA studies non-renewable energy use and land use are typically considered as environmental impacts [5,6,8,32] representing 'use of resources' ([32], p.72). Use of resources can be seen both as resources to be used more efficiently and as environmental impacts resulting from the use of non-renewable resources, for
instance CO₂ emissions from combustion of fossil fuels [32]. That is why in this study non-renewable energy use and land use were considered as environmental impacts.

**DEA and eco-efficiency scores for French specialized dairy farms.** DEA was developed by Charnes et al. [11], originating from Farrell's [33] work. It is a non-stochastic, non-parametric technique that benchmarks different decision-making units (DMUs) performing the same task in terms of their capacity to convert inputs into outputs. DEA calculates dimensionless and aggregated efficiency indices without requiring *a priori* assumptions on the importance of each variable for the DMUs' performance, i.e. the variables' weights are obtained from the data themselves. DEA constructs an efficient frontier, that is, a piece-wise linear surface over observed data points against which all DMUs are benchmarked. This frontier comprises of the best performers and the performance of all other DMUs is evaluated by deviations from the frontier line [34]. This is a fundamental difference between DEA and methods such as regression as the latter reflects 'average' or 'central tendency' behaviour [34,35].

The characteristics of DEA make it a particularly attractive tool for the calculation of aggregate eco-efficiency 'scores' as an alternative to multiple eco-efficiency ratios. This is easily demonstrated with the study of Soteriades et al. [10]. These authors fed a DEA model with the five aforementioned dairy farm LCA indicators (non-renewable energy use, land use, eutrophication, acidification, global warming potential) generated by, and three outputs (milk, meat and crop production) produced by, each French specialized dairy farm in the sample. The DEA model (known as the range adjusted measure of inefficiency [36], described in S1 Appendix) then calculated the 'distance' of each farm from the efficient frontier by determining by how much each LCA environmental impact should be reduced, and each output should be increased, for each farm to reach the frontier. These 'distances' are called inefficiencies, as they indicate that a farm is over-generating environmental impacts and/or under-producing outputs relatively to other farms in the sample. The DEA model averaged all inefficiencies for each farm to produce a score between 0 and 1. A score less than 1 meant that a farm was inefficient as it had to eliminate its inefficiencies to reach the frontier. A score of 1 meant that the farm was efficient, as its distance from the frontier- and thus its inefficiencies- were 0. The aggregation was done by the DEA model itself, which was able to weight and sum the inefficiencies altogether. Weighting the inefficiencies cancelled out the variables' different measurement units (e.g. land use in ha, non-renewable energy use in 10⁶ J, milk in kg protein), making the aggregation meaningful. The DEA model calculated the weights from the data themselves and thus no subjective weighting choices were required. Finally, note that the range adjusted measure carries with it a ranking property allowing for the ranking of farms by their eco-efficiency scores (see S1 Appendix).

Using the aforementioned DEA model and LCA and output data, Soteriades et al. [10] calculated eco-efficiency scores for 185 French specialized dairy farms. Note that the eco-efficiency measure did not include operational inputs (e.g. labour, capital, on-farm electricity use etc.) and 'undesirable' outputs (e.g. kg CO₂-equivalents, wastewater etc.) because the idea was to aggregate altogether the two elements used in LCA ratios: environmental impacts and outputs. In other words, Soteriades et al. [10] were concerned with the environmental impacts rather than the amount of operational inputs and undesirable outputs of DMUs (see [16], p.715). This is standard practice in dairy farm LCA exercises [5,6,8]. Because eco-efficiency measures an organization's capability to produce goods and services with minimal environmental impacts [37], Kuosmanen and Kortelainen [37] argue that it represents a societal rather than a managerial perspective. Therefore, they argue that although inputs such as labour and capital are expenditures for the owners of the firm, these inputs represent income (wages and rents) for the society. Thus, such inputs are irrelevant to the context of an eco-efficiency
indicator. See also sub-subsection 'Choice of DEA variables' in Soteriades et al. [10] and Jan et al. [16] p.715.

After obtaining the scores from the whole sample, Soteriades et al. [10] split the sample into two regions, West and Continental France, to examine if the differences in bio-physical conditions between them [38] would favour a specific region in terms of eco-efficiency performance. 'West' was defined by regions with oceanic climate (Basse-Normandie, Bretagne, Haute-Normandie, Pays de la Loire, Poitou-Charentes) and 'Continental' by regions with continental climate (Alsace, Centre, Champagne-Ardennes, Franche-Comté, Lorraine, Rhône-Alpes). Farms on the West were called Oceanic Specialized Systems (OSS, n = 126) as opposed to Continental Specialized Systems (CSS, n = 59) in Continental France. The farms were also split into three feeding strategies (regardless of region), defined by the proportion of maize silage in the total forage area of each farm.

Results showed that OSS systems ranked higher, on average, than CSS in terms of eco-efficiency scores (mean ranks of the DEA scores were 84 for CSS and 97 for OSS). Moreover, farms with <10% and with 10-30% maize in the total forage area ranked, on average, higher than farms with >30% maize in the total forage area in terms of eco-efficiency scores (mean ranks of DEA scores were, respectively, 114, 99 and 78). As commented earlier, the latter results justified a more holistic analysis of the relationships between dairy farm eco-efficiency, intensification and self-sufficiency, leading to the current study. The former results warranted the inclusion of regional heterogeneity effects in the PLS-SEM exercise of this article.

Exploring the relationship of eco-efficiency with intensification and self-sufficiency using PLS-SEM

PLS-SEM is a structural equation modelling (SEM) approach, the latter being a general term for methods used to study the relationships among latent variables indicated by multiple manifest variables [39]. The setting of the PLS-SEM model built in the current study involved four latent variables: for eco-efficiency (ECO), animal-level intensification (INTENS-A), farm-level intensification (INTENS-F) and self-sufficiency (SELF), see Table 1 and Fig 1. ECO was represented by one manifest variable, the DEA eco-efficiency measure of Soteriades et al. [10]. INTENS-A was represented by three manifest variables, milk yield (in kg of raw milk) per cow/year ('milk/cow'); meat production (kg carcass weight produced) per livestock unit (LU) per year ('meat/LU'); and concentrate fed (10³ kg concentrate) per LU per year ('concentrate/LU'). INTENS-F was represented by four manifest variables, kg mineral and organic nitrogen (N) per ha total on-farm area ('N/on-farm ha'); kg mineral and organic phosphorous (P) per ha total on-farm area ('P/on-farm ha'); stocking density (LU/ha main forage area); and the proportion of maize silage in the total forage area ('maize/forage ha').

SELF was represented by four manifest variables of self-sufficiency: 'economic' (gross operating profit/turnover); 'energy' (on-farm energy use/total energy use); 'feed' (feedstuff produced on-farm/total use of feedstuff); and 'land' (on-farm land use/total land use). These four variables captured different dimensions of a farm's capacity to produce goods from its own resources (i.e. to require little purchased inputs) and are complementary, providing a more holistic view of farm self-sufficiency than earlier studies (see [21]). For instance, economic self-sufficiency is an indicator of the degree of economic dependence on external inputs, e.g. pesticides (a similar indicator is used by Lebacq et al. [21]). On the other hand, land self-sufficiency indicates the extent to which a farm depends on 'imported' land but is unable to do so for other inputs such as pesticides and manure. Energy self-sufficiency accounts for, among others, energy required for manure for crop production, with the latter being either imported or recycled on-farm. At the same time, feed self-sufficiency alone fails to consider that farms
<table><caption>Table 1. Statistics for intensification, self-sufficiency and eco-efficiency variables per system.</caption><thead><tr><th rowspan="2"></th><th colspan="3">CSS (n = 59)</th><th colspan="3">OSS (n = 126)</th></tr><tr><th>Mean</th><th>SD</th><th>Median</th><th>Mean</th><th>SD</th><th>Median</th></tr></thead><tbody><tr><td colspan="7"><strong>INTENS-A</strong></td></tr><tr><td>Milk/cow (kg raw milk/cow/year)</td><td>6986</td><td>1560</td><td>7016</td><td>7015</td><td>1196</td><td>7095</td></tr><tr><td>Meat/LU (kg carcass weight produced/LU/year)</td><td>166</td><td>48</td><td>162</td><td>172</td><td>50</td><td>174</td></tr><tr><td>Concentrate/LU (10<sup>3</sup> kg concentrate fed/LU/year)</td><td>1.00</td><td>0.41</td><td>0.98</td><td>0.95</td><td>0.38</td><td>0.98</td></tr><tr><td colspan="7"><strong>INTENS-F</strong></td></tr><tr><td>N/on-farm ha (kg mineral and organic N/ha total on-farm area)</td><td>149</td><td>57</td><td>141</td><td>164</td><td>48</td><td>165</td></tr><tr><td>P/on-farm ha (kg mineral and organic P/ha total on-farm area)</td><td>54</td><td>23</td><td>51</td><td>60</td><td>19</td><td>57</td></tr><tr><td>Stocking density (LU/ha main forage area)</td><td>1.42</td><td>0.54</td><td>1.26</td><td>1.49</td><td>0.34</td><td>1.45</td></tr><tr><td>Maize/forage ha (% maize silage area in the total forage area)</td><td>0.20</td><td>0.20</td><td>0.20</td><td>0.28</td><td>0.16</td><td>0.30</td></tr><tr><td colspan="7"><strong>SELF</strong></td></tr><tr><td>Economic (gross operating profit/turnover)</td><td>0.41</td><td>0.07</td><td>0.40</td><td>0.40</td><td>0.08</td><td>0.40</td></tr><tr><td>Energy (on-farm energy use/total energy use)</td><td>0.60</td><td>0.17</td><td>0.57</td><td>0.63</td><td>0.15</td><td>0.59</td></tr><tr><td>Feed (feedstuff produced on-farm/total use of feedstuff);</td><td>0.90</td><td>0.06</td><td>0.89</td><td>0.90</td><td>0.07</td><td>0.90</td></tr><tr><td>Land (on-farm land use/total land use)</td><td>0.88</td><td>0.08</td><td>0.91</td><td>0.89</td><td>0.07</td><td>0.91</td></tr><tr><td colspan="7"><strong>ECO</strong></td></tr><tr><td>DEA eco-efficiency</td><td>0.94</td><td>0.05</td><td>0.93</td><td>0.95</td><td>0.05</td><td>0.95</td></tr></tbody></table>

CSS: continental specialized systems; OSS: oceanic specialized systems; LU: livestock unit; N: nitrogen; P: phosphorous; DEA: data envelopment analysis. The multiple measures building the latent constructs of animal and farm-level intensification (INTENS-F and INTENS-A) and self-sufficiency (SELF) were obtained from a comprehensive Life Cycle Analysis (LCA) exercise [40] resulting from a partnership involving voluntary participation of farmers ('Inosys Réseaux d'Elevage'), the Chambers of Agriculture (France) and the French Livestock Institute. The LCA exercise is described in [10]. See same reference, text above and S1 Appendix for more details on the DEA eco-efficiency scores building the latent construct of eco-efficiency (ECO).

doi:10.1371/journal.pone.0166445.t001

also produce crops in addition to animal products. From this example it is evident that the four variables agree with the self-sufficiency definition adopted in this study and are complementary.

The dataset used for the PLS-SEM exercise is summarized on Table 1. See S1 File for the whole dataset. Fig 1 is a graphical representation of the PLS-SEM model built in the current study. The model is further explained below.

PLS-SEM comprises of two models (see [41]), the structural model and the measurement model (Fig 1). The structural model assumes a linear relationship between latent variables (also called constructs) and uses linear regression to estimate the path coefficients, representing the strength and direction of the relationships between the response latent variable, or target (endogenous) construct (ECO), and the predictor latent variables, or predictor constructs (INTENS-A, INTENS-F and SELF). In Fig 1 this relationship is represented by single-headed arrows from the predictor latent variables towards the response latent variable. In a similar manner, the measurement model uses linear regression to estimate the loadings, i.e. the correlations between a latent variable and its manifest variables. In Fig 1 this relationship is represented by single-headed arrows from the latent variables towards their manifest variables. The arrows' directions imply that the constructs (INTENS-A, INTENS-F, SELF, ECO) cause the measurement (more precisely, the covariation) of their corresponding manifest variables [30]. The mathematical equations representing the PLS-SEM model in Fig 1 are presented in S2 Appendix.

PLS-SEM was deemed the appropriate SEM method to examine the relationship of ECO with INTENS-A, INTENS-F and SELF for the following reason. The procedure in PLS-SEM
<figure><img src="image_4.png" /><figcaption>**Fig 1. Partial Least Squares Structural Equation Model for eco-efficiency, animal and farm-level intensification, and self-sufficiency.** LU: livestock unit; N: nitrogen; P: phosphorous; DEA: data envelopment analysis; INTENS-A: animal-level intensification; INTENS-F: farm-level intensification; SELF: self-sufficiency; ECO: eco-efficiency.</figcaption></figure>

doi:10.1371/journal.pone.0166445.g001

is an ordinary least squares regression-based estimation and thus is especially useful for developing theories in exploratory research, that is, to search for latent patterns in the data when there is no or only little prior knowledge on how the variables are related (as in this study). It does so by focusing on explaining the variance in the dependent variables (ECO) when examining the model, i.e. it places emphasis on prediction and theory development. This is by contrast with other SEM methods such as covariance-based SEM, which follows a maximum likelihood procedure and is thus better-suited to confirming (or rejecting) theories, i.e. it focuses on causation. See [30,42,43].
## Putting all methods together

Fig 2 summarizes graphically how the different elements of this study's analysis integrate. Soteriades et al. [10] used LCA impacts and farm outputs to produce the DEA eco-efficiency scores for the sample of French specialized dairy farms. The results suggested a possible relationship between eco-efficiency, farming intensity and farm self-sufficiency, which could differ between regions. Further evaluation of this relationship in this study required additional data capturing several aspects of intensity and self-sufficiency. Also, modelling this relationship required a model not based on an *ad hoc* assumption of how the three elements are related, thus PLS-SEM was chosen. The eco-efficiency scores of Soteriades et al. [10] were used to build the manifest variable for *ECO* in the current study. Additional data were used in this study to create manifest variables for *INTENS-A*, *INTENS-F* and *SELF* so as to build the PLS-SEM model. Spatial heterogeneity was accounted for in the PLS-SEM model to examine if regional differences could lead to different findings in West and East France.

All calculations were performed in the R language [44]. The PLS-SEM exercise and evaluation were performed with the R package 'plspm' [41] and, where necessary, the first author's own R functions.

## Results

This section summarizes the PLS-SEM results for CSS and OSS. The PLS-SEM models for both CSS and OSS complied with the minimum sample size requirements (see [30], p.20-22). Before running the exercise, the raw data were standardized (mean = 0, variance = 1) because doing so allowed the PLS-SEM model to calculate standardized coefficients between -1 and +1 for every relationship in the measurement and structural models [30].

Because PLS-SEM does not have a single goodness-of-fit criterion the measurement and structural models had to be evaluated independently and step-by-step using several assessment criteria and guidelines outlined in [30,41,42,45,46]. Given the multitude of criteria and guidelines, it is standard practice in many PLS-SEM studies (e.g. [43,47,48]) to briefly describe criteria, guidelines and their results altogether, and step-by-step. The same practice was adopted here. It should be noted that many of the criteria did not apply for single-item and/or endogenous constructs such as *ECO*, unless explicitly stated below.

The results are presented in four phases. Phase 1 is the measurement model evaluation for CSS and OSS. Phase 2 describes an issue related with the measurement model in Phase 1 and proposes alternative PLS-SEM models for CSS and OSS. Phase 3 is the measurement model evaluation of the new models. Phase 4 is the structural model evaluation of the new models.

### Phase 1: step-by-step evaluation of the measurement models for CSS and OSS

The measurement models were first evaluated for indicator reliability, which requires that the constructs should explain over 50% of their manifest variables' variance [30]. Thus, manifest variables with loadings less than 0.70 should be removed because their construct explains less than $0.70^2 = 0.49 \approx 50\%$ of their variance. The manifest variables meat/LU, P/on-farm ha, stocking density and economic and land self-sufficiency had loadings less than 0.70 for CSS. For OSS, variables with loadings less than 0.70 were meat/LU, N/on-farm ha, P/on-farm ha, stocking density and economic self-sufficiency.

Two things should be noted here. First, leaving the PLS-SEM model for OSS with just one manifest variable (maize/forage ha) representing *INTENS-F* is at odds with conventional measurement theory, according to which constructs should be typically represented by several
<figure><img src="image_5.png" /><figcaption>**Fig 2.** Graphical summary of how the different elements of this study's analysis integrate. LCA: life cycle analysis; DEA: data envelopment analysis; PLS-SEM: partial least squares structural equation modelling; INTENS-A: animal-level intensification; INTENS-F: farm-level intensification; SELF: self-sufficiency; ECO: eco-efficiency.</figcaption></figure>

doi:10.1371/journal.pone.0166445.g002
(reflective) manifest variables (see [49], p.436). (This was not the case with ECO as its single manifest variable was a comprehensive measure of multiple LCA impacts and outputs. Hence, representing ECO by several manifest variables was unnecessary.) For that reason, N/on-farm ha, which had a loading of 0.51, was not deleted from the OSS model. Keeping manifest variables with loadings less than 0.70- but at least 0.50- is usual when situations like the present one require it (see [50], p.198). Second, the results indicated that land self-sufficiency should be retained for OSS but deleted for CSS. Retaining it in OSS would render impossible any comparisons between the PLS-SEM models for CSS and OSS as the manifest variables must be identical among models [51]. Therefore, it was removed from both models.

The final set of variables of the three exogenous constructs was, for both CSS and OSS, milk/cow, concentrate/LU, N/ on-farm ha, maize/forage ha and energy and feed self-sufficiency. The two models were then evaluated for internal consistency reliability, that is, manifest variables in the same construct should be highly correlated since they measure the same construct [30]. The criteria for internal consistency reliability used [52] were Cronbach's *alpha* and Dillon-Goldstein's *rho*, which should be at least 0.70, and the eigenvalues of the correlation matrix, where the first and second eigenvalues should be greater than and smaller than 1 respectively. Both models complied with these criteria, except for OSS with Cronbach's *alpha* equal to 0.65. This value was close to 0.70 and was considered acceptable, given the aforementioned issue with N/on-farm ha but also that the criteria for Dillon-Goldstein's *rho* and the eigenvalues were fulfilled.

The next step was to evaluate the two models in terms of convergent validity, that is, the extent to which a construct converges in its manifest variables by explaining the items' variance [42]. Convergent validity is assessed by the average variance extracted, which equals the mean of the squared loadings for all manifest variables in a construct [42]. An average variance extracted value of 0.50 or higher is acceptable as it indicates that on average, a construct explains over 50% of the variance of its manifest variables [42]. Convergent validity was achieved in both models.

The final step was to evaluate the two models' discriminant validity, that is, the extent to which a construct is truly distinct from other constructs in terms of how much it correlates with them, as well as how much manifest variables represent only a single construct [30]. Discriminant validity was not achieved for either model, because they failed to comply with the HTMT criterion [46]. The HTMT criterion is an estimate of the correlations between two constructs so an absolute value over 0.85 can be interpreted as a violation of discriminant validity [46,53]. All pairwise HTMT correlations between INTENS-A, INTENS-F and SELF had absolute values over 0.85, with their bias-corrected confidence intervals (95%, 5000 samples with replacement of farms within regions; see [30]) containing 1. The only exception was the correlation between INTENS-A and INTENS-F for OSS, with a HTMT value of 0.72 and a bias-corrected confidence interval of (0.55, 0.88). HTMT signs were positive for correlations between INTENS-A and INTENS-F and negative between INTENS-A and SELF and between INTENS-F and SELF.

The HTMT results provided strong evidence that INTENS-A, INTENS-F and SELF were not truly distinct from each other in either model. This required a complete reconsideration of the PLS-SEM model displayed in Fig 1.

## Phase 2: re-considering the PLS-SEM model

The HTMT results suggested that the three exogenous latent variables *INTENS-A*, *INTENS-F* and *SELF* measured 'the same thing'. In fact, the negative HTMT correlations of *SELF* with *INTENS-A* and *INTENS-F* suggested that *SELF* could be better suited as a measure of
<figure><img src="image_6.png" /><figcaption>**Fig 3. Partial Least Squares Structural Equation Model for eco-efficiency and intensification (including inverse-coded self-sufficiency variables).** LU: livestock unit; IC: inverse-coded; DEA: data envelopment analysis; INTENS-ALL: animal and farm-level intensification, as well as inverse-coded self-sufficiency; ECO: eco-efficiency.</figcaption></figure>

doi:10.1371/journal.pone.0166445.g003

intensification; that is, the lower the self-sufficiency levels the higher the intensification levels. Therefore, the PLS-SEM model displayed in Fig 1 was replaced by the PLS-SEM model in Fig 3, where the manifest variables for SELF were inverse-coded to reflect intensification, because that way the higher the self-‘insufficiency’ levels the higher the intensification levels. (Inverse-coding is standard practice in PLS-SEM when facing similar issues. It is done by multiplying the manifest variables of interest by -1, see [30,41].) The new exogenous latent variable, comprising of INTENS-A, INTENS-F and inverse-coded SELF is denoted as INTENS-ALL. It should be made clear though that by no means can this result be generalized to suggest that more intensive farms are less self-sufficient. This result was rather specific to the available dataset and is further discussed in the Discussion section. Therefore, a second PLS-SEM model was considered, where the manifest variables for SELF were completely removed (Fig 4). In this case, the new exogenous latent variable, comprising of INTENS-A and INTENS-F is denoted as INTENS-AF. The two new PLS-SEM models were named accordingly, that is, PLS-SEM-ALL (Fig 3) and PLS-SEM-AF (Fig 4). The mathematical equations representing these two models are presented in S3 Appendix.

PLS-SEM-ALL and PLS-SEM-AF were run for CSS and OSS. The results for the new measurement models for CSS and OSS are reported in Phase 3 below.

## Phase 3: step-by-step evaluation of the new measurement models for CSS and OSS

The final set of manifest variables kept in PLS-SEM-ALL for both CSS and OSS was milk/cow, concentrate/LU, maize/forage ha, and energy and feed self-‘insufficiency’. For PLS-SEM-AF, the final set of manifest variables was milk/cow, concentrate/LU and maize/forage ha for both CSS and OSS. Skewness and kurtosis values for all aforementioned variables, as well as for
<figure><img src="image_7.png" /><figcaption>**Fig 4.** Partial Least Squares Structural Equation Model for eco-efficiency and intensification (self-sufficiency variables removed). LU: livestock unit; DEA: data envelopment analysis; INTENS-AF: animal and farm-level intensification; ECO: eco-efficiency.</figcaption></figure>

doi:10.1371/journal.pone.0166445.g004

DEA eco-efficiency, were between -1 and +1 and thus complied with the requirement that with PLS-SEM data can only slightly depart from normality (this minimizes the chances of obtaining unreliable results, see [30]). The only exception was milk/cow for CSS with kurtosis of -1.29, which was not ‘too far’ from -1 and was not considered a problem [30,54].

The final measurement models complied with all criteria for indicator reliability, internal consistency reliability and convergent validity (S1, S2, S3 and S4 Tables). Bootstrapped confidence intervals (95%, 5000 samples with replacement of farms within regions; see [30]) showed that all the relationships between the manifest variables and their construct were significant (S1, S2, S3 and S4 Tables).

## Phase 4: evaluating the structural models for CSS and OSS

The first step of the structural model assessment was to determine the structural models' predictive accuracy and relevance by means of R², R²-adjusted and Stone-Geisser's cross-validated Q² (blindfolding) [30].

The R² is a measure of the proportion of the endogenous construct's (ECO) variance that is explained by the exogenous constructs (INTENS-ALL or INTENS-AF). For PLS-SEM-ALL the R² values were 0.13 (CSS) and 0.33 (OSS). The respective R²-adjusted values, that is, R² adjusted for sample size to allow for comparisons between models, were 0.11 for CSS and 0.32 for OSS. For PLS-SEM-AF, the R² (R²-adjusted) values were 0.09 (0.08) for CSS and 0.31 (0.30) for OSS.

The Q² is a measure of the structural models' predictive relevance and values above 0 are considered acceptable. For PLS-SEM-ALL Q² ranged between 0.08 and 0.11 for CSS and between 0.29 and 0.32 for OSS, indicating, respectively, small and large predictive relevance for CSS and OSS. For PLS-SEM-AF Q² ranged between 0.05 and 0.08 for CSS and between 0.27 and 0.29 for OSS, again indicating, respectively, small and large predictive relevance for CSS and OSS.

The second and final step was to evaluate the strength and significance of the structural models' path coefficients. The results and their bootstrapped estimates (95%, 5000 samples
<table><caption>Table 2. Path coefficients and their bootstrapped confidence intervals, standard errors, *t*-values and *p*-values for the final structural models for CSS and OSS.</caption><thead><tr><th>Path</th><th>Path coefficient (95% CI)</th><th>SE</th><th>t-value</th><th>p-value</th></tr></thead><tbody><tr><td>INTENS-ALL → ECO (CSS)</td><td>-0.36<sup>a</sup> (-0.55, -0.21)</td><td>0.124</td><td>-2.860</td><td>0.001</td></tr><tr><td>INTENS-ALL → ECO (OSS)</td><td>-0.57<sup>b</sup> (-0.67, -0.47)</td><td>0.074</td><td>-7.770</td><td>&lt; 0.001</td></tr><tr><td>INTENS-AF → ECO (CSS)</td><td>-0.31<sup>a</sup> (-0.51, -0.15)</td><td>0.126</td><td>-2.416</td><td>0.019</td></tr><tr><td>INTENS-AF → ECO (OSS)</td><td>-0.55<sup>b</sup> (-0.65, -0.45)</td><td>0.075</td><td>-7.369</td><td>&lt; 0.001</td></tr></tbody></table>

INTENS-ALL: exogenous construct in model PLS-SEM-ALL, consisting of manifest variables for animal and farm-level intensification, and inverse-coded self-sufficiency; INTENS-AF: exogenous construct in model PLS-SEM-AF, consisting of manifest variables for animal and farm-level intensification; ECO: eco-efficiency, endogenous construct in both PLS-SEM-ALL and PLS-SEM-AF; CSS: continental specialized systems; OSS: oceanic specialized systems; CI: confidence interval; SE: standard error. Values within a column with different superscripts differ significantly at *p* < 0.05. Note that superscripts are not comparable between PLS-SEM-ALL and PLS-SEM-AF.

doi:10.1371/journal.pone.0166445.t002

with replacement of farms within regions; see [30]) are presented in Table 2. The path coefficients INTENS-ALL → ECO and INTENS-AF → ECO were significant and negative for both CSS and OSS. By means of a permutation test (5000 repetitions) [55], significant differences at *p* < 0.05 were found between the path coefficients of CSS and OSS for both PLS-SEM-ALL (*p* = 0.036) and PLS-SEM-AF (*p* = 0.024) (Table 2).

## Discussion

This paper studied the relationship between dairy farm eco-efficiency, self-sufficiency, and farm and animal-level intensification so as to quantify the extent to which agricultural management practices linked to animal production and land use affect output and environmental outcomes at a larger scale. Importantly, this paper accounted for spatial heterogeneity in dairy farms, an increasingly recognized potentially differentiating factor of dairy farm performance [10,18,56–59]. From a policy perspective, this study provides information for the identification of dairy farming systems that can help mitigate dairy farm impacts while ensuring food security.

## Eco-efficiency, intensification and self-sufficiency: relationships, overlaps and recommendations

As noted earlier, a possible confusion could arise from the HTMT results in that they suggested that the higher the (animal and farm-level) intensification, the lower the self-sufficiency. This finding resulted from the available dataset, but it is not always the case. Counterexamples include (i) the ‘average’ New Zealand dairy farm that is intensive per ha on-farm land, however imports less than 10% of total feed offered [5]; and (ii) intensive grazing systems in Belgium with high milk yields per cow and ha [22]. Both cases suggest a synergy between intensification and self-sufficiency, contrary to what was found in this study.

Nevertheless, in the case of this study’s 185 French specialized dairy farms, the negative relationship between intensification and self-sufficiency indicated by the HTMT results is fairly intuitive. Indeed, half the country’s milk production comes from more intensive systems making heavy use of maize silage (especially in the West), which requires supplementation with protein-rich complements, mainly soybean meal imported from Brazil [26,27,60], consequently reducing the farms’ self-sufficiency levels.

The heavy use of maize silage also has implications for the farms’ environmental performance. We confirmed earlier studies [5,10,20] arguing that heavy reliance on maize silage and/or concentrates can reduce dairy farm environmental performance. This is because
importing soybean meal to complement maize silage-based diets of more intensive systems
implies major environmental impacts at the point of production from the expansion of soy
area and forest clearing [26,61]. Note that these impacts were fully accounted for in the LCA
indicators used to derive the eco-efficiency scores thus revealing the true implications of on-
farm management decisions for the environment.

The results also confirmed that although increasing milk yield per cow can improve envi-
ronmental performance at the cow-level [8,9], this does not necessarily hold at whole farm-
level [62,63]. The role of DEA enhanced the validity of this finding because the method does
not express each impact per some functional unit (e.g. unit of milk or land), choice of which
may lead to radically different conclusions [64]; DEA rather aggregates farm impacts and out-
puts altogether, providing a single eco-efficiency measure for the whole farm. It is noteworthy,
however, that two recent DEA studies found that intensification at cow-level actually resulted
in better environmental performance at farm-level, regardless of feed self-sufficiency levels
as reflected by two different feeding practices (lower versus higher reliance on imported con-
centrates) [13,19]. These findings demonstrate the potential for intensive farms to improve
environmental performance [9], perhaps by adoption of highly self-sufficient management
practices, in line with of agro-ecological concepts.

The concept of agro-ecology considers agro-ecosystems as a whole in terms of several
aspects (biological, technical, economic, etc.) and aims at eliminating the disconnectedness of
livestock farming from the land [23]. More self-sufficient systems are more in line with agro-
ecology because they can better regulate biogeochemical cycles and environmental fluxes to
the atmosphere and hydrosphere through interactions among different farm units in space
and time [23]. When stocking densities are kept at lower levels, adequate cropland on-farm
offers a high potential for the recycling of manure-related nutrients [65]. By coupling agro-
ecological concepts with other practices that reduce reliance on external nutrient sources (e.g.
by improving the efficiency of nutrient utilization by the animals; [23]), self-sufficiency then
has the potential for high environmental performance without necessarily lowering a farm's
intensification level [5,22] or, in some cases, by just moderately reducing productivity/ha (see
RAD, 2010 in [23]). Given that in our sample more grass-based farms exhibited higher eco-
efficiency than more maize silage-based farms (see [10]), intensive yet self-sufficient farms
with higher reliance on grass than maize silage may then be preferable because (i) land in
France is relatively cheap [66]; (ii) farms could benefit from agri-environment schemes aimed
at maintaining grasslands, for example France's 'prime à l'herbe' grassland support scheme
[2].

In summary, contrasting results between studies leave less room for a definitive answer as to whether or not intensification is beneficial for the environment; it depends on particular farming systems and circumstances. However, increasing self-sufficiency could offer a way to improve the eco-efficiency of intensive dairy farms, in line with agro-ecological concepts and the consensus that food production should be increased with the least possible environmental impacts [67,68]. In any case, PLS-SEM is very advantageous for exploring relationships between potentially overlapping elements, especially when these relationships are theory-skeletal, which is the case with eco-efficiency, intensification and self-sufficiency. Finally, feeding PLS-SEM with data derived by LCA and DEA enhanced the credibility of the findings.

## The pros and cons of PLS-SEM to study farm sustainability

A wealth of methods has been employed to aggregate different farm sustainability and/or
efficiency indicators and/or to explore the relationships between them. Methods include
bivariate and multivariate statistics, linear programming, DEA, LCA, simulation, and monitoring tools [4,8,9,16,18,69–73]. In this study, PLS-SEM was deemed as the appropriate exploratory tool for two main reasons. First, PLS-SEM allows for the aggregation of manifest variables into latent variables, as opposed to simply feeding the manifest variables altogether into, for example, a regression [9] or principal component analysis [8] model. Second, with PLS-SEM it is possible to simultaneously explore the relationships between more than two latent variables at a time, contrary to bivariate statistics that can only handle pairs of variables. Actually, the first and second aforementioned advantages derive from the fact that PLS-SEM can analyse the whole model as a unit, rather than dividing it into pieces (adopted from [54]).

PLS-SEM combined with LCA and DEA has potential as a guiding tool for the identification of more sustainable dairy farming practices and the formation of environmental regulations. Recent studies argue that environmental externalities are unlikely to decrease solely as a result of market-based instruments and the free market-oriented CAP reforms, hence some policy intervention is necessary [7,74]. Consider, for example, dairy farms in West France. These farms are generally highly industrialized and competitive and have demonstrated a greater ability to respond to CAP's shift from milk support prices to direct payments and towards trade liberalization than farms in other areas [60]. Yet, this study demonstrated that such farms might perform poorly in terms of environmental performance when considering market and non-market goods, as well as local and 'imported' impacts and spatial heterogeneity (recall results in Table 2). Hence, our holistic, 'global' framework could help guide the formation of local environmental regulations that typically ignore 'imported' impacts of farming, potentially resulting in self-sufficient farms experiencing greater enforcement as they tend to generate higher local impacts than less self-sufficient farms.

On the downside, a widely recognized problem with the assessment of PLS-SEM is that, unlike other SEM models, it does not have a standard goodness-of-fit statistic [75]. For instance, there is no universal agreement as to which value of R² is considered acceptable and in some cases a value as low as 0.10 is satisfactory [42]. A general guideline is to interpret R² in the context of the study at hand by considering R² values from related studies [42]. This was impossible to do in this study for the following reason. From the six applications of SEM to dairy farming identified in the literature [47,76–80], only Gyau et al. [47] employed PLS-SEM as the preferred SEM method and the objective of their study had no relevance to that of the current study. Consequently, it was hard to draw any conclusions on the structural model's performance based on the R² values obtained in this study. It should be admitted though that the R² values for CSS were probably too low and a future step would be to develop a better PLS-SEM model for CSS with more data. At least, the positive Q² values indicated the predictive relevance of the structural models, especially for OSS where the values were 'well above zero' ([42], p.111). It is noteworthy that in other studies R² and Q² values as low as for CSS were considered acceptable (e.g. [43]). Other indices to judge the overall model fit in PLS-SEM models have been suggested, such as the global and relative goodness-of-fit indices (see [75]). However, these indices have proven unsuitable for model validation [75] and were not used here.

In summary, the advantages of PLS-SEM as an exploratory tool to study farm sustainability lie on the fact that it can analyse the whole model as a unit, rather than dividing it into pieces. When the PLS-SEM analysis is enriched by holistic LCA and DEA data, the method has potential as a guiding tool for the identification of more sustainable dairy farming practices and the formation of environmental regulations. A disadvantage of PLS-SEM is that it has no standard goodness-of-fit statistic and so its performance can only be assessed with general guidelines and by comparisons with other studies.
## Conclusions

The findings of the current study suggested that whether intensive dairy farming is ‘good’ or ‘bad’ for the environment depends on particular farming systems and circumstances. In this study, on-farm management practices such as increased reliance on maize silage and bought-in concentrates reduced eco-efficiency when the latter was assessed at a ‘global’ level with the LCA-based DEA eco-efficiency scores. The same was true for the effect of increasing milk yield/cow on eco-efficiency because other products (meat and crops) and their associated impacts were also accounted for in the aggregate DEA eco-efficiency scores. However, these results might not necessarily hold for intensive farms with higher self-sufficiency levels and thus better ability to recycle elements on-farm. Intensive yet self-sufficient farms with higher reliance on grass than maize silage may be preferable in cases where they may benefit from relatively low land prices and agri-environment schemes aimed at maintaining grasslands. Finally, our holistic, ‘global’ framework proved to be a valuable tool for quantifying the extent to which agricultural management practices linked to animal production and land use affect output and environmental outcomes, and could help guide the formation of local environmental regulations that typically ignore ‘imported’ impacts of farming.

## Supporting Information

**S1 Appendix. RAM: range adjusted measure.**
(DOCX)

**S2 Appendix. Equations for the PLS-SEM model in Fig 1.**
(DOCX)

**S3 Appendix. Equations for PLS-SEM-ALL and PLS-SEM-AF.**
(DOCX)

**S1 File. The dataset used in this study.** CSS: continental specialized systems; OSS: oceanic specialized systems; INTENS-A: animal-level intensification; INTENS-F: farm-level intensification; SELF: self-sufficiency; ECO: eco-efficiency. Note that the data for ECO were derived by [10].
(XLSX)

**S1 Table. Indicator reliability, internal consistency reliability and convergent validity of the measurement model for PLS-SEM-ALL for Continental Specialized Systems (CSS).**
INTENS-ALL: animal and farm-level intensification, and inverse-coded self-sufficiency indicators; ECO: eco-efficiency; IC: inverse-coded; CI: confidence interval; AVE: average variance extracted; LU: livestock unit; DEA: data envelopment analysis; N/A: not applicable. Rules of thumb [30,42,45,52]: loadings, Cronbach's *alpha* and Dillon-Goldstein's *rho* should be at least 0.70. Communalities (squared loadings) and AVE should be at least 0.50. First and second eigenvalues should be above and below 1 respectively. Note: *ECO* is a single-item construct so assessment criteria above do not apply.
(DOCX)

**S2 Table. Indicator reliability, internal consistency reliability and convergent validity of the measurement model for PLS-SEM-ALL for Oceanic Specialized Systems (OSS).** INTENS-ALL: animal and farm-level intensification, and inverse-coded self-sufficiency indicators; ECO: eco-efficiency; IC: inverse-coded; CI: confidence interval; AVE: average variance extracted; LU: livestock unit; DEA: data envelopment analysis; N/A: not applicable. Rules of thumb [30,42,45,52]: loadings, Cronbach's *alpha* and Dillon-Goldstein's *rho* should be at least 0.70. Communalities (squared loadings) and AVE should be at least 0.50. First and second
eigenvalues should be above and below 1 respectively. Note: *ECO* is a single-item construct so
assessment criteria above do not apply.
(DOCX)

**S3 Table. Indicator reliability, internal consistency reliability and convergent validity of the measurement model for PLS-SEM-AF for Continental Specialized Systems (CSS).**

*INTENS-AF*: animal and farm-level intensification; *ECO*: eco-efficiency; *CI*: confidence interval; *AVE*: average variance extracted; *LU*: livestock unit; *DEA*: data envelopment analysis; *N/A*: not applicable. Rules of thumb [30,42,45,52]: loadings, Cronbach's *alpha* and Dillon-Goldstein's *rho* should be at least 0.70. Communalities (squared loadings) and *AVE* should be at least 0.50. First and second eigenvalues should be above and below 1 respectively. Note: *ECO* is a single-item construct so assessment criteria above do not apply.

(DOCX)

**S4 Table. Indicator reliability, internal consistency reliability and convergent validity of the measurement model for PLS-SEM-AF for Oceanic Specialized Systems (OSS).** *INTENS-AF*: animal and farm-level intensification; *ECO*: eco-efficiency; *CI*: confidence interval; *AVE*: average variance extracted; *LU*: livestock unit; *DEA*: data envelopment analysis; *N/A*: not applicable. Rules of thumb [30,42,45,52]: loadings, Cronbach's *alpha* and Dillon-Goldstein's *rho* should be at least 0.70. Communalities (squared loadings) and *AVE* should be at least 0.50. First and second eigenvalues should be above and below 1 respectively. Note: *ECO* is a single-item construct so assessment criteria above do not apply.

(DOCX)

**Acknowledgments**

A warm thank you to two anonymous reviewers whose constructive comments helped signifi-
cantly improve an earlier version of this manuscript. The authors take responsibility for any
errors that may remain. Special thanks to Professors Marko Sarstedt and Joe F. Hair for their
tips on PLS-SEM, as well as to Prof. David Elston for reviewing an earlier version of this
manuscript.

**Author Contributions**

**Conceptualization:** ADS AWS PF.

**Data curation:** MB SM TC.

**Formal analysis:** ADS.

**Funding acquisition:** AWS PF.

**Methodology:** ADS AWS JL PF.

**Project administration:** AWS PF.

**Resources:** AWS PF.

**Software:** ADS.

**Supervision:** AWS PF.

**Validation:** ADS AWS JL PF.

**Visualization:** ADS.
Writing – original draft: ADS.

Writing – review & editing: ADS AWS JL PF.

## References

1. Foresight. The Future of Food and Farming. London: The Government Office for Science; 2011.

2. Matthews A. Greening agricultural payments in the EU's Common Agricultural Policy. Bio-based and Applied Economics 2013; 2: 1–27.

3. van Berkum S, Helming J. European dairy policy in the years to come: impact of quota abolition on the dairy sector. The Hague: Agricultural Economics Research Institute; 2006.

4. Acosta-Alba I, Lopez-Ridaura S, van der Werf HMG, Leterme P, Corson MS. Exploring sustainable farming scenarios at a regional scale: an application to dairy farms in Brittany. J Clean Prod. 2012; 28: 160–167.

5. Basset-Mens C, Ledgard S, Boyes M. Eco-efficiency of intensification scenarios for milk production in New Zealand. Ecol Econ 2009; 68: 1615–1625.

6. Guerci M, Bava L, Zucali M, Sandrucci A, Penati C, Tamburini A. Effect of farming strategies on environmental impact of intensive dairy farms in Italy. J Dairy Res. 2013; 80: 300–308. doi: 10.1017/S0022029913000277 PMID: 23806128

7. Santarossa JM, Stott AW, Woolliams JA, Brotherstone S, Wall E, Coffey MP. An economic evaluation of long-term sustainability in the dairy sector. Anim Sci. 2004; 79: 315–325.

8. Bava L, Sandrucci A, Zucali M, Guerci M, Tamburini A. How can farming intensification affect the environmental impact of milk production? J Dairy Sci. 2014; 97: 4579–4593. doi: 10.3168/jds.2013-7530 PMID: 24792806

9. Casey JW, Holden NM. The relationship between greenhouse gas emissions and the intensity of milk production in Ireland. J Environ Qual. 2005; 34: 429–436. PMID: 15758094

10. Soteriades AD, Faverdin P, Moreau S, Charroin T, Blanchard M, Stott AW. An approach to holistically assess (dairy) farm eco-efficiency by combining life cycle analysis with data envelopment analysis models and methodologies. Animal 2016; 10: 1899–1910. doi: 10.1017/S1751731116000707 PMID: 27126890

11. Charnes A, Cooper WW, Rhodes E. Measuring the efficiency of decision-making units. Eur J Oper Res. 1978; 2: 429–444.

12. Berre D, Boussemart JP, Leleu H, Tillard E. Economic value of greenhouse gases and nitrogen surpluses: Society vs farmers’ valuation. Eur J Oper Res. 2013; 226: 325–331.

13. Soteriades AD, Faverdin P, March M, Stott AW. Improving efficiency assessments using additive data envelopment analysis models: an application to contrasting dairy farming systems. Agr Food Sci. 2015; 24: 235–248.

14. Berre D, Blancard S, Boussemart JP, Leleu H, Tillard E. Finding the right compromise between productivity and environmental efficiency on high input tropical dairy farms: a case study. J Environ Manage. 2014; 146: 235–244. doi: 10.1016/j.jenvman.2014.07.008 PMID: 25178529

15. Iribarren D, Hospido A, Moreira MT, Feijoo G. Benchmarking environmental and operational parameters through eco-efficiency criteria for dairy farms. Sci Total Environ 2011; 409: 1786–1798. doi: 10.1016/j.scitotenv.2011.02.013 PMID: 21371738

16. Jan P, Dux D, Lips M, Alig M, Dumondel M. On the link between economic and environmental performance of Swiss dairy farms of the alpine area. Int J Life Cycle Ass. 2012; 17: 706–719.

17. Pérez Urdiales M, Oude Lansink A, Wall A. Eco-efficiency among dairy farmers: the importance of socio-economic characteristics and farmer attitudes. Environ Resour Econ. 2016; 64: 559–574.

18. Shortall OK, Barnes AP. Greenhouse gas emissions and the technical efficiency of dairy farmers. Ecol Indic. 2013; 29: 478–488.

19. Toma L, March M, Stott AW, Roberts DJ. Environmental efficiency of alternative dairy systems: a productive efficiency approach. J Dairy Sci. 2013; 96: 7014–7031. doi: 10.3168/jds.2013-6911 PMID: 24054279

20. Arsenault N, Tyedmers P, Fredeen A. Comparing the environmental impacts of pasture-based and confinement-based dairy systems in Nova Scotia (Canada) using life cycle assessment. Int J Agr Sustain. 2009; 7: 19–41.

21. Lebacq T, Baret PV, Stilmant D. Role of input self-sufficiency in the economic and environmental sustainability of specialised dairy farms. Animal 2015; 9: 544–552. doi: 10.1017/S1751731114002845 PMID: 25431140
22. Meul M, Van Passel S, Nevens F, Dessein J, Rogge E, Mulier A, et al. MOTIFS: a monitoring tool for integrated farm sustainability. Agron Sustain Dev. 2008; 28: 321–332.

23. Dumont B, Fortun-Lamothe L, Jouven M, Thomas M, Tichit M. Prospects from agroecology and industrial ecology for animal production in the 21st century. Animal 2013; 7: 1028–1043. doi: 10.1017/S1751731112002418 PMID: 23257276

24. Maulfair D, Heinrichs J, Ishler V. Feed efficiency for lactating dairy cows and its relationship to income over feed costs. Available: http://extension.psu.edu/animals/dairy/nutrition/nutrition-and-feeding/diet-formulation-and-evaluation/feed-efficiency-in-lactating-dairy-cows-and-its-relationship-to-income-over-feed-costs; 2011.

25. Acosta-Alba I, Corson MS, van der Werf HMG, Leterme P. Using reference values to assess environmental sustainability of dairy farms. Renew Agr Food Syst. 2012; 27: 217–227.

26. Lehuger S, Gabrielle B, Gagnaire N. Environmental impact of the substitution of imported soybean meal with locally-produced rapeseed meal in dairy cow feed. J Clean Prod. 2009; 17: 616–624.

27. Nguyen TTH, Corson MS, Doreau M, Eugene M, van der Werf HMG. Consequential LCA of switching from maize silage-based to grass-based dairy systems. Int J Life Cycle Ass. 2013; 18: 1470–1484.

28. Lohmöller J-B. Latent variable path modeling with partial least squares. Heidelberg: Physica-Verlag; 1989.

29. Wold H. Path models with latent variables: the NIPALS approach. In: Blalock HM, Aganbegian A, Bordodkin FM, Boudon R, Capecchi V, editors. New York: Academic; 1975.

30. Hair JF, Hult GT, Ringle CM, Sarstedt M. A primer on partial least squares structural equation modelling (PLS-SEM). SAGE Publications Inc; 2014.

31. Tomich TP, Brodt S, Ferris H, Galt R, Horwath WR, Kebreab E, et al. Agroecology: a review from a Global-Change perspective. Annu Rev Env Resour. 2011; 36: 193–222.

32. de Boer IJM. Environmental impact assessment of conventional and organic milk production. Livest Prod Sci. 2003. 80: 69–77.

33. Farrell MJ. The measurement of productive efficiency. J Roy Stat Soc A. 1957; 120: 253–290.

34. Cooper WW, Seiford LM, Tone K. Data envelopment analysis: a comprehensive text with models, applications, references and DEA-Solver software. 2nd ed. Springer Science+Business Media, LLC; 2007.

35. Cooper WW. Origins, Uses of, and relations between goal programming and data envelopment analysis. Journal of Multi-Criteria Decision Analysis 2005; 13: 3–11.

36. Cooper WW, Park KS, Pastor JT. RAM: a range adjusted measure of inefficiency for use with additive models, and relations to other models and measures in DEA. J Prod Anal. 1999; 11: 5–42.

37. Kuosmanen T, Kortelainen M. Measuring eco-efficiency of production with data envelopment analysis. Journal of Industrial Ecology 2005; 9: 59–72.

38. Gac A, Manneville V, Raison C, Charroin T, Ferrand M. L'empreinte carbone des élevages d'herbivores: présentation de la méthodologie d'évaluation appliquée à des élevages spécialisés lait et viande. Renc Rech Rum. 2010; 17: 335–342.

39. Lei P-W, Wu Q. Introduction to structural equation modeling: issues and practical considerations. Educ Meas. 2007; 26: 33–43.

40. Charroin T, Palazon R, Madeline Y, Guillaumin A, Tchakerian E. Le système d'information des Réseaux d'Elevage français sur l'approche globale de l'exploitation. Intérêt et enjeux dans une perspective de prise en compte de la durabilité. Renc Rech Rum. 2005; 12: 335–338

41. Sanchez G. PLS path modeling with R. Berkeley: Trowchez Editions; 2013.

42. Sarstedt M, Ringle CM, Smith D, Reams R, Hair JF. Partial least squares structural equation modeling (PLS-SEM): a useful tool for family business researchers. Journal of Family Business Strategy 2014; 5: 105–115.

43. Rodriguez-Entrena M, Salazar-Ordonez M, Sayadi S. Applying partial least squares to model genetically modified food purchase intentions in southern Spain consumers. Food Policy 2013; 40: 44–53.

44. R Core Team. R: a language and environment for statistical computing. Version 3.1.3. Vienna, Austria: R Foundation for Statistical Computing; 2015. Available: http://www.r-project.org

45. Hair JF, Sarstedt M, Ringle CM, Mena JA. An assessment of the use of partial least squares structural equation modeling in marketing research. J Acad Market Sci. 2012; 40: 414–433.

46. Henseler J, Ringle CM, Sarstedt M. A new criterion for assessing discriminant validity in variance-based structural equation modeling. J Acad Market Sci. 2015; 43: 115–135.

47. Gyau A, Spiller A, Wocken C. Price or relational behaviours? Supplier relationship management in the German dairy industry. Brit Food J. 2011; 113: 838–852.
48. Ringle CM, Sarstedt M, Zimmermann L. Customer satisfaction with commercial airlines: the role of perceived safety and purpose of travel. Journal of Marketing Theory and Practice 2011; 19: 459-472.

49. Diamantopoulos A, Sarstedt M, Fuchs C, Wilczynski P, Kaiser S. Guidelines for choosing between multi-item and single-item scales for construct measurement: a predictive validity perspective. J Acad Market Sci. 2012; 40: 434-449.

50. Hulland J. Use of partial least squares (PLS) in strategic management research: a review of four recent studies. Strategic Manage J. 1999; 20: 195-204.

51. Henseler J, Ringle CM, Sarstedt M. Testing measurement invariance of composites using partial least squares. Int Market Rev. 2016; 33: 405-431.

52. Tenenhaus M, Vinzi VE, Chatelin YM, Lauro C. PLS path modeling. Comput Stat Data An. 2005; 48: 159-205.

53. Voorhees CM, Brady MK, Calantone R, Ramirez E. Discriminant validity testing in marketing: an analysis, causes for concern, and proposed remedies. J Acad Market Sci. 2016; 44: 119-134.

54. Goodhue DL, Lewis W, Thompson R. Does PLS have advantages for small sample size or non-normal data? MIS Quarterly 2012; 36: 981-1001.

55. Chin WW, Dibbern J. An introduction to a permutation based procedure for multi-group PLS analysis: results of tests of differences on simulated data and a cross cultural analysis of the sourcing of information system services between Germany and the USA. In: Vinzi VE, Chin WW, Henseler J, Wang H, editors. Handbook of Partial Least Squares Concepts, Methods and Applications. Berlin, Heidelberg: Springer-Verlag; 2010. pp. 171-193.

56. Barnes AP. Does multi-functionality affect technical efficiency? A non-parametric analysis of the Scottish dairy industry. J Environ Manage. 2006; 80: 287-294. doi: 10.1016/j.jenvman.2005.09.020 PMID: 16549233

57. Hoch I. Estimation of production function parameters combining time-series and cross-section data. Econometrica 1962; 30: 34-53.

58. Novo AM, Slingerland M, Jansen K, Kanellopoulos A, Giller KE. Feasibility and competitiveness of intensive smallholder dairy farming in Brazil in comparison with soya and sugarcane: Case study of the Balde Cheio Programme. Agr Syst. 2013; 121: 63-72.

59. Olesen JE, Schelde K, Weiske A, Weisbjerg MR, Asman WAH, Djurhuus J. Modelling greenhouse gas emissions from European conventional and organic dairy farms. Agr Ecosyst Environ. 2006; 112: 207-220.

60. Ben Arfa N, Rodriguez C, Daniel K, Shonkwiler JS. Spatial structure of French dairy sector: a spatial HAC estimation. In: 3èmes journées de recherches en sciences sociales, INRA-SFER-CIRAD, 2009 Dec 9-Dec 11, Montpellier, France. Available: http://www.sfer.asso.fr/content/download/3014/27370/file/

61. Bickel U, Dros JM. The impacts of soybean cultivation on Brazilian ecosystems. Three case studies. Amsterdam: WWF Forest Conversion Initiative; 2003.

62. Puillet L, Agabriel J, Peyraud JL, Faverdin P. Modelling cattle population as lifetime trajectories driven by management options: a way to better integrate beef and milk production in emissions assessment. Livest Sci. 2014; 165: 167-180.

63. Zehetmeier M, Baudracco J, Hoffmann H, Heissenhuber A. Does increasing milk yield per cow reduce greenhouse gas emissions? A system approach. Animal 2012; 6: 154-66. doi: 10.1017/S1751731111001467 PMID: 22436163

64. Salou T, Le Mouël C, van der Werf HMG. Environmental impacts of dairy system intensification: the functional unit matters! J Clean Prod. (in press). doi: 10.1016/j.jclepro.2016.05.019

65. Saam H, Powell JM, Jackson-Smith DB, Bland WL, Posner JL. Use of animal density to estimate manure nutrient recycling ability of Wisconsin dairy farms. Agr Syst. 2005; 84: 343-357.

66. Chatellier V, Pflimlin A. Dairy systems in the European regions of the Atlantic area. In: Proceedings of the final seminar of the Green Dairy Project, 2006 Dec 13-Dec 14 2006, Rennes, France. Available: http://ciam.gal/pdf/Green_Dairy.pdf

67. Godfray HCJ, Garnett T. Food security and sustainable intensification. Phil T Roy Soc B 2014. doi: 10.1098/rstb.2012.0273 PMID: 24535385

68. Eisler MC, Lee MRF, Tarlton JF, Martin GB, Beddington J, Dungait JAJ, et al. Steps to sustainable livestock. Nature 2014; 507: 32-34. PMID: 24605375

69. Reig-Martinez E, Gomez-Limon JA, Picazo-Tadeo AJ. Ranking farms with a composite indicator of sustainability. Agr Econ. 2011; 42: 561-575.

70. van Calker KJ, Berentsen PBM, Giesen GWJ, Huirne RBM. Maximising sustainability of Dutch dairy farming systems for different stakeholders: a modelling approach. Ecol Econ. 2008; 65: 407-419.
71. Del Prado A, Scholefield D. Use of SIMS(DAIRY) modelling framework system to compare the scope on the sustainability of a dairy farm of animal and plant genetic-based improvements with management-based changes. J Agr Sci. 2008; 146: 195-211.

72. Meul M, Van Passel S, Nevens F, Dessein J, Rogge E, Mulier A, et al. MOTIFS: a monitoring tool for integrated farm sustainability. Agron Sustain Dev. 2008; 28: 321-332.

73. Barnes AP, Thomson SG. Measuring progress towards sustainable intensification: how far can secondary data go? Ecol Indic. 2014; 36: 213-220.

74. Zimmermann A. Optimization of sustainable dairy-cow feeding systems with an economic-ecological LP farm model using various optimization processes. J Sustain Agr. 2008; 32: 77-94.

75. Henseler J, Sarstedt M. Goodness-of-fit indices for partial least squares path modeling. Computation Stat. 2013; 28: 565-580.

76. de los Campos G, Gianola D, Boettcher P, Moroni P. A structural equation model for describing relationships between somatic cell score and milk yield in dairy goats. J Anim Sci. 2006; 84: 2934-2941. doi: 10.2527/jas.2006-016 PMID: 17032786

77. Detilleux J, Theron L, Beduin JM, Hanzen C. A structural equation model to evaluate direct and indirect factors associated with a latent measure of mastitis in Belgian dairy herds. Prev Vet Med. 2012; 107: 170-179. doi: 10.1016/j.prevetmed.2012.06.005 PMID: 22770804

78. Hansson H, Ferguson R. Factors influencing the strategic decision to further develop dairy production— a study of farmers in central Sweden. Livest Sci. 2011; 135: 110-123.

79. Ondersteijn CJM, Giesen GWJ, Huirne RBM. Identification of farmer characteristics and farm strategies explaining changes in environmental management and environmental and economic performance of dairy farms. Agr Syst. 2003; 78: 31-55.

80. Richards TJ, Jeffrey SR. Efficiency and economic performance: an application of the MIMIC model. J Agr Resour Econ. 2000; 25: 232-51.