Author Archives: cannabisgrow

The effect of conflict on discriminatory workplace behavior does not decay in the nine months after conflict ended

So far we have seen that output in factory production in Kenya is lower when individuals of different ethnic backgrounds work together, and that the reason appears to be that biased upstream workers under supply downstream workers of other ethnic groups and misallocate intermediate goods across coethnic and non-coethnic downstream workers. We have also seen that distortionary workplace discrimination is greater durings times of conflict, and that firms introduce policies in response in order to reduce workers’ incentive to discriminate. By studying how discriminatory preferences are shaped, and how firms choose their response to distortionary discrimination, researchers can go beyond identifying a source of ethnic diversity effects in production and begin to address why those effects vary across space and time and how profit motives in the private sector can reduce the aggregate effect of ethnic diversity. In the model of taste-based discrimination above, the impact of conflict on output in diverse teams should persist for as long as attitudes towards workers of other ethnic groups are affected. Periods of increased antagonism may entail significant hidden economic costs if “mean reversion” in taste for discrimination is slow . The evolution of output in teams of different ethnicity configurations across the three sample periods was depicted in figure 2. After the introduction of team pay, average output in both homogeneous and mixed teams was steady for the remainder of the sample period, cannabis grow setup suggesting that the impact of conflict on social preferences was long-lived.How did the response to conflict of distortionary discrimination at work vary across individuals?

Modeling θC and θNC as parameter values shared by all workers is a simplification: in reality some workers will have a higher taste for discrimination than others. Figure 9 plots the distribution, across individual suppliers, of the difference in output between homogeneous and mixed teams supplied, before and after conflict began. It appears that most suppliers discriminate against non-coethnic processors during the pre-conflict period. Conflict led to an increase in the output gap between homogeneous and mixed teams supplied for most upstream workers, but also to a notable widening of the distribution of the output gap. The figure indicates that some upstream workers respond more to conflict than others, differentially increasing the extent to which they discriminate against non-coethnics downstream. Some workers in the sample were more exposed to the conflict period of early 2008 than others. Though the workers at the plant and their co-habitating family-members were not themselves directly affected, 22 percent of workers report to have “lost a relative” during the conflict. The decrease in output in mixed teams when conflict began was significantly greater in teams supplied by such workers, as seen in columns 1 and 2 of table 10. These results indicate that personal grievances exacerbate individuals’ workplace response to conflict. Younger individuals may have more malleable social preferences. In columns 3 and 4 of table 10 we see that, although output in homogeneous teams led by old and young suppliers was similar, output in mixed teams with young suppliers was significantly higher during the first year of the sample period. Young suppliers were less discriminatory towards noncoethnic co-workers than old suppliers before conflict began, it appears. This finding is consistent with an expectation expressed by many Kenya commentators before 2008.

It was argued that the young coming of age at the time would be the country’s first “post-tribal” generation . The results of table 10 also show that the decrease in output in mixed teams when conflict began was significantly greater in teams with young suppliers, however. Output in mixed teams with young suppliers was no higher than in mixed teams with older suppliers during the conflict period. These results suggest that youth start out relatively tolerant, but that the attitudes of the young towards non-coethnics respond more negatively to conflict. The results discussed in this section paint a consistent picture of how distortionary attitudes towards workers of other ethnic groups respond to ethnic conflict. It appears that conflict may entail significant hidden economic costs because distortionary social preferences are updated in a “Bayesian” fashion when conflict occurs, at least in the Kenyan context. A serious episode of violent, political conflict between the Kikuyu and Luo blocs led to a significant shift in the average weight attached to the well-being of non-coethnics, a shift that did not decay in the nine months after conflict ended. The negative response was greater among those more affected and among those likely to have a less cemented “prior”.Segregating workers of different ethnic groups would appear to be the profit-maximizing response to distortionary discrimination, from the viewpoint of the econometrician. The results in tables 4 and 8 suggest that segregation would have increased plant productivity by four percent before conflict and by eight percent after conflict began, relative to the status quo of arbitrary assignment to teams. Are these expected benefits of a magnitude that is likely to be salient to supervisors? Consider the output increase expected from optimally assigning workers to teams and positions based on ethnicity, productivity or both. If we view a worker as having three characteristics – the tercile to which she belongs in the distribution of processor productivity, the tercile to which she belongs in the distribution of supplier productivity, and her ethnicity – then an average output will be associated with teams of each of 3 ethnicity configurations, 18 productivity configurations and 63 ethnicity productivity configurations. 

In theory, supervisors can then solve the linear programming problem of maximizing total output subject to the expected output associated with a given type of team and the “budget set” of workers available . The optimal assignments and associated expected output gains are shown in table 11. Throughout the period observed, the output gains expected from assigning workers to teams based on ethnicity were larger than those expected from assigning workers based on productivity – twice as large during the conflict period. In fact segregation achieves about half the output gains of the “complete” solution. The complete solution assigns workers optimally to fully specified teams and thus takes into account interactions between the three workers’ ethnicities and productivities – a complicated “general equilibrium” problem that is likely infeasible for supervisors to solve. It thus appears that the expected productivity gain of segregation is sizable relative to the expected effect of changing other comparable factors under supervisors’ control.It is possible that a similar effect occurs in a Kenyan workplace, although in a situation in which mixed teams are characterized by discriminatory behavior it is also possible that interaction increases tensions and exacerbates ethnic biases. To investigate, I compare the behavior of suppliers with greater versus lower experience working with non-coethnics, in table 12. Focusing on output during the second half of 2007 and the first six weeks of 2008, I contrast teams with suppliers with above-average versus below-average time spent in mixed teams during the first half of 2007. Because most workers at the farm had already spent significant time working with non-coethnics before 2007, columns 3 and 4 restrict the sample to those with below-average tenure. The results show no significant effect of time spent working with non-coethnics on the output gap between mixed and homogeneous teams supplied, vertical grow system neither before nor after conflict began. Workers who have interacted more with individuals of other ethnic groups thus appear no less discriminatory in production. The results in table 12 do not rule out the possibility that complete segregation between the two ethnic groups over time would have a negative influence on attitudes or behavior towards non-coethnics, however. Carrell, Sacerdote, and West find that implementing an estimated optimal assignment can have unintended consequences due to unforeseen responses on the part of individuals to out-of-sample assignments. In the context of the sample farm, in a country that has experienced periodical violent clashes between ethnic groups, and where workers of different ethnic groups reside in the same quarters, complete segregation at the plant could for example lead to increased social tensions on the farm.Nevertheless, it is arguably surprising that a supposedly profit-maximizing firm chose to leave large productivity gains “on the table” by not segregating workers of different ethnicities. Ethical considerations add complexity to the issue of team assignment in Kenya, but we would perhaps expect longer-term costs of segregation to be incurred primarily by society, rather than the firm itself, in which case a case can be made for government intervention to enforce integration within firms. Becker pointed out that discriminatory employers should go out of business as their profits suffer.

A priori, the same argument should hold for flower farms that allow workplace discrimination to influence productivity. However, the floriculture business is not particularly competitive, as evidenced by high profit margins . Moreover, as the literature in macroeconomics on across-firm misallocation has highlighted, it is not necessarily the most productive firms that survive in poor countries’ economies . Further, plant managers did respond to the increase in distortionary discrimination when conflict began, as we have seen. The introduction of team pay for processors was likely motivated by the decrease in productivity in diverse teams in early 2008. It is unsurprising that the dramatic differential decrease in mixed teams’ output when conflict began led managers to respond, even though the lower output observed in diverse teams during the first year of the sample period did not. A doubling of the output gap of diverse teams during a short period of time is likely more salient to managers than potential foregone productivity gains from arbitrary assignment to teams. It appears that managers considered an adjustment to contractual incentives a more desirable response to distortionary discrimination than segregating workers. But note that it is likely not possible to eliminate discrimination through contractual incentives, without entirely breaking the link between workers’ output and pay. At the sample plant, vertical discrimination continued to significantly affect output after the introduction of team pay.Evidence suggests that ethnic diversity negatively affects public goods provision and the quality of macroeconomic policies. While the possibility of an additional, direct effect on micro-level productivity has long been recognized, corresponding evidence is largely absent. In this paper, I begin by identifying a sizable, negative productivity effect of ethnic diversity in teams in Kenya. I do so using two years of daily output data for 924 workers, almost equally drawn from two rival tribes, at a flower-packing plant. The packing process takes place in triangular production units, one upstream “supplier” supplying two downstream “processors” who finalize bunches of flowers. I show that an arbitrary position rotation system led to quasi-random variation in teams’ ethnicity configuration. As predicted by a model in which different weight is attached to coethnic and non-coethnic downstream workers’ utility, suppliers discriminate both “vertically” – undersupplying downstream noncoethnics – and “horizontally” – shifting flowers from non-coethnic to coethnics downstream workers. By doing so, upstream workers lower their own pay and total output. I show that less distortionary, non-taste-based ethnic diversity effects are unlikely to explain this paper’s results. As Becker points out, significant aggregate effects “could easily result from the manner in which individual tastes for discrimination allocate resources within a free-enterprise framework” . Discrimination should lead to misallocation of resources in most joint production situations in which individuals influence the output and income of others. I take advantage of two natural experiments during the time period observed to begin to explore how the productivity effects of ethnic diversity are likely to vary across time and space. When contentious presidential election results led to political conflict and violent clashes between the two ethnic groups represented in the sample in early 2008, a dramatic, differential decrease in the output of mixed teams followed, as predicted by themodel. The reason appears to be that workers’ taste for discrimination against non-coethnic co-workers increased. I estimate a decrease in the weight attached to non-coethnics’ utility of approximately 35 percent in early 2008, through a reduced form approach. A back-of the-envelope calculation suggests that the increase in distortionary workplace discrimination may have cost the plant half a million dollars in annual profit, had it not responded. Six weeks into the conflict period, the plant implemented a new pay system in which downstream workers were paid for their combined output .

The effects of ethnic divisions are of particular importance in the Kenyan context

The best method for determining incorporation timing is to walk the field with a shovel and dig numerous holes and “feel” soil moisture at various depths throughout the field. In medium- and heavy textured soils you want to be able to form a ball of soil in your hand and then break it apart easily. For this exercise it is important to get soil from at least 8 inches deep. If the soil “ribbons” easily when squeezed between your thumb and index finger it is probably still too wet to work . Optimum soil moisture is critical for good incorporation and breakdown; in average rainfall years early April is commonly the best time for incorporation in the Central Coast region. After flail mowing the residue needs to be mixed with the soil to enhance microbial breakdown and facilitate seedbed formation. The best tool for this is a mechanical spader. Spaders are ideal for cover crop incorporation for many reasons. When operated in optimal soil moisture conditions spaders have minimal impact on soil aggregation and create almost no compaction compared to other primary tillage tools. Spaders are capable of uniformly mixing the cover crop residue into the tilled zone while at the same time leaving the soil lofted and well aerated, allowing for ideal conditions for microbial breakdown of the residue. Spaders also have two major drawbacks: they are expensive, cannabis grow equipment and they require very slow gearing and high horse power to operate; 10 horse power per working foot of spader is the basic requirement depending on soil conditions and depth of operation. They run at a very slow ground speed, often in the range of 0.6 to 0.8 mph. Thus a 7-foot wide spader requires 70 HP and takes between 3 and 4 hours to spade an acre.

Although time consuming, the results are impossible to replicate with any other tillage options now available. If a mechanical spader is not available the next best and probably most commonly used tool for cover crop incorporation is a heavy offset wheel disc. Depending on the size and weight of the disc multiple passes are often required for adequate incorporation. Chiseling after the first several passes will facilitate the disc’s ability to turn soil and will also help break up compaction from the disc.Rapid urbanization, compounded by globalization, has had lasting effects on the agricultural sector and on both urban and rural communities: As urban populations increase, they place more demands on a shrinking group of rurally-based food suppliers. And as the movement for locally based food systems grows to address this and other food system problems, urban agriculture has become a focal point for discussion, creativity, and progress. Indeed, the production of food in and around densely populated cities bears much promise as part of any solution to food supply and access issues for urban populations. Growing Power, Inc. , a Milwaukee-based organization, exemplifies the way that urban agriculture can address some of the needs of rapidly growing urban communities, particularly those with poor, minority populations. Will Allen, Growing Power’s founder—born to former sharecroppers in 1949—was drawn back to agriculture after a career in professional basketball. Aside from his love for growing food, he saw that the mostly poor, black community near his roadside stand in North Milwaukee had limited access to fresh vegetables or to vegetables they preferred. Confirming his observation, a 2006 study found that more diverse food options exist in wealthy and white neighborhoods than in poor and minority neighborhoods.2 Allen decided that he would serve this unmet demand by growing fresh food in the neighborhood where his customers lived and involve the community in the process. In 1993, long before urban agriculture bloomed into the movement it is today, Growing Power began. The importance of equal access to fresh food cannot be overestimated.

