Mathematical Modelling of Hydrophilic Ionic Fertiliser Diffusion in Plant Cuticles: Lipophilic Surfactant Effects

The global agricultural industry requires improved efficacy of sprays being applied to weeds and crops to increase financial returns and reduce environmental impact. Enhancing foliar penetration is one way to improve efficacy. Within the plant leaf, the cuticle is the most significant barrier to agrochemical diffusion. It has been noted that a comprehensive set of mechanisms for ionic active ingredient (AI) penetration through plant leaves with surfactants is not well defined, and oils that enhance penetration have been given little attention. The importance of a mechanistic mathematical model has been noted previously in the literature. Two mechanistic mathematical models have been previously developed by the authors, focusing on plant cuticle penetration of calcium chloride through tomato fruit cuticles. The models included ion binding and evaporation with hygroscopic water absorption, along with the ability to vary the AI concentration and type, relative humidity, and plant species. Here, we further develop these models to include lipophilic adjuvant effects, as well as the adsorption and desorption, of compounds on the cuticle surface with a novel Adaptive Competitive Langmuir model. These modifications to a penetration model provide a novel addition to the literature. We validate our theoretical model results against appropriate experimental data, discuss key sensitivities, and relate theoretical predictions to physical mechanisms. The results indicate the addition of the desorption mechanism may be one way to predict increased penetration at late times, and the sensitivity of model parameters compares well to those present in the literature.


