Leaf trait co-ordination in relation to construction cost, carbon gain and resource-use efficiency in exotic invasive and native woody vine species

Background and Aims: Success of invasive plant species is thought to be linked with their higher leaf carbon fixation strategy, enabling them to capture and utilize resources better than native species, and thus pre-empt and maintain space. However, these traits are not well-defined for invasive woody vines. Methods: In a glass house setting, experiments were conducted to examine how leaf carbon gain strategies differ between non-indigenous invasive and native woody vines of south-eastern Australia, by investigating their biomass gain, leaf structural, nutrient and physiological traits under changing light and moisture regimes. Key Results: Leaf construction cost (CC), calorific value and carbon : nitrogen (C : N) ratio were lower in the invasive group, while ash content, N, maximum photosynthesis, light-use efficiency, photosynthetic energyuse efficiency (PEUE) and specific leaf area (SLA) were higher in this group relative to the native group. Trait plasticity, relative growth rate (RGR), photosynthetic nitrogen-use efficiency and water-use efficiency did not differ significantly between the groups. However, across light resource, regression analyses indicated that at a common (same) leaf CC and PEUE, a higher biomass RGR resulted for the invasive group; also at a common SLA, a lower CC but higher N resulted for the invasive group. Overall, trait co-ordination (using pair-wise correlation analyses) was better in the invasive group. Ordination using 16 leaf traits indicated that the major axis of invasive-native dichotomy is primarily driven by SLA and CC (including its components and/or derivative of PEUE) and was significantly linked with RGR. Conclusions: These results demonstrated that while not all measures of leaf resource traits may differ between the two groups, the higher level of trait correlation and higher revenue returned (RGR) per unit of major resource need (CC) and use (PEUE) in the invasive group is in line with their rapid spread where introduced.

Responses of Super Rice (Oryza sativa L.) to Different Planting Methods for Grain Yield and Nitrogen-Use Efficiency in the Single Cropping Season

To break the yield ceiling of rice production, a super rice project was developed in 1996 to breed rice varieties with super high yield. A two-year experiment was conducted to evaluate yield and nitrogen (N)-use response of super rice to different planting methods in the single cropping season. A total of 17 rice varieties, including 13 super rice and four non-super checks (CK), were grown under three N levels [0 (N0), 150 (N150), and 225 (N225) kg ha−1] and two planting methods [transplanting (TP) and direct-seeding in wet conditions (WDS)]. Grain yield under WDS (7.69 t ha−1) was generally lower than TP (8.58 t ha−1). However, grain yield under different planting methods was affected by N rates as well as variety groups. In both years, there was no difference in grain yield between super and CK varieties at N150, irrespective of planting methods. However, grain yield difference was dramatic in japonica groups at N225, that is, there was an 11.3% and 14.1% average increase in super rice than in CK varieties in WDS and TP, respectively. This suggests that high N input contributes to narrowing the yield gap in super rice varieties, which also indicates that super rice was bred for high fertility conditions. In the japonica group, more N was accumulated in super rice than in CK at N225, but no difference was found between super and CK varieties at N0 and N150. Similar results were also found for N agronomic efficiency. The results suggest that super rice varieties have an advantage for N-use efficiency when high N is applied. The response of super rice was greater under TP than under WDS. The results suggest that the need to further improve agronomic and other management practices to achieve high yield and N-use efficiency for super rice varieties in WDS.

Nitrogen use efficiency in summer sorghum grown on clay soils