While ever more exotic fruits and vegetables from around the world stock health and natural foods stores in wealthy and predominantly white neighborhoods, poor and minority communities face fewer and less healthy food choices in the form of convenience stores, fast food restaurants, and disappearing supermarkets. This lack of access can lead to higher rates of diet related illnesses . Growing Power has been working to create an alternative food system based on intensive fruit and vegetable production, fish raising, and composting, in order to make healthy food available and affordable to the surrounding community, and to provide community members with some control over their food choices. But as anyone who has initiated an urban agricultural project knows, fertile, uncontaminated land is often difficult to find in a city. Even if land with soil is available, most empty lots are in former industrial areas where toxic contamination often renders land unusable . In Milwaukee, Growing Power sat on a lot with no soil and five abandoned greenhouses. Compost became the foundation for all of Growing Power’s activities. The raw materials needed to produce it were in abundant and cheap supply in the city— food waste, brewery grains, coffee grounds, newspaper waste, grass clippings, and leaf mold are all by-products of urban life destined, in most places, for the landfill. Businesses will often donate these materials to urban agriculture projects, saving the cost of garbage hauling services. Compost, and vermicompost in particular, also provides a renewable source of fertilizer that doesn’t rely on fossil-fuel inputs and can itself be used as a growing media. With a healthy compost-based system, Growing Power discovered a low-cost, renewable, and easy to-duplicate solution to one of the biggest hurdles people face when growing food in cities. Since 1993, Growing Power has grown in size and scope, starting gardens in Chicago as well as Milwaukee, and training centers in 15 cities, and including youth training, outreach and education, and policy initiatives in its mission. Interest in urban agriculture has also blossomed into a movement that includes commercial urban farms, scores of community farms and gardens, and educational gardens and training programs growing food and flowers and raising chickens, bees, goats, and other livestock for local consumption.

Urban agriculture has grown so rapidly in the last two decades that in 2012, the USDA granted $453,000 to Penn State University and New York University for a nationwide survey of the “State of Urban Agriculture”3 with an eye toward providing technical assistance, evaluating risk management, and removing barriers for urban farmers. The federal government’s interest in urban agriculture comes on the heels of state and local initiatives to encourage urban agriculture in numerous cities, including Milwaukee, Chicago, New York and San Francisco. The driving force behind these initiatives and the urban agriculture movement as a whole has always been groups of committed individuals in urban communities in search of food, community, opportunity, security, and access. What makes the Growing Power model work is not just its innovative techniques and creative use of urban spaces, but the partnership with its neighbors who not only receive the program’s services, vertical grow rack but contribute significantly to its success.There is evidence to suggest that ethnic heterogeneity may impede economic growth. A negative influence on decision-making in the public sphere has been documented: public goods provision is lower and macroeconomic policies of lower quality in ethnically fragmented societies . The possibility of an additional direct effect on productivity in the private sector has long been recognized, however. Individuals of different ethnicities may have different skill-sets and therefore complement each other in production, but it is also possible that workers of the same ethnic background collaborate more effectively . Evidence from poor countries on the productivity effects of ethnic diversity is largely absent. This paper provides novel microeconometric evidence on the productivity effects of ethnic divisions. I identify a negative effect of ethnic diversity on output in the context of joint production at a large plant in Kenya where workers were quasi-randomly assigned to teams. I then begin to address how output responds to increased conflict between ethnic groups, how firms respond to lower productivity in diverse teams, and how workplace behavior responds to policies implemented by firms to limit ethnic diversity distortions. A model of taste-based discrimination at work explains my findings across these dimensions. I study a sample of 924 workers working in teams at a plant in Kenya. The workers package flowers and prepare them for shipping: productivity is observed and measured by daily individual output. Tribal competition for political power and economic resources has been a defining character of Kenyan society since independence . Workers at the flower plant are almost equally drawn from two historically antagonistic ethnic blocs – the Kikuyu and the Luo . Production takes place in triangular packing units. One upstream “supplier” supplies and arranges roses that are then passed on to two downstream “processors” who assemble the flowers into bunches, as illustrated in figure 1a. The output of each of the two processors is observed. During the first period of the sample, processors were paid a piece rate based on own output and suppliers a piece rate based on total team output. Inefficiently low supply of roses to downstream workers of the rival ethnic group was thus costly for suppliers. I show that the plant’s system of assigning workers to positions through a rotation process generates quasi-random variation in team composition.

A worker’s past productivity and observable characteristics are orthogonal to those of other workers in her assigned team. The productivity effect of ethnic diversity can thus be identified by comparing the output of teams of different compositions. Two natural experiments during the time period for which I have data allow me to go further. During the second period of the sample, in early 2008, contentious presidential election results led to political and violent conflict between the Kikuyu and Luo ethnic groups, but production at the plant continued as usual. In the third period of the sample, starting six weeks after conflict began, the plant implemented a new pay system in which processors were paid for their combined output . By taking advantage of the three periods observed, I identify the source of productivity effects of ethnic diversity in the context of plant production in Kenya; how the economic costs of ethnic diversity vary with the political and social environment; and how managers responded to ethnic diversity distortions at the plant, and how workplace behavior changed as a consequence of the policies implemented in response. I model ethnic diversity effects as arising from a “taste for discrimination” among upstream workers: suppliers attach a potentially differential weight to coethnics’ and noncoethnics’ utility, a formulation that follows Becker , Charness and Rabin and others. The model predicts that discriminatory suppliers in mixed teams will “misallocate” flowers both vertically – under supplying downstream workers of the other ethnic group – and horizontally – shifting flowers from non-coethnic to coethnic downstream workers.1 The impact of horizontal misallocation on total output will depend on the relative productivity of favored and non-favored downstream workers. If conflict led to a decrease in non-coethnics’ utility-weight, a differential fall in mixed teams’ output in early 2008 is predicted. Under team pay, a positive output effect of a reduction in horizontal misallocation is expected to offset negative free riding effects, in teams in which the two processors are of different ethnic groups. The reason is that suppliers can no longer influence the relative pay of the two processors through relative supply under team pay. Quasi-random assignment led to teams of three different ethnicity configurations. About a quarter of observed teams are ethnically homogeneous, another quarter are “vertically mixed” teams in which both processors are of a different ethnic group than the supplier, and about half are “horizontally mixed” teams in which one processor is of a different ethnic group than the supplier. The ethnicity configurations are displayed in figure 1b. I test the model’s predictions by comparing the average output of teams of different ethnicity configurations within and across the three sample periods. In the first main result of the paper, I find that vertically mixed teams were eight percent less productive and horizontally mixed teams five percent less productive than homogeneous teams during the first period of the sample.

There are many commercial mixes available that come close to meeting most of the above criteria

Synthetic nitrogen-based fertilizers were made possible because of the Haber-Bosch process, which converts stable, inert nitrogen gas unavailable to plants into the reactive ammonia molecule readily available for plant uptake. Once the process was commercialized, synthetic fertilizer use skyrocketed, as farmers were no longer dependent only on their soil organic matter, compost, cover crops, and livestock manure for nitrogen. Fertilizer use in the United States increased from about 7.5 million tons in 1960 to 21 million tons in 2010. In 2007, California farmers applied 740,00 tons of nitrogen in fertilizers to 6.7 million acres of irrigated farmland. With cheap sources of nitrogen and water available, our current agricultural system is based on the liberal application of synthetic fertilizers and irrigation water to ensure high yields, often at the expense of environmental and public health. California’s Central Valley is home to some of the most heavily fertilized cropland and some of the most polluted water in the United States. Communities there are particularly vulnerable to public health effects of nitrate contamination because groundwater provides drinking water for the majority of residents. Additionally, rural communities in the valley are generally poor and populated by immigrants and minorities least able to afford treatment costs and most vulnerable to discriminatory decision-making. Tulare County, the second most productive agricultural county in California, includes many of these communities. Though it generates nearly $5 billion in revenue from agriculture each year, it has the highest poverty rate in California and is populated mainly by minorities , most of whomare Latino. The average per capita income in the county is $18,021. Here, indoor cannabis growing one in five small public water systems and two in five private domestic wells surpass the maximum contaminant level for nitrates.

As a result, residents of towns like Seville, East Orosi, and Tooleville are paying $60 per month for nitrate-contaminated water they can’t safely use, and must spend an additional $60 to purchase bottled water for drinking and bathing. In contrast, San Francisco water customers pay $26 per month for pristine water from the Hetch Hetchy water system in Yosemite. The economic cost of nitrate contamination in drinking water is not the only cost to these communities. Farm workers make up a significant segment of the population of small towns throughout the Central Valley and are both directly exposed to the hazards of heavy fertilizer use in the fields and in the air, and through excess nitrogen leached into groundwater drinking supplies. Scientists estimate that 50–80% of nitrogen applied in fertilizer is unused by plants. Of that, about 25% volatilizes into the atmosphere . As a result, approximately 30–50% of nitrogen applied in fertilizer—about 80 pounds per acre in California—leaches into groundwater beneath irrigated lands and into public and private water supplies.13 High nitrate levels in water can cause a number of health problems, including skin rashes, eye irritation, and hair loss. More severe is “Blue Baby Syndrome” , a potentially fatal blood disorder in infants caused by consumption of nitrate-contaminated water. Direct ingestion, intake through juices from concentrate, and bottle-fed infant formula are all potential threats to children. Nitrate contamination has also been linked to thyroid cancer in women. Widespread contamination of groundwater through leached fertilizer has rendered drinking water in rural communities across the country not only unusable, but dangerously so. While nitrate contamination is an acute problem in California, it exists across the country. The EPA estimates that over half of all community and domestic water wells have detectable levels of nitrates. Rural communities that rely on private wells , or lack access to adequate water treatment facilities, have the most insecure water supplies. In the short term, municipalities must devise a plan to reduce the disproportionately high cost of water to these communities.

One potential solution is a fee attached to the purchase of fertilizer used to subsidize water costs for communities with contaminated water. Communities with contaminated water could also be added to a nearby water district with access to clean water. In the longer term, the obvious solution is to substantially reduce synthetic fertilizer and water use in agriculture. Treatment, while effective on a small scale, cannot keep up with the vast quantities of nitrates continually entering groundwater supplies through fertilizer application. Similarly, reduced irrigation on farms, drawn mostly from uncontaminated sources, frees up new sources of drinking water for nearby communities. Lastly, to reach a truly sustainable and equitable system of water distribution, residents of rural communities must be included in the planning and decision-making process as members of local water boards, irrigation districts, and planning commissions to establish and safeguard their right to uncontaminated water.Currently, there is mounting evidence that suggests sustainable agriculture practices, exemplified by those used in agroecological systems, provide an opportunity to achieve the dual goals of feeding a growing population and shrinking agriculture’s carbon footprint, in addition to the social benefits of increased food security and stronger rural economies. This is in contrast with industrial-scale conventional systems that rely on fossil fuel-based fertilizers, pesticides, and heavy tillage and look to genetic engineering to help plants cope with climate change, e.g. by developing drought-resistant crop varieties, which themselves require high inputs of fertilizers and pesticides to produce optimally. Agroecological systems, on the other hand, can mitigate climate change by reducing fossil fuel use, and employing farming techniques that reduce GHG emissions by sequestering carbon in the soil. Of the range of practices in an agroecological system that address issues related to climate change, cover cropping is perhaps the most effective. As climate change continues to affect weather patterns and cause more frequent and severe weather events, protecting against soil erosion will become increasingly important. Cover crops provide an effective mitigation strategy by protecting soil against wateror wind-driven erosion. Cover cropping also provides other climaterelated benefits, including: an on-farm source of fertility, less dependence on fossil fuels and their derived products, and adaptability and resilience. Most of all, while the specific species, timing, and primary purpose of a cover crop vary geographically, the principles behind their cultivation are universally applicable and their benefits universally available. The use of a leguminous cover crop to fix nitrogen in the soil over the wet season for the next season’s crop is widely recognized as an effective fertility management tool. According to an FAO report on agriculture in developing countries, using cover crops in a maize/pigeon pea rotation led to increased yields and required less labor for weeding than continuous maize cropping systems with conventional fertilizer use.3 Nitrogen-fixing cover crops also greatly reduce, and in some cases eliminate, reliance on off-farm sources of fertility, thus reducing the overall carbon footprint of the farm while maintaining high fertility levels in the soil. Note that even organic fertilizers have a high embedded energy cost as they are mostly derived from manure from animals raised in confined feedlots, so the ability to grow one’s fertility needs on farm is important across different agricultural systems. Cover crops are not only a mitigation strategy for climate change, but also a cost-saving measure. Synthetic fertilizer costs have steadily increased over the last half-century, cannabis growing supplies causing hardship for farmers in developing countries especially where fertilizer prices are already two to three times the world price. Organic farmers are less vulnerable to price shifts in fertilizer, but can equally benefit from the reduced need for compost as a result of cover cropping. By saving seeds from their cover crops, farmers can close the loop in their cover crop management, save on annually purchased seed, and develop strains well-adapted to local conditions.

