Low temperature is one of the critical factors for determining agricultural production. Therefore, the main objective was to develop the durum wheat lines, which are tolerant to cold particularly in Central Anatolia. The recombinant inbred lines (RILs) (141) derived from F6 progeny of durum wheat, developed from Kunduru-1149 (female parent) and Cham-1 (male parent) cross, were characterized in terms of their cold tolerance both in the field and under controlled (in growth chamber) conditions. A 0–9 scale was used for the visual evaluation of cold damage under field conditions, and a quantitative measure of photochemical efficiency (Fv/Fm = variable fluorescence/maximum fluorescence) method was used for the evaluation of cold damage under controlled conditions. Seventeen RILs displayed the scores within the range of 1–2.5 in 0–9 scale evaluation, while 41 RILs showed the scores equal to or higher than 0.70 approaching to optimal (~0.83) in Fv/Fm evaluation. Some quality parameters tested and some gliadin fragments and random amplified polymorphic DNA loci showed meaningful correlations with the data of 0–9 scale and Fv/Fm ratios from lower to the medium level. The RILs considered as cold resistant; might have the potential to be used as parental lines to develop cold-resistant durum wheat varieties in the future durum wheat-breeding programmes.

Heat stress during the reproductive and grain-filling stages leads to severe yield losses in rice, especially in irrigated areas during the dry season in Thailand. Thus, breeding for heat-tolerant rice is one of the strategies for developing rice crops adapted to the effects of high temperature in major growing regions. M9962, a fast neutron-induced mutant, was identified as a heat-tolerant line that was then used as a donor parent and crossed with the Phisanulok 2 (PSL2) cultivar (the recurrent parent) in a backcross breeding programme. Heat tolerance scoring under controlled greenhouse and field conditions was based on the percentage of spikelet fertility. Yield trials for the BC2F7 and BC3F6 lines were subsequently conducted under field conditions at three locations and under controlled greenhouse conditions during the dry season of 2019 (DS2019). In the early generation, 210 plants of the BC1F1 population and 440 plants of BC2F1 population were evaluated for their heat tolerance, and the spikelet fertility distribution within the BC1F1, BC2F1 and BC3F1 populations showed a skewed-right pattern. In addition, most rice plants in every generation were heat sensitive. Finally, compared with PSL2, the four-candidate heat-tolerant lines exhibited a high percentage of spikelet fertility and high yield. In addition, the one candidate line, BC2F7-6-5-4-1-1-21 was most similar to PSL2 in terms of genetic background, plant type and grain quality. Thus, backcross breeding programmes via the spikelet fertility can be used as an indirect trait to select for heat tolerance.

Owing to its contribution to the maintenance of carbon stocks, soil nitrogen and nutrient cycling for subsequent crops, the integrated systems become increasingly important for agricultural conservation. Thus, the objective of this study was to evaluate the biomass production of and total nutrient in Brachiaria spp. and Panicum maximum forage grasses used as mulch and soybean yields in an integrated crop–livestock system and second-crop maize succession system. The treatments consisted of the following cropping systems: Xaraes palisadegrass intercropped with soybean, Congo grass intercropped with soybean, Mombaça guinea grass intercropped with soybean, Tamani guinea grass intercropped with soybean and a soybean/maize succession system. The forage grasses were established during the soybean R6–R7 stage. Compared with Congo grass, Xaraes palisadegrass, Mombaça guinea grass and Tamani guinea grass produced more biomass and equivalent amounts of fertilizer returned to the soil and resulted in greater nutrient cycling, indicating the benefits of these grasses for use as mulch in integrated production systems. Maize had a greater C/N ratio, but the forage grasses also exhibited high potential by protecting the soil until the end of the soybean development cycle. The use of an integrated crop–livestock system combined with a forage cropping system provided greater soil nutrient cycling than the maize cropping system did, which resulted in increased soybean yields, thus contributing to the sustainability of agricultural systems.

Water consumption in agriculture is impossible without considering relations between water, soil and plant. In this regard, there are various models and developed software in order to evaluate relation between soil, water and crop growth stages. These models can be used for irrigation planning if properly optimized and applied. AquaCrop is one of the known crop models, which was developed by the Food and Agriculture Organization of the United Nations. In order to optimize this model for crop production and irrigation management, an experiment was developed in a hydroponic cucumber greenhouse. Various parameters including water consumption volume, crop yield and leaf area index were measured during a season. A fuzzy control system was utilized for controlling temperature, relative humidity, planting bed moisture, light intensity and carbon dioxide values. The main purpose of designing a control system in the greenhouse is to achieve the desired values of temperature and relative humidity. In this model, evapotranspiration, irrigation requirements and crop yield were simulated. The results show that the AquaCrop model can estimate evapotranspiration with the least error in the greenhouse environment, which is controlled by a fuzzy controller. Also the system has estimated the crop yield and biomass of the product with a good degree of precision and it may support crop production in a greenhouse, including crop management and environmental control.