Introduction
The global agricultural industry requires improved efficacy of sprays applied to crops and weeds [1]. While spray application of agrochemicals is known to be effective, it is often inefficient [2]. Enhancing the efficacy of agrochemicals has many benefits [3][4][5]. For all foliar-applied agrochemicals, inefficiencies arise during deposition (sprays are not reaching the target) and retention (sprays are reaching the target but are not retained). In the case of systemic pesticides and foliar fertilisers, additional inefficiencies arise during penetration (not all active ingredient (AI) penetrates the plant) and translocation (not all AI within the target plant is transported to the site of biological activity).
The plant cuticle is considered the rate-limiting barrier to foliar penetration of agrochemicals [6]. Aqueous pores are dynamic nanopores within the cuticle that form only in the presence of water [4,7,8].
weakly affected by wax extraction and less molecular size-selective compared to the lipophilic pathway [11,41,42]. It has been noted that for lipophilic compounds, when the temperature increases, voids appear and disappear more frequently, which leads to diffusion rates greatly increasing [4,35]. Lipophilic compounds mobility is significantly influenced by lipophilic accelerators, which dissolve in the cutin and wax domains and act as plasticisers [11,35,43,44].
Surfactants have the ability to increase penetration of ionic AIs, especially at late times. Kraemer et al. [22] measures penetration of ionic CaCl 2 , which has been formulated with the lipophilic surfactant RSO 5, through isolated astomatous tomato fruit (Solanum lycopersicum L., cultivar 'Panovy') cuticles. If we focus on the penetration at late times, between 24 and 48 h, when CaCl 2 is formulated with RSO 5, penetration increases by 5.7% more than without RSO 5. This increase in penetration at late times is confirmed [45], with glyphosate penetration and ethoxylated RSO, along with other adjuvants with a range of HLBs, through barley (Hordeum vulgare L., cv. Plaisant) and ryegrass (Lolium multiflorum Lam.) leaves. They found a similar penetration profile, where penetration of glyphosate significantly increased to 72 h, and in some cases, a 40% increase was obtained between 24 and 72 h. Coret and Chamel [46] also confirmed this trend with glyphosate and non-ionic adjuvant penetration through isolated tomato fruit cuticles, where a mean of 20% additional penetration was obtained between 20 and 100 h.
The mathematical models for ionic penetration in plant leaves and cuticles have been reviewed elsewhere [12,13,47,48]. The authors Tredenick et al. [12] and Tredenick et al. [13] previously introduced two mechanistic models to simulate ionic penetration in plant cuticles and incorporated pore swelling, relative humidity, droplet evaporation, and hygroscopic water absorption that is influenced by the point of deliquescence.
Mechanistic models for adsorption and desorption, with a moving droplet radius due to evaporation, are limited in the well-established literature. We previously presented a model in Tredenick et al. [13], where ions bind to the leaf surface under the moving droplet radius, but do not desorb. Several adsorption-desorption models are available in the well-established literature, such as the competitive Langmuir model [49,50] and an adsorption-desorption model for chemicals in soil [51]. However, the models in [49][50][51] do not consider desorption with a moving droplet radius that is undergoing evaporation and Tredenick et al. [13] does not include desorption. We wish to consider this case for the adsorption and desorption of agrochemicals, which includes the moving droplet radius due to evaporation on the cuticle surface and create a novel mechanistic model.
There are several models on penetration of lipophilic compounds through plant tissues present in the literature. Previous modelling reviews are published elsewhere [47,48]. Some models incorporate diffusion [52][53][54][55], while others employ empirical expressions [56]. Whole plant leaf [52][53][54] and isolated cuticle models [55,56] have been studied. The empirical model [56] assumes equilibrium has been reached and relies on many experimental parameters, such as partitioning coefficients. Penetration in plant leaves is related to the physiochemical properties of the lipophilic compound, especially molecular size and lipophilicity, but penetration cannot be predicted by either property [25]. Therefore, a mechanistic model would be an advantage over an empirical model that relies solely on such parameters.
Keymeulen et al. [52] employed a model commonly applied in biological applications, based on Fick's first law of diffusion and Trapp and co-workers have modified this formulation to account for whole plant transport in several publications [53,57]. These models rely on many experimental parameters and cannot be applied to the experimental setup considered in Kraemer et al. [22]. Three experimental systems have been developed [54] to simplify plant penetration. These experimental systems are a novel way to test the governing mechanisms as they simplify the complex process of plant leaf penetration. The systems were simulated using finite-element techniques with Fick's first law of diffusion with convection of vapour, incorporating wind. However, they focused on leaf and stomatal transport and applied different experimental techniques to those in Kraemer et al. [22], so the models have limited applicability here.
We wished to consider the transport of lipophilic adjuvants by creating a novel mechanistic model that accounts for adsorption, desorption, and diffusion mechanisms and allows for any plant species, lipophilic adjuvant and concentration, temperature, and relative humidity to be used. We aimed to simulate the complex governing mechanisms involved in surfactant enhanced hydrophilic ionic AI penetration through isolated astomatous plant cuticles by utilising a predictive mathematical model. We will expand upon the models in previous works [12,13] by including adsorption and desorption of the ionic AI and lipophilic surfactant from the cuticle surface. In this work, we will investigate the mechanism, where the lipophilic adjuvant causes the ionic AI to desorb from the cuticle surface, and this may increase penetration at late times. Without lipophilic adjuvant, desorption of ionic AI is not possible. We are motivated to describe the increased penetration seen in the Kraemer et al. [22] experimental data at late times, and this was the key objective of this work. The modelling results were validated against experimental data [22] and a sensitivity analysis was conducted.