Fertility management systems based on cover crops insulate conventional farmers from increasingly frequent spikes in fertilizer prices and provide organic farmers with a cheap and renewable source of fertility. Adaptation and resilience are also crucial to farmers’ long-term success in the face of unpredictable and disruptive effects from a changing climate because so much of agriculture depends on constantly changing climatic conditions. Added to climate change are increasing input prices and a growing demand for food that put pressure on farmers to maintain high yields while paring down on costs. Cover crops can provide farmers with the flexibility they need by protecting topsoil from wind and water erosion, storing a reliable supply of nutrients to the soil, and—if managed correctly— minimizing costly weeding requirements. For many resource poor farmers who maintain livestock, cover crops provide a path to financial independence and food security as they can be grown both for soil fertility and livestock feed. Cover crops as part of a climate mitigation strategy also make sense at every scale of agriculture. Large conventional farms require consistently high yields to stay profitable as they often operate on razor thin margins. To achieve this goal, these farms rely heavily on fossil fuel-based sources of energy and fertility. Whether used on conventional or organic farms, cover cropping not only reduces farm emissions, but also contributes to the biological health of the farm’s aggressively cultivated soils. Many organic farms at all scales already use cover crops as part of their fertility management program, contributing to the sustainability of the overall system. Subsistence and small-scale farmers in developing countries who do not already practice cover cropping can benefit greatly in production and climate-related sustainability from adopting locally relevant techniques. And finally, low-cost, locally available sources of fertility are vital to the viability and success of urban agriculture projects that rely on cost minimization and closed-loop systems since external resources are not as readily available or economical in cities.In areas of the Central Coast where winter rainfall typically exceeds 25 inches per year, and especially on sloped ground, cover cropping in annual vegetable cropping systems is highly advisable to protect non-cropped soil from both erosion and nutrient leaching. Based on numerous studies, the optimum time for planting winter cover crops on the Central Coast is mid October. In our mild winter climate we can plant cover crops as late as January, however the best results in terms of weed suppression, stand uniformity, and biomass production are from cover crops planted in mid to late October or early November. Depending on rainfall patterns it is often critical to get winter cover crops planted prior to the onset of heavy winter rainfall. Cover crop ground preparation and planting are best accomplished when soil is dry enough to work without the risk of compaction, which can result in poor drainage and clod formation. This is especially important on heavier soils. Because timing is critical, growers need accurate long-range weather forecasts to help determine when to prepare ground and plant fall cover crops. Timing these operations is directly related to soil type and rainfall amounts, so each farm will have a different set of criteria on which to base ground preparation and planting schedules: the heavier the soil and the greater the rainfall, the tighter the window for fall planted cover crops. There is often a very tight window between cover crop planting and harvest of fall crops which, coupled with the potential for significant rain events, can add considerably to the excitement.Selecting optimum cool season cover crop mixes is challenging since there are so many factors involved. The optimum mix provides early and uniform stand establishment, good weed competition, and minimal pest and disease pressure. It “catches” potentially leachable nutrients, does not lodge or fall over in high wind and heavy rainfall events, does not set viable seed prior to incorporation, fixes nitrogen, does not get too carbonaceous prior to incorporation, and is relatively easy to incorporate and quick to break down once incorporated. The ideal mix also improves overall soil health and helps form stable soil aggregates by providing adequate amounts of carbon as a food source for the soil microbial communities. A good standard mix that has proven successful at the Center for Agroecology & Sustainable Food Systems Farm on the UC Santa Cruz campus over the past 20 years is a 50/50 mix of bell beans and lana vetch with no more than 7% cayuse oats, planted at a rate of about 175 lbs per acre with a no-till drill.There are many options available for mid- and late summer cover crops in the Central Coast region. Water use and “land out of production” are the two biggest challenges with summer cover crops, but in a diverse system they can provide good weed suppression and nutrient cycling, and can significantly improve soil tilth and aggregation when planted in rotation with mixed vegetables.

It is also interesting to note just how little temporal heterogeneity is found among the nonactive and rumination axes

The resulting data was formatted as a three-dimensional tensor, with cows on the first axis, time on the second axis, and mutually exclusive behaviors on the third axis . In order to more intuitively represent this data as the proportion of time devoted to each behavior, the total minutes a given cow was recorded as engaging in a given behavior on a given day was normalized by the total minutes that a given cow was recorded on a given day of observation.An ensemble of data mimicries was created for this observed data tensor using the simTimeBudget utility previously developed for the LIT package. Observational error attributable to the precision of the sensor was again simulated by stochastically resampling each observed hourly time budget using the joint Dirichlet-Multinomial sampling strategy . The resulting mimicry of the raw sensor data was then conditionally aggregated by date of observation to create a four-dimensional tensor formatted identical to the observed data tensor but with simulation number on the fourth axis to comply with LIT package formatting standards. Ensemble variance estimates were calculated over the simulation axis for each combination of cow, day, and behavior indices. As with overall time budget, these ensemble variance estimates could subsequently be used to scale dissimilarity estimates for the observed data to account for heterogeneity in multinomial-formatted data both within and between behavioral axes. In previous analysis of overall time budgets, trimming tray for weed jackknife resampling could be performed within-animal to nonparametrically estimate the reliability of the underlying behavioral signal by leveraging information in the temporal subsamples used to create the aggregate record .

With daily time budget records, however, this strategy could not be employed, ashomogeneity of hourly time budgets could not be assumed due to fluctuations in behaviors imposed by the management schedule and the circadian rhythms of the animals themselves . Thus, while the behavioral axis could be rescaled using the precision penalized ensemble weighted distances, systematic differences between days of observation was here accommodated by the empirically-driven iterative re-weighting of the time axis at the heart of the data mechanics algorithm .In order to encode daily time budgets, the basic data mechanics algorithm was extended to a three-dimensional tensor. As before, column clusters were used to reweight dissimilarity matrices calculated for row observations, and row clusters were used to re-weight dissimilarity matrices calculated for column observations, allowing structural information to be shared between the two axes . To accommodate a multivariate response, dissimilarity values were aggregated over the third axis at each matrix index prior to aggregation over the remaining matrix index for which the dissimilarity matrix was calculated. The efficacy of three dissimilarity estimators were explored. The first was a simple unweighted Euclidean distance , which is functionally equivalent to the standard two-dimensional implementation of the data mechanics algorithm, except that in the tensor implementation all behaviors for a given time budget observation are forced into the same cluster. The second dissimilarity estimator considered was the Kullback-Leibler Divergence Distance, which is the sum of the asymmetric relative entropy estimates for any two probability distributions vectors that sum to one. Here KLD-distance was calculated at each index for the first two axes and prior to aggregation over the axis for which the dissimilarity matrix was computed. The third and final dissimilarity estimator explored wasthe ensemble-weighted Euclidean distance, wherein the squared distance between observed values were normalized by the sum of the ensemble variance estimates calculated over all simulated datasets at the corresponding tensor indices.

The dissimilarity matrices calculated over the subset of row or column indices were reweighted iteratively using the cluster results of the opposite matrix axes until either cluster sets became stable or a maximum of ten iterations were reached . For all three dissimilarity estimators, clustering results were here calculated on a grid of metaparameter values: from one to ten cow clusters and one to six day clusters. Full details on the implementation of the tensorMechanics algorithm are provided in Supplemental Materials. To visualize the final dendrograms produce from each tensor mechanics optimization, the pheatmap package was used to create heatmaps wherein cows were arranged along the row axis, days along the column axis, and cells colored to represent the proportion of time that a given cow dedicated to a given behavior on a given day . To help identify temporal patterns captured by these clustering results, the column axis was annotated in purple with the day on trial that each set of records were recorded. In order to visualize patterns recovered across all five behavioral axes simultaneously, the ggpubr package was used to create a composite image of the final heatmaps created using equivalent 0 to 1 scales to facilitate direct visual comparisons of behavioral investments .Visualizations for final clustering results for all three dissimilarity estimators for all candidate metaparameter values are provided in Supplemental Materials. Comparisons of dissimilarity estimators for encodings of daily time budgets largely mirrored the clustering dynamics found with encodings of overall time budgets. Tensor mechanics encodings using thestandard unweighted Euclidean norm, which employed no rescaling of behavioral axes, overemphasized patterns in high frequency behaviors such as eating and rumination, recovering very little systematic differences in any of the activity axes across days or animals. KLD distance, on the other hand, produced a more balanced encoding across the five behavioral axes, but again was prone to over-estimate dissimilarity at the extremes of the time budget distribution, which resulted in a number of animals with extremely low or high eating times being classified as outlier in clusters containing only one animal, which would effectively remove these animals from consideration in downstream analyses of bivariate associations using these clustering results.

The encoding created using the ensemble-weighted Euclidean distance is visualized in Figure 1. This dissimilarity estimator provided arguably the most balanced representation of heterogeneity in daily time budgets across all five behavioral axes without over-pruning at the extremes of the distribution. Perhaps the most striking feature of this visualization, however, is the remarkable consistency in daily time budget records across this observational window. As with overall time budget, eating time remains the primary driver of difference between animals in this encoding. At the extremes of these eating time budgets, there are few systematic differences in clusters across the temporal axis, nor even much variability in observations within clusters, suggesting that neither transient environmental fluctuations or more persistent changes in management or the biology of these animals had much influence on the time investments of these individuals with more extreme behavioral strategies. Among cows with more moderate eating times, there is certainly more variability within clusters, but systematic differences in cells across the temporal axis are still surprisingly subtle. In contrasting what systematic patterns are apparent against the time indices of these records, it appears that time spent eating was slightly elevated during roughly the first 30-60 days on trial relative to the remainder of the observation window, which would have encompassed the transition period for many of these animals and the earliest stages of lactation for nearly all cows enrolled on this trial. This result is counter intuitive, as we would expect appetites to be suppressed immediately following calving and gradually increase throughout this observational window; however, it has been shown previously that time spent eating is not always well correlated with feed intake, weed trim tray as cows can compensate for reduced time invested in mastication at initial ingestion with increases in time masticating during subsequent rumination . Therefore, it is also possible that this early surge in time spent at the feed bunk might represent increase in feed sorting behaviors, which might represent latent behavioral strategies cows employ to cope with this period of negative energy balance, or it could simply reflect palatability issues related to the nutritional supplement or some other element of the total mixed ration. It is also interesting to note that, amongst animals with moderate eating times, the drop off in eating times seems to correspond with a slight increase in time spent highly active observed across all cow clusters. While this temporal subperiod would likely correspond to the return to estrus for some animals in this herd , the pervasiveness of this shift would seem more easily attributable to some latent biological or managerial shift, such as an increase in appetite during peak milk that led to a greater number of trips to the feed bunk between milkings . As these cows were under stocked with respect to bunk space during this observation window, this may indicate that conditions in this pen support good lying time irrespective of fluctuations in the management environment .In comparing the results of the tensor mechanics clustering with the results for data aggregated to an overall time budget, which are visualized in Figure 2, we see that the two encodings are in close agreement. The contingency table reveals that cluster assignments are nearly identical between the two data sets for the coarser branches of the dendrogram nearer the trunk, which reflect differences among more extreme time budgets, and differ only slightly in the cut off values established between branches representing more subtle differences in time budgets. Given the temporal consistency in this daily time budget data, this result is not necessarily surprising, as there is little additional information or complexity to be recovered from disaggregating this information across days for this herd.

