Hybrid metaheuristic techniques for optimising sugarcane rail operations

Mixed integer programming and parallel-machine job shop scheduling are used to solve the sugarcane rail transport scheduling problem. Constructive heuristics and metaheuristics were developed to produce a more efficient scheduling system and so reduce operating costs. The solutions were tested on small and large size problems. High-quality solutions and improved CPU time are the result of developing new hybrid techniques which consist of different ways of integrating simulated annealing and Tabu search techniques.

Balanço de gases de efeito estufa e emissões de CO2 do solo nos sistemas de colheita da cana-de-açúcar manual queimada e mecanizada crua

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Soil and crop residue CO2-C emission under tillage systems in sugarcane-producing areas of southern Brazil

Appropriate management of agricultural crop residues could result in increases on soil organic carbon (SOC) and help to mitigate gas effect. To distinguish the contributions of SOC and sugarcane (Saccharum spp.) residues to the short-term CO2-C loss, we studied the infl uence of several tillage systems: heavy offset disk harrow (HO), chisel plow (CP), rotary tiller (RT), and sugarcane mill tiller (SM) in 2008, and CP, RT, SM, moldboard (MP), and subsoiler (SUB) in 2009, with and without sugarcane residues relative to no-till (NT) in the sugarcane producing region of Brazil. Soil CO2-C emissions were measured daily for two weeks after tillage using portable soil respiration systems. Daily CO2-C emissions declined after tillage regardless of tillage system. In 2008, total CO2-C from SOC and/or residue decomposition was greater for RT and lowest for CP. In 2009, emission was greatest for MP and CP with residues, and smallest for NT. SOC and residue contributed 47% and 41%, respectively, to total CO2-C emissions. Regarding the estimated emissions from sugarcane residue and SOC decomposition within the measurement period, CO2-C factor was similar to sugarcane residue and soil organic carbon decomposition, depending on the tillage system applied. Our approach may define new emission factors that are associated to tillage operations on bare or sugarcane-residue-covered soils to estimate the total carbon loss.

The impact of fallowing and green manuring on soil conditions and the growth of sugarcane

There are currently concerns within some sugar industries that long-term monoculture has led to soil degradation and consequent yield decline. An investigation was conducted in Swaziland to assess the effects of fallowing and green manuring practices, over a seven-month period, on sugarcane yields and the physical properties of a poorly draining clay soil. In the subsequent first sugarcane crop after planting, yields were improved from 129 t ha(-1) under continuous sugarcane to 141-144 t ha(-1) after fallowing and green manuring, but there were no significant responses in the first and second ratoon crops. Also, in the first crop after planting, root length index increased from 3.5 km m(-2) under continuous sugarcane to 5.2-6.8 km m(-2) after fallowing, and improved rooting was still evident in the first ratoon crop where there had been soil drying during the fallow period. Soil bulk density, total porosity and water-holding capacity were not affected by the fallowing practices. However, air-filled porosity increased from 11% under continuous sugarcane to 16% after fallowing, and steady state ponded infiltration rates were increased from 0.61 mm h(-1) to 1.34 mm h(-1), but these improvements were no longer evident after a year back under sugarcane. Levels of soil organic matter were reduced in all cases, probably as a result of the tillage operations involved. In the plant crop, root length was well correlated with air-filled porosity, indicating the importance of improving belowground air supply for crop production on poorly draining clay soils.

Greenhouse gas balance due to the conversion of sugarcane areas from burned to green harvest in Brazil

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Scheduling techniques to optimise rail operations

In Australia, railway systems play a vital role in transporting the sugarcane crop from farms to mills. The sugarcane transport system is very complex and uses daily schedules, consisting of a set of locomotives runs, to satisfy the requirements of the mill and harvesters. The total cost of sugarcane transport operations is very high; over 35% of the total cost of sugarcane production in Australia is incurred in cane transport. Efficient schedules for sugarcane transport can reduce the cost and limit the negative effects that this system can have on the raw sugar production system. There are several benefits to formulating the train scheduling problem as a blocking parallel-machine job shop scheduling (BPMJSS) problem, namely to prevent two trains passing in one section at the same time; to keep the train activities (operations) in sequence during each run (trip) by applying precedence constraints; to pass the trains on one section in the correct order (priorities of passing trains) by applying disjunctive constraints; and, to ease passing trains by solving rail conflicts by applying blocking constraints and Parallel Machine Scheduling. Therefore, the sugarcane rail operations are formulated as BPMJSS problem. A mixed integer programming and constraint programming approaches are used to describe the BPMJSS problem. The model is solved by the integration of constraint programming, mixed integer programming and search techniques. The optimality performance is tested by Optimization Programming Language (OPL) and CPLEX software on small and large size instances based on specific criteria. A real life problem is used to verify and validate the approach. Constructive heuristics and new metaheuristics including simulated annealing and tabu search are proposed to solve this complex and NP-hard scheduling problem and produce a more efficient scheduling system. Innovative hybrid and hyper metaheuristic techniques are developed and coded using C# language to improve the solutions quality and CPU time. Hybrid techniques depend on integrating heuristic and metaheuristic techniques consecutively, while hyper techniques are the complete integration between different metaheuristic techniques, heuristic techniques, or both.