Model Framework
The model takes the form of a quasi-one-dimensional, diffusion model. We will briefly describe the modelling formulation, that is based on the authors previous works [12,13] and we refer the reader to Tredenick et al. [13] for a full description of the auxiliary equations governing evaporation in Appendix A. We account for three components, hydrophilic ionic active ingredient (AI), water (H 2 O), and lipophilic adjuvant (ADJ). The ionic AI can adsorb and desorb from the cuticle surface and the free ionic AI in the droplet solution can diffuse through the aqueous pores. The lipophilic adjuvant can adsorb and desorb from the cuticle surface, then the adsorbed adjuvant can diffuse through the cuticle via the lipophilic pathway. Water can diffuse through both the aqueous pores and the lipophilic pathway [4,34] within the cuticle. The authors previous model in Tredenick et al. [13] is simplified here by assuming the porosity does not significantly change in time, to focus and expand on the adsorption mechanism. Pore swelling may be influential under other environmental, initial and boundary conditions, as discussed Tredenick et al. [12] and Tredenick et al. [13], and this will be the subject of future works. All three components change primarily along the cuticle membrane thickness, x (0 ≤ x ≤ b). Initially, a droplet with known contact angle, volume, radius, concentration of ionic AI, and adjuvant is placed on the outer cuticle surface, at x = 0. A well stirred water bath exists at the inner cuticle surface, at x = b. The model, including the variables and parameters, as described in Table 1, governing partial differential equations, initial conditions (ICs), boundary conditions (BCs), and parameters, is as follows: BC -AI (bath) : BC -H 2 O (bath) : BC -ADJ (bath) : Parameters : CCR mode : CCA mode : dV % Ca penetration(t) = m(t) 100%   AI can travel through plant cuticle aqueous pores via Fickian diffusion and this is described by Equation (1). Water diffuses through aqueous pores and the lipophilic pathway (Equation (2)). Therefore, two diffusion coefficients are required, one for the diffusion of water through aqueous pores, D H 2 O,AqP (Equation (18)) and the lipophilic pathway, D H 2 O,L (Equation (17)). The adsorbed lipophilic adjuvant can diffuse through the cuticle via the lipophilic pathway (Equation (3)). We note the adsorbed adjuvant diffuses in terms of Γ ads ADJ in mol/m 2 , unlike ionic AI and water, which diffuse in terms of a concentration, c, in mol/m 3 .
Initially, the lipophilic adjuvant has not been adsorbed to the cuticle surface or through the cuticle (Equation (9)). An initially applied concentration of adjuvant is present on the cuticle surface (Equation (6)). We note the parameters described in Equation (15)- (19) are constants, as the porosities are constant. The diffusion coefficient for water in the lipophilic pathway, D H 2 O,L (Equation (17)), is the same formulation as the diffusion coefficient for water in the aqueous pathway, D H 2 O,AqP (Equation (18)), but is based on the porosity of the lipophilic pathway, ε L . Following from this, the diffusion coefficient for adjuvant in the lipophilic pathway, D ADJ (Equation (19)), is also a similar formulation to Equation (18). Droplet evaporation is governed by Equation (20)- (22) in constant contact radius mode (CCR) then constant contact angle mode (CCA), which includes hygroscopic water absorption, and further details are available in Tredenick et al. [13]. Water travels along the two pathways so the two diffusivities (D H 2 O,AqP and D H 2 O,L ) are included in the term for evaporation on the third line of Equation (22), which accounts for the possibility of water diffusing towards the droplet from the cuticle.
We introduce a mechanism for the droplet boundary condition on the cuticle surface, as depicted in Figure 1, where ionic AI and lipophilic adjuvant can be adsorbed and desorbed. Focusing on ionic AI, it first adsorbs, then the adsorbed ionic AI can be desorbed by the adjuvant that is adsorbing onto the cuticle surface. More adjuvant adsorbs as the concentration of adjuvant in the droplet increases due to evaporation. The same process occurs for adjuvant and ionic AI. The process is competitive, where ionic AI and adjuvant compete for the same sites on the cuticle surface. The adsorbed ionic AI can only be desorbed by adjuvant in the droplet solution, not ionic AI in solution. The most important feature of this model is that the adjuvant will desorb ionic AI that is adsorbed, and this desorbed ionic AI can then diffuse into the cuticle at later times. The adsorbed ionic AI would otherwise not diffuse through the cuticle.
Cuticle Figure 1. The Adaptive Competitive Langmuir model that includes adsorption and desorption. The initial condition is at t = 0 and t > 0 is a short time later after some droplet evaporation has occurred. Two generic molecules A (orange circles) and B (purple circles) sit on the cuticle surface.
Initially, A and B are in a droplet solution. Then as time progresses, A and B adsorb and desorb in a competitive process. As the droplet has moved due to evaporation, some of the molecules have been left behind. We note A in this case is ionic AI and B is lipophilic adjuvant. This figure is a simplification and describes the Adaptive Competitive Langmuir model and moving droplet radius, and we note other processes are involved in droplet evaporation, as described in text (not to scale).
The ionic AI and adjuvant are contained within a droplet that evaporates and the radius changes with time. This evaporation causes the droplet area to become smaller over time. Only the adsorbed ionic AI or adjuvant that is under the droplet can desorb into the droplet solution. Therefore, a portion of ionic AI and adjuvant that are adsorbed onto the cuticle surface are left behind by the moving droplet radius and are not available for diffusion. As discussed in Tredenick et al. [13], Schönherr [42] and Schönherr and Luber [75], hygroscopic water absorption in the context of the Kraemer et al. [76] experimental setup, where relative humidity is 70% and the point of deliquescence of CaCl 2 is 32%, will cause the total evaporation and penetration time to increase. If ionic AI is present (as free molecules) in the droplet, the droplet will not completely evaporate. We have assumed that adsorbed ions (that are not free) do not contribute to the hygroscopic water absorption, both inside (adsorbed under the droplet) and outside the droplet (left behind by the moving droplet radius).
The adjuvant adsorbs as a monolayer [33], as does the ionic AI. It is known that adjuvants, including the surfactant RSO 5, adsorb to various surfaces [29][30][31][32][33], lipophilic compounds adsorb and desorb on the cuticle surface [40], and ionic AIs adsorb to the cuticle surface [18]. However, we assume that ionic AIs can desorb and the process is competitive.
We modify the well-known competitive Langmuir formulation [49,50], to describe a model for ionic AI and adjuvant adsorption (ads) and desorption on the cuticle surface. The well-known competitive Langmuir formulation is described elsewhere for comparison ([77] p. 188). This novel model is known henceforth as the "Adaptive Competitive Langmuir model" and adapts to a moving droplet radius that is undergoing evaporation, which influences the amount of compound adsorbed and desorbed.
In Figure 1, we see the scenario of the Adaptive Competitive Langmuir model. A droplet sits on the cuticle surface that contains ionic AI (A-orange circles) and lipophilic adjuvant (B-purple circles), initially, at t = 0 and a short time later, at t > 0. Initially, ionic AI and adjuvant are in a droplet solution. Then after some time, the droplet has partially evaporated and some adsorption of ionic AI and adjuvant has occurred. Ionic AI and adjuvant can desorb each other, competing for the same sites, shown with the arrows. As the droplet has moved due to evaporation, some of the molecules have been left behind. We note that Figure 1 demonstrates the Adaptive Competitive Langmuir model and moving droplet radius, and is an oversimplification, as the ionic AI and adjuvant molecules may not be homogeneously distributed in the droplet solution [21,29,31,33].
Our novel Adaptive Competitive Langmuir model for ionic AI is as follows: Sites unoccupied by AI Sites occupied by AI where Ω AI is the fraction of sites occupied by ionic AI on the cuticle surface, k 1 is the rate constant for adsorption for ionic AI, c AI is the concentration of ionic AI in the drop, k 2 is the rate constant for desorption of AI and c ADJ is the concentration of adjuvant in the drop. The term (1 − Ω AI ) is the fraction of sites that are either vacant or unoccupied by ionic AI on the cuticle surface. Therefore the first term in Equation (26) describes the rate of adsorption, based on the sites available (vacant or no ionic AI) and the concentration of ionic AI in the drop, c AI . The second term in Equation (26) describes the rate of desorption, based on the fraction of sites occupied by ionic AI on the cuticle surface, Ω AI , and the concentration of adjuvant on the cuticle surface, c ADJ . The ionic AI desorption term of Equation (26) only allows adsorbed ionic AI to be desorbed by adjuvant in the droplet (c ADJ ), not ionic AI in the droplet. We then define the fraction of sites occupied by ionic AI under the droplet on the cuticle surface, Ω AI , which evolves in time with the changing droplet area during evaporation, as: where A drop,0 is the initial droplet contact area, Γ ads AI is the concentration of adsorbed ionic AI per unit area of cuticle surface under the droplet, Γ ads ADJ is the concentration of adsorbed adjuvant per unit area of cuticle surface under the droplet, Γ S is the saturated adsorbed molecules per droplet area and A drop is the droplet area on the cuticle surface.
The equation describing the concentration of adsorbed ionic AI per unit area of cuticle surface under the droplet, Γ ads AI , is found by substituting Equation (27) into Equation (26) and simplifying, producing Equation (28) and an analogous equation for Γ ads ADJ in Equation (29).
Change in droplet area (28) dΓ ads Change in droplet area where k 3 is the adsorption rate constant for the adjuvant and k 4 is the desorption rate constant for the adjuvant. The first term in Equation (28) accounts for the adsorption of the ionic AI to the cuticle surface based on the concentration of ionic AI in the droplet, c AI , and the sites that are unoccupied by ionic AI under the droplet, (Γ S A drop /A drop,0 − Γ ads AI ). The second term accounts for the desorption of ionic AI from the cuticle surface based on the concentration of adjuvant on the droplet surface, c ADJ , and the amount of ionic AI adsorbed per area, Γ ads AI . The third term accounts for the changing droplet area, A drop . Equation (29) is similar to Equation (28), except the last term is added to account for the adsorbed adjuvant, Γ ads ADJ , being transported away from the cuticle surface and into the lipophilic pathway of the cuticle, based on the porosity of the lipophilic pathway, ε L , the changing area of the droplet, A drop , the diffusion coefficient for the adjuvant, D ADJ , and the flux of the adsorbed adjuvant at the droplet boundary.
The equations for the concentration of ionic AI and adjuvant in the drop, c AI and c ADJ , are as follows: where k 5 is the kinetic rate constant for the ionic AI and k 6 is the kinetic rate constant for the adjuvant in mol. In Equation (30), the droplet becomes more concentrated due to evaporation, governed by the volume of water in the drop, V drop H 2 O , and the first term is formulated to incorporate the changing droplet area, A drop . The first term in Equation (30) is found by substituting Equation (27) into k 5 dΩ AI /dt and simplifying. The second term in Equation (30) describes how ionic AI is transported into the cuticle via diffusion, based on the porosity of the aqueous pores, ε AqP , the density of aqueous pores in cuticle, η pore , the control volume area, A Π , the changing area of the droplet, A drop , the diffusion coefficient for ionic AI, D AI , and the flux of the ionic AI at the droplet boundary. Equation (31) does not include the flux term for the concentration of adjuvant, c ADJ , transported into the cuticle, as the adsorbed adjuvant, Γ ads ADJ , is transported into the cuticle in Equation (29). The k 5 and k 6 rate constants in Equations (30)-(31) are as follows: where k 5 and k 6 are the kinetic rate constants in m 3 and c drop AI,0 and c drop ADJ,0 are the initial concentration of ionic AI and adjuvant in the drop in mol/m 3 . The parameters k 5 and k 6 were found to be a function of the initial concentrations, as changing the initial concentration alters the initial droplet contract angle θ 0 (see Tredenick et al. [13]).
In summary, the novel Adaptive Competitive Langmuir formulation consists of Equations (28) and (33). The two points of difference from the well-known competitive Langmuir model and the adaptive model described here, are that ionic AI desorbs adjuvant (and vice-versa) and the droplet contact area changes in time, causing some of the adsorbed ionic AI and adjuvant to be left behind by the moving droplet radius and are not available for diffusion. The Adaptive Competitive Langmuir model assumes the following: ionic AI and lipophilic adjuvant adsorb as a monolayer, that the process is competitive, all sites are equivalent, there are no interactions between the adsorbed ionic AI and adjuvant molecules on the cuticle surface and the adsorption-desorption process is not in equilibrium.
As this model assumes a competitive process, when one compound is adsorbing and desorbing, then the other must also, so k 5 or k 6 can only be zero or non-zero pairwise.