Accommodation of temporal heterogeneity in the tensormechanics encoding appears to have largely only served to more finely distinguish between animals with low eating times. Closer inspection reveals that these distinctions appear to have been made based on differences in tradeoffs between the nonactive and highly active behavioral axes across the early and later phases of this observation window. As a result, the tensor mechanics encoding produces a coarser encoding of animals with more moderate time budgets. This has caused some of the more moderate overall time budget clusters to be consolidated in the tensor mechanics encoding, with some ambiguity in the resulting cutoffs that appears again to be driven by greater weight being placed on how the eating-nonactivity tradeoff and the eating-highly active tradeoff shifts over the observation window.To determine if the temporal heterogeneity found in daily time budget records would modify bivariate associations found in previous analyses between overall time budgets and other farm data streams, bivariate tree tests were conducted using the tensor mechanics results for all three dissimilarity estimators. The first set of tests conducted explored if patterns in daily time budgets differed between animals fed control diets and those whose TMR rations were amended with the Organilac fat supplement. In prior analyses with overall time budgets, no significant bivariate associations were recovered using the bivariate tree test framework ±a result that was not necessarily surprising given that control and treatment animals were housed together throughout the duration of the trial except while head locked for the morning feeding and herd check. Significant bivariate associations were, however, recovered between treatment group and all three dissimilarity estimators used to create tensor mechanics encodings. Visual characterizations of these relationships, which can be found in Supplemental Materials, were created using the compare Encoding utility, wherein contingency table cells were colored by point wise mutual information estimates that were deemed statistically significant by simulation against the null using multinomial resampling . The relationship recovered using the ensemble-weighted Euclidean distance is visualized in Figure 3. Organilac cows were significantly over represented among daily time budgets that were consistent throughout the observation window and characterized by relatively high time spent eating, moderate rates of rumination, and low non-activity. Cows in the control group, on the other hand, were over-represented among daily time budgets that were characterized by time spent eating that varied between moderately-low to moderately-high , rumination rates that were consistently relatively low, elevated rates of nonactivity during the first half of the observation window, and elevated rates of high activity during the second half of the observation period. Collectively, these results might suggest that, amongst cows with more moderate time budgets, control animals may have struggled in the early stages of this study that encompassed much of the transition milk period, whereas fat supplemented cows were able to maintain a robust and well balanced time budget throughout this early lactation period . Alternatively, Organilac cows might simply have spent more time eating throughout the first phase of this research trial in order to sort through their TMR ration in order to avoid the fat supplement.

Animals enrolled in this feed trail were also fitted with a CowManager ear tag accelerometer

Observations collected on a single animal over extended observation windows at high sampling frequencies can, however, contain a range of complex temporal patterns ± cyclicity, non-stationarity, autocorrelation, etc . Further, when sensors are applied to large heterogenous groups of animals housed socially in spatially restricted environments, recorded behaviors may also contain complex interdependencies between animals at the dyadic, triadic, clique, and herd levels . Failing to accommodate all these complex structural and stochastic features in a conventional model-based approach tostatistical inference risks returning spurious insights into the underlying behavioral dynamics. Developing such a model with a single PLF data stream can be challenging. Provided multiple data streams, however, the logistical challenges presented by model-based analytical frameworks can rapidly compound, creating significant barriers to cross-sensor inferences and thereby impeding researchers from extracting more holistic behavioral inferences from increasingly data-rich farm environments. Unsupervised Machine Learning tools may provide a more flexible and forgiving approach to knowledge discovery in the context of large sensor datasets . Such algorithms excel at identifying and characterizing complex nonrandom behavioral patterns lying beneath the stochastic surface of a dataset, while often employing relatively few structural assumptions about the data . Hierarchical clustering-based techniques offer an intuitive and highly adaptable approach to visualizing high dimensional datasets that is particularly well-suited to exploratory data analysis . Indeed, trim tray with screen by reducing the complex behavioral signals present in a sensor dataset into a series of discrete clusters, such algorithms may be viewed as an empirical extension of classical ethological techniques.

Discrete data, however, can be challenging to work with in most frequentist and even many Bayesian frameworks. Estimators based on information entropy, on the other hand, are purpose-made to quantify uncertainty in discretely encoded data without knowledge of the underlying distribution, andthus naturally complement hierarchical clustering-based algorithms . In these analyses, data mechanics algorithms were able to recover complex nonstationarity in the order in which cows entered the milking parlor. Some of these changes in queuing patterns could attributed to the shift to spring pasture access, but other transient and persistent shifts in entry order recovered in these encodings that may have been driven by environmental factors not experimentally recorded . Entropy-based non-parametric permutation tests were also successful in recovering preliminary evidence of significant nonlinear associations between encodings of entry order patterns and activity patterns recorded using ear tag accelerometers. In this paper we will explore how novel ensemble simulation techniques that emulate and adjust for the complex sources of error in PLF data streams may be used to produced more balanced encodings of multi-dimensional behavioral data. We also introduce a new dendrogram pruning algorithm that is able to efficiently repurpose these same ensemble simulations to ensure that that the power of hierarchical clustering tools do not exceed there solution of the sensor. Finally, we demonstrate the utility of information decomposition techniques within our existing non-parametric mutual information testing framework to better facilitate visual characterization of complex behavioral patterns across sensor data sets that might be overlooked in more conventional model-based analyses.

To demonstrate the efficacy of our analytical approach, data was repurposed from a feed trial assessing the impact of an organic fat supplement on cow health and productivity through the first 150 days of lactation. All animal handling and experimental protocols were approved by the Colorado State University Institution of Animal Care and Use Committee . The study ran from January through July in 2017 on a USDA Certified Organic dairy in Northern Colorado, enrolling a total of 200 cows over a 1.5 month period into a mixedparity herd of animals with predominantly Holstein genetics. Cows were maintained in a closed herd in an open-sided free stall barn, stocked at roughly half capacity with respect to both feed bunk spaces and stalls. Cows had free access to an adjacent outdoor dry lot while in their home pen, and beginning in April were moved onto pasture at night to comply with Organic grazing standards. Cows were milked three times a day, with free access to TMR between milkings, and were head locked each morning to facilitate data collection and daily health checks. For more details on feed trial protocols see Manriquez et al. and Manriquez et al. .In addition to standard production and health assessments, behavioral data was also obtained from several PLF data streams . Milking order, or the sequence in which cows enter the parlor to be milked, is automatically recorded as metadata in all modern RFID-equipped milking systems. Study cows were here milked in a DelProTM rotary parlor . At each morning milking, raw milking logs were exported from the parlor software, and the data processed to extract the single-file order that cows entered the rotary . A total of 80 milk order records ±26 recorded while cows remained overnight in a freestall barn, and 54 following the transition to overnight access to spring pasture ±were used to create discrete encodings for parlor entry patterns via data mechanics clustering .

The dendrograms summarizing the distribution of cow entry order patterns and subsequent heatmap visualizations will here be subjected to further analysis without modifications to the previously reported encodings. This commercial sensor platform, while designed and optimized for disease and heat detection, also provides hourly time budget estimates for total time engaged in five mutually exclusive discrete behaviors – eating, rumination, non-activity, activity, and high activity . Time budget data was collected on all animals for a contiguous period of 65 days , with the observation window beginning shortly after trial enrollment was completed on February 17th, and ending on April 23rd when the grazing season commenced and cows were moved overnight beyond the range of the receiver antennae. After dropping cows that were removed prematurely from the observation herd due to acute clinical illness, as well as several cows with persistent receiver failure, complete sensor records were available for 179 animals. In order to more fully focus on the logistical challenges in encoding and characterizing the complex multivariate dynamics of this system, we have here chosen to compress this data over the time axis to consider only the overall time budgets of these cows, and will leave explorations of the longitudinal and cyclical complexity of this dataset for future work.Domain constraints are not, however, the only stochastic feature that need be accommodated when working with time budget data. There is also the measurement error attributable to the sensor itself. Returning to the previous example, suppose that we also know that our rumination records are only accurate to r1hr. Is it then still appropriate to give more weight to the one-hour difference between Daisy and Delilah than between Betty and Betsy? Since both observations are within the bounds of error, attempting to enhance the underlying biological signal may instead only succeed in amplifying measurement noise. A closed form estimator, however, may not be readily generalizable to the wide range of measurement error models encountered with PLF sensors. We therefore propose that a simulation-based approach may offer a more flexible means of accounting for measurement error in dissimilarity estimates . The LIT package provides a built-in simulation utility for time budget data that seeks to mimic the stochastic error structure of the original data while still preserving the underlying behavioral signal . Data is provided as a tensor, with cow indexed on the first axis, time indexed on the second, and the component behaviors on the final axis. The count data at each cow-by-time index is then used to redraw a simulated data point from one of three optional distributions . In the first, the user may sample directly from a multi-nomial distribution centered around the normalized observed count vector. This model assumes that measurement error should shrink as a cow dedicates larger proportions of an observation window to specific behaviors, and intrinsically prevents estimates from being generated outside the domain of support. Variance can be under-estimated at the extremes of the domain, however, cannabis trimming trays if the probability for a behavior is non-negligible but the observed count is zero due to under-sampling. This issue may be addressed in sampling option two, where samples are redrawn from a Multivariate Beta Distribution , also known as a Dirichlet distribution, again parameterized using the normalized observed count. While this sampling strategy slightly biases the simulation towards the center of the distribution, it prevents undersampling at the extremes of the domain. Finally, users may combine these sampling strategies in sampling option three, wherein the probability vector used to parameterize the multinomial is drawn first from the Dirichlet, in order to further increase the uncertainty in the simulated data.

After simulation has been completed by redrawing samples at the finest level of temporal granularity supported by the sensor, the data can then be conditionally or fully aggregated along temporal axis as required for downstream analysis as a time budget. This simulation routine was used to create an ensemble of B = 500 simulated overall time budget matrices that mimicked the stochasticity attributable to a reasonable approximation of the measurement error of the sensor. Stored as a tensor with replication on the last axis, the variance of the ensemble of simulations could then be easily calculated for each combination of cow index and behavioral axis. If the underlying simulation strategy is a reasonable representation of the noise in the sensor, then these variance terms will then serve as a sufficient approximation of the relative uncertainty in each data point. We propose that that this information can then be incorporated into the calculation of dissimilarity estimates by servingas penalty terms in the calculation of an ensemble-weighted distance estimator defined in Equation 3.The rescaling strategy employed in our proposed dissimilarity estimator is strongly inspired by traditional Analysis of Variance techniques, thereby providing several insights into its anticipated behavior. First, because the simulations were generated using the multinomial or one of its analogs, we can infer that these penalty terms will not be homogenous across the domain of support, but should shrink as observations approach the boundary. This will allow the ensemble weighted distance estimator to emulate the rescaling dynamic achieved with the KL distance, but here rescaling at the extremes of the domain will ultimately be bounded by our simulated measurement error so as not to exceed the precision of the sensor. Second, because we have here emulated measurement error in our simulation using sampling uncertainty, the central limit theorem will apply . Thus, we can anticipate that as the number of observations per animals increases, the impact of measurement error on our inferences will shrink, allowing progressively more subtle differences between animals to come into resolution. Taking this property to its limit, however, can it be said that with enough observation minutes the differences between cows can be inferred with near certainty? That intuition, of course, is at odds with our characterization of a dairy herd as a complex system, and highlights an additional stochastic element that must be accommodated ±the behavioral plasticity of the cows themselves in response to changes in the production environment . Given the extended observation window of this particular data set, it would be possible to recalculate time budget conditional on day of observation, and then use the variance in daily time budget along each behavioral axis as a penalty term. Such estimates would collectively reflect heterogeneity in variance attributed domain constraints, measurement error, and behavioral plasticity. Such an approach would not, however, be feasible for datasets collected over shorter time intervals with fewer replications or in applications with behavioral responses where there is no clear hierarchy in the temporal structure of the same. We therefore propose that our stochastic simulation model can be extended to also provide a generalizable means to approximate the uncertainty of the underlying behavioral signal. As before, the measurement error was simulated by redrawing samples at the finest temporal granularity provided by the sensor. Prior to compression along the temporal axis, however, a random subsample of observations days was selected across all cows, and only these values used to calculate overall time budget. If all cows demonstrated comparable levels of consistency in their daily time budgets, then reducing the effective sample size of our simulated data sets through a subsampling routine would increase the ensemble variance estimates. This in turn would make our approximation of measurement error hyper conservative, but this increase would be uniform across all cows.

