1. Introduction
Huanglongbing (HLB), also known as "Yellow Dragon" or "Citrus Greening," is a bacterial disease that affects citrus trees. The causative agent is a gram-negative, persistent, and propagative bacteria (Aidoo et al., 2023; Darolt et al., 2021; Nguyen et al., 2023; Stover et al., 2018). In the Americas, the Candidatus Liberibacter americanus strain is prevalent (Bové et al., 2006; Clarke & Brown, 2018; De Chaves et al., 2023; Hansen et al., 2022). It is transmitted by the psyllid insect Diaphorina citri (Aubert, 1987; Garnier & Bové, 1996; Li et al., 2022; Liu et al., 2022; Pan et al., 2023). This psyllid primarily prefers grapefruit trees (Citrus paradisi) (Phuyal et al., 2020; Santiago et al., 2023) and the curry leaf tree (Murraya paniculata) (Hernández-Landa et al., 2018; Kalile et al., 2023; Meng et al., 2022; SENASICA, 2017) as hosts. Psyllid adults are not strong fliers, and their dispersion occurs through wind or human movement (Aidoo et al., 2023; Chen et al., 2022; Kalile et al., 2023). When the vector feeds on an infected plant, it acquires the pathogen and subsequently spreads the disease (Grafton-Cardwell et al., 2013; Oke et al., 2020; Setamou et al., 2022). HLB leads to the development of small, deformed fruits with bitter juice, fruit drop, and shrinkage. Symptoms may not become apparent until two years after infection, and the affected tree is no longer economically viable between five and eight years (Bové et al., 2006; Lin et al., 2017; Pandey et al., 2022; Spreen & Baldwin, 2013; Thakuria et al., 2023).
The disease was first reported in Brazil, specifically in the State of São Paulo in 2004 (De Carvalho et al., 2019; Zorzenon et al., 2021), followed by cases in Minas Gerais in 2005 (Alves et al., 2022; Sulzbach et al., 2018) and Paraná in 2006 (Cruz et al., 2021; Sauer et al., 2015). In the United States, it emerged in Florida in 2005 (Lin et al., 2023; Lopez & Durborow, 2014; Zapien-Macias et al., 2022) and California in 2008 (McRoberts et al., 2019; Nguyen et al., 2023). Cuba reported cases in 2007 (Bertaccini et al., 2019; Paredes-Tomás et al., 2023), and it entered Yucatán, Mexico, in 2009 (López-Collado et al., 2013; Pérez-Valencia et al., 2019). The Dominican Republic experienced its first cases in 2010 (Luis-Pantoja et al., 2021), followed by Belize in 2011 (Petri & Jorge, 2015) and Argentina in 2012 (Badaracco et al., 2017, 2022; Petri & Jorge, 2015). Paraguay documented occurrences in 2013 (COSAVE, 2017), while Colombia faced its first cases in 2015 (Ramírez-Godoy et al., 2018; Wang et al., 2021). In Ecuador, the presence of the vector has been noted since 2013 (Chavez et al., 2019).
The mapping analysis reveals that the United States, China, and Brazil stand out as prominent countries in these research areas. Scientists have established connections between HLB and citrus cultivation, insect vectors, and disease control. Additionally, research has delved into the relationship between HLB and parasitology, as well as microorganisms. Notably, the themes of "prevention and control" have gained increasing attention since 2018 (Supplementary Material 1).
China, Brazil, and the United States, which are the major global citrus growers, have all grappled with the impacts of HLB (Cui et al., 2022; Da Silva et al., 2019; Ding et al., 2020; Gao et al., 2022; Hallman et al., 2023). This disease has had far-reaching consequences on the economy, the society, and the environment. For instance, Durborow & Lopez (2012) explored the economic repercussions of HLB on orange production in California. Their research reveals staggering losses totaling US$ 2.2 billion over a 20-year span. In Florida, the economic toll of HLB was estimated at US$ 9.3 billion over a three-year period (Lopez & Durborow, 2014). In Brazil, Miranda et al. (2012) observed that in the absence of HLB, the projected production after 20 years would reach 6.7 billion boxes of 40.8 kg each. However, without effective HLB control measures, this figure saw a significant 45.3% decline. Meanwhile, in the state of Paraná, Da Costa et al. (2021) estimated HLB-related losses amounting to US$ 11.8 million dollars for the period 2011-2013, coupled with additional losses of US 39.2 million attributable to reduced citrus plantation productivity. Salcedo et al. (2011) conducted an assessment in Mexico that indicated that total production would decrease by 25% in a scenario with low HLB impact and by 41% in a scenario with a high HLB impact.
Citrus farming plays a significant role in Peru, with the industry experiencing substantial growth, now ranking as the 19th largest global grower. It contributed 3.3% to the gross production value of the agricultural gross domestic product (GDP). Notably, from 2018 to 2022, the Region Junín alone contributed a remarkable 56% to the country's total orange production (SIEA, 2023).
In 2018, the HLB vector Diaphorina citri was detected in Tumbes (SENASA, 2018) due to the informal trade of ornamental plants, particularly Murraya paniculata. Subsequently, three more detections occurred in 2019. In response, the Servicio Nacional de Sanidad Agraria (SENASA) and Agrocalidad of Ecuador initiated collaborative surveillance efforts along the northern border of Peru (SENASA, 2020). Given the presence of HLB in neighboring countries like Colombia and Brazil, coupled with the vector's presence in Ecuador and Bolivia (Giancola et al., 2019), the entry of the disease in Peru is highly probable. Such an incursion could result in significant economic losses unless growers prepare to combat it. Without a national phytosanitary program (PNF) dedicated to prevention, control, and monitoring in place to address the contingency, the importance of prevention cannot be overstated, especially since there is no known cure for HLB (Bové et al., 2006; Mora-Aguilera et al., 2019; Ortiz-Isabel & Manzo-Ramos, 2019). Furthermore, in Latin American countries, there has been a lack of ex ante studies examining the economic impact of HLB prior to its introduction. This study aims to simulate scenarios of economic losses, considering the possibility of HLB entering the Central Jungle of Peru (Junín).
The Junín Region is the leading grower of oranges in Peru. As per the SIEA (2023), it boasts 13,897 hectares dedicated to orange production, with impressive yields of 21.8 tons per hectare. This thriving industry involves the participation of 2,433 small and medium-sized growers. The region's growth is primarily driven by substantial investments in orchards to meet market demands.
The aim of this study was to assess potential economic losses within the direct orange production chain in Peru's Central Jungle, specifically the Region Junín. This assessment involved simulating the introduction and spread of HLB through a range of plausible future scenarios spanning from 2026 to 2045. The core objective was to estimate the cost-benefit of implementing a PNF as a proactive measure to mitigate damages caused by HLB. The study's ultimate purpose was to provide empirical evidence that could inform and support evidence-based public policy decisions.
2. Methodology
The study focused on Agricultural Units (AU) engaged in orange production within the Region Junín. To gather data, questionnaires were administered using a stratified sampling approach, targeting a total of 345 growers from a population of 3,366 AU (Equation 1). This population was distributed across two provinces: Chanchamayo (comprising the districts of Perené with 99 questionnaires, Pichanaqui with 41, Chanchamayo with 25, San Ramón with 16, and San Luis de Shuaro with 14) and Satipo (encompassing the districts of Satipo with 61 questionnaires, Pangoa with 49, Mazamari with 20, and Río Negro with 20) as illustrated in Figure 1. Details are given in. the questionnaire (Annex 1).
where n is sample size, N is the population, W h is the weight of the stratum in the population, Z is the normal distribution, and e is the error; while ph, and qh are respectively the proportion of success and failure in stratum h.
Three scenarios were developed to assess the economic impact of HLB introduction. The business as usual (BAU) scenario represents the trend production of oranges in the absence of HLB, essentially the status quo. The epidemiological scenario assumed that HLB present without the implementation of a PNF for disease prevention, control, and monitoring. In the HLB entry with variable compliance scenario, HLB is introduced and various levels of compliance and adoption by orange growers are considered. Compliance rates include 25%, 40%, 50%, 75%, and 100%, which represent the extent to which growers adhere to the measures outlined in the PNF for HLB treatment and management.
The first scenario (BAU) was built with the autoregressive integrated moving average (ARIMA) model that considers information on the value of orange production 1980-2022 (Xn), where α is the autoregression coefficient, εt is the error in time n and, θ is the coefficient of moving average part of the model (Equation 2).
In addition, growers were classified according to risk factors in the spread of HLB using the Kprototypes algorithm, and the cluster was segmented using the ELBOW method (Equation 3), considering the place of acquisition of seedlings, knowledge of symptoms, if there was pest control, the distance between the UA and the nearby populated center, associativity, training, movement of fruits in the UA and willingness to eliminate trees with HLB.
Figure 1. Region Junín and its main orange producing provinces and districts The Chanchamayo province in Region Junín and five of its districts producing oranges: (A) San Ramón, (B) Chanchamayo, (C) San Luis de Shuaro, (D) Perené, (E) Pichanaqui. The Satipo province including four districts producing oranges: (F) Río Negro, (G) Satipo, (H) Mazamari, (I) Pangoa. Source: INEI (2023).
The epidemiological model of Bassanezi & Bassanezi (2008) (Equation 4) allowed to estimate the severity curves, by considering st (proportion of tree crown area with HLB), age ranges of the trees (from 0 to 2 years, 3 to 5 years, 6 to 10 years, and more than 10 years), and the annual rate of severity progress by tree age (rL) with values of 3.68, 1.84, 0.92 and 0.69; For the proportion of tree crown area with HLB in year zero (s0), values of 0.2, 0.1, 0.05 and 0.025 were used per age range of the plantation; the values of rL and S0, the time in years (n) start at zero. Because there is no HLB disease in Peru, the parameters were field observations validated in orange production in São Paulo (Brazil), and the age of trees in plantations in Region Junín (Peru).
Subsequently, disease incidence progress curves (represented by yt, which is the proportion of trees displaying HLB symptoms in the orchard) were estimated using the Gompertz model, with annual incidence progress rate (rG) values set at 1,300, 0.6500, 0.325, and 0.244, as per the methodology outlined by Bassanezi & Bassanezi (2008) (Equation 5). The initial values of y0, indicating the proportion of trees with symptoms at year zero within different age ranges, were calculated using the inverse of the number of trees per hectare (1/number of trees). This method was adapted from Oliveira et al. (2013) and applied to the State of Bahia in Brazil. In this context, n represents the number of trees within each age category in the year 2024, with values of 38, 56, 98, and 63 corresponding to the categories spanning from 0 to 2 years, 3 to 5 years, 6 to 10 years, and more than 10 years, respectively.
With the values obtained from equations 4 and 5, the total severity (Sn) of HLB in the orchard over time for each age range was calculated (Equation 6).
Figure 2. Flow of the epidemiological model for calculating the production of oranges with HLB in scenarios with different degrees of adoption of the growers to the PNF. Source: After Bassanezi & Bassanezi (2008). Note. In the HLB-affected zone, the total severity (Sn) is a function of both the severity (sn) and incidence (yn), which then determines the performance under HLB conditions (RYn). Production calculations involving HLB vary per hectare in terms of the number of trees, harvested area, and overall production. Meanwhile, yields remain consistent across all scenarios, as they are contingent on the age of the plantations rather than the number of trees.
Table 1. Number of trees in a hectare by age category over time
The production yields in the presence of HLB at time 'n' (RYn) were derived concerning both healthy and diseased trees, accounting for the severity of HLB symptoms (Sn) within each age range. These estimates were obtained through fitting to the negative exponential model as outlined in Bassanezi & Bassanezi (2008) (Equation 7).
Given that the epidemiological scenario involving HLB in the studies spans a 20-year period, the analysis assumes the detection of HLB in 2026, without the presence of a PNF. Consequently, the projected impacts extend until the year 2045. Insights were sought from experts affiliated with the Asociación de Productores de Cítricos del Perú (PROCITRUS) and the Universidad Nacional Agraria La Molina (UNALM), both actively engaged in the Junín region. Additionally, data from the administered questionnaires in Region Junín were incorporated. Figure 2 provides an overview of the impact within various scenarios as modeled epidemiologically.
Epidemiological scenario (HLB without PNF)
The production scenario in 2045, under the influence of HLB without the implementation of a National Phytosanitary Program (PNF), was based on several assumptions. It assumed that tree renewal occurred every 25 years, involving the removal of seven trees per hectare while planting the same number. A constant density of 255 trees per hectare was kept, with trees aged 0-2 years old not bearing fruit. Additionally, it was considered that HLB comes to Peru in 2026. The harvested area considered variation in the number of trees within each age category (as detailed in Table 1). Hence, the number of trees was determined by the initial values at time zero, with 34, 58, 114, and 49 trees respectively allocated to categories ranging from 0 to 2 years, 3 to 5 years, 6 to 10 years, and more than 10 years.
To calculate the harvested area, the last three categories were considered, and being their sum kn at time n. Furthermore, an, bn, and cn denote the respective proportions of harvested area within each category (Table 2).
The calculation of harvested area per hectare considered an average yield of 22 t/ha (averaged over 2018-2022), with time (n) incrementing by 1 for each subsequent year. This calculation involved dividing the total production for year n (Pn) by 22 t, thereby resulting in the total harvested area over time (Spcn). Subsequently, this value was multiplied by an, bn, and cn to determine the area associated with each age category (for details, Annex 2, part a).
The calculations over time were resolved using Excel 2019 Solver. In this optimization process, the objective was to minimize α to achieve a value of zero. The variables considered include Pn (total production), PAn (production from 3 to 5 years), PBn (production from 6 to 10 years), and PCn (production over 10 years), as expressed in Equation 8 (for details, Annex 2, part b):
The total production under HLB conditions without the presence of a PNF was computed. PAn, PBn, and PCn, which represent the production by age range (without HLB), were multiplied by RYn, thus corresponding to yields under HLB for their respective categories. As a result, the total summation of the outcomes (δn + εn + θn) over time constituted the total production (Pn) under HLB conditions and in the absence of a PNF.
Table 2. Share of harvested area in one ha by age category over time
Full adoption scenario (with HLB and 100% of growers applying the PNF)
This scenario was simulated using the methodology developed by Bassanezi & Bassanezi (2008), wherein all orange growers adopt the National Phytosanitary Program (PNF). To estimate the trend production without HLB and under strict PNF implementation, several assumptions were made. Tree renewal was every 25 years, involving the removal of seven trees and the planting of seven new seedlings per hectare. A constant tree density of 255 trees per hectare was kept, with trees aged 0-2 years not bearing fruit. An annual removal rate of 2% of trees per hectare was used. Consequently, the harvested area varied, thereby leading to a change in the sequence presented in Table 1 (Annex 2, part c). However, the underlying variables remained consistent, and the number of trees within each category retained values of w0 for 0 to 2 years, x0 for 3 to 5 years, y0 for 6 to 10 years, and z0 for more than 10 years. To determine the share of harvested area by age in the epidemiological scenario, the analysis calculated the new area by age category (PAn, PBn, and PCn). These areas were then multiplied by the yields without HLB, which were computed using the Solver, to obtain the total production. The total new production was derived by summing these three categories.
Partial adoption scenarios (with HLB and a percentage of growers complying with the PNF protocol)
For these cases, the harvested area was divided between adopters and non-adopters. The scenarios of 25%, 40%, 50% and 75% of growers adopting the PNF caused the removal of 3% of trees annually due to HLB. For growers who did not carry out any control and where the disease spread, the epidemiological scenario methodology was used. Finally, the production of both types of growers was added.
To estimate the production in the partial scenarios, the following assumptions were made. Tree renewal occurred every 25 years, involving the removal of 7 trees and the planting of 7 new seedlings per hectare. A constant tree density of 255 trees per hectare was maintained, with trees aged 0-2 years not bearing fruit. Additionally, 3% of trees per year per hectare were removed from growers who adopted the National Phytosanitary Program (PNF). The harvested area varied over time, reflecting the total adoption scenario as presented in Table 1 (Annex 2, part d). Consequently, the number of trees within each category changed over time, with initial values set at w0 for 0 to 2 years, x0 for 3 to 5 years, y0 for 6 to 10 years, and z0 for more than 10 years. At time zero, these values were 34, 58, 98, and 77 trees, respectively. The percentage of growers who did not adopt the PNF was equivalent to the harvested area not covered by the program, following the methodology of the epidemiological scenario. Finally, the sum of production from both adopting and non-adopting growers equated to the total production within these scenarios.
Cost-benefit analysis (CBA)
The costs and benefits associated with each proposed scenario were quantified in monetary terms using the formula employed by Miranda et al. (2012), Dahlui et al. (2021), and Kotchen & Levinson (2023), which is represented by Equation 9.
where, for the epidemiological scenario, Rj is benefits without adopting the PNF in the year j, and Cj is the costs without implementing PNF in year j; while for the scenarios with PNF Rj is benefits after adopting PNF in year j, Cj: is the cost of implementing PNF in year j. Likewise, relying on the direct orange production chain in the Central Jungle and the production simulations spanning from 2026 to 2045, we estimated the income losses attributable to HLB in each of the proposed scenarios. We also assessed their impacts on wages and profits, in conjunction with the application of CBA. Consequently, the disparities in production between the baseline scenario and scenarios with HLB allowed for a comparison of the disease's impact, thereby enabling the estimation of potential future economic losses. This assessment considered the net present value (NPV) in 2026 and the social discount rate (SDR) of 8%. The SDR signifies the opportunity cost incurred by the Government of Peru when allocating resources to finance its projects, as outlined in government documents (MEF, 2019; Seminario de Marzi, 2017).
3. Results and discussion
Cluster analysis
Table 3 presents the classification of orange grow ers in the Junín region, emphasizing the risk factors associated with the dissemination of HLB. Cluster one comprises 190 growers, constituting 55% of the total, with the majority from the province of Chanchamayo (97%). Meanwhile, cluster two consists of 155 growers (45%), with 93% of them originating from the province of Satipo.
The identified risks upon the entry of HLB into the region include: (i) limited acquisition of seedlings from certified nurseries; (ii) less than 50% of farmers possessing the knowledge to identify HLB symptoms, (iii) infrequent pest control efforts, with most growers dedicating less than one day per year to this task.; (iv) inadequate levels of grower association, involving less than a third of the growers; and (v) proximity of Agricultural Units (UAs) to populated areas or main roads, which is a known factor in HLB transmission, as the disease is primarily transported by humans. This transmission occurs through infected fruits, cuttings, seedlings, or ornamental plants that contain the insect vector (Bassanezi et al., 2020; Zhang et al., 2020).
Table 3. Clusters of orange growers and risk attributes before the entry of Huanglongbing (HLB)
Source: After a questionnaire survey of 345 growers in Junín (2021).
These risks closely align with those noted by Sulzbach et al. (2018) in the State of Rio Grande do Sul, Brazil, where HLB has not yet been detected. In this region, orange cultivation is predominantly done by small growers with limited phytosanitary knowledge. Moreover, the proximity of the disease to neighboring Paraguay and Argentina poses an additional threat. These authors emphasize that the primary preventive measure against the disease is quarantine, aiming to prevent the introduction of citrus trees infected with HLB and Murraya paniculata carrying Diaphorina citri, which is the vector of this disease. Rigorous inspection is crucial, as the vector can also be transported via infected fruit (Halbert et al., 2010; Hall & McCollum, 2011).
Scenarios with HLB
Following the epidemiological model (Equations 4, 5, 6, and 7), Figure 3a reveals that the severity of HLB increases more rapidly in young trees. Consequently, by the year 2026, following the entry of HLB, the crown infection rate in 0-2-year-old trees may reach 20%, while for other age categories, it may remain at 10% or lower.
Figure 3. Percentage of severity, incidence, total severity and production yields with Huanglongbing (HLB) depending on the age of the trees in Region Junín. These curves represent the impact of HLB from 2026 to 2045.
Figure 3b shows the proportion of trees displaying HLB symptoms per hectare (incidence). In 2026, all age categories may exhibit infection rates of at least 1%. By 2027, the infection rate for 0-2-year-old trees could rise to 38%, while for 3-5-year-old trees, it may be 12%. In comparison, the last two categories may maintain infection rates below 3%. Figure 3c provides curves depicting the total severity per hectare, considering both the production of diseased trees and the ratio of healthy ones. It also shows the severity of HLB in relation to the age of the trees when the first symptom appears. The contagion and spread of the disease may occur most rapidly in 0-2-year-old trees, reaching 100% infection by 2033. Subsequently, 3-5-year-old trees may attain full infection by 2039, while the last categories may reach 96% and 83% infection rates in 2045, respectively. Figure 3d presents the returns with HLB. The reduction in returns appears to be most pronounced in the 0-2-year age category, quickly plateauing at 17% by 2031. For 3-5-year-old trees, returns may also stabilize at 17% by 2036, whereas 6-10-year-old trees could reach 18% in 2044. In the case of trees older than 10 years, returns may stand at 22% by 2045.
Figure 4 displays a range of prospective scenarios for orange production in Region Junín up to the year 2045. In the baseline scenario (BAU, without HLB), production is projected to increase at an average annual rate of 2%. Conversely, in the epidemiological scenario, the production of oranges affected by HLB may initially experiences minimal growth until 2025, after which it could undergo a significant annual decline of -6% until 2045, given the uncontrolled nature of the disease. The full adoption scenario, where production follows an annual growth rate of 2%, shows an increase that is somewhat lower compared to the trend production (BAU, without HLB). By 2045, this scenario anticipates a production loss of 4% in comparison to the BAU scenario (without HLB). For the partial adoption scenarios involving HLB, with varying percentages of growers adopting the PNF program (25%, 40%, 50%, and 75% of the total), the results demonstrate that a higher number of growers embracing the program leads to a reduced production loss over time. It's worth noting that in Region Junín, in termediate scenarios are considered more plausible because most orange growers (60% of respondents) are characterized as small and unorganized.
The introduction of the HLB disease would result in losses spanning the entire direct production chain. These losses encompass not only reduced income for growers but also a decrease in wages across various sectors, including production, transporttation, processing, and both internal and external trade. These losses are quantified in terms of decreased production, the reduction in the num ber of wages, and unrealized monetary gains. This quantitative data is integral to the prospective eval uation of the cost-benefit analysis (CBA), Detailed information is given in Annex 3 and Annex 4.
By the year 2045, in the presence of HLB and without the implementation of the PNF as per the epidemiological scenario, the citrus production would only reach 20% of what is expected in the BAU trend scenario, with a production of 89,839 tons compared to 457,022 tons. However, if 25% of producers were to adopt the PNF, production would increase to 38% of the expected output. Remarkably, with 100% adoption of the PNF, the losses would be limited to just 4%, as depicted in Figure 4. These findings align with research by Vera-Villagrán et al. (2016) on lime production in Colima, Mexico, in the presence of HLB, and Miranda et al. (2012) on citrus production in the State of São Paulo, Brazil. Both investigations underscore that without the implementation of a comprehensive phytosanitary program to control the disease, production experiences significant reductions. For Vera-Villagrán et al. (2016), the scenario with effective HLB control results in production totaling 227,446 thousand tons, whereas the scenario limited to vector control alone only achieves 123,578 thousand tons. Similarly, Miranda et al. (2012) show that a phytosanitary program significantly curbs production losses compared to scenarios without this critical element. These outcomes collectively emphasize the effectiveness of phytosanitary programs in curbing the spread of HLB disease and mitigating losses in citrus production.
Cost-benefit analysis (CBA)
The CBA for the scenario involving HLB and partial adoption of PNF by 25% producers, along with the HLB control costs within the direct productive chain of Junín orange, cumulatively assessed from 2026 to 2045, yields a NPV of S/ 1,224,531.6 thousand by the year 2026. This NPV is attributed to production losses, equivalent to US$ 368,345.1 thousand using the exchange rate of the year 2023. Additionally, there is an unrealized gain of S/ 224,713.2 thousand soles (Table 4). Considering the incremental costs incurred by various entities, the government's expenditure for implementing the PNF stands at S/ 25,781 thousand, while PROCITRUS contributes S/ 5,024.2 thousand. Furthermore, factoring in the additional expenses that producers will need to bear due to HLB, amounting to S/ 44,458.8 thousand (which include costs related to the removal of diseased plants, replanting, enhanced sanitary control, increased use of chemical inputs, among others), the simulation accounts for the impact of high HLB-related costs as referenced in Miranda et al. (2012). Consequently, the total incremental costs across these three entities sum up to S/ 75,264 thousand. These figures translate into a benefit-to-cost ratio (B/C) of 6.25 for growers, 10.78 for the government, and 9.02 for the combined efforts of the Government of Peru and PROCITRUS.
For every S/ 1 invested by growers to prevent and manage HLB in their crops, they would avoid the loss of S/ 5.25. In the case of the Government of Peru, for every S/ 1 it invests it would avoid losses of S/ 9.78 in the productive chain. direct production of oranges from Junín, and the three agents together (Government of Peru, PROCITRUS and growers) would avoid a loss of S/ 2.69 (Table 4), with the public intervention ratios being positive.
The results of the CBA and associated ratios are presented across a range of scenarios involving the adoption of the PNF by growers, including adoption rates of 25%, 40%, 50%, 75%, and 100%. It is evident that as the level of acceptance and adherence to the disease control guidelines by producers in creases, losses diminish, and the B/C ratios improve. This underscores the potential for greater success in preventing, managing, and controlling the dis ease, ultimately resulting in reduced losses across production, income, jobs, and the broader eco nomic activity. These outcomes emphasize that, from the government's perspective, the PNF re mains viable even when only 25% of producers adopt it, as the B/C ratios consistently surpass the 8% threshold (TSD) typically used to assess the fea sibility of financing public projects. Consequently, to enhance the effectiveness of public policy, it is im perative to concurrently promote the acceptance of the PNF among growers. Furthermore, to assess the feasibility of implementing the PNF, variations in producer prices and production costs are consid ered, with fluctuations of +/- 10%. Even under these dynamic conditions, the B/C ratio remains viable in all scenarios, consistently exceeding the 8% TSD (Annex 5).
Figure 4. Orange production projection (2023-2045) in Region Junín (t), trend scenario (without Huanglongbing [HLB]), epidemiological scenario (with HLB), total adoption (with HLB and national phytosanitary program [PNF]), and scenarios with HLB and partial adoption of PNF growers (25%, 40%, 50% and 75%). Note. The orange production projection from 2023 to 2045 (in t) in Junín considers the trend scenario (BAU, without HLB), the epidemiological scenario (with HLB), total adoption scenario (with HLB and PNF), and various scenarios with existence of the HLB and partial adoption of the PNF by growers (at 25%, 40%, 50% and 75% adoption). To estimate the effects of HLB on production, it is considered that the Diaphorina citri vector enters in 2024 or 2025 and in 2026 the existence of HLB is formally declared in the country (although the disease may not be declared in 2026, the analysis It is valid from the year HLB is detected and a PNF is implemented).
Table 4. Economic evaluation of Region Junín orange production under various cost-benefit scenarios considering wage losses and foregone profits due to HLB-induced production decline, simulated costs (in thousands of S/ and thousands of US$) Accumulated from 2026 to 2045, and net present value (NPV) as of 2026
(*) The conversion from soles (S/) to US dollars (US$) is with reference to the exchange rate of S/ 3.68 per US$ 1 dated 8th Nov. 2023. Source: BCRP (2023).
The findings from Oliveira et al. (2013) in Brazil indicate that the implementation of disease control measures would yield a significant net benefit, totaling 1,002 billion reais. Over a 20-year horizon, this translates into a net present value gain of 685 million reais. Similarly, Miranda et al. (2012) provide B/C indicators for 20-year timeframe. In their study, they report a B/C ratio of 4.6 for the Government of Brazil, Fundecitrus, and growers, which closely aligns with the results observed at a 40% adoption rate of PNF in Peru (B/C of 4.86). However, it is noteworthy that in Brazil, the Government and Fundecitrus achieve an exceptionally high B/C ratio of 57.3, which significantly surpasses all scenarios examined in the case of Peru, with the highest reaching 32.23.
The findings of Vera-Villagrán et al. (2016) regarding HLB management in diverse citrus settings in Colima, Mexico, present various production scenarios. These encompass a scenario devoid of HLB (BAU), another involving an epidemiological model, and several intermediate scenarios, thereby reflecting varying degrees of HLB control and tree eradication. Notably, one of these intermediate scenarios -when government intervention is factored in yields a benefit-cost (B/C) ratio of 10.9 monetary units, akin to the B/C return observed in Junín with a 25% adoption rate among producers (10.78). Moreover, the study by IICA et al. (2019), which investigated the impact of phytosanitary strategies against HLB in Brazil, revealed a B/C ratio of 2.28 for society over a decade. This figure is lower than the one obtained for Region Junín in Peru (B/C of 3.69 with a 25% adoption rate). Furthermore, the production scenarios and their associated B/C indicators consider the possibility that, in response to the incursion of HLB into the Central Jungle, investors may opt to halt the expansion of production and instead aim to maintain production levels by 2030. The representation of this simulation and its corresponding ratios are presented in Supplementary Material 2.
4. Conclusions
The presence of HLB disease poses a substantial threat to the citrus agribusiness, potentially rendering it economically unviable. In Peru, while swift vector control efforts have successfully contained detections on the North Coast, the risk of HLB introduction still looms over orange growers in the Central Jungle region. These growers remain susceptible to significant losses, given their lack of preparedness to combat the disease. As a result, prospective scenarios for the year 2045 consider the growth of orange production under different circumstances: the baseline scenario without HLB (BAU), scenarios with HLB alongside full compliance with the PNF by all growers, and the epidemiological scenario without the PNF. In the latter scenario, HLB is projected to lead to an annual production decline of -6%.
The benefit-cost (B/C) ratio for scenarios involving HLB and varying levels of PNF adoption by growers, accumulated up to 2045, reflects the benefits arising from the prevention of production losses. However, it is important to note that the incremental costs associated with these scenarios are expected to be borne jointly by the Government of Peru, PROCITRUS, and the growers themselves. From the government's perspective, the B/C ratio ranges from 10.8 with 25% adoption among growers to 37.3 when 100% of growers adopt the PNF, thus making the PNF economically viable in all cases. For growers, a positive B/C ratio under scenarios with HLB and the PNF is beneficial because it not only prevents production losses but also safeguards income and jobs within the direct productive chain. Moreover, when considering potential interactions with other sectors of the economy, the benefits become even more significant.
The implementation of effective public policies is strongly recommended to mitigate economic losses for the country's citrus growers. Failure to do so could result in losses throughout the entire production chain. Consequently, public policy should encompass not only the management of the disease itself but also the preparation of growers to effectively confront it.