Numerical Solution Procedure
The model as described in Equation (1)- (33), is solved numerically. This is achieved using a finite volume method and discretising the model's partial differential equations with second order central differences to approximate the spatial derivatives. The resulting system of ordinary differential equations is solved using 'ode15i' [78] within MATLAB ® [79]. The fitted parameter F s in Table 1 was found by using the value in Tredenick et al. [13], then the parameter was adjusted by focusing on penetration for only the initial applied mass of 25 µg (5 g/L) of calcium using trial-and-error. The other fitted parameters in Table 1 were found using trial-and-error for only the initial applied mass of 25 µg of calcium. The experimental data [22] includes penetration of five different initial applied masses of Ca of 5, 25, 50, 75, and 150 µg. All the fitted parameters were then kept constant and used to solve the penetration of the other four initial applied masses.

Results and Discussion
The plant cuticle model, as described in Equations (1)- (33), is solved numerically, with parameters described in Table 1. In Figure 2, we can see the numerical solution results as percentage of ionic AI penetration with time (purple line) compared to experimental data (orange crosses- [22]) for CaCl 2 at 5 g/L (25 µg) with RSO 5 (1 g/L). The results of the model compare well to the experimental data, especially at late times, between 4 and 48 h. Penetration increases rapidly over the first 4 h, then begins to level out and the penetration can then continue to increase between 4 and 48 h by 15%, due to ionic AI adsorption with desorption. Penetration at 48 h has not yet reached equilibrium and can continue to increase at much later times, beyond 48 h.   Table 1.
In Figure 2, the fit of the model in the current work to the experimental data is reasonable; however, the fit is reduced at 24 h. We surmise this could be due to the Kraemer et al. [22] experimental setup, where each data point represents 10 repeats. Other authors [41,42,75] have noted they conducted 50-100 repeats as variability among individual cuticles can be large [80].
In Figures 3 and 4, we can see the validation results of the model, compared with experimental data for CaCl 2 with RSO 5, across 5 initial applied concentrations of CaCl 2 . Percentage Ca penetration is compared to time in hours in Figure 3 and as a mass in µg in Figure 4. The parameters, as described in Table 1, are the same across all 5 graphs, except the initial applied concentration of CaCl 2 ; 5 µg (A), 25 µg (B), 50 µg (C), 75 µg (D), 150 µg (E). Figure 3 is calculated according to Equation (25) and Figure 4 according to Equation (24). We have included these two figures as Figure 4 shows that increasing the initial concentration of ionic AI increases the mass penetration, while Figure 3 shows the penetration trend clearly.  The results in Figures 3 and 4 show a reasonable fit across all five plots, considering the complex mechanisms involved, producing an R-squared value of 97.95%, including the error bars. This indicates the model can predict changes in the initial applied concentration. In the figure, the desired trend of increased penetration at late times, after penetration has begun to level out, can be seen in all five graphs.
In Figure 2-4, we can see that the penetration timescale of the current work is significantly extended to more than 48 h. This extended timescale is caused by two mechanisms. The first mechanism is the desorption of ionic AI. At late times, ionic AI is continuously released back into the solution. The second mechanism is the hygroscopic water absorption. As the concentration of ionic AI increases in the drop, due to desorption, hygroscopic water absorption is increased. This in turn extends the evaporation timescale and the penetration timescale. Both the ionic AI desorption and the hygroscopic water absorption promote penetration at late times and both mechanisms must be included in a plant cuticle model for extended penetration to occur at late times.
The increase in penetration in Figure 2 in the current work, between 4 and 48 h, confirms the addition of the desorption mechanism may be one way to achieve increased penetration at late times. It has been well established that surfactants including RSO increase the penetration of ionic AI, especially at late times between 24 and 100 h, to a significant degree [22,45,46]. Earlier work by [13], which focused on penetration and droplet evaporation, found that penetration ceased after 9 h, as shown by the black line in Figure 2. Penetration ceased after 9 h as the ionic AI was either adsorbed on the cuticle surface to a significant degree or transported through the cuticle. Desorption was not included in the previous work. The current work produces a penetration profile that is closer to the experimental data.
In future works, the model presented here can be adapted to consider a range of constant temperatures. Equation (17)- (19) for the diffusivities for adjuvant and water in the lipophilic pathway could be modified to be a function of temperature by utilising an Arrhenius equation [81]. It has been noted that for lipophilic compounds, when the temperature increases, voids appear and disappear more frequently, which leads to diffusion rates greatly increasing [4,35]. Likewise, an Arrhenius equation could be used to formulate the k i terms if temperature varied. For lipophilic compounds, it has been found that desorption hysteresis occurs [40] at the inner water bath boundary (in real leaves the boundary condition would represent the apoplast of the epidermal cells), where the adsorbed lipophilic compound may desorb into the bath at a different rate and a portion of compound may be retained within the cuticle. This can be incorporated into the model by utilising a time varying function for Γ ads ADJ at the bath boundary condition in Equation (14).
The inner cuticle boundary condition for ionic AI is based on a well stirred water bath, so the concentration of ionic AI is instantaneously removed and is zero in Equation (11). In real plants, the inner cuticle surface would be adjacent to the apoplast of the epidermal cells, and transport through the apoplast is thought to be governed by diffusion [82,83]. A flux condition based on diffusion could replace Equation (11) for the ionic AI, and a similar equation for the lipophilic adjuvant in Equation (14).
In this work, we have added the mechanism, where the lipophilic adjuvant causes the ionic AI to desorb from the cuticle surface, and this increases penetration at late times. Without lipophilic adjuvant, desorption of ionic AI is not included, and the model can be applied when no adjuvant is included in the formulation by setting the initial applied concentration of adjuvant in Equation (10) to zero.
The fitted model parameters are shown in Table 1. Some of these parameters have physiological significance, such as the tortuosity or fractal scaling dimension, F s . A low value for F s of 1.15 was found with the fitting exercise. A low value for F s indicates the aqueous pores in tomato fruit cuticles are in the lower range of tortuosity, tending to straight. Penetration is quite rapid, even without adjuvant, so this value is logical. Fitting F s is reasonable as it facilitates the diffusion path length calculation, which currently cannot be found by a physical measurement [84]. The value for the partial molar volume of the adjuvant,v ADJ , was not available in the literature and was found by fitting with trial-and-error. Likewise, the value for the bulk diffusion coefficient of the adjuvant, D bulk ADJ , was not available in the literature and was found by fitting with trial-and-error. The value for the porosity of the lipophilic pathway, ε L , was found to be similar to the value for the porosity of the aqueous pathway, ε AqP . The values in Tredenick et al. [13] for the aqueous pore density, η pore , and logistic decay evaporation constant, χ, were utilised. We note the model has many fitted parameters. However, the penetration of ionic AI showed little sensitivity to D bulk ADJ and ε L , when fitting with trial-and-error. The parameters D bulk ADJ ,v ADJ and Γ S may be found with further experimental works, reducing the number of fitted parameters. The parameters D bulk ADJ andv ADJ for RSO 5 can then be used in any material, not just plant cuticles.