Each data point is then projected into the resulting low dimension linear space

Any days where less than 75% of the herd was successfully recorded in the parlor were also dropped. This left a total of 80 days of milk order observations ± 26 recorded while cows remained overnight in their pen, and 54 after the transition to overnight pasture. Finally, cows that were not present in at least 50% of the remaining milkings were excluded from further analysis. Of the 177 cows with sufficient records, 114 had no recorded health events.With this metric, the more consistently a smaller set of cows are observed in a given segment of the queue, the smaller the entropy values becomes to reflect less stochasticity in the system. In standard statistical models, the nominal value of estimators such as log likelihood and AIC scale with the size of the data set, and must be interpreted relative the value of equivalent terms assessed against a null model. Analogously, the nominal value of the entropy estimates scales with the number of discrete categories used. The maximum theoretical value occurs when no underlying deterministic structures are present and all categories are equally likely to occur, which algebraically simplifies to the log of the number of discrete categories used . Here the maximum theoretical entropy value would be log2 = 6.83. To visually contrast differences in stochasticity across the queue, the observed entropy values were plotted against the median entry quantile of the corresponding queue segment using the ggplot2 package, vertical grow racks with maximum theoretical entropy added as a horizontal reference line . Nonrandom patterns in queue formation could also be explored by tracking the entry position of individual cows over time. As entry quantile has a numerical value, we can now also use variance to quantify and contrast stochasticity between animals.

As with all analytical approaches reviewed in this paper, there are both strengths and shortcomings to either approach . In this system there are two potential drawbacks to this conventional summary statistic. The first is that variance estimates are quite sensitive to outliers, making it difficult to empirically distinguish between cows that occupy a wider range of queue positions and animals who typically occupy a narrower range but might have gotten jostled far from their normal position on one or several occasions. The second drawback is that, because variance quantifies dispersion about a central value, it cannot distinguish between cows that demonstrate little consistency in entry position and multimodal queuing patterns. For example, if a cow always entered the parlor either first or last, we would intuitively determine that this pattern is nonrandom, but the corresponding variance estimate would be the largest in the herd. Having recovered evidence of nonrandom patterns, the next step was to begin characterizing the behavioral mechanisms driving this heterogeneity. The most fundamental question that need be answered to inform further analysis was the degree to which queueing patterns were driven by individual or collective behaviors. Because cows jockey for position with one another in the crowd pen, where they are pushed up to enter the parlor, we know intuitively that entry quantile records cannot be considered truly independent observations. If cows move through this melee as independent agents, such that their position within the queue is determined by individual attributes ±preferences, dominance, etc ±then a linear model may still provide a reasonable approximation of the underlying system. Early observational work on milking order, however, has suggested that cows may form consistent associations when entering the milking parlor, particularly when heifers are reared together .

If cows move into the parlor in cohesive units, such that queue position is more determined by clique-level than individual attributes, then network analyses may be a more appropriate. Principal Component Analysis is commonly employed to visualize relationships between observational units in high dimensional datasets. In this approach, redundancy between variables, here each milking record, is captured using either covariance or correlation assessed across all data points, here all animals. An eigenvector decomposition is then used to linearly compress the information contained in the data via rotation of the orthogonal axes. New axes are added iteratively such that each new dimension is pointed in the direction of greatest remaining variability until only noise remains . PCA was here performed only on animals with no recorded health in order to prevent any anomalous queuing behaviors recorded from acutely or chronically ill animals from obscuring the queuing patterns of the broader herd. The correlation matrix was constructed using all pairwise complete observations, and a scree plot was used to determine the dimensionality of the resulting space . The plotly package was then used to visualize the final embedding. While PCA provides a computationally expedient means of visualizing high dimensional data, the underlying assumption of linearity is not always appropriate . In some data sets complex geometric constraints, such as those commonly found with images or raw accelerometer data, and other latent deterministic features may project data points onto high dimensional geometric surfaces collectively called manifolds .

When these topologies are nonlinear , the spatial relationships between data points cannot always be reliably maintained when projected directly into a linear space, which can lead to incorrect inferences . Imagine, for example, you had a round globe of the world and wanted instead a flat map. Applying PCA to this task would be analogous to smooshing the globe flat on a table. Some of the original geographic relationships would be discernable, but some locations would appear erroneously close, and some landscapes would be entirely obscured. Modern manifold learning algorithms strive to more reliably project the complex geometric relationships between observational units into a standard Some geographic features will still be lost, particularly over sparsely sampled regions like the oceans, but the spatial relationships between landmarks would collectively prove more representative of the original topography. To further explore the underlying structure of this data absent assumptions of linearity, and thereby potentially accommodate any complex geometric constraints imposed on milk order records by latent social structures within the herd, a diffusion map algorithm was implemented using functions provided in base R . This was done here by first calculating the Euclidean distance between temporally aligned vectors of parlor entry quantiles for each pairwise combination of cows, scaled to adjust for missing records, and then inverting these values to create a similarity matrix. From this similarity matrix a weighted network was created by progressively adding links for the k = 10 nearest neighbors surrounding each data point. A spectral value decomposition was then performed on the corresponding graph Laplacian matrix . The resulting eigenvalues were used to select the appropriate number of dimensions, and the corresponding eigenvectors visualized using the 3D scatter tools from the plotly package . Finally, as a means of comparing geometric structures identified in the observed dataset with those of a completely randomized queuing process, the permutated dataset generated in the previous section was also embedded and visualized using plotly graphics .Having determined from the previous visualizations that a linear model might be a reasonable representation of the underlying deterministic structures of this system, the next step was to explore the temporal dynamics of this dataset. In a standard repeated measures model, multiple observations from the same animal are assumed to be identically and independently sampled, mobile shelving system implying that sampling order should not affect the observed value. If the observation period is sufficiently long to allow the underlying process to shift or evolve over time, however, stationarity cannot be assumed. Failure to statistically accommodate a temporal trend can not only lead to spurious inferences due to incorrect estimation of error variance, but also risks overlooking dynamic features of the behaviors under consideration . In practice temporal trends are often assessed by first fitting a stationary model and analyzing the resulting residuals. This may suffice when the temporal trend is uniform across animals, but risks overlooking more complex nonhomogeneous temporal affects. This could occur if only a subset of the larger group displays a non-stationary pattern, a risk that is likely heightened in large socially heterogeneous groups.

In this physically constrained system, where we know that every cow moving forwards in the queue must force other cows backwards, compensatory trends could also be easily overlooked in collective assessment of residuals. We first assessed temporal trend using two conventional EDA techniques. First, the ggplot2 package was used to generate scatter plots of entry quantile values against the corresponding observation date for each individual cow, with pasture access annotated with a verticalline. Plots were visually inspected for non-stationary, and are provided in supplementary materials. Next, to further explore the impact of the shift from pen to overnight pasture access on morning queueing patterns, median queue positions from the two subperiods were plotted against using the ggplot2 package , and Pearson correlation and Kendall Tau were computed using the stats package . While these preliminary visualizations were easy to both generate and interpret, both treat cows as independent and somewhat isolated units. With such a large number of animals to consider, the capacity for human pattern detection is quickly overwhelmed, making it difficult to contextualize trends within the broader herd. Further, this approach fails to leverage non-independence between animals entering the parlor, and thus risks overlooking subtler collective responses. Data mechanics visualizations were implemented to simultaneously explore systematic heterogeneity in milk entry quantiles both between animals and across the temporal axis. This was done by first using entry quantile values to compute two Euclidean distance matrices: one quantifying the similarity between pairwise combinations of cows, the second quantifying similarity between pairwise combinations of daily milking sequences. These distance matrices were then used to generate two independent hierarchical clustering trees using the Ward D2 method . By cutting both trees at a fixed number of clusters, observation days and cows were both partitioned into empirically defined categories, and a contingency table was then formed with cow clusters as the row variable and day clusters as the column variable. The original distance matrices were then updated, using the clustering structure between cows to create a weighted distance matrix between days and vice versa, thereby allowing mutual information to be shared between the temporal and social axes of the dataset . After several iterations of this algorithm, clusterings converged towards a contingency table with minimal entropy, wherein the entry quantile values within each cell were as homogenous as possible. When the entry quantile values were subsequently visualized using a heat map, this highly generalizable entropy minimization technique served to visually enhance heterogeneity within the data driven by nonrandom patterns along either axis. Further, by facilitating the transfer of information between axes, interaction effects between the social and temporal dimensions of this system were magnified, which here provided a means to explore nonhomogeneous temporal non-stationary between subgroups within the herd . The data mechanics pipeline was used to analyze the temporal dynamics present in both the complete milking order dataset and the subset of animals with no recorded health events. Instead this algorithm was applied on a grid from 1 to 10 clusters for either axis. The resulting 100 heat maps scanned visually to determine the clustering granularity required to bring into resolution any interactions between social and temporal mechanisms. While this process may be computationally cumbersome, it is empirically analogous to systematically varying the focus of a light microscope to bring into resolution microbes of unknown size ±a tedious but effective means of identifying all relevant structures within a sample . Finally, the RColorBrewer package was used to add color annotations to the column margin, to clarify temporal patterns, and to the row margins, which served to visualize potential relationships between queue position, a selection of individual cow attribute variables, and the onset of recorded health complications.Having thoroughly characterized the stochastic structures present in this dataset, the insights gleaned from the preceding visualizations were incorporated into a linear model to evaluate the relationship between queue position and several cow attributes. The four days identified as outliers by the data mechanics visualizations were first removed and the dataset converted to long format to be analyzed as a repeated measures model using the nlme package . Cow was fit as a random intercept via maximum likelihood method. Guided by the results of entropy and data mechanics visualizations, VarIdent was used to estimate separate error variance terms for each cow, and the necessity of this data-hungry heterogeneous variance model confirmed via likelihood ratio test against the null model with homogenous variance .

We multiplied temperature and soil scores to make a preliminary suitability map

As California’s temperatures get hotter and precipitation becomes increasingly variable with climate change , we expect a further systematic overestimation of suitable areas identified based on the past 30 years of weather data. For the suitability analysis we assigned temperature and soil texture to three categories that were each associated with a score: good , tolerable , and intolerable , while precipitation was divided into ranges that were suitable with no additional irrigation, suitable with additional irrigation, and unsuitable.For temperature, we considered the average maximum temperature in the three hottest months of the growing season , categorizing them separately with the scores described above . We then multiplied these three categorized scores together and took the cube root to get temperature suitability scores for the state, also excluding any areas whose monthly 30-year minimum temperature was above 59o F. We followed a similar procedure for soil texture, using SSURGO estimates of clay content averaged across soil horizons at a 90m resolution . Because farmers did not give numeric estimates of how much clay was needed in dry farm soils, we made sure our defined ‘tolerable’ range encompassed the full range of clay content observed in participating farms’ soils . To define the ‘good’ range , we excluded the farm with the lowest clay content, which was also the only farm where farmers stated that they could not grow tomatoes of a high enough quality to consistently market them as “dry farm.” This multiplication reflects the interaction between temperature and soil texture, grow rack in which good texture can compensate for higher temperatures by increasing soil water holding capacity, and lower temperatures can lessen the evapotranspirative demand that would be particularly problematic for plants growing in sandier soils with a lower soil water holding capacity.

We then separated the dataset into three areas based off of farmers’ understandings of where tomato dry farming could occur with no added irrigation and where it could occur with supplemental irrigation , and excluding areas that would not get enough winter rain to grow a suitable winter cover crop . The final map shows suitability scores in all areas that are categorized a ‘cropland’ in the 2019 National Land Cover Database . These areas are superimposed onto groundwater basins categorized as high priority in California’s Sustainable Groundwater Management Act . Crop totals on land that was deemed suitable for tomato dry farm management in these areas were calculated using the 2021 Cropland Data Layer .By focusing on the characteristics that limited water can give a tomato, these farmers highlight a recurring theme in understanding the functional definition of dry farming tomatoes. As the Central Coast faces increasingly limited water availability, the idea of dry farming has gained traction among policymakers purely by virtue of offering a means to continue farming while maintaining a restricted water budget. However, these farmers are quick to recognize that dry farming is only a management style that they can afford to choose for their operations insofar as it can excite customers and return a reasonable profit. In this way, the product that dry farming creates, which is valuable enough to consumers that they are willing to pay a significant premium for it, is the outcome that defines the management approaches farmers can use. Farmers know that they could alter the schedule for the minimal irrigation they do put on their dry farm tomatoes to increase yields . However, while defining the practice by some maximum threshold of water application, and then choosing to allocate irrigation water to maximize yields, may be appealing from a water savings perspective, farmers recognize that they must define the practice in terms of outcomes and not inputs.