Nitrogen fertilizer inputs dominate the fertilizer budget of grain sorghum growers in northern Australia, so optimizing use efficiency and minimizing losses are a primary agronomic objective. We report results from three experiments in southern Queensland sown on contrasting soil types and with contrasting rotation histories in the 2012-2013 summer season. Experiments were designed to quantify the response of grain sorghum to rates of N fertilizer applied as urea. Labelled 15N fertilizer was applied in microplots to determine the fate of applied N, while nitrous oxide (N2O) emissions were continuously monitored at Kingaroy (grass or legume ley histories) and Kingsthorpe (continuous grain cropping). Nitrous oxide is a useful indicator of gaseous N losses. Crops at all sites responded strongly to fertilizer N applications, with yields of unfertilized treatments ranging from 17% to 52% of N-unlimited potential. Maximum yields ranged from 4500 (Kupunn) to 5450 (Kingaroy) and 8010 (Kingsthorpe) kg/ha. Agronomic efficiency (kg additional grain produced/kg fertilizer N applied) at the optimum N rate on the Vertosol sites was 23 (80 N, Kupunn) to 25 (160N, Kingsthorpe), but 40-42 on the Ferrosols at Kingaroy (70-100N). Cumulative N2O emissions ranged from 0.44% (Kingaroy legume) to 0.93% (Kingsthorpe) and 1.15% (Kingaroy grass) of the optimum fertilizer N rate at each site, with greatest emissions from the Vertosol at Kingsthorpe. The similarity in N2O emissions factors between Kingaroy and Kingsthorpe contrasted markedly with the recovery of applied fertilizer N in plant and soil. Apparent losses of fertilizer N ranged from 0-5% (Ferrosols at Kingaroy) to 40-48% (Vertosols at Kupunn and Kingsthorpe). The greater losses on the Vertosols were attributed to denitrification losses and illustrate the greater risks of N losses in these soils in wet seasonal conditions.

Effects of weed control and fertilization at early establishment on tree nitrogen and water use in an exotic F1 hybrid pine of subtropical Australia

Purpose We investigated the effects of weed control and fertilization at early establishment on foliar stable carbon (δ13C) and nitrogen (N) isotope (δ15N) compositions, foliar N concentration, tree growth and biomass, relative weed cover and other physiological traits in a 2-year old F1 hybrid (Pinus elliottii var. elliottii (Engelm) × Pinus caribaea var. hondurensis (Barr. ex Golf.)) plantation grown on a yellow earth in southeast Queensland of subtropical Australia. Materials and methods Treatments included routine weed control, luxury weed control, intermediate weed control, mechanical weed control, nil weed control, and routine and luxury fertilization in a randomised complete block design. Initial soil nutrition and soil fertility parameters included (hot water extractable organic carbon (C) and total nitrogen (N), total C and N, C/N ratio, labile N pools (nitrate (NO3 −) and ammonium (NH4 +)), extractable potassium (K+)), soil δ15N and δ13C. Relative weed cover, foliar N concentrations, tree growth rate and physiological parameters including photosynthesis, stomatal conductance, photosynthetic nitrogen use efficiency, foliar δ15N and foliar δ13C were also measured at early establishment. Results and discussion Foliar N concentration at 1.25 years was significantly different amongst the weed control treatments and was negatively correlated to the relative weed cover at 1.1 years. Foliar N concentration was also positively correlated to foliar δ15N and foliar δ13C, tree height, height growth rates and tree biomass. Foliar δ15N was negatively correlated to the relative weed cover at 0.8 and 1.1 years. The physiological measurements indicated that luxury fertilization and increasing weed competition on these soils decreased leaf xylem pressure potential (Ψxpp) when compared to the other treatments. Conclusions These results indicate how increasing N resources and weed competition have implications for tree N and water use at establishment in F1 hybrid plantations of southeast Queensland, Australia. These results suggest the desirability of weed control, in the inter-planting row, in the first year to maximise site N and water resources available for seedling growth. It also showed the need to avoid over-fertilisation, which interfered with the balance between available N and water on these soils.

Modelling the nitrogen-driven trade-off between nitrogen utilisation efficiency and water use efficiency of wheat in eastern Australia.