Sensitivity Analysis
A sensitivity analysis was performed with the results from the model, using parameters in Table 1 at an initial concentration of CaCl 2 of 5 g/L (or 25 µg) with RSO 5 (1 g/L). We have investigated the adsorption-desorption parameters, with the method described in Tredenick et al. [12] and Tredenick et al. [13]. The results of the sensitivity analysis are seen in Figure 5A Table 1 at an initial concentration of CaCl 2 of 5 g/L and RSO 5 of 1 g/L. The kinetic rate constants for the ionic AI are k 1 and k 2 (A) and k 5 (C). The kinetic rate constants for the lipophilic adjuvant are k 3 and k 4 (B) and k 6 (D).
In Figure 5A, we investigate the sensitivity of the models results for percentage penetration with the adsorption and desorption rate constants for ionic AI; k 1 and k 2 . In the figure, the ionic AI adsorption rate constant, k 1 , controls the percent penetration where the penetration starts to level out and is inversely proportional to penetration. A larger value for k 1 indicates more ions are being adsorbed so penetration is reduced. The ionic AI desorption rate constant, k 2 , is directly proportional to penetration and influences the amount of desorption of ionic AI. If k 2 is zero, penetration levels out earlier and no desorption occurs. As k 2 increases, more ionic AI is released at later times, which increases penetration and the penetration timescale is extended. Only when no adsorption occurs (k i = 0) can penetration reach 100%. The sensitivity in Figure 5A aligns to a similar adsorption-desorption model in the well-established literature [51].
In Figure 5B, we investigate the sensitivity of the model's results for ionic AI percentage penetration with the adjuvant adsorption and desorption rate constants; k 3 and k 4 . In the figure, the adjuvant adsorption rate constant, k 3 , has minimal influence over AI penetration and is inversely proportional to penetration. As k 3 increases, ionic AI penetration decreases as more adjuvant is adsorbed, which decreases the concentration of adjuvant in the drop and this causes less ionic AI to desorb and penetration ceases earlier. The adjuvant desorption rate constant, k 4 , has minimal influence over ionic AI penetration and is directly proportional to penetration. As k 4 increases, ionic AI penetration increases at late times as more adjuvant is being desorbed, which increases the concentration of adjuvant in the drop and increases ionic AI desorption. Figure 5C shows the sensitivity of the model to the kinetic rate constant for ionic AI, k 5 , which controls how influential the adsorption and desorption of ionic AI is, when included into the droplet boundary condition in Equation (30). Penetration of ionic AI is very sensitive to k 5 and penetration is inversely proportional to k 5 . As the value of k 5 increases, both adsorption and desorption increase. As adsorption increases, penetration levels out at a lower value. Increasing k 5 also increases the amount of desorption, which can be seen as the gradient for k 5 at 1.1 × 10 −9 is larger than at 2.2 × 10 −10 . A balance needs to be obtained between adsorption and desorption to explain the experimental data trends. Figure 5D shows the sensitivity of the model results to the kinetic rate constant for adjuvant, k 6 , which controls the level of adsorption and desorption of adjuvant. Penetration is inversely proportional to k 6 . When k 6 is small, this reduces both the adsorption and desorption of adjuvant, resulting in the concentration of adjuvant in the droplet increasing. This causes an increased concentration of adjuvant in the droplet solution that increases ionic AI desorption, which increases ionic AI penetration at late times. Conversely, when k 6 is large, both the adsorption and desorption of adjuvant are increased and the increase of adjuvant adsorption causes the concentration of adjuvant in the drop to significantly reduce and thus decreases ionic AI desorption. Therefore, a balance needs to be maintained.
The parameters discussed in the sensitivity analysis show appropriate sensitivities within the model and it is clear that a balance needs to be maintained between adsorption and desorption of ionic AI. We note that there are few experimental studies that measure the adsorption and desorption kinetic rate constants for adjuvants. Future experimental works, similar to Yamada et al. [18], could be conducted to directly measure the k i rate constants for RSO 5 and identify the mechanisms that are involved in adjuvant transport through plant leaves and cuticles.
The Adaptive Competitive Langmuir model applies to sets of two substances, which are adsorbing and desorbing in a competitive manner, under a dynamic droplet contact area. The plant cuticle diffusion model can apply to the penetration of most hydrophilic ionic AIs including lipophilic adjuvants in the formulation. It can apply to isolated astomatous plant leaf or fruit cuticle, where the aqueous pores are sufficiently large to allow ionic AI to be transported through the cuticle via Fickian diffusion. The model for lipophilic transport of adjuvant through the cuticle can be applied to both lipophilic adjuvants and lipophilic AIs. Adaptations can be made for hydrophilic adjuvants as transport is thought to be similar to hydrophilic ionic AIs. The model cannot be applied, in its current form, to whole leaf penetration, hydrophilic uncharged compounds, nor compounds, such as Fe chelates [85,86], that dehydrate aqueous pores. Adaptations could be made to this model to account for these compounds and this is the subject of future work. We note that if the type of hydrophilic ionic AI, lipophilic adjuvant, or plant species was changed, the model and mechanisms would not change. However, new parameter values would need to be found, and this may present a challenge depending on the extent of the available experimental literature.

Conclusions
We created a novel Adaptive Competitive Langmuir model that allows a dynamic droplet contact area, which influences the amount of adsorption and desorption of molecules. A mechanistic model has, therefore, been developed to simulate diffusion of hydrophilic ionic AIs, including lipophilic adjuvants through plant cuticles. This model makes novel additions to a simple diffusion model by incorporating the important governing mechanisms of the adsorption and desorption of ionic AI and adjuvant at the plant cuticle surface, along with the adsorption and subsequent diffusion of the adsorbed lipophilic adjuvant through the cuticle. The model has been solved numerically, producing results that show reasonable agreement with the experimental data. The addition of the desorption mechanism at the cuticle surface enables penetration of AI to increase at later times and further explains the experimental data trends. The sensitivity analysis indicated the parameters that controlled the adsorption and desorption had a significant impact on penetration. The results indicate the chemistry within the droplet on the cuticle surface is significant. Overall, the results of the validation and sensitivity analysis imply the model has included many important governing mechanisms simulating ionic AI penetration through plant cuticles. Λ (1 + cos(θ)) 2 f (θ) dθ.