The primary sources of carbon loss during the crop-establishment period in a subtropical Oxisol under contrasting tillage systems

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Soil fragmentation study applying different tillage systems

Runoff generation depends on rainfall, infiltration, interception, and surface depressional storage. Surface depressional storage depends on surface microtopography, usually quantified trough soil surface roughness (SSR). SSR is subject to spatial and temporal changes that create a high variability. In an agricultural environment, tillage operations produce abrupt changes in roughness. Subsequent rainfall gradually decreases roughness. Beside it, local variation in soil properties and hydrology cause its SSR to vary spatially at different scales. The methods commonly used to measure it involve collecting point elevations in regular grids using laser profilers or scanners, digital close range stereo-photogrammetry and terrestrial laser scanning or LIDAR systems. In this case, a laser-scanning instrument was used to obtain representative digital elevation models (DEMs) at a grid resolution of 7.2x7.2mm that cover an area of 0.9x0.9m. The DEMs were obtained from two study sites with different soils. The first study site was an experimental field on which five conventional tillage methods were applied. The second study site was a large olive orchard with trees planted at 7.5x5.0m and bare soils between rows. Here, three tillage treatments were applied. In this work we have evaluated the spatial variability of SSR at several scales studying differences in height calculated from points separated by incremental distances h were raised to power values q (from 0 to 4 in steps of 0.1). The q = 2 data were studied as a semivariogram model. The logarithm of average differences plotted vs. log h were characterized by their slope, ?(q). Structure functions [?(q) vs. q] were fitted showing that data had nonlinear structure functions typical of multiscale phenomena. Comparisson of the two types of soil in their respective structure functions are shown.

Spatial variability patterns of some Vertisol properties at a field scale using standardized data

Spatial variability of Vertisol properties is relevant for identifying those zones with physical degradation. In this sense, one has to face the problem of identifying the origin and distribution of spatial variability patterns. The objectives of the present work were (i) to quantify the spatial structure of different physical properties collected from a Vertisol, (ii) to search for potential correlations between different spatial patterns and (iii) to identify relevant components through multivariate spatial analysis. The study was conducted on a Vertisol (Typic Hapludert) dedicated to sugarcane (Saccharum officinarum L.) production during the last sixty years. We used six soil properties collected from a squared grid (225 points) (penetrometer resistance (PR), total porosity, fragmentation dimension (Df), vertical electrical conductivity (ECv), horizontal electrical conductivity (ECh) and soil water content (WC)). All the original data sets were z-transformed before geostatistical analysis. Three different types of semivariogram models were necessary for fitting individual experimental semivariograms. This suggests the different natures of spatial variability patterns. Soil water content rendered the largest nugget effect (C0 = 0.933) while soil total porosity showed the largest range of spatial correlation (A = 43.92 m). The bivariate geostatistical analysis also rendered significant cross-semivariance between different paired soil properties. However, four different semivariogram models were required in that case. This indicates an underlying co-regionalization between different soil properties, which is of interest for delineating management zones within sugarcane fields. Cross-semivariograms showed larger correlation ranges than individual, univariate, semivariograms (A ≥ 29 m). All the findings were supported by multivariate spatial analysis, which showed the influence of soil tillage operations, harvesting machinery and irrigation water distribution on the status of the investigated area.

A near Optimal Cane Rail Scheduler under Limited and Unlimited Capacity Constraints