The nitrogen-driven trade-off between nitrogen utilisation efficiency (yield per unit nitrogen uptake) and water use efficiency (yield per unit evapotranspiration) is widespread and results from well established, multiple effects of nitrogen availability on the water, carbon and nitrogen economy of crops. Here we used a crop model (APSIM) to simulate the yield, evapotranspiration, soil evaporation and nitrogen uptake of wheat, and analysed yield responses to water, nitrogen and climate using a framework analogous to the rate-duration model of determinate growth. The relationship between modelled grain yield (Y) and evapotranspiration (ET) was fitted to a linear-plateau function to derive three parameters: maximum yield (Ymax), the ET break-point when yield reaches its maximum (ET#), and the rate of yield response in the linear phase ([Delta]Y/[Delta]ET). Against this framework, we tested the hypothesis that nitrogen deficit reduces maximum yield by reducing both the rate ([Delta]Y/[Delta]ET) and the range of yield response to evapotranspiration, i.e. ET# - Es, where Es is modelled median soil evaporation. Modelled data reproduced the nitrogen-driven trade-off between nitrogen utilisation efficiency and water use efficiency in a transect from Horsham (36°S) to Emerald (23°S) in eastern Australia. Increasing nitrogen supply from 50 to 250 kg N ha-1 reduced yield per unit nitrogen uptake from 29 to 12 kg grain kg-1 N and increased yield per unit evapotranspiration from 6 to 15 kg grain ha-1 mm-1 at Emerald. The same increment in nitrogen supply reduced yield per unit nitrogen uptake from 30 to 25 kg grain kg-1 N and increased yield per unit evapotranspiration from 6 to 25 kg grain ha-1 mm-1 at Horsham. Maximum yield ranged from 0.9 to 6.4 t ha-1. Consistent with our working hypothesis, reductions in maximum yield with nitrogen deficit were associated with both reduction in the rate of yield response to ET and compression of the range of yield response to ET. Against the notion of managing crops to maximise water use efficiency in low rainfall environments, we emphasise the trade-off between water use efficiency and nitrogen utilisation efficiency, particularly under conditions of high nitrogen-to-grain price ratio. The rate-range framework to characterise the relationship between yield and evapotranspiration is useful to capture this trade-off as the parameters were responsive to both nitrogen supply and climatic factors.

Relative contributions of light interception and radiation use efficiency to the reduction of maize productivity under cold temperatures

Maize (Zea mays L.) is a chill-susceptible crop cultivated in northern latitude environments. The detrimental effects of cold on growth and photosynthetic activity have long been established. However, a general overview of how important these processes are with respect to the reduction of productivity reported in the field is still lacking. In this study, a model-assisted approach was used to dissect variations in productivity under suboptimal temperatures and quantify the relative contributions of light interception (PARc) and radiation use efficiency (RUE) from emergence to flowering. A combination of architectural and light transfer models was used to calculate light interception in three field experiments with two cold-tolerant lines and at two sowing dates. Model assessment confirmed that the approach was suitable to infer light interception. Biomass production was strongly affected by early sowings. RUE was identified as the main cause of biomass reduction during cold events. Furthermore, PARc explained most of the variability observed at flowering, its relative contributions being more or less important according to the climate experienced. Cold temperatures resulted in lower PARc, mainly because final leaf length and width were significantly reduced for all leaves emerging after the first cold occurrence. These results confirm that virtual plants can be useful as fine phenotyping tools. A scheme of action of cold on leaf expansion, light interception and radiation use efficiency is discussed with a view towards helping breeders define relevant selection criteria. This paper originates from a presentation at the 5th International Workshop on Functional–Structural Plant Models, Napier, New Zealand, November 2007.

Stretching water - Queensland’s water use efficiency cotton and grains adoption program

The Cotton and Grain Adoption Program of the Queensland Rural Water Use Efficiency Initiative is targeting five major irrigation regions in the state with the objective to develop better irrigation water use efficiency (WUE) through the adoption of best management practices in irrigation. The major beneficiaries of the program will be industries, irrigators and local communities. The benefits will flow via two avenues: increased production and profit resulting from improved WUE and improved environmental health as a consequence of greatly reduced runoff of irrigation tailwater into rivers and streams. This in turn will reduce the risk of nutrient and pesticide contamination of waterways. As a side effect, the work is likely to contribute to an improved public image of the cotton and grain industries. In each of the five regions, WUE officers have established grower groups to assist in providing local input into the specific objectives of extension and demonstration activities. The groups also assist in developing growers' perceptions of ownership of the work. Activities are based around four on-farm demonstration sites in each region where irrigation management techniques and hardware are showcased. A key theme of the program is monitoring water use. This is applied both to on-farm storage and distribution as well as to application methods and in-field management. This paper describes the project, its activities and successes.

Effect of irrigation amounts applied with subsurface drip irrigation on corn evapotranspiration, yield, water use efficiency, and dry matter production in a semiarid climate

Quantifying the local crop response to irrigation is important for establishing adequate irrigation management strategies. This study evaluated the effect of irrigation applied with subsurface drip irrigation on field corn (Zea mays L.) evapotranspiration (ETc), yield, water use efficiencies (WUE = yield/ETc, and IWUE = yield/irrigation), and dry matter production in the semiarid climate of west central Nebraska. Eight treatments were imposed with irrigation amounts ranging from 53 to 356 mm in 2005 and from 22 to 226 mm in 2006. A soil water balance approach (based on FAO-56) was used to estimate daily soil water and ETc. Treatments resulted in seasonal ETc of 580-663 mm and 466-656 mm in 2005 and 2006, respectively. Yields among treatments differed by as much as 22% in 2005 and 52% in 2006. In both seasons, irrigation significantly affected yields, which increased with irrigation up to a point where irrigation became excessive. Distinct relationships were obtained each season. Yields increased linearly with seasonal ETc (R 2 = 0.89) and ETc/ETp (R 2 = 0.87) (ETp = ETc with no water stress). The yield response factor (ky), which indicates the relative reduction in yield to relative reduction in ETc, averaged 1.58 over the two seasons. WUE increased non-linearly with seasonal ETc and with yield. WUE was more sensitive to irrigation during the drier 2006 season, compared with 2005. Both seasons, IWUE decreased sharply with irrigation. Irrigation significantly affected dry matter production and partitioning into the different plant components (grain, cob, and stover). On average, the grain accounted for the majority of the above-ground plant dry mass (≈59%), followed by the stover (≈33%) and the cob (≈8%). The dry mass of the plant and that of each plant component tended to increase with seasonal ETc. The good relationships obtained in the study between crop performance indicators and seasonal ETc demonstrate that accurate estimates of ETc on a daily and seasonal basis can be valuable for making tactical in-season irrigation management decisions and for strategic irrigation planning and management.

The limit to wheat water-use efficiency in eastern Australia. I. Gradients in the radiation environment and atmospheric demand

In the wheatbelt of eastern Australia, rainfall shifts from winter dominated in the south (South Australia, Victoria) to summer dominated in the north (northern New South Wales, southern Queensland). The seasonality of rainfall, together with frost risk, drives the choice of cultivar and sowing date, resulting in a flowering time between October in the south and August in the north. In eastern Australia, crops are therefore exposed to contrasting climatic conditions during the critical period around flowering, which may affect yield potential, and the efficiency in the use of water (WUE) and radiation (RUE). In this work we analysed empirical and simulated data, to identify key climatic drivers of potential water- and radiation-use efficiency, derive a simple climatic index of environmental potentiality, and provide an example of how a simple climatic index could be used to quantify the spatial and temporal variability in resource-use efficiency and potential yield in eastern Australia. Around anthesis, from Horsham to Emerald, median vapour pressure deficit (VPD) increased from 0.92 to 1.28 kPa, average temperature increased from 12.9 to 15.2°C, and the fraction of diffuse radiation (FDR) decreased from 0.61 to 0.41. These spatial gradients in climatic drivers accounted for significant gradients in modelled efficiencies: median transpiration WUE (WUEB/T) increased southwards at a rate of 2.6% per degree latitude and median RUE increased southwards at a rate of 1.1% per degree latitude. Modelled and empirical data confirmed previously established relationships between WUEB/T and VPD, and between RUE and photosynthetically active radiation (PAR) and FDR. Our analysis also revealed a non-causal inverse relationship between VPD and radiation-use efficiency, and a previously unnoticed causal positive relationship between FDR and water-use efficiency. Grain yield (range 1-7 t/ha) measured in field experiments across South Australia, New South Wales, and Queensland (n = 55) was unrelated to the photothermal quotient (Pq = PAR/T) around anthesis, but was significantly associated (r2 = 0.41, P < 0.0001) with newly developed climatic index: a normalised photothermal quotient (NPq = Pq . FDR/VPD). This highlights the importance of diffuse radiation and vapour pressure deficit as sources of variation in yield in eastern Australia. Specific experiments designed to uncouple VPD and FDR and more mechanistic crop models might be required to further disentangle the relationships between efficiencies and climate drivers.

The limit to wheat water-use efficiency in eastern Australia. II. Influence of rainfall patterns

We investigated the influence of rainfall patterns on the water-use efficiency of wheat in a transect between Horsham (36°S) and Emerald (23°S) in eastern Australia. Water-use efficiency was defined in terms of biomass and transpiration, WUEB/T, and grain yield and evapotranspiration, WUEY/ET. Our working hypothesis is that latitudinal trends in WUEY/ET of water-limited crops are the complex result of southward increasing WUEB/T and soil evaporation, and season-dependent trends in harvest index. Our approach included: (a) analysis of long-term records to establish latitudinal gradients of amount, seasonality, and size-structure of rainfall; and (b) modelling wheat development, growth, yield, water budget components, and derived variables including WUEB/T and WUEY/ET. Annual median rainfall declined from around 600 mm in northern locations to 380 mm in the south. Median seasonal rain (from sowing to harvest) doubled between Emerald and Horsham, whereas median off-season rainfall (harvest to sowing) ranged from 460 mm at Emerald to 156 mm at Horsham. The contribution of small events (≤ 5 mm) to seasonal rainfall was negligible at Emerald (median 15 mm) and substantial at Horsham (105 mm). Power law coefficients (τ), i.e. the slopes of the regression between size and number of events in a log-log scale, captured the latitudinal gradient characterised by an increasing dominance of small events from north to south during the growing season. Median modelled WUEB/T increased from 46 kg/ha.mm at Emerald to 73 kg/ha.mm at Horsham, in response to decreasing atmospheric demand. Median modelled soil evaporation during the growing season increased from 70 mm at Emerald to 172 mm at Horsham. This was explained by the size-structure of rainfall characterised with parameter τ, rather than by the total amount of rainfall. Median modelled harvest index ranged from 0.25 to 0.34 across locations, and had a season-dependent latitudinal pattern, i.e. it was greater in northern locations in dry seasons in association with wetter soil profiles at sowing. There was a season-dependent latitudinal pattern in modelled WUEY/ET. In drier seasons, high soil evaporation driven by a very strong dominance of small events, and lower harvest index override the putative advantage of low atmospheric demand and associated higher WUEB/T in southern locations, hence the significant southwards decrease in WUEY/ET. In wetter seasons, when large events contribute a significant proportion of seasonal rain, higher WUEB/T in southern locations may translate into high WUEY/ET. Linear boundary functions (French-Schultz type models) accounting for latitudinal gradients in its parameters, slope, and x-intercept, were fitted to scatter-plots of modelled yield v. evapotranspiration. The x-intercept of the model is re-interpreted in terms of rainfall size structure, and the slope or efficiency multiplier is described in terms of the radiation, temperature, and air humidity properties of the environment. Implications for crop management and breeding are discussed.

Effect of timing of a deficit-irrigation allocation on corn evapotranspiration, yield, water use efficiency and dry mass

Water regulations have decreased irrigation water supplies in Nebraska and some other areas of the USA Great Plains. When available water is not enough to meet crop water requirements during the entire growing cycle, it becomes critical to know the proper irrigation timing that would maximize yields and profits. This study evaluated the effect of timing of a deficit-irrigation allocation (150 mm) on crop evapotranspiration (ETc), yield, water use efficiency (WUE = yield/ETc), irrigation water use efficiency (IWUE = yield/irrigation), and dry mass (DM) of corn (Zea mays L.) irrigated with subsurface drip irrigation in the semiarid climate of North Platte, NE. During 2005 and 2006, a total of sixteen irrigation treatments (eight each year) were evaluated, which received different percentages of the water allocation during July, August, and September. During both years, all treatments resulted in no crop stress during the vegetative period and stress during the reproductive stages, which affected ETc, DM, yield, WUE and IWUE. Among treatments, ETc varied by 7.2 and 18.8%; yield by 17 and 33%; WUE by 12 and 22%, and IWUE by 18 and 33% in 2005 and 2006, respectively. Yield and WUE both increased linearly with ETc and with ETc/ETp (ETp = seasonal ETc with no water stress), and WUE increased linearly with yield. The yield response factor (ky) averaged 1.50 over the two seasons. Irrigation timing affected the DM of the plant, grain, and cob, but not that of the stover. It also affected the percent of DM partitioned to the grain (harvest index), which increased linearly with ETc and averaged 56.2% over the two seasons, but did not affect the percent allocated to the cob or stover. Irrigation applied in July had the highest positive coefficient of determination (R2) with yield. This high positive correlation decreased considerably for irrigation applied in August, and became negative for irrigation applied in September. The best positive correlation between the soil water deficit factor (Ks) and yield occurred during weeks 12-14 from crop emergence, during the "milk" and "dough" growth stages. Yield was poorly correlated to stress during weeks 15 and 16, and the correlation became negative after week 17. Dividing the 150 mm allocation about evenly among July, August and September was a good strategy resulting in the highest yields in 2005, but not in 2006. Applying a larger proportion of the allocation in July was a good strategy during both years, and the opposite resulted when applying a large proportion of the allocation in September. The different results obtained between years indicate that flexible irrigation scheduling techniques should be adopted, rather than relying on fixed timing strategies.

Measuring and predicting canopy nitrogen nutrition in wheat using a spectral index—The canopy chlorophyll content index (CCCI).

Varying the spatial distribution of applied nitrogen (N) fertilizer to match demand in crops has been shown to increase profits in Australia. Better matching the timing of N inputs to plant requirements has been shown to improve nitrogen use efficiency and crop yields and could reduce nitrous oxide emissions from broad acre grains. Farmers in the wheat production area of south eastern Australia are increasingly splitting N application with the second timing applied at stem elongation (Zadoks 30). Spectral indices have shown the ability to detect crop canopy N status but a robust method using a consistent calibration that functions across seasons has been lacking. One spectral index, the canopy chlorophyll content index (CCCI) designed to detect canopy N using three wavebands along the "red edge" of the spectrum was combined with the canopy nitrogen index (CNI), which was developed to normalize for crop biomass and correct for the N dilution effect of crop canopies. The CCCI-CNI index approach was applied to a 3-year study to develop a single calibration derived from a wheat crop sown in research plots near Horsham, Victoria, Australia. The index was able to predict canopy N (g m-2) from Zadoks 14-37 with an r2 of 0.97 and RMSE of 0.65 g N m-2 when dry weight biomass by area was also considered. We suggest that measures of N estimated from remote methods use N per unit area as the metric and that reference directly to canopy %N is not an appropriate method for estimating plant concentration without first accounting for the N dilution effect. This approach provides a link to crop development rather than creating a purely numerical relationship. The sole biophysical input, biomass, is challenging to quantify robustly via spectral methods. Combining remote sensing with crop modelling could provide a robust method for estimating biomass and therefore a method to estimate canopy N remotely. Future research will explore this and the use of active and passive sensor technologies for use in precision farming for targeted N management.

Nitrogen use for high productivity and sustainability in cashew.

Interest in cashew production in Australia has been stimulated by domestic and export market opportunities and suitability of large areas of tropical Australia. Economic models indicate that cashew production is profitable at 2.8 t ha-1 nut-in-shell (NIS). Balanced plant nutrition is essential to achieve economic yields in Australia, with nitrogen (N) of particular importance because of its capacity to modify growth, affect nut yield and cause environmental degradation through soil acidification and off-site contamination. The study on a commercial cashew plantation at Dimbulah, Australia, investigated the effect of N rate and timing on cashew growth, nut production, N leaching and soil chemical properties over five growth cycles (1995-1999). Nitrogen was applied during the main periods of vegetative (December-April) and reproductive (June-October) growth. Commercial NIS yields (up to 4.4 t ha-1 from individual trees) that exceeded the economic threshold of 2.8 t ha-1 were achieved. The yield response was mainly determined by canopy size as mean nut weight, panicle density and nuts per panicle were largely unaffected by N treatments. Nitrogen application confined to the main period of vegetative growth (December-April) produced a seasonal growth pattern that corresponded most consistently with highest NIS yield. This N timing also reduced late season flowering and undesirable post-November nut drop. Higher yields were not produced at N rates greater than 17 g m-2 of canopy surface area (equating to 210 kg N ha-1 for mature size trees). High yields were attained when N concentrations in Mveg leaves in May-June were about 2%, but this assessment occurs at a time when it is not feasible to correct N deficiency. The Mflor leaf of the preceding November, used in conjunction with the Mveg leaf, was proposed as a diagnostic tool to guide N rate decisions. Leaching of nitrate-N and acidification of the soil profile was recorded to 0.9 m. This is an environmental and sustainability hazard, and demonstrates that improved methods of N management are required.

Increase water use efficiency and productivity, and reduce energy and water usage and costs, of dairy and fodder enterprises, to reduce costs of milk production

Increase water use efficiency and productivity, and reduce energy and water usage and costs, of dairy and fodder enterprises, to reduce costs of milk production.