Farmers must produce what consumers have come to expect from a dry farm tomato if they are going to make dry farming an economically viable choice for their operation.To better understand where tomatoes might conceivably be farmed in California given the environmental constraints identified above, we modeled dry farm suitability on California cropland as a function of precipitation, temperature, and percent clay in soil. The resulting map shows what lands could potentially support a dry farm crop, with and without supplemental irrigation, using constraints that are relaxed to encompass the least restrictive farmer-elicited constraints . The map therefore errs on the side of including land that is not an ideal candidate for dry farming, rather than leaving off land that may potentially be a good fit. With rising temperatures and less reliable rainfall, this map, which is based off of 30-year normals, likely also systematically overestimates what areas might fall into these thresholds when projecting into future climatic conditions. All areas in blue indicate land that meets a threshold where dry farming could be considered in a non-drought year without adding any irrigation. Areas in orange indicate that, while there is likely enough rain to sustain a winter cover crop, some amount of irrigation would often be needed to grow a successful dry farm crop. Areas in darker colors connote land that falls in conditions that are closer to ideal, whereas lighter colors indicate that more conditions are tolerable, rather than ideal, for dry farming.It is crucial to note that areas that show up as “suitable” on the map–including the most ideal locations–will likely require years of diversified management for soils to build the water holding capacity and fertility that allow for peak dry farm performance. These areas should therefore be considered candidates for long-term dry farm management, rather than ready-to-go dry farm fields. Because the constraints used to build the model were elicited specifically with regard to tomatoes, this of course is not a comprehensive map of everywhere that might be considered for dry farming non-tomato crops.

Particularly when it comes to grains and perennials , the range of possible locations is likely much broader. In the case of grains, winter varietals can be planted that take advantage of rain in winter months, while tree crops have far more extensive root systems that can reach water well beyond that which might be available to a tomato, in both cases relaxing the temperature and precipitation constraints that tomatoes need to survive without irrigation. Tomatoes are likely a better proxy for other vegetable crops , though each will have its unique requirements . As we imagine a shift towards dry farm agriculture in California, it is also important to consider how land that is suitable for dry farming is currently being used. Combining areas that are suitable for tomato dry farming with and without irrigation, we compiled a list of the top ten crops by area that are currently grown on these lands . Some of them are currently being dry farmed with some regularity in the state and could signal particularly easy targets for a shift to low-water practices. Others are dry farmed in other Mediterranean climates and suggest an important opportunity for management exploration in lands that might be particularly forgiving to experimentation. The remaining crops are some of the most water intensive in the state and would therefore lead to substantial water savings if the land could be repurposed. While unrealistic in the near future, calculating potential water savings from a complete conversion of suitable lands to dry farming allows for comparison with other water saving strategies. Even assuming that an acre-foot of irrigation is added to each acre of dry farm crops every year , vertical racks if all the land listed in Table 3 were converted to dry farming and irrigated to the statewide averages listed in the table , California would save 700 billion gallons of water per year, or nearly half the volume of Shasta Lake, the largest reservoir in the state. Given the overlap between suitable dry farm areas and high priority groundwater basins, these potential water savings are especially valuable as water districts scramble to balance their water budgets in light of SGMA. Perhaps the largest caveat to these potential water savings–and any analysis of dry farm suitability that relies solely on environmental constraints–is the economic reality in which conversions to dry farming currently occur. As discussed above, while a dramatic reduction in irrigation inputs might be feasible from a crop physiological perspective, whether farms can remain profitable through such a transition is an entirely different question. Given a dramatically increased supply of dry farm tomatoes, the profits that current dry farmers rely on could easily crumble. When considering other, less charismatic crops that could be good candidates for dry farming , customers’ likely hesitance to pay as steep a premium for high quality produce as they do for tomatoes also casts doubt on the viability of a large-scale dry farm transition given current profit structures for farmers.Our suitability map shows potential for vegetable dry farming to be practiced on California croplands that are currently irrigated, though its expansion is inherently limited.

Even if markets could be adapted to support an influx of dry farmed vegetables, our map indicates that climatic constraints will largely require dry farming to be practiced in coastal regions or other microclimates that can provide cool temperatures and sufficient rainfall. However, the Central Coast’s tomato dry farming offers principles–but not a blueprint–for low water agriculture in other regions. Based on themes from our interviews, these principles show a cycle of water savings that connect reduced inputs, management diversification, and market development . The cycle begins with lower irrigation , which can be accomplished in concert with soil health practices that build soil water holding capacity and increase long-term fertility. Reduced weed pressure and lower biomass production can then lead to reducing other inputs, such as labor and fertilizers, while also allowing for further water savings. The combination of reduced inputs and soil health practices then gives rise to a product that is unique in its water saving potential, and may also be of unusually high quality. By encouraging consumers to appreciate the products, or through novel policy support, farmers can develop markets that will provide a premium for these low-water products–or payment for the practice itself–which in turn creates an opportunity to expand the practice, further lowering inputs.As we ask how policies may impact dry farm production systems, we find a forking path in what types of expansion may result from different policies. An increase in production can be accomplished through both scaling size and scaling number . Both options can tap into the water saving cycle to decrease water usage; however, the search for just, agroecological transitions has pointed time and again to the need for scaling number . On the Central Coast, small, diversified farms have used this water saving cycle to both cut water use and develop a specialty product that allows growers to farm in areas with high land values by increasing their land access, profits, and resilience to local water shortages. Through these principles, small-scale operations have differentiated their management from both industrial farms and even other small farms in the region by creating a system based in localized knowledge, soil health practices, and thought-intensive management. However, it cannot be taken as a given that this water saving cycle will continue to uplift the small scale operations on which it started. Recent work highlights the potential for biophysical and sociopolitical conditions to combine to shrink–rather than grow–the use and viability of agroecological systems . In the case of dry farm tomatoes, socio-political attention is already beginning to target the biophysical need to decrease water consumption. If well-intentioned policy interventions designed to decrease irrigation water use build markets that value the fact of dry farming, rather than the high quality fruits it produces , growers will be able to scale the size of dry farm operations without needing to rely on the highly localized knowledge required to produce high quality fruits. As large grocers scale up dry farm produce sales without worrying about quality-based markets that may quickly saturate at industrial scales, the agroecological systems that originally produced dry farm tomatoes may be edged out of the market. On the other hand, if policies build guaranteed markets for small farms growing dry farm produce, dry farming may grow by scaling out to more small-scale operations. Policies focused on water savings may then favor industrial or small-scale farms, depending on how interventions shape the “Market Development” aspect of the cycle.

It is also possible that Sordariomycetes themselves are improving dry farm outcomes

The ability to analyze rotations at the field scale across the entire Midwestern US allows us to ask how farmers optimize their rotations in complex economic and biophysical landscapes that include pressures to both simplify and diversify. Several biophysical and policy variables show statistically clear relationships with rotational complexity: high land capability, high rainfall during the growing season, and proximity to bio-fuel plants are all associated with rotational simplification. Given policy incentives, farmers often find that “corn on corn on dark dirt usually pencil out to be the way to go,” with farmers growing corn year after year when high quality soil is available. However, when that proverbial “dark dirt” is not available, calculations are not so simple. If growing conditions are sufficiently poor , these intensive corn systems may not be profitable, and farmers will have to rely more heavily on non-corn crops to maintain crop health and profitability in their fields.We see this dynamic at play with land capability in the present analysis. Despite—or rather because of— the fact that more diverse rotations improve soils, the most degrading cropping systems counter intuitively tend to occur on the highest quality land. Highly capable lands can be farmed intensively without dipping into a production “danger zone” in years with weather that is historically typical for the region, creating a pattern of land use that is likely to degrade these high quality lands in the long term and potentially jeopardize future yields, particularly in the face of climate change. Recent analyses show that enhanced drought tolerance and resilience for crops is one of the key benefits of diverse crop rotations. In the present analysis, mobile rolling shelving mean rainfall during the growing season correlates positively with rotational simplification. Farmers may therefore be employing crop rotation in areas of low rainfall to achieve production levels that will keep a farm solvent, as was seen with rotational complexity increases in Nebraska during a drought period from 1999 to 200773.

This trend is further accentuated by the negative interaction between land capability and rainfall variance in our analysis, where higher rainfall variability leads to even more diverse rotations on marginal lands.Proximity to bio-fuel plants, the main policy indicator in our model, showed a statistically clear trend towards rotational simplification, likely due to increased economic profits. Local corn prices increase by $0.06 – $0.12/bushel in the vicinity of a bio-fuel plant, amplifying incentives to grow corn more frequently. Wang and Ortiz Bobea were surprised not to find an impact of bio-fuel plant proximity on county-level frequencies of corn cropping in their own analysis, and the present analysis—done at a field rather than county scale—shows exactly this expected effect: corn-based rotations are simplified when in closer proximity to a bio-fuel plant. In the current economic and policy landscape, farmers are pushed to simplify rotations through more frequent corn cropping, especially in proximity to bio-fuel plants, while marginal soils and low rainfall pull fields towards more diverse rotations.RCI’s ability to classify rotational complexity across large regions at the field scale and with low computational cost opens doors to future analyses that explore the interplay between localized landscape conditions, management choices, and agricultural, environmental, and economic outcomes. We see a strong potential to employ this metric not only in new regions, but in analyses that address how results from field experiments with crop rotation may scale up to regional levels. We also note that the metric should be used with caution. For example, because RCI cannot recognize functional groups in crop sequences , it cannot capture the added benefits that diverse functional groups often add to a rotation. In addition, though RCI includes a perennial correction that avoids penalizing multiple consecutive years of perennials the metric likely still underestimates the benefits of perennials in rotations. RCI is neutral to the soil benefits of annuals vs. perennials, while in practice the year-round cover and crop species mixes that often accompany perennials may boost soil benefits beyond those of annuals. Consecutive years of perennials are uncommon in our study area , and we encourage caution before applying the metric to regions with a more substantial perennial presence.

We therefore recommend using RCI in studies that explore a wide range of cropping sequences where large differences in RCI are very likely to be meaningful, rather than as a tool to rank sequences that give similar scores. It is also important to note that, though the index can be applied to data of any sequence length, RCI values from different sequence lengths cannot be compared to each other; a rotation that results in a 2.2 from examining a six-year sequence will not be a 2.2 when examining a five or seven-year sequence. We also note that in using crop sequence as a proxy for crop rotation, RCI cannot fully capture the cyclical nature of true crop rotations. Because RCI examines a fixed number of years, it may “split up” identical rotations in ways that give slightly different scores or ABBAAB in a six-year sequence). As these discrepancies will decrease when longer sequences are considered, we recommend applying RCI to sequences that are as long or longer than the longest expected rotation in the study region.We hope to see RCI used in future analyses that extend beyond the Midwest; however, regional and historic patterns of crop production likely influence farmers’ rotational decisions and may render RCI scores calculated from disparate geographical regions difficult to interpret when called into direct comparison. We therefore see great promise in RCI as a rotational metric, and caution against applications that are overly narrow and overly broad .The time period chosen in this study, 2012 – 2017, coincides with the introduction of the Renewable Fuel Standard, or “bio-fuel mandate,” which took full effect in 2012. This policy mandates that 7.5 billion gallons of bio-fuel be blended with gasoline annually, and caused bio-fuel plants to open and local corn prices to soar across the Midwestern US. Now in 2021, there is significant political pressure both to maintain the bio-fuel mandate in its current state and to relax the standards, and new exemptions to the mandate have already caused several bio-fuel plants to close in the region. Given the link between bio-fuel plant proximity and rotational complexity, our analysis suggests that these closures, if continued, would likely be associated with an increase in mean RCI in the Midwestern US.