Increasing crop diversity has been highly recommended because of its environmental and economic benefits. However, the impacts of crop diversity on soil properties are not well documented. Thus, the present study was conducted to assess the impacts of crop diversity on selected soil quality indicators. The cropping systems investigated here included wheat (Triticum aestivum L.) grown continuously for 5 years as mono-cropping (MC), and a 5-year cropping sequence [(wheat–cover crop (CC)–corn (Zea mays L.)–pea (Pisum sativum L.) and barley (Hordeum vulgare L.)–sunflower (Helianthus annuus L.)]. Each crop was present every year. This study was conducted in the northern Great Plains of North America, and soil quality data were collected for 2016 and 2017. Selected soil quality indicators that include: soil pH, organic carbon (SOC), cold water-extractable C (CWC) and N (CWN), hot water-extractable C (HWC) and N (HWN), microbial biomass carbon (MBC), bulk density (BD), water retention (SWR), wet soil aggregate stability (WAS), and urease and β-glucoside enzyme activity were measured after the completion of 5-year rotation cycle and the following year. Crop diversity did not affect soil pH, CWC, CWN, HWC, HWN and SWR. Cropping systems that contained CC increased SOC at shallow depths compared to the systems that did not have CC. Crop diversity increased WAS, MBC, and urease and β-glucoside enzyme activity compared with the MC. Comparison of electrical conductivity (EC) measured in this study to the baseline values at the research site prior to the establishment of treatments revealed that crop rotation decreased EC over time. Results indicate that crop diversity can improve soil quality, thus promoting sustainable agriculture.

Shaded coffee systems can mitigate climate change by fixation of atmospheric carbon dioxide (CO2) in soil. Understanding soil organic carbon (SOC) storage and the factors influencing SOC in coffee plantations are necessary for the development of sound land management practices to prevent land degradation and minimize SOC losses. This study was conducted in the main coffee-growing regions of Yunnan; SOC concentrations and storage of shaded and unshaded coffee systems were assessed in the top 40 cm of soil. Relationships between SOC concentration and factors affecting SOC were analysed using multiple linear regression based on the forward and backward stepwise regression method. Factors analysed were soil bulk density (ρb), soil pH, total nitrogen of soil (N), mean annual temperature (MAT), mean annual moisture (MAM), mean annual precipitation (MAP) and elevations (E). Akaike's information criterion (AIC), coefficient of determination (R2), root mean square error (RMSE) and residual sum of squares (RSS) were used to describe the accuracy of multiple linear regression models. Results showed that mean SOC concentration and storage decreased significantly with depth under unshaded coffee systems. Mean SOC concentration and storage were higher in shaded than unshaded coffee systems at 20–40 cm depth. The correlations between SOC concentration and ρb, pH and N were significant. Evidence from the multiple linear regression model showed that soil bulk density (ρb), soil pH, total nitrogen of soil (N) and climatic variables had the greatest impact on soil carbon storage in the coffee system.

Perception of freshwater use varies between nations and has led to concerns of how to evaluate water use for sustainable food production. The water footprint of beef cattle (WFB) is an important metric to determine current levels of freshwater use and to set sustainability goals. However, current WFB publications provide broad WF values with inconsistent units preventing direct comparison of WFB models. The water footprint assessment (WFA) methodologies use static physio-enviro-managerial equations, rather than dynamic, which limits their ability to estimate cattle water use. This study aimed to advance current WFA methods for WFB estimation by formulating the WFA into a system dynamics methodology to adequately characterize the major phases of the beef cattle industry and provide a tool to identify high-leverage solutions for complex water use systems. Texas is one of the largest cattle producing areas in the United States, a significant water user. This geolocation is an ideal template for WFB estimation in other regions due to its diverse geography, management-cultures, climate and natural resources. The Texas Beef Water Footprint model comprised seven submodels (cattle population, growth, nutrition, forage, WFB, supply chain and regional water use; 1432 state variables). Calibration of our model replicated initial WFB values from an independent study by Chapagain and Hoekstra in 2003 (CH2003). This CH2003 v. Texas production scenarios evaluated model parameters and assumptions and estimated a 41–66% WFB variability. The current model provides an insightful tool to improve complex, unsustainable and inefficient water use systems.

The successful ensiling of lucerne (Medicago sativa L.) depends on a rapid acidification in the silo and consequently relies on a sufficient proliferation of, particularly homofermentative, lactic acid bacteria. Similarly, growth of spoilage bacteria, such as enterobacteria and clostridia, must be suppressed and silage additives are therefore frequently applied to promote favourable conditions during ensiling. Three silage additives or soil were applied during lucerne ensiling and investigated for their effects on silage quality characteristics and abundances of total bacteria as well as the bacterial key players Lactobacillus spp., homofermentative Lact. plantarum, heterofermentative Lact. buchneri, Clostridium spp. and Enterobacteriaceae after 30 days of storage. Inoculation with viable Lact. plantarum resulted in highest concentration of this species and excellent silage quality, i.e. high lactic acid concentration coupled with low acetic acid and ammonia-nitrogen concentrations. A sodium nitrite and hexamine-based additive did not support growth of lactic acid bacteria, which was also apparent by higher pH and low lactic acid concentration. No effect of treatments was found on spoilage-related enterobacteria and clostridia, even not when adding soil to lucerne to increase initial clostridial contamination. However, soil treatment resulted in increased ammonia-nitrogen and acetic acid concentrations. Consequently, among the bacterial key players, lactic acid bacteria concentrations were related to silage quality. Regarding spoilage bacteria, however, alterations in silage quality characteristics were not reflected in the abundances of enterobacteria and clostridia. Future investigations should underpin the present findings and help to understand how silage additives affect microbial key players and silage fermentation.

This study aimed to evaluate the fermentation profile and nutritional quality of silages composed of cactus pear and maniçoba. Two experiments were performed: the first evaluated the fermentation characteristics, chemical composition and determined the organic acids in cactus pear silages with the inclusion of five levels of maniçoba (0, 25, 50, 75 and 100%) in six opening times (1, 7, 15, 30, 60 and 90 days). The second experiment determined the nutrient intake, digestibility, water balance and nitrogen balance in Canindé goats fed diets based on cactus pear silage with the inclusion of four levels of maniçoba (25, 50, 75 and 100%), with six animals per treatment. The increase in maniçoba levels in cactus pear silage provided a linear increase in the butyric acid, dry matter (DM), ether extract, crude protein, neutral detergent fibre, acid detergent fibre, lignin, cellulose, hemicellulose, water intake via drinking fountain and metabolic water, and reduced the pH, lactic acid, acetic acid, mineral matter, total carbohydrates, non-fibrous carbohydrates, water intake via food, total water intake, water excreted in the faeces, water excreted in the urine, total water excretion and water balance. Quadratic behaviour was observed for N-NH3, DM recovery and propionic acid, with an increase in maniçoba levels in cactus pear silages. Regarding the different opening times, there was a significant effect in pH, N-NH3, acetic acid, lactic acid and butyric acid (P < 0.050). The inclusion of maniçoba in cactus pear silage improved the fermentation characteristics and nutritional quality to be used in diets for goats.

The paper presents results from a study examining the relationship between large-scale modes of climate variability with the fluctuations in the yield of barley, durum wheat, olives and sunflower crops in Tuscany, Italy. In particular, the blocking circulation over the growing season, with associated hot and dry conditions, decreased yield for olive crops, barley and durum wheat. The teleconnections analysed in this study are the winter North Atlantic Oscillation (NAO) and the Summer North Atlantic Oscillation (SNAO); the West African Monsoon (WAM) and the Intertropical Front (ITF); and although NAO, SNAO, ITF and WAM are not strictly related to each other, the values of these indices are strongly related to the atmospheric circulation regimes and related weather types. Thus, they have an impact on precipitation and temperature patterns in Italy and on yields of important crops in Tuscany. Results show that the large-scale temperate and tropical variability directly influences the crop yield through three main circulation regimes. These patterns illustrate the importance of the large-scale modes, which, together with the associated weather types, have an impact directly on Tuscan crop yields; both barley and olive yields decline significantly when the ITF is further north with warmer and drier conditions in Italy.

For ensuring food demand of the fast growing population in developing countries, quantification of crop yield gaps and exploring production constraints are very crucial. Sorghum is one of the most important climate change resilient crops in the rainfed farming systems of semi-arid tropics. However, there is little information about yield gaps and production constraints. This study aimed at analysing existing yield gaps and exploring major constraints of sorghum production in Southwest Ethiopia. A crop simulation model approach using AquaCrop and DSSAT was used to estimate potential yield and analyse the yield gaps. Model calibration and evaluation was performed using data from field experiments conducted in 2018 and 2019. Sorghum production constraints were assessed using a survey. The actual and water-limited yield of sorghum ranged from 0.58 to 2.51 and 3.6 to 6.47 t/ha, respectively for the period 2003–17. The regional yield gaps of sorghum for the targeted period were 3.02–3.95 t/ha with a mean value of 3.51 t/ha. Majority of respondent farmers considered seasonal rainfall risk (98%), poor soil fertility (86%), lack of improved varieties (78%) and inadequate weed management (56%) as major factors responsible for the existing yield gaps. The mean exploitable yield gap (2.5 t/ha) between water-limited and actual yield showed the level of existing opportunity for improvement in the actual productivity of sorghum. The gaps also call for introduction of proper interventions such as adoption of improved varieties, planting date adjustment, conservation tillage, fertilizer application and on time weed management.

The determination of livestock growth patterns is important for meat or milk production systems, and nonlinear models are used to summarize and interpret the information. The aim of this study was to more accurately estimate growth curve parameters in buffalo cows by evaluating and selecting nonlinear mixed models that employ different types of residuals and include or not contemporary groups (CG) as a covariate. Weight records from 720 animals obtained over a period of 60 months were used. The growth curves were fit using nonlinear mixed-effects models. The Bertalanffy, Gompertz and Logistic models were evaluated. Modelling residuals using four structures (constant, combined, exponential and proportional) and the inclusion or not of CG in the models were also evaluated. The Akaike information criterion (AIC) and Bayesian information criterion (BIC) were used to select the model. In addition to estimating the parameters of the nonlinear growth models and their correlations, the instantaneous growth rate and inflection point were obtained. The Bertalanffy model with a combined residual structure and CG exhibited the lowest AIC and BIC values. Asymptotic weight (A) estimates ranged from 621.8 to 742.1 kg, and the maturity rate (k) ranged from 0.068 to 0.115 kg/month. The correlation between A and k ranged from −0.32 to −0.82 among the models evaluated. The selection criteria indicated that the Bertalanffy model was the most suitable for growth curve analysis in buffaloes.

Tropical soils tend to harden during drying due to the generally low content of free-iron and organic carbon, combined with high fine sand and silt proportions. It was hypothesized that the change in soil physical condition induced by the addition of a leguminous mulch in cohesive tropical soil enriched with calcium may mitigate soil hardening through wetting and drying cycles by rain or irrigation, thereby improving the soil rootability. A leguminous mulch was added in different concentrations to a structurally fragile tropical soil enriched with calcium, which then had different irrigation intervals. The treatments were with or without mulch (10 t/ha), with or without added nitrogen (100 kg/ha at 2 intervals) and two irrigation intervals. In 2015 the irrigation intervals were either 4 or 8 days, and in 2016 they were either 6 or 9 days. Two years were used in the attempt to achieve greater differences, as for tested variables, between treatments. Maize planted in these soil treatments was measured for physiological performance, water use efficiency and yield. Mulch used on structurally fragile tropical soil enriched with calcium was found to delay increased penetration resistance from hardening by wet/dry cycles. In this context, an improved soil rootability led to an enlargement of the leaf area index, greater nitrogen uptake and increased CO2 assimilation. This had important physiological consequences due to the positive effect on increased dry matter production and maize yield. In addition, these results suggested that mulch, used with urea, can delay the water supply for 3 or 4 days due to improvements in soil rootability caused by calcium and organic matter interactions. This may be crucial to a region where small intervals without rain are increasingly common due to global climate change. Therefore, due to a greater water use efficiency, this strategy may be a profitable way to increase crop productivity in tropical conditions rather than increasing water and nutrient application alone.

Oat (Avena spp. L.) intercropped with common vetch (Vicia sativa L.) has shown high forage yield and crude protein content (CP) in rainfed Mediterranean environments of Greece. The objective was to study the oat genotypic interactions for yield, quality and agronomic traits of four oat accessions (cv. Flega and lines B2, B3 and B4) grown as intercrops with a common vetch cultivar (Pegasus) and as monocrops without or with nitrogen (N) supplementation (100 kg N/ha). In all mixtures, the land equivalent ratio exceeded the unit indicating that there was a yield advantage. Intercrops produced, on average, 5.0% less dry matter (DM) and 38.9% more crude protein yield (CPY) compared to monocropped oats with N fertilization and 3.4% more DM and 75.5% more CPY compared to monocropped oats without N fertilization. The DM, CP and quality differences between the intercropped and monocropped oats were attributed to oats by vetch cultivar interactions. The oat B2 intercropped was the shortest, the latest flowering and the less vigorous early in the season, a trait that allowed common vetch to dominate in the intercrop (vetch participation 770 g/kg). The high vetch fragment in the intercrop contributed to high CP, CPY, acid detergent lignin, total digestible nutrients, relative feeding value, relative forage quality, net energy for lactation and to low neutral detergent fibre, acid detergent fibre and hemicellulose. The replacement value index indicated that the intercrop was 21% more economically viable in comparison to the fertilized monocrop.

Sheep are grazers and goats are intermediate feeders. By employing O2 consumption and heart rate measurements, resting metabolic rate (RMR) and field metabolic rate (FMR) were determined in four male fat-tailed Awassi sheep (44.0 ± 3.94) and four male Baladi goats (35.5 ± 5.42 kg) that were co-grazing natural pasture in the Negev Desert. There were 67.7 ± 3.75 g DM/m2 of herbaceous vegetation biomass, which was rapidly becoming senescent and more fibrous. We hypothesized that FMR of these desert-adapted ruminants would be relatively low when compared to other sheep and goat breeds, as animals in arid areas tend to have low metabolic rates. Both sheep (n = 6) and goats (n = 6) foraged 71% of the allotted 11 h free-pasture period; however, sheep grazed more than goats (P < 0.001); whereas goats browsed more than sheep (P < 0.001). RMR was higher (P = 0.007) in sheep than in goats (529 ± 23.5 v. 474 ± 25.4 kJ/kg0.75 BW/d), but FMR did not differ between species (618 ± 55.7 v. 613 ± 115.2 kJ/kg0.75 BW/d). In addition, the cost of activities, as a proportion of FMR, did not differ between sheep and goats; FMR increased by 89 kJ/kg0.75 BW/d or 17% in sheep and by 138 kJ/kg0.75 BW/d or 29% in goats. In comparing FMRs of sheep and goats in this study with these species in other studies, differences were inconsistent and, therefore, our hypothesis was not supported.

Data of 2780 Markhoz kids originated from 1216 dams and 211 sires during 1993–2016 in Markhoz Goat Breeding Station, located in Sanandaj, Iran, were used. Traits investigated were body weights at birth, weaning, six-month age [six months weight (6MW)], nine-month age and yearling age [yearling weight (YW)]. Two considered multivariate models including standard multivariate model (SMM) and fully recursive multivariate model (FRM) were compared using deviance information criterion (DIC) and predictive ability measures including mean square of error (MSE) and Pearson's correlation coefficient between the observed and predicted values $(r(y,\hat{y}))$ of records. Spearman's rank correlation coefficients between posterior means of direct genetic effects of the studied traits of kids under SMM and FRM were also calculated across all, 50, 10 and 1% top-ranked animals. In general, FRM performed better than SMM in terms of lower DIC and MSE and also higher $r\lpar y\comma \;\hat{y}\rpar$. For all traits, the lowest MSE and the highest $r\lpar y\comma \;\hat{y}\rpar$ were obtained under FRM. All structural coefficients estimated under FRM were statistically significant except for that of 6MW on YW. Comparisons of Spearman's rank correlations between posterior means of direct genetic effects of kids for growth traits under SMM and FRM revealed that taking the causal relationships among the studied growth traits of Markhoz goat into account may cause considerable re-ranking for the animals in terms of estimated breeding values, especially for the top-ranked animals. It may be concluded that FRM had more plausibility over SMM for genetic evaluation of the studied growth traits in Markhoz goat.

Feeding fodder beet (FB) to dairy cows in early lactation has recently been adopted by New Zealand dairy producers despite limited definition of feeding and grazing management practices that may prevent acute and sub-acute ruminal acidosis (SARA). This modelling study aimed to characterize changes of rumen pH, milk production and total discomfort from FB and define practical feeding strategies of a mixed herbage and FB diet. The deterministic, dynamic and mechanistic model MINDY was used to compare a factorial arrangement of FB allowance (FBA), herbage allowance (HA) and time of allocation. The FBA were 0, 2, 4 or 7 kg dry matter (DM)/cow/day (0FB, 2FB, 4FB and 7FB, respectively) and HA were 18, 24 or 48 kg DM/cow/day above ground. All combinations were offered either in the morning or afternoon or split across two equal meals. Milk production from 2FB diets was similar to 0FB but declined by 4 and 16% when FB increased to 4 and 7 kg DM, respectively. MINDY predicted that 7FB would result in SARA and that rumen conditions were sub-optimal even at moderate FBA (pH < 5.6 for 160 and 90 min/day, 7FB and 4FB respectively). Pareto front analysis identified the best compromise between high milk production and low total discomfort was achieved by splitting the 2FB diet into two equal meals fed each day with 48 kg DM herbage. However, due to low milk response and high risk of acidosis, it is concluded that FB is a poor supplement for lactating dairy cows.

It was aimed to simultaneously study standardized ileal digestible (SID) tryptophan (Trp) and lysine (Lys) for gilts. A digestibility assay was previously conducted to determine the SID amino acid in the basal diet (low levels of SID Trp and Lys). Sixty-four gilts (15.04 ± 1.44 kg) were allotted to 16 diets in a 4 × 4 factorial arrangement (1.55, 1.85, 2.15 and 2.45 g/kg SID Trp and 9.72, 11.12, 12.52 and 13.92 g/kg SID Lys) with four replicates per treatment. Performance, longissimus muscle (LM), backfat thickness (BF) and blood variables were evaluated. An interaction was observed for G:F, and by response surface model, the optimum Trp level was achieved at 2.15 g/kg (0.159 g/MJ of ME). A quadratic effect of Trp was observed on body weight (BW) and average daily gain (ADG); the daily feed intake increased linearly as Trp increased. The optimum Trp levels of 2.25 and 2.24 g/kg were estimated for BW and ADG, respectively. The BF increased with increasing levels of Trp. There was a quadratic and linear effect of Trp and Lys, respectively, on the LM, in which the optimum Trp level was determined as 2.05 g/kg in the diet. Plasma urea nitrogen decreased as Trp and Lys levels increased. Using estimates provided by response surface, maximized G:F ratio was obtained at 2.15 g SID Trp/kg of diet and at least 13.92 g SID Lys/kg of diet is necessary to optimize the G:F for 15–30 kg gilts, providing a Trp:Lys ratio of 15.4:100.

The development of resilient cropping systems with high yield stability is becoming increasingly important due to future climatic and agronomic challenges. Consequently, it is essential to compare the effects of different agronomic management practices, such as cropping sequences and nutrient supply, on the stability of crop yields. Long-term experiments are a valuable resource for investigating these effects, as they provide enough time to accurately estimate stability parameters. The objective of the current study was to compare the effects of different cropping sequencing (#1: continuous v. rotational), fertilization (#2: mineral v. organic) and straw management techniques (in the case of continuous wheat; #3: removal v. incorporation) on the yield stability of winter wheat; yield risk (the probability of yield falling below a threshold yield level) and inter-annual yield variability were used as stability indicators of the effects. Long-term yield data from the Broadbalk Wheat Experiment (Rothamsted, UK) were analysed using a mixed model. Overall, the results showed that rotational cropping combined with sufficient mineral N fertilizer, with or without organic manure, ensured stable wheat yields while reducing yield risk. In contrast, higher yield risks and inter-annual yield variabilities were found in continuous wheat sections with less mineral N fertilizer or with organic manure only.

The aim of this study was to investigate whether social rank (SR) and stability in the order of entrance to the milking parlour are associated with production traits. The study was conducted on a dairy farm where cows (n = 215) were managed in three groups according to lactation stage (group 1: 78 cows, 0–100 days in milk (DIM); group 2: 65 cows, 101–200 DIM and group 3: 72 cows, >200 DIM). SR was calculated from observations made from agonistic behaviour performed at the water troughs and feed bunks of each pen (n = 3). The animals were classified into three levels of dominance based on at least five clear interactions, resulting in: 61 dominant, 75 intermediate and 69 subordinate cows based on SR. Stability in the order of entry was estimated as the standard deviation of the entry position. SR was weakly correlated with milk yield, urea and protein content in milk. The results showed that stable cows had higher milk production and entered the milking parlour after the non-stable animals. Stability in the order of entry to the milking parlour was not affected by SR. Overall, the use of milking facilities appears to be associated with production traits rather than SR.