Implications of within-farm transmission for network dynamics: Consequences for the spread of avian influenza

Highlights • Cross-scale dynamics were investigated for avian influenza in British poultry.• Transmission risk is dependent on the assumed within-flock transmission mode.• Transmission risk may not scale with transmissibility or flock size.• Transmission risk corresponds with between-farm impact for 28% of farms.• These results have implications for targeted disease control at the farm-level.


Introduction
In social network analysis, when applied to epidemiology, it has typically been assumed that links between network nodes are static and therefore represent fixed, persistent contacts during an individual's infectious period (Huerta and Tsimring, 2002;Keeling, 1999;Meyers et al., 2005). Comparatively few models have allowed links to be dynamic -reflecting that, in reality, an individual's contacts are likely to change over time (Eames and Keeling, 2004; Kao the network dynamics (Bajardi et al., 2011;Vernon and Keeling, 2009). The impact of "super-shedding" cattle on persistence provides one example where the role of within-herd (node-level) dynamics in population transmission has been studied (Liu et al., 2007). In general, however, relatively little attention has been given to the interaction between node-level dynamics and network links.
Whether it is important to consider these factors in unison depends on the relative rates of change at the farm and network levels (Kao et al., 2007;Ochab and Gora, 2011;Volz and Meyers, 2009). For diseases such as highly pathogenic avian influenza (HPAI), which spreads rapidly at the farm level (i.e. node dynamics) (Bos et al., 2010;Elbers et al., 2004;Yoon et al., 2005), the opportunity for onward spread via epidemiologically relevant industry movements will depend on the timing of these movements.
HPAI is a public health threat with the potential to cause large economic losses; therefore, all countries require good contingency plans to limit the impact of outbreaks in commercial poultry. As the control of HPAI hinges on rapid detection and notification, improved understanding of within-flock transmission dynamics could help design measures to reduce transmission from infected premises. These dynamics are known to be influenced by management factors, such as flock size, and influence the time-to-detection of an outbreak (Savill et al., 2008). Although the implications of the within-farm transmission of HPAI for indirect farm-to-farm transmission via movements of people, vehicles and equipment has been considered (Dorea et al., 2010), previous studies have not considered in detail the dynamics of within-flock transmission (such as the impact of flock size and assumed transmission mode) or the impact of the timing of on-farm movement activities.
Catching team personnel (through clothing and equipment), forklift trucks, and slaughterhouse vehicles are considered to be important risk factors for the spread of disease within the British commercial poultry industry (Anon, 2006;Gittins and Canning, 2006). Using temporally explicit data describing the daily on-farm visit schedules for a major poultry catching company located in England (Dent et al., 2011), we explored the interaction between the within-flock transmission dynamics of HPAI, measured in continuous time, and temporally explicit catching-team visits, measured in discrete time. These on-farm visits made by catching teams represent the risk of fomites (i.e. people, vehicles or equipment) being exposed to HPAI and subsequently spreading infection through their movements, and not the movement of the birds themselves.
We present an explicit exploration of the interaction across these time-scales by combining a deterministic model of disease spread (given the typically large flock sizes for commercial poultry) with the empirical data describing on-farm catching-team visits. As these on-farm visits represent 'impulses' enabling the spread of HPAI between-farms, we gain insights that have implications for the dynamics of network links between premises. Our results highlight the impact of the characteristics of bird-to-bird transmission on the potential for spread from infected premises. Additionally, we identify the farm-level factors that are predictive of high relative transmission risk (see below) and discuss how this can be used to inform the risk of HPAI spread between farms.

Within-flock transmission
A S-E-I-R model was used to track the number of birds in susceptible (S), exposed (E), infectious (I) and removed (R) classes. The 'R' class represents the end of a birds' infectious period as a result of HPAI-induced mortality and therefore these birds no longer contribute to the force of infection within the flock. In addition, a 'F' class was incorporated to represent the environmental buildup of infectious faecal material within the poultry house. This 'F' class determined the potential for the exposure of fomites such as catching team personnel, forklift trucks, slaughterhouse vehicles and equipment. In order to gain a better understanding of the relationships amongst several flock-level factors under different transmission mode assumptions, we captured average effects using a deterministic framework and thus do not account for stochastic variability in the within-farm dynamics. However, we do consider uncertainty and between farm variability in the parameter values of this model (see below).
The within-flock model was developed to reflect floor-reared broiler poultry, which was the predominant production type captured by our catching-company movement data. These data reflect the daily schedules of the catching teams that are responsible for transporting poultry to the slaughterhouse vehicle; in the case of broiler chicken production, this either relates to flock thinning, partial bird depopulation or total bird depopulation at the end of the production cycle.
Each bird was assumed to excrete ε grams of faeces per hour; εI therefore gave the quantity of infectious faeces excreted by the flock each hour. The amount of infectious faecal material was assumed to decay at 5% per hour, consistent with available experimental data (Shortridge et al., 1998) and a published within-flock model for HPAI in poultry (Savill et al., 2006(Savill et al., , 2008. As the relative contribution of different transmission mechanisms is poorly understood (Spekreijse et al., 2011a), the total infection pressure, ˇ, was modelled as the sum of the infection pressure due to two independent routes; namely transmission via: (i) aerosol (ˇa), which represented direct bird-to-bird contact (i.e. dependent on 'I'); and (ii) via infectious faeces (ˇf), which represented indirect contact via faecal contamination of dust and drinking/feeding equipment (i.e. dependent on 'F'). As the mechanisms of bird-to-bird contact (direct or indirect) and the likely rates of contact are poorly understood, both frequency-dependent (FD) and density-dependent (DD) modes of transmission were considered. For FD transmission the ODEs are given by The expressions for DD transmission are identical except that the force of infection is replaced witȟ where the transmission rates are scaled by the mean flock size (n) to facilitate comparison with the FD models. Table 1 contains a full list of parameter values and ranges.
We chose a range of values for the basic reproduction number, R 0 , from 2 to 38 (corresponding to ˇ values ranging 0.01-10 h −1 ) to explore a spectrum of scenarios that are covered by experimental and modelling studies. For example, estimates of R 0 are as low as 1.2 and 2.2 in experimental transmission studies and field data respectively (Bouma et al., 2009;Spekreijse et al., 2011b;Tiensin et al., 2007), while estimates as high as 22-66 have been assumed in mathematical models due to the greater flock size and infection pressure expected within a commercial flock (Savill et al., 2006;Sharkey et al., 2008;Truscott et al., 2007). Supplementary Material Section 1 includes details of model sensitivity to the R 0 range. R 0 was calculated as the dominant eigenvalue of the Next Gen- and where x a = ˇa/N and x f = ˇf/N for FD transmission, and x a = ˇa/n and x f = ˇf /n for DD transmission. The average flock size per farm (N, the ratio of number of poultry to number of poultry houses) was divided into quintiles, using the distribution obtained from the Catching Company Database (mean = 22,465 birds, range = 4571-45,667 birds). Twenty values each of ˇa and ˇf were chosen to cover the possible range likely on farms (0.01-10 h −1 ). All flock size and ˇ combinations (400 pairwise combinations of ˇa and ˇf) were run for both FD and DD transmission modes and for each quintile. In total, 4000 outbreaks (400 ˇ combinations × 5 flock sizes × 2 transmission scenarios) were simulated using Matlab v.7.8.0 (The MathWorks, Inc., Natick, MA, USA). Each simulation represented transmission within a single flock and one single infected bird initiated each outbreak at time zero.

Relative transmission risk
Model simulations were used to explore the effect of withinflock transmission characteristics on a farms' relative transmission risk (TR) (i.e. relative across farms). Empirical data was obtained from a major catching company in England which provided the daily schedules of catching-teams recorded on an hourly basis. This Catching Company Database represented 68 catching teams visiting a total of 415 farms over a total period of 950 days (see Dent et al., 2011 and Supplementary Material Section 2). Each farms' recorded catching-team visit days were given a value determined by the total number of slaughterhouse vehicle loads (ϕ) associated. The TR of an infected farm was then computed by matching each of their daily on-farm visits to the corresponding day of all the simulated within-flock outbreaks that matched the relevant flock-size.
Specifically, for farm i given a specified day of incursion z (the day within the movement data that was identified as day zero of the simulated outbreak), the relative transmission risk TR i,z was given by the product of the amount of infectious faecal material (F i,z,t ) and the number of vehicle loads (ϕ i,z,t ), summed over each day, t, of the outbreak: To account for the intra-farm heterogeneity in movement pattern (i.e. both frequency of catching-team visits and bird transportation capacity of slaughterhouse vehicles), iterations over all possible incursion days (z) were run. As a result of this intra-farm heterogeneity in the movement patterns, the number of incursion day iterations varied per farm. An overall relative transmission risk at the farm-level, TR i,sim , was calculated by averaging across all incursion day iterations (z): where n i represents the total number of iterations (z) corresponding to farm i. We consider two types of infection incursion. First, it was assumed that an incursion event would occur only on movement days, corresponding to farm-to-farm transmission via a catching team whose previous visit was to an infectious premises. Second, incursions could occur on any day thereby representing other sources of introductions such as a wildlife reservoir (for details see Supplementary Material Section 3). It was assumed that once a threshold-level of dead birds (known as the mortality threshold, MT) had been reached, the outbreak would be detected and any further risk of transmission prevented. Following an outbreak of HPAI H7N7 in the Netherlands in 2003, a MT of 0.5% (of the initial susceptible flock) per day for two consecutive days was recommended for Dutch broiler producers . However, MTs are recognised to vary according to several factors such as production type, management practices and bird age and a wider likely range of 0.03-3.33% was found for poultry producers in Georgia, USA (Vieira et al., 2009). As the likely MT triggering HPAI detection in British poultry farms is not known we considered a range of MTs and present results for an intermediate threshold of 0.5%, corresponding to the Dutch recommendation.

Identifying predictors of high relative transmission risk
To identify factors that have the greatest effect on the probability of transmission from a farm, we produced a general linear regression model to compare the model inputs with the predicted relative transmission risk. To remove the effect of uncertainty in transmission parameters, the farm-level transmission risk (TR i ) corresponding to a mid-range transmissibility scenario (ˇa + ˇf ∼ 10) was chosen. We wished to identify the most important drivers of transmission risk, therefore only main effects in The median dispersal index (variance-to-mean ratio) for the farm-level time-interval between consecutive visit days was 26 days (range: 11-348 days). b Based on 100 records for which integration status was known. c 75.9%, 12.96% and 9.26% of farms were located within the North and West of England and Wales, respectively; note that results for the East of England were excluded due to an extremely small sample size (1.85% of farms). the statistical model were considered. For full details of predictor variables see Table 2.
A square root transformation was applied to TR i as this improved the normality of residuals and decreased heteroscedasticity. The models were built using a backwards stepwise method and AIC to assess the fit at each stage in the model development. The inclusion of 'company integration' as a random effect was considered but not included in the final models due to the few companies with an adequate group size (only two companies had more than two associated farm premises).
The most influential data points, as identified by their Cook's statistic, were assessed for their impact on the model coefficients and their significance levels. To assess the overall effect of covariates in the model, the relative transmission risk was predicted using the linear model with all parameters set at the median (covariates) or modal (factors) value. The effect of varying a single covariate or factor across the range observed within the dataset was then computed. Those covariates or factors with the greatest effect were recorded. All statistical analyses were carried out in R software v.2.14.2.

Correlation between relative transmission risk and network connectivity
To assess the potential impact on infection propagation between farms, TR i was cross-classified with a measure of between-farm association frequency, informed by the Poultry Network Database (PND). This database consisted of surveys of: (i) single-site and (ii) multi-site farm premises, (iii) slaughterhouses and (iv) catching companies; these data were used to infer potential associations between farms arising through shared industry associations. For example, farms that used the same slaughterhouse, catching company or that were integrated as part of a larger company were assumed to be epidemiologically linked (Dent et al., 2008;Nickbakhsh et al., 2011).
Using these associations informed by the PND, a between-farm association matrix was generated for the farms captured by the Catching Company Database that were used in the computation of TR i (n = 108). The total number of associations per farm was assumed to represent a maximum potential for propagation via the poultry network. In this way, correlation between the farmlevel risk of exposing personnel, vehicles and equipment during farm visits (i.e. as measured through TR i ), and the potential risk of propagation occurring through the industry network, was determined.

Movement data
A descriptive analysis of the input data showed that the farmlevel mean number of caught birds and average flock size both increased linearly with the mean daily number of vehicle loads (Fig. 1a). In contrast, the mean total number of catching days (across the entire dataset) and flock size decreased with the mean timeinterval between consecutive visit days, T b (Fig. 1b). Overall, larger flocks tended to have a greater number of daily vehicle loads, a greater total number of visits, and a lower T b . Further inspection of T b showed it to have a variance-to-mean ratio exceeding 1, indicating temporal clustering of visits. This temporal clustering at the farm-level reflects two important activities related to these catching teams -flock thinning mid-production cycle and partial depopulation at the end of the production cycle (see Supplementary Material Section 3 for an example data snapshot highlighting this temporal pattern of catching team visits). For FD transmission, the time-to-detection increased with increasing flock size (solid red arrow compared to dashed red arrow, Fig. 2a), as previously highlighted (Savill et al., 2008). However, for DD transmission, time-to-detection decreased with increasing flock size (solid red arrow compared to dashed red arrow, Fig. 2b). In this case, greater infection pressure and greater movement activity, but more rapid outbreak detection, interacts to influence the potential for transmission.

Cross-scale interactions at the farm-level
The model simulations, combining the deterministic withinflock transmission model for HPAI with the catching-team movement data, showed that for a given flock size, transmission mode and mortality threshold, the average TR (across all farms) increased with transmissibility (approximately scaling with the build-up of infectious faeces), peaked at high transmissibilities (R 0 ∼25-30) and then dropped sharply for even higher transmissibilities -a consequence of the more rapid accumulation of dead birds. This occurred even for relatively small increases to ˇa oř up to a level which depended on the assumed transmission mode (discussed further below).
For frequency-dependent (FD) transmission, an increase to flock size increased TR ( Fig. 3a and b). However, for density-dependent (DD) transmission, although TR scaled with flock size for low-midrange transmission rates (i.e. ˇa + ˇf < 10), for mid-to-high range transmission rates (i.e. ˇa + ˇf > 10) earlier disease detection caused TR to peak at lower transmission rates as flock size increased. This narrowing of the region of parameter space between the lowest and highest TR values with an increase to flock size is illustrated by the relative sizes of the white arrows in Fig. 3c and d. As the number of dead birds approached the mortality threshold, TR became highly sensitive to the balance between infectious faeces and timeto-detection under DD transmission (see bottom right red region of Fig. 3d); the greater environmental build-up of infectious faecal material counteracted the shortened time-to-detection for small regions of parameter space causing TR to oscillate as transmission parameters varied.
The results presented assume a mortality threshold of 0.5%. Given the likely variability and lack of recent experience of HPAI by farmers in Great Britain, we compare this to MTs of 0.3% and 0.7% (Supplementary Material Section 4). We note that the range of infection pressures corresponding with high TR increases with increasing MT, whilst the rapid drop in TR phenomenon occurs later.   (c) and (d) demonstrate the relative difference in the region of parameter space between the lowest and highest TR as flock size increases under DD transmission. These analyses assumed outbreak detection occurred when at least 0.5% mortality was reached for two consecutive days. All other parameters were set to default values (see Table 2).
Section 5 for details of the multivariable model diagnostics. Overall, the time between consecutive catching-team visit days, T b , had the largest average effect on TR i for both transmission modes (linear coefficients = −0.033, p ≤ 0.008; for both model scenarios). However, when considering the range of possible effect size across all farms for each predictor in the multivariable model, 'mean daily vehicle loads' had the potential to have the largest impact under FD transmission (see Supplementary Material Section 6).
The relationship of TR i with flock size varied depending on the transmission mode, as shown in Fig. 4; a significant effect was found only in the univariable model under FD transmission (results not shown, linear coefficient = 4 × 10 −5 , 95% CI: 2 × 10 −5 to 6 × 10 −5 , p < 0.0001), whilst under DD transmission the average flock size had a significant negative impact in the multivariable model, an effect which was further modified by 'mean daily vehicle loads' (interaction parameter between 'flock size' and 'mean daily vehicle loads' = −1.21 × 10 −7 , p < 0.0001). Therefore, for DD transmission, the relationship with 'mean daily vehicle loads' was seen most clearly for small and medium sized flocks, with TR peaking at midrange flock sizes (∼25,000-30,000 birds).  Table 2 for variable definitions and value ranges. b Results obtained for square root transformed values of TR i . c Numeric variables were centred and therefore the intercept corresponds to their average values. d Between-farm association frequency was estimated as a measure of between-farm network connectivity. e West and North refer to geographical regions of England. Predictors in bold had the largest possible range of effect sizes across all farms -see Supplementary Material Section 6. f The most influential data points (based on their Cook's statistic) were identified and removed from these analyses (n = 2 and n = 4 data points were removed from FD and DD analyses, respectively).  Table 2 for a full description of the model variables. The transmission risk corresponds with a mid-range level of transmissibility where ˇa ∼ 5 and ˇf ∼ 5.
Generating between-farm profiles for relative transmission risk Fig. 5 shows the distribution of between-farm association frequencies for slaughterhouse, catching company, integrated company and combined layers of the poultry network, for all farms recorded in the Catching Company Database (n = 108). However, it should be noted that the median between-farm association frequency overall combining all network layers, when based on the full Poultry Network Database, was estimated to be 343 farms (range = 105-1453 farms).
Farms with high estimates for both between-farm associations and TR i would produce the highest risk of a widespread epidemic. When cross-classifying the between-farm association frequency with TR i , a large proportion of farms had a below median value for both factors (21% of farms, see bottom left quadrant of Fig. 6). However, a further 21% of farms had a relatively low estimate of TR i and a relatively high estimate of between-farm association frequency (see top left quadrant of Fig. 6). The largest fraction of farms (28%) had the highest risk combination, with above median estimates for both factors (see top right quadrant of Fig. 6). No difference was found in the distribution of farms between FD and DD transmission modes when dichotomising TR i by median values.

Cross-population scale interactions
The integration of infection transmission dynamics at the within-group and population level is increasingly important as greater demands are made of predictive mathematical models (Haydon and Matthews, 2007;Kao et al., 2007). For pathogens that are likely to spread rapidly within a farm, the risk of onward transmission via fomites depends on the opportunity for an on-to-farm movement to coincide with an outbreak. We have explored the importance of these cross-scale interactions for highly pathogenic avian influenza (HPAI) within the British commercial broiler poultry industry. By considering the interaction between the farm-level infectious period and explicit temporal pattern of catching-team visits, our results have implications for the between-farm spread, or the network dynamics, of HPAI.
Poultry flock size may be indicative of both the time to detection of an outbreak, and the expected amount of on-to-farm movement activity. However, this does not account for the trade-off between the build-up of faecal virus within the poultry house and the opportunity for virus exposure. Our results suggest that the overall effect of flock size on the relative transmission risk (TR) is sensitive to the bird-to-bird transmission mode. For frequency-dependent (FD) transmission, the time-to-detection increased with flock size (Fig. 2a), as highlighted previously by Savill et al. (2008). However, for density-dependent (DD) transmission, greater infection pressure and greater movement activity, but more rapid outbreak detection, interacts to influence the opportunity for onward spread from larger flocks (Fig. 2b). This interaction had not previously been considered with respect to HPAI within commercial poultry.
As the infection pressure (i.e. ˇa + ˇf) increased, so did the amount of infectious faeces thus increasing TR; this counteracted the shortened time-to-detection for DD transmission for some regions of transmissibility but TR was highly sensitive to the balance between these factors causing an oscillatory effect (Fig. 3d). Overall, mid-range sized flocks (∼25,000-35,000 birds) were found to have the greatest TR under DD transmission, in contrast to larger flocks (∼35,000-45,600 birds) under FD transmission.
Faeces are not removed mid production cycle in British commercial broiler chicken farms, therefore it was assumed that the TR depended on the product of available infectious faecal material and the number of vehicle loads. However, a nonlinear dependence (e.g. an upper contamination threshold) could change the sensitivity to flock size and the impact of transmission mode assumptions. Given the lack of empirical data available in relation to viral survival times we assumed a constant rate of viral decay in line with previous published studies (Savill et al., 2006(Savill et al., , 2008. Whilst we cannot postulate the appropriateness of alternative assumptions, we anticipate that this would not impact on our qualitative results.

Generating farm-level risk profiles
The identification of predominantly movement-related predictors that depend on the bird-to-bird transmission mode has important implications for disease control strategies. For example, the multivariable regression analyses suggest that reducing the overall frequency of on-to-farm visits would be more effective for FD transmission, in contrast to limiting the clustering of consecutive visit days for DD transmission (Tables 3a and 3b). As the temporal clustering of visits at the farm-level was strong, activities requiring on-to-farm visits mid-production cycle, such as the practice of flock-thinning, might be particularly detrimental for DD transmission characteristics. Such practices are recognised risk fac- of the likely bird-to-bird contact behaviour, and therefore the predominant mechanism of transmission within a poultry flock, is needed to better characterise this risk. In our model, the highest potential for spread via catchingteam visits was associated with outbreaks of HPAI virus with high transmissibility (R 0 ∼25-30), in large broiler flocks, under FD transmission. Delayed detection of the outbreak by a farmer (as demonstrated by increases to the mortality threshold triggering outbreak detection, see Supplementary Material Section 4) would further increase this risk. Delayed notification might be expected for less virulent poultry pathogens such as low pathogenic avian influenza (LPAI), which despite sometimes significant mortality (Bano et al., 2003), typically has less severe morbidity and lower mortality (Defra, 2006). This observation has important implications for the spread of HPAI, as the undetected circulation of LPAI virus precursors may increase the likelihood of a subsequent HPAI outbreak (Mannelli et al., 2006). Though particular HPAI subtypes, such as H5N1, are expected to be detected relatively rapidly, with some analyses indicating no longer than one week (Yoon et al., 2005), evidence from experimental work using the H7N7 subtype, implicated in an outbreak in the Netherlands in 2003, suggests that HPAI may have been left to circulate within a flock undetected for 11-15 days (Bos et al., 2007). Due to the difficulties in directly extrapolating from experimental studies typically involving very few birds (i.e. typically a maximum of five susceptible in-contact birds, in the absence of concurrent infection and environmental factors that may enhance transmission under field conditions), the upper limit for transmission rates in relation to commercial poultry flocks may be greater than estimated by these studies.
In view of this likely variability between virus subtypes, and the possible under-estimation of transmission rates by experimental studies, we explored a range of time to detections (2-10 days), representing a range of plausible HPAI outbreak characteristics, reflecting in particular the most virulent HPAI viruses. Furthermore, in view of the impact of farmer behaviour, we also considered a range of mortality thresholds, MTs, and find that as the MT increases, TR increases for a greater range of infection pressures; for large flocks and thresholds beyond 0.7% we might not observe the "rapid drop in TR" phenomenon as the infection pressure increases within the range considered (see Supplementary Material Section 4).

Generating between-farm risk profiles
If the most highly connected farms also have a high TR, the overall risk of a widespread outbreak could be high; this was the case for 28% of farms in these analyses. However, a large proportion of farms (42%) had discordant network-connectivity and relative transmission risk characteristics (Fig. 6); approximately half of these (21%) had a low estimate of TR (corresponding to a mid-range transmissibility scenario, TR i ) but a high between-farm association frequency. In this case a high impact at the between-farm level but with a low associated risk of occurrence could be construed as representing a medium-level risk.
A more detailed understanding of the relative importance of farm-level transmission "risk" and between-farm "impact" would allow us to better characterise farms and determine the most appropriate control measures. For example, categorising farms based on combinations of exposure "risk" and overall "impact" may enable the effective targeting of control measures -via a relatively greater focus on preventive measures such as farm-level biosecurity for farms with high estimated transmission risk, and a relatively greater focus on contact-tracing during an outbreak for farms with high estimated between-farm association frequency. Farms with high estimates for both factors warrant the greatest focus from prevention and control measures.
The most crucial factor for minimising the risk of spread between farms will be the rapidity with which farmers notify and report a suspected outbreak, as discussed for farm-level risk profiles in relation to assumed mortality thresholds. In this regard, further consideration must be given to the potential negative consequences associated with notifiable HPAI and the risk of false-alarms (Savill et al., 2008) that may be perceived by individual farmers. Ultimately, when developing a risk profile different aspects of riskinitial incursion vs. the potential for spread between farms, as well as the scale of observation -each have different implications for implementing control measures.

Conclusions
We have shown how particular assumptions for within-farm dynamics result in counter-intuitive implications for pathogen spread at the population-level; increased transmissibility at the within-farm level sometimes presents a reduced transmission risk to other farms. This is the first demonstration of the potential importance of these cross-scale interactions for HPAI in a commercial poultry industry. Together with consideration to the relative importance of network connectivity, these findings have implications for the targeting of HPAI control measures at the farm-level.