Using our current model, simulations of randomly closing 20 of the 198 bio-fuel plants in the region lead to an increase of 0.003 in average RCI in the region, driven by greater distance to the nearest bio-fuel plant. In turn, increasing average RCI by 0.003 represents, for instance, the equivalent of 41,000 ha of cropland switching from corn-soy rotations to the most diverse rotation possible . Rotational simplification near bio-fuel plants is a pertinent example of the influence that policy can have on farm management decisions and its landscape repercussions. Bio-fuel mandates are one of several policies, including crop insurance and research funding priorities, that currently maintain the profitability of corn production; however, these policies need not be the ones that define rotational landscapes, and increased funding for policies such as the Conservation Stewardship Program could better align farmers’ economic incentives with improved environmental health80 . When strong economic incentives encourage rotational simplification, our analysis suggests that it is more likely to occur on land with favorable biophysical conditions for corn growth. With our current policy structure, the highest quality lands in the Midwestern US therefore become the most prone to degradation through intensive management.As rainfall becomes more variable with changing climates, farmers around the world are contending with droughts that are increasing in both intensity and duration. For many farmers, restricted water use has become a constant and looming threat, forcing the agricultural sector to confront a key question: how can we adapt to water scarcity without jeopardizing farmer livelihoods? This question is particularly salient for California’s agricultural system, heavy duty industrial pallet racks which has become increasingly fragile in recent decades due to its dependence on a shifting and shrinking water supply. Changing climates have caused droughts that not only result in massive financial losses , but also raise major concerns for farmers’ ability to maintain continuity in their farming operations. Because California’s waters are over-allocated even in years of typical rainfall, the Sustainable Groundwater Management Act,which requires sustainable groundwater use by 2040, implies that irrigation will need to be discontinued on hundreds of thousands of cropland acres. In this backdrop, dry farming, a practice in which farmers grow crops with little to no irrigation, has quickly garnered interest from farmers and policy-makers around the state. While dry farming is an ancient practice with rich histories in many regions, perhaps most notably the Hopi people in Northeast Arizona, dry farming emerged more recently in California, with growers first marketing dry farm tomatoes as such in the Central Coast region in the early 1980’s. In a lineage that likely traces back to Italian and Spanish growers, dry farming on the Central Coast relies on winter rains to store water in soils that plants can then access throughout California’s rain-free summers, allowing farmers to grow produce with little to no external water inputs. As water-awareness gains public attention, dry farming has been increasingly mentioned as an important piece in California’s water resiliency puzzle, however, while some extension articles exist, no peer-reviewed research has been published to date on vegetable dry farming in California. We therefore assembled a group of six dry farming operations on the Central Coast to collaboratively identify and answer key management questions in the dry farm community.

Growers identified three main management questions that would benefit from further research: 1. Which depths of nutrients are most influential in determining fruit yield and quality? 2. Are AMF inoculants effective in this low-water system, and more broadly. 3. How can farmers best support high-functioning soil fungal communities to improve harvest outcomes? Growers were primarily concerned with fruit yield and quality, with blossom end rot prevention and overall fruit quality being of particular interest due to the water stress and high market value inherent to this system. Managing for yields and quality present a unique challenge in dry farm systems, as the surface soils that farmers typically target for fertility management in irrigated systems dry down quickly to a point where roots will likely have difficulty accessing nutrients and water. Because plants are likely to invest heavily in deeper roots as compared to irrigated crops, we hypothesized that nutrients deeper in the soil profile would be more instrumental in determining fruit yields and quality. As deficit irrigation and drought change microbial community composition in other agricultural and natural systems, we hypothesized that dry farm management would cause shifts in fungal communities in response to dry farm management, which could in turn improve tomato harvest outcomes. Beyond general shifts in fungal communities, farmers were specifically interested in arbuscular mycorrhizal fungi inoculants, which are increasingly available from commercial sellers. Recent research has shown that AMF can help plants tolerate water stress, and we therefore hypothesized that commercial AMF inoculants might be beneficial. We organized a season-long field experiment from early spring to late fall of 2021 to answer these questions, sampling soils and collecting harvest data from plots on seven dry farm tomato fields on the Central Coast. Each farmer managed the fields exactly as they normally would, with AMF inoculation being the only experimental manipulation. We sampled soils for nutrients and water content at four depths down to one meter throughout the growing season to determine which nutrient depths influenced harvest outcomes. We also took DNA samples from soils and roots in surface and subsurface dry farm soils, as well as nearby irrigated and non-cultivated soils, sequencing the ITS2region to analyze the fungal community to verify inoculation establishment and more broadly characterize soil fungal communities to see how fungal communities changed under dry farm management and determine whether these changes or the introduction of an inoculant influenced harvest outcomes. We then used Bayesian generalized linear mixed models to estimate the effects of nutrient depths and fungal community metrics on yield and fruit quality data from 10-20 weekly harvests on each field. Our results highlight a tension between managing nutrients for fruit yield and quality, while fungal community metrics show promise for increasing fruit quality.

Other experiments offered restricted credit that can only be used to purchase agricultural inputs

Once the innovation is proven profitable and is locally available, its adoption may still be hampered by constraints facing SHF in accessing liquidity, risk-reducing instruments, information, and markets. These four categories of constraints have been extensively analyzed using in particular randomized control trials to identify their causal relations to adoption . These studies typically seek to identify ways of overcoming these constraints that could be implemented by governments, international organizations, NGOs, and benevolent agents such as philanthropic foundations and corporate social responsibility initiatives.Due to seasonality, especially under rainfed farming conditions which is where most of the lag in modernization currently prevails , there is a lack of correspondence between the timing of agricultural incomes and that of expenditures. As a consequence, the inter-temporal displacement of liquidity through credit and savings appears to be important for farmers to invest in new technologies, purchase inputs, optimize the timing of sales, buy consumption goods, and cover timely expenditures such as school fees. Financial services for SHFs appear to frequently be ill-designed for their purpose, expensive, excessively risky, and not easily available. Even when they have formal land titles, SHF are typically unwilling to put their land at risk as collateral with a commercial bank, thus acting as “risk constrained” . Microfinance products that effectively circumvent the collateral problem by relying on group lending and joint liability tend to be too expensive for the long agricultural cycles and have repayment conditions that are typically ill adapted to the timing of farmers’ capacity to pay . Availability of credit from formal sources, both commercial and non-profit, is consequently limited, weed dry rack and SHFs must either self-finance or rely on informal lenders with prohibitive interest rates. Hence, there would appear to exist a largely unresolved liquidity constraint on adoption originating on the supply side of the financial market. Yet, this is often not the main reason for low adoption which may be on the demand side.

Recent field experiments are providing evaluations of interventions aiming at relaxing the liquidity constraint on SHFs, with fertilizer the most commonly used indicator of technology adoption because of its ubiquitous recognition and yet massive underuse. While contexts and interventions vary for these experiments, they surprisingly tend to show that a liquidity constraint is not the reason why a majority of SHFs are under-investing in fertilizers. The main constraint may be instead lack of profitability in adopting fertilizers.A first category of experiments consists in providing unrestricted access to credit to a defined eligible population, as was done in Morocco , Mali , and Ethiopia . While interest rates in these studies were variously subsidized , uptake remained low: only 17% of eligible farmers took a loan in Morocco, 21% in Mali, and 36% in Ethiopia. Furthermore, farmers that did take a loan only used a small fraction of the liquidity to increase their expenditures on fertilizer or other agricultural inputs . Such credit displaces the equilibrium allocation of liquidity in favor of the targeted inputs, similarly to what a price discount would do. And yet, uptake remained low. In Malawi input credit for high-yielding maize and groundnuts was taken by 33% of the farmers . This low demand for credit thus seems to be reflective of a low demand for the inputs themselves. Low demand for fertilizer is exemplified in two rather extreme experiments. In Mali, Beaman et al. provided to another group of farmers a pure cash grant, rather than the credit described above. This only increased expenditures on fertilizer by 15%, in comparison with 11% with a credit that had to be paid for, showing that credit is not the major constraint to adoption. In Mozambique, a group of progressive and well positioned farmers, with good access to extension services and to input and output markets, were offered vouchers with a 75% discount on fertilizer price, and yet only 28% of the farmers redeemed their vouchers .

There are a few experiments that show the importance of well-tailored credit and savings products to support the adoption of profitable innovations by SHF. Two cases exhibit the importance of accounting for the seasonal distribution of farmer income. In Kenya, One Acre Fund offered harvest-time loan at 10% interest rate with repayment expected 9 months after harvest, collateralized with stored maize. The objective was to allow farmers to avoid selling their harvest at the time of the year where prices are lowest, postponing sale to the period of high prices. 63% of eligible farmers took the loan . A similar savings scheme through group-based grain storage was introduced in Kenya. Fifty eight percent of the farmers took-up the product, and were twice as likely to sell maize on the market . While these financial products contribute to the households’ welfare, and indirectly increase the return to agricultural production, there is no evidence that they induced higher adoption of fertilizer or other modern technologies. These experiments point to the existence of other constraints to fertilizer use. This can be lack of complementary inputs , excessively high risk, or high transaction costs in reaching markets that all make fertilizer use not profitable. These other constraints need to be jointly addressed with credit availability. An example is new financial products such as Risk Contingent Credit–where repayment is insured by an index insurance, where insurance serves as collateral for the loan, and where the insurance premium is paid with loan repayment at the end of the season–that have promise and are under experimentation . Another example is precision farming where soil testing allows to customize fertilizer recommendations to heterogenous local conditions and to design comprehensive technological packages . Conclusion is thus that, in spite of presumptions, credit availability is not the main constraint to fertilizer use for a large majority of SHF in SSA and SA. For them, low fertilizer use is mainly due to low profitability associated with physical, market, and institutional conditions.

Interventions to jointly secure the profitability of fertilizers and the availability of well-designed financial products for their own particular circumstances are necessary for large scale adoption to occur.Smallholder farmers are exposed to many risks that can put their livelihoods and assets in jeopardy and deter investment. Shocks include weather, plagues, prices, and health. As a consequence, SHF engage in shock-coping adjustments after an adverse event has occurred, and in risk-management strategies in anticipation of shocks difficult to cope with. Both responses are costly. Shock-coping includes dis-saving, emergency borrowing, sale of productive assets, emergency migration, use of child labor by taking children out of school, and postponement of consumption expenditures. Some of these responses can have long-term consequences, particularly when they imply decapitalization of assets, including child human capital and health . SHF also engage in risk-management practices. This includes preferring to use low return-low risk traditional technologies, holding large precautionary savings and productive assets that are biased toward liquidity , and engaging in income diversification at an efficiency cost . These costly strategies to deal with risk reduce the resources available for investment and technology adoption. Risk exposure also has a direct consequence on the adoption of technologies in reducing their expected return. Thinking of fertilizer, for example, dry racks for weed the possibly that drought or flood wipes out the harvest and hence the return to fertilizer discourages its application in the first place The obviously missing piece in the panoply of risk management and risk coping instruments used by smallholder farmers is insurance. Agricultural insurance is common in developed countries, although usually heavily subsidized. For SHF in developing countries, the administrative and implementation costs for agricultural insurance that requires verification of losses by an adjust or are too high to make it cost effective. Index insurance, where payments are triggered by a verifiable local index of rainfall or small area average yield has promise . Yet, take up has been very low, typically not exceeding 6 to 18% at market price for multiple reasons including basis risk, lack of trust, liquidity constraints, and limited salience of benefits The question for this paper, however, is whether insurance, when taken, induces farmers to adopt technology. Only a few studies in which the insurance uptake was sufficient permit this analysis. Experimental results for Ghana have shown that farmers that purchased insurance increased their use of chemical inputs by 24% . Mobarak and Rosenzweig offering a insurance service to farmers in India, show that it induced them to replace their use of traditional risk tolerant rice varieties by higher yielding varieties. Cai shows that a weather insurance policy in China induced tobacco farmers to increase their production of this risky but highly profitable crop by 16 percent and their borrowing by 29 percent. These experimental results suggest that if one could improve the design and marketing of the insurance product so that uptake could increase, technology adoption may follow. Promising avenues include products with reduced basis risk , financial training to help farmers better understand the value of insurance , group insurance on the presumption that managers have a better understanding of the product than farmers , and combining index insurance with other risk-reducing instruments, including social assistance for large shocks . Also promising is to reduce risk through resilient technology such as drought and flood tolerant seed varieties. SwarnaSub1, a superior rice technology with flood resilience properties, is appealing to farmers in flood-prone areas of India. Emerick et al. show that adoption of Swarna-Sub1 enhances agricultural productivity by crowding in modern inputs and cultivation practices , and increasing credit demand.