Australia is the world’s third largest exporter of raw sugar after Brazil and Thailand, with around $2.0 billion in export earnings. Transport systems play a vital role in the raw sugar production process by transporting the sugarcane crop between farms and mills. In 2013, 87 per cent of sugarcane was transported to mills by cane railway. The total cost of sugarcane transport operations is very high. Over 35% of the total cost of sugarcane production in Australia is incurred in cane transport. A cane railway network mainly involves single track sections and multiple track sections used as passing loops or sidings. The cane railway system performs two main tasks: delivering empty bins from the mill to the sidings for filling by harvesters; and collecting the full bins of cane from the sidings and transporting them to the mill. A typical locomotive run involves an empty train (locomotive and empty bins) departing from the mill, traversing some track sections and delivering bins at specified sidings. The locomotive then, returns to the mill, traversing the same track sections in reverse order, collecting full bins along the way. In practice, a single track section can be occupied by only one train at a time, while more than one train can use a passing loop (parallel sections) at a time. The sugarcane transport system is a complex system that includes a large number of variables and elements. These elements work together to achieve the main system objectives of satisfying both mill and harvester requirements and improving the efficiency of the system in terms of low overall costs. These costs include delay, congestion, operating and maintenance costs. An effective cane rail scheduler will assist the traffic officers at the mill to keep a continuous supply of empty bins to harvesters and full bins to the mill with a minimum cost. This paper addresses the cane rail scheduling problem under rail siding capacity constraints where limited and unlimited siding capacities were investigated with different numbers of trains and different train speeds. The total operating time as a function of the number of trains, train shifts and a limited number of cane bins have been calculated for the different siding capacity constraints. A mathematical programming approach has been used to develop a new scheduler for the cane rail transport system under limited and unlimited constraints. The new scheduler aims to reduce the total costs associated with the cane rail transport system that are a function of the number of bins and total operating costs. The proposed metaheuristic techniques have been used to find near optimal solutions of the cane rail scheduling problem and provide different possible solutions to avoid being stuck in local optima. A numerical investigation and sensitivity analysis study is presented to demonstrate that high quality solutions for large scale cane rail scheduling problems are obtainable in a reasonable time. Keywords: Cane railway, mathematical programming, capacity, metaheuristics

Agricultural intensification and biodiversity partitioning in European landscapes comparing plants, carabids, and birds

Effects of agricultural intensification (AI) on biodiversity are often assessed on the plot scale, although processes determining diversity also operate on larger spatial scales. Here, we analyzed the diversity of vascular plants, carabid beetles, and birds in agricultural landscapes in cereal crop fields at the field (n = 1350), farm (n = 270), and European-region (n = 9) scale. We partitioned diversity into its additive components alpha, beta, and gamma, and assessed the relative contribution of beta diversity to total species richness at each spatial scale. AI was determined using pesticide and fertilizer inputs, as well as tillage operations and categorized into low, medium, and high levels. As AI was not significantly related to landscape complexity, we could disentangle potential AI effects on local vs. landscape community homogenization. AI negatively affected the species richness of plants and birds, but not carabid beetles, at all spatial scales. Hence, local AI was closely correlated to beta diversity on larger scales up to the farm and region level, and thereby was an indicator of farm-and region-wide biodiversity losses. At the scale of farms (12.83-20.52%) and regions (68.34-80.18%), beta diversity accounted for the major part of the total species richness for all three taxa, indicating great dissimilarity in environmental conditions on larger spatial scales. For plants, relative importance of alpha diversity decreased with AI, while relative importance of beta diversity on the farm scale increased with AI for carabids and birds. Hence, and in contrast to our expectations, AI does not necessarily homogenize local communities, presumably due to the heterogeneity of farming practices. In conclusion, a more detailed understanding of AI effects on diversity patterns of various taxa and at multiple spatial scales would contribute to more efficient agri-environmental schemes in agroecosystems.

Assessment of soil particle redistribution on two contrasting cultivated hillslopes

The particle size distributions of surface soils from two cultivated silty fields (Moorfield and Railway South) in Herefordshire, UK, were assessed by sampling on 20-m grids across the fields. Moorfield (8 ha) had a uniform landscape sloping mainly in a North-South direction while Railway South (12 ha) had complex undulating landscape characteristics. Samples from 3 surficial layers were also taken at 3 landscape positions at Moorfield to investigate recent (within-season) soil particle redistribution. Size fractions were determined using chemical dispersion, wet sieving (to separate the sand fractions) and laser gramilometry (for the finer fractions). The distribution of various fractions and the relationships between elevation and the various fractions suggest preferential detachment and movement of coarse to very coarse silt fractions (16-63 mu m), which were found mostly at downslope or depositional areas. Upper slope samples had higher clay to fine silt (< 16 mu m) contents than bottom slope samples. The upslope-downslope patterns of size fractions, particularly on uniformly sloping areas, of the 2 fields were similar and their deposited sediments were dominated by coarse silt fractions. Samples from 3 landscape positions at Moorfield became coarser from the less eroded summit, through the eroding side-slope to the bottom-slope depositional area. Within each of these landscape positions the top 0-2.5 cm layers were more enriched in coarse silt fractions than the bottom layers. The spatial patterns of soil particle size distributions in the 2 fields may be a result of sediment detachment and deposition caused by water erosion and tillage operations. (c) 2005 Elsevier B.V. All rights reserved.