Finally, an innovative financial product introduced by BRAC in Bangladesh is contingent credit lines indexed on events such as flooding . This experiment shows that households given pre-approval to take a loan if they experienced flooding in their local area increased investment in risky production practices as part of their risk-management response. Since offering contingent credit lines has little ex-ante cost, the behavioral effect can be very large.A farmer’s decision to adopt a technology relies on his assessment of its value for himself. Hence beyond being aware of the existence of the technology, the farmer needs fairly complete information on the specificity of the technology, how to use it or adapt it, and how it would perform in his own context. To take some examples, the adoption of a new variety or management practice requires being aware and informed on their characteristics, associated best practices, and benefits. But even the adoption of fertilizers which are broadly familiar to most farmers, requires reliable information on their quality and the very specific quantity and timing of application that depends on local conditions. This learning process is particularly difficult due to the heterogeneity of contexts across farmers and across years . The traditional model of public or private sector extension agents faces multiple limitations. In addition to customization of the advice that should be delivered, the sheer number of smallholder farmers and their geographical dispersion limit what a public or private extension service can possibly achieve. In India, for example, fewer than 6% of the agricultural population reports having ever received information from extension services . In Uganda and Malawi, these numbers are higher but still imply receiving a service less than once a year . Given the multitude of smallholder farmers and the limited number of extension agents, extension services have typically focused their efforts on training “contact farmers” and expecting that these entry points in the farming community will circulate information in their social networks, inducing other farmers to adopt . The idea is to reinforce the process of learning and diffusion of technology from farmers experiencing for themselves to watching others experience the technology . Recent research has attempted to improve on this model by identifying optimum entry point farmers to maximize the subsequent diffusion of information and adoption. The theory of social networks gives useful clues as to which farmers to potentially select as contact farmers based on their position in the network . Other less theoretically based experiments compare the diffusion of the technology through contact farmers selected as “peer” farmers , large farmers, extension officer-designated farmers, community-designated farmers, members of women’s groups, etc. . Results do not converge to a general finding, suggesting that who is a better contact farmer is context specific. To be good entry points, contact farmers also need to be good experimenters, demonstrating the benefits of the new technology. Finding out who those good demonstrators are in not easy.

It normally takes a longer time to reach equilibrium in a temperate soil than in a tropical soil

Around 50% of the land area is used for agricultural purposes and is characterized by tropical, dry, and temperate climates along with diverse ecosystems, land uses, and management practices. The region is densely populated, and per capita land availability in some countries is less than 0.1 ha and is continuously decreasing. The possibility of increasing crop area is limited. The region is undergoing rapid industrialization contributing to greater emission of GHGs. In addition, there is rapid degradation of soil quality with low SOM content due to fertility-mining practices . Lal reported a C sequestration potential of 7–10 Tg C yr1 and 18–35 Tg C yr1 from restoration of degraded land in India and South Asia, respectively. With the adoption of recommended management practices on the cropland of South Asia, SOC potential was estimated to be 11–22 Tg C yr1 . The underlying assumptions included were the implementation of appropriate policies to promote recommended management practices such as conservation agriculture , mulch farming, cover crops, integrated nutrient management with manuring and biological nitrogen fixation, weed dryers and water conservation and harvesting. Lal also reported a soil inorganic C sequestration potential of 19–27 Tg C yr 1 of secondary carbonates and 26–38 Tg C yr 1 of leaching of carbonates in the arid regions of South Asia. Using International Soil Reference and Information Centre Soil Grids 250 m and FAO GLC Share Land Cover database, Zomer et al. reported a C sequestration potential of 0.11–0.23 Pg C hr 1 in South Asia. Assuming that C sequestration continues for 20 years, the current soil C stock of 7.68 Pg is likely to increase to 9.87 or 12.18 Pg for medium and high sequestration scenarios, respectively.

Grace et al. , using IPCC methodology together with local data, calculated a sequestration potential of 44.1 Mt C over 20 years from the implementation of zero tillage practices in rice-wheat systems of India.SOC sequestration is a dynamic process, and the amount and duration of C storage depends on the pools and their cycling , the form of stabilization , and the physical location of the C in the soil . Rates of turnover of organic matter depend on soil properties such as clay content and nutrient status. Clay is one of the key carbon-capture materials and tends to bind organic matter in soil and helps to protect it from microbial breakdown . Yang et al. also showed that the quasi-irreversible sorption of high molecular-weight sugars within clay aggregates, inaccessible by the microbes is responsible for clay-C protection. In addition, temperature plays a crucial role, which is complex because of variation in the temperature sensitivities of different SOM fractions . The impact of temperature becomes more crucial with a rise in ambient temperature due to climate forcing, resulting in microbially-driven increases in decomposition. Therefore, there are limits to C sequestration which are not only biophysical but also include technical and economic barriers.Over time, SOC reaches a steady-state equilibrium, balancing C gains and losses. Since organic inputs vary in quality, quantity, and subsequent interactions with soil constituents and environment, the ability of a soil to retain C is not unlimited. Carbon saturation is often used to describe the maximum capacity of a soil to retain C as a stabilized fraction based on soil properties . Sanderman et al. opined that while the term ‘soil C saturation’ is conceptually and theoretically appealing, the results from some of the long-term experiments may not support it. For example, Blair et al. found that total C stocks increased linearly with input levels of up to 200 Mg dry weight ha 1 for 15 years, without showing any signs of saturating behavior. However, Stewart et al. found some evidence of saturation.

Likewise, Johnston et al. reported that the annual addition of farmyard manure in the Broadbalk long-term experiment at Rothamsted increased C over the 160-year period, but the higher increase in early years was followed by a slower increase in later years, arriving at a new equilibrium. The time taken for soil to reach a new equilibrium tends to vary not only between soils within a temperate or tropical environment but also between the environments. It has been suggested that SOC saturation depends on clay and silt content and that there is a critical C concentration below which a soil’s function is reduced . Nevertheless, most current SOC models assume first-order kinetics for the decomposition of various conceptual pools of organic matter , which means that equilibrium C stocks are linearly proportional to C inputs .Carbon stored in soils is non-permanent. With changes in land use and land management, soil loses C, which can only be maintained or increased with the continuous addition of C input. By changing agricultural management or land use, soil C is lost more rapidly than it accumulates . Soil clay plays an important role in retaining C. Agricultural soil with a 50% clay content requires >2.2 Mg C ha 1 annually to maintain a given C level, while agricultural soil containing a 30% clay content requires more than 6.5 Mg C ha 1 annually. In addition, the rate of C input must be higher at existing soil C levels to maintain a level stock of C in the soil . Microbial decomposition and mineralization to CO2 is the major outcome of organic C. Approximately 1–2% of crop residues are stabilized as humified SOM for a period that are composed of large complex macromolecules, carbohydrates, proteinaceous materials, and lipids. This could be 60–85% of the total SOM . However, this notion, which was based on chemical analysis of the extracted materials has been challenged, and recent understanding suggests that humic substances are marginally important . Based on direct high-resolution in situ observations with non-destructive techniques, it has been established that humic substances are rather simple, smaller biomolecules . Although Hayes and Swift however, strongly disagreed with these views. They presented a detailed account of decomposition processes leading to the formation of a range of products including soil humic substances with a degree of resistance to microbial degradation. The new thinking in SOM research suggests that the molecular structure of plant inputs and organic matter has a secondary role in determining C residence times over decades to millennia, and that C stability depends mainly on the biotic and abiotic environment . The biotic and abiotic factors along with dynamics of labile C pools are required to evaluate management, land use, and climate change effects on SOC changes and soil functionalities .

New findings suggest that microbial decomposition actually facilitates long-term C sequestration by maintaining C flow through the soil profile , and that infrequent tillage may not cause sufficient disruption of soil aggregates leading to C loss . Schmidt et al. proposed that a new generation of experiments and soil C models will be needed to make advances in our understanding of SOM and our responses to global warming.While many land and crop management practices are known to enhance SOC sequestration, benefit accrual is constrained by the existence of numerous adverse forces on the ground. Table 5 provides key adoption constraints to an effective SOC sequestration strategy, the existing practice, and their implications. There are major barriers for farmers to adopt SOC sequestration practices because of the trade-offs involved. For example, the removal of crop residues from the field for other uses such as fodder, fuel, and fencing are traditional practices for managing residues. Not only is this an economical option for farmers, but there is also a lack of knowledge and capacity which discourages the adoption of practices promoting SOC sequestration. Likewise, shifting to zero- or reduced-tillage requires altering farm implements/equipment and the substitution of conventional crop and weed-control methods. The adoption of practices to enhance SOC also involves additional costs and the risk of getting lower yields in the short term. Much remains unknown about SOC storage, so it is difficult to estimate total benefits and to know which soil management practices offer the most potential for a given soil type, climate, and crop.Not only does SOC sequestration involve economic and biological costs but there can also be environmental cost. When mismanaged, some management practices that are known to result in C sequestration and GHG mitigation risk losing SOC and/or enhancing GHG emissions. Notably, N fertilization, either from organic or inorganic source, drying cannabis has negative consequences when applied sub-optimally–used either insufficiently or excessively. On one hand, when applied in inadequate amounts over time, for example in Africa, then there is no or negligible soil C build up . On the other hand, when applied in excess, for example in China and India, then soil C decreases from enhanced decomposition, which increases N2O emission, NH3 volatilization, and/or NO3 leaching. No-till compared to conventional tillage is another example of a practice that is reported to result in higher N2O emissions . No-till adoption may also increase the use of herbicides and pesticides, potentially affecting the environment . Sub-optimal or excess organic amendment to soil can also have an adverse effect on grain yield from nutrient immobilization. A growing interest in bio-fuel, resulting in a competition for fixed C, could also be a threat to SOC sequestration , as the use of bio-fuel involves burning of C which originated recently from photosynthetic activity.Lal et al. proposed six soil C management strategies to increase SOC: minimum disturbance of soil, maintenance of permanent ground cover, intensification of nutrient recycling mechanisms, creation of a positive nutrient balance, enhancement of biodiversity, and reduction in losses of water and nutrients. These strategies are generally applicable in South Asia and could be achieved notably through conversion of degraded land to perennial vegetation, increasing the NPP of agricultural ecosystems, and converting conventional tillage to no-till farming opined that a C-management strategy should not only be able to increase SOC content, but also should have some potential for reducing GHG emissions. Carbon management practices are aimed at increasing the ecosystem C balance by adding more C into the soil , increasing below- and above-ground biomass , sequestering SOC , and also reducing C losses from the soil . In the eastern Indo-Gangetic Plains of India, CA management practices like zero tillage with partial residue retention in rice-wheat systems could increase SOC content by 4.7 Mg C ha 1 after seven years of practice . Avoidance of adverse land use, management strategies, and restoration of degraded land can help in maintaining SOC stocks in soil . Table 6 provides details of various management options for increasing soil SOC stocks and reducing GHG emissions.Proper land leveling is known to enhance input use efficiency, crop growth, and yield . In South Asia, the majority of agricultural lands are poorly leveled by traditional land-leveling practices . Precision land leveling is laser-assisted, and very fine leveling of land is achieved with the desired grade within 2 cm of its average micro elevation . PLL is known to lower GHG emission by improving water and N use efficiency . Under Indian conditions, PLL could reduce almost0.15 Mg of CO2-e ha 1 year 1 of GHG emissions due to less time spent for pumping irrigation water and decreased cultivation time . PLL is critical for efficient water use and for increasing water productivity, and improves crop productivity through better crop establishment practices . There has been 6–11% and 10–25% increases in wheat yields in Punjab, India due to PLL. The associated increase in NPP in terms of crop residues and below ground biomass can be a source of soil C if further managed properly.Loss of SOC is often attributed to the practice of tilling the soil. Adoption of zero or reduced tillage will enable SOC sequestration, and is believed to be one of the key global mitigation strategies of climate change. Zero tillage has been widely reported as a viable option in increasing the C storage in soils , although few have reported no change . Most of the cases where zero tillage showed SOC increase, were mostly sampled to a depth of 30 cm or less, thereby not-revealing changes down the profile. In limited studies, where soil sampling was beyond 30 cm, no apparent difference in SOC between conservation and conventional tillage was recorded . Soil aggregates are stabilized under reduced and zero tillage practice, which physically protect C from mineralization , however, the effect is realized over the long run .