Optimization of N, P, K for soilless production of Cannabis sativa in the flowering stage using response surface analysis

Following legalization, cannabis has quickly become an important horticultural crop in Canada and increasingly so in other parts of the world. However, due to previous legal restrictions on cannabis research there are limited scientific data on the relationship between nitrogen (N), phosphorus (P), and potassium (K) supply (collectively: NPK) and the crop yield and quality. This study examined the response of a high delta-9-tetrahydrocannabinol (THC) Cannabis sativa cultivar grown in deep-water culture with different nutrient solution treatments varying in their concentrations (mg L -1 ) of N (70, 120, 180, 250, 290), P (20, 40, 60, 80, 100) and K (60, 120, 200, 280, 340) according to a central composite design. Results demonstrated that inflorescence yield responded quadratically to N and P, with the optimal concentrations predicted to be 194 and 59 mg L -1 , respectively. Inflorescence yield did not respond to K in the tested range. These results can provide guidance to cultivators when formulating nutrient solutions for soilless cannabis production and demonstrates the utility of surface response design for efficient multi-nutrient optimization.


Introduction
Drug-type cannabis (Cannabis sativa) is an important horticultural crop grown for medicinal and recreational purposes. Historically, many countries have prohibited the cultivation of drug-type cannabis which consequently provided a significant barrier into research on this crop. However, change in social attitudes toward consumption of cannabis has led to the repeal of cannabis prohibition in several countries/regions around the world. Following the 2018 repeal of cannabis prohibition in Canada, production of cannabis has quickly become an important part of the Canadian horticulture industry worth billions of dollars annually (Zheng, 2021). However, cannabis cultivators still lack scientific information about optimal growing conditions, such as supply of mineral nutrients, to help maximize crop yields, quality, and profits.
Proper supply of mineral nutrients is essential for the efficient and sustainable cultivation of any crop. Among the most important mineral nutrients for plants are nitrogen (N), phosphorus (P), and potassium (K). However, few studies have investigated the response of cannabis to these nutrients. As a result, cannabis cultivators currently rely on fertilizer nutrient recipes developed by fertilizer companies, or by community consensus based on previously clandestine production. This poses a problem because deficient or excessive supply of nutrients may reduce yield (Caplan et al., 2017b(Caplan et al., , 2017a or lead to environmental pollution from runoff of excess nutrients (Beerling et al., 2014;Zheng, 2018). Nutrient runoff is an issue in many agricultural areas of the world because excess nutrients, specifically P, can lead to the eutrophication of water bodies (Schindler et al., 2016). In Ontario (the Canadian province in which this study was conducted) disposal of waste greenhouse nutrient solution, including from cannabis production facilities, is regulated by law at considerable cost to the cultivators (Ontario Ministry of Agriculture Food and Rural Affairs, 2019). An understanding of cannabis' mineral nutrient requirements can help us better synchronize nutrient supply and demand to maximize production while reducing nutrient waste and resulting environmental impacts.
Recent peer-reviewed studies have started to examine the response of Cannabis to mineral nutrients, but this area of research remains largely unexplored. These studies indicate the optimal N supply for both vegetative and flowering stages of cannabis production using conventional fertilizers is approximately 160 mg L -1 Bernstein, 2020, 2021). Plants supplied with N below 160 mg L -1 during the vegetative stage saw reduced photosynthetic capacity and plant growth, and during the flowering stage saw reduced inflorescence yield, though cannabinoid concentrations (not total production) were greater at extremely low N rates. The optimal N supply for plants grown with liquid organic fertilizers seems to be higher, with the highest yields being achieved with an organic N supply of approximately 390 and 260 mg L -1 for the vegetative and flowering stages, respectively (Caplan et al., 2017b(Caplan et al., , 2017a. Given the limited number of studies and the relative importance of N on plant growth and development, collecting more information about cannabis response to N are needed to establish more accurate recommendations. Phosphorus nutrition has long been a focus in cannabis cultivation. Growers often supply plants with relatively high P concentrations (up to 200 mg L -1 ) during the flowering stage based on a belief that this promotes flower development. However, there is little evidence to support this practice. A recent study found that cannabis plants in the vegetative stage supplied with 100 mg L -1 P performed similar to those supplied with 30 mg L -1 P (Shiponi and Bernstein, 2021). High P concentration in the nutrient solution creates a situation where environmental pollution from excess P is more likely. Clearly, the practice of supplying cannabis with very high concentrations of P needs to be evaluated.
While there are no published studies examining effect of K on inflorescence quality, some recent studies have looked at how K impacts inflorescence yield. Yield of aquaponically-grown cannabis (g/plant) increased linearly with increasing nutrient solution K concentration in the range of 15-150 mg L -1 (Yep and Zheng, 2020). The nitrogen concentration (75 mg L -1 ) used by Yep and Zheng (2020) reflects that of a typical aquaponic solution, but this N concentration is fairly low compared to conventional hydroponic nutrient solutions and may have been a limiting factor for plant growth and yield (Yep et al., 2020b). For the vegetative stage, cannabis plants supplied with 15 mg L -1 K had reduced growth and displayed foliar symptoms characteristic of K deficiency, while plants that received 60-240 mg L -1 K produced substantially more biomass and did not display K deficiency symptoms (Saloner et al., 2019). While there is a lack of recommendation based on scientific research, some fertilizer companies are recommending 300-400 mg L -1 K. More research is needed to determine the optimal nutrient solution K concentration during cannabis flowering in soilless production systems when other nutrient elements are not limiting.
A challenge in developing fertilizer recommendations is that the number of combinations of nutrient concentrations that can be empirically tested is limited due to logistical and statistical considerations. As a result, most nutrient studies have a limited range of nutrient compositions that can overlook potential nutrient interactions across a broad range of nutrient compositions. Studies on cannabis response to nutrients so far have either investigated different concentrations of one nutrient while holding the others constant (Saloner et al., 2019;Bernstein, 2020, 2021;Shiponi and Bernstein, 2021), or provided different concentrations of NPK in a set ratio (Caplan et al., 2017b(Caplan et al., , 2017aBernstein et al., 2019). Neither of these approaches can evaluate nutrient interactions, which could have substantial impacts on the recommendations of optimum application rates.
Response surface methodology (RSM) is an alternative experimental design capable of concurrently optimizing multiple factors over a wide range of levels using fewer experimental units compared to traditional designs (Myers et al., 2016). The efficiency of this design is achieved through using the fewest number of experimental units which conserves space, time, and resources. Nutrient solution optimization has been approached by some researchers as a 'mixture system' which is a type of multifactor optimization similar to response surface analysis Schrevens, 1998, 1999a). However, the experimental design of a mixture system only optimizes the nutrient composition of the solution but not the nutrient concentration as the design maintains a constant total concentration of nutrients in the solution. RSM allows the optimization of both the nutrient solution composition and the concentrations of individual components without this limitation. Given the high cost of cannabis and growing space being limited to government-approved production facilities, the reduced number of experimental units required for a RSM approach is an advantage over conventional experimental designs.
The objective of this study was to determine the optimal concentrations of NPK for the flowering stage of cannabis in a soilless production system using the RSM approach.

Plant material and growing conditions
The experiment was conducted in a controlled-environment growth room at a Health Canada approved cannabis production facility located in Southern Ontario. A clonal selection of a high delta-9-tetrahydrocannabinol (THC), low cannabidiol (CBD) Cannabis sativa cultivar 'Gelato' was used for this trial. Plants were grown in deep-water culture (DWC) systems. Each DWC unit used a 19 L white plastic bucket (36 cm height x 30.5 cm top outside diameter x 26.4 cm bottom outside diameter) as the nutrient solution reservoir. DWC units were placed on the floor in five double rows of ten DWCs each (i.e., 100 DWCs total), each with one plant, spaced ten cm between adjacent units, 15 cm within the rows, and a one metre aisle-space between rows (Figure 1). Uniform two-week-old cuttings (~15cm tall, 5-6 nodes trimmed to 3-4 leaves) rooted in rockwool cubes were transplanted into each DWC unit using a mesh pot (FHD Plastics, 0.62 L, 10.3 cm height x 12.5 cm diameter) filled with 8-16 mm expanded clay pebbles (Liapor, Hallerndorf, Germany) and inserted flush to the top of the bucket lids, with the bottom three cm of the mesh pot submerged in the nutrient solution. Each DWC bucket was supplied with nutrient solution and had an air-stone (Pawfly ASC030, 30 mm height x 18 mm diameter) providing 1.5 litres of air per minute to continuously mix and aerate the solution. The nutrient solutions in all DWC units were drained and replaced with 17 L of fresh nutrient solution weekly. The initial pH of the nutrient solutions was adjusted to 5.6 with 1 M sulphuric acid or 1 M sodium hydroxide, as needed. DWC units were topped up with pH-adjusted (5.6) rainwater three to four days after each weekly nutrient solution replacement to replace water lost due to evapotranspiration. Nutrient solution pH and electrical conductivity (EC, mS cm -1 ) were measured using a hand-held meter (BLU2300E Combo Meter, Bluelab Corporation, New Zealand). EC and pH of treatment feed solution and of the final leachate are listed in Table 1. 1.9 1.9 5.6 5.9 16 1.9 1.8 5.6 6.0 17 1.9 1.8 5.6 5.6 18 1.9 1.8 5.6 6.1 19 1.9 1.8 5.6 5.9 20 1.9 1.9 5.6 5.9

Figure 1. Rows of deep-water culture units containing trial plants at the end of the three-week vegetative stage.
Plants were grown in the DWC systems vegetatively, under 18/6-h light/dark conditions, for three weeks before switching to a 12/12-h light/dark (i.e., short-day) photoperiod, to induce flowering. Plants were grown under short-day conditions for seven weeks before being harvested. Light was provided by 1000 W metal halide bulbs at an average canopy-level photosynthetically photon flux density of 570 μmol m -2 s -1 . The air temperature and relative humidity were set at 25 C and 65%, respectively. There was no CO2 supplementation.
The experimental unit was one plant grown in its individual DWC unit. There were five replicates per treatment. Plants were randomly assigned to each nutrient solution treatment by generating a random sequence of numbers from 1 to 100 arranged in ten columns and ten rows (matching DWC unit arrangement). For the first three weeks following transplant, all plants received the same nutrient solution containing (mg L -1 ): 112.8 N-NO3, 7.2 N-NH4, 40 P, 180 K, 110 Ca, 45 Mg, and 60 S. Once switched to short-day conditions, plants received their respective treatment nutrient solutions for the remainder of the experiment. The major cation and anion compositions of the treatment nutrient solutions are detailed in Table 3. All treatments were formulated to have the same N-NH4/N-NO3 ratio (1:16

Aboveground growth
Plant height and spread of the first three plants in each treatment were measured during the fifth week of the flowering stage. Plant height (cm) was measured from the lid of the DWC unit to the top of the apical inflorescence, and plant spread (cm) was measured at the widest point on the plant and then perpendicular to this measurement. Growth index (GI) was then calculated using the formula (GI = (height x width1 x width2) / 300) (Caplan et al., 2017a). Plants were destructively harvested during the eighth week of flowering. To assess aboveground (including inflorescence) fresh weight (FW), plants were cut at substrate level and individually weighed on a digital balance.

Root weight
During harvest, roots from the first three replicates of each treatment were cut from around the outer surface of the mesh pot and air dried for several days and then oven-dried at 92C for 72 hours and weighed (EG2200-2NM, KERN & SOHN, Balingen, Germany) to obtain root dry weight (DW).

Inflorescence yield
Inflorescence material was trimmed of leaf tissue, removed from the stem, and then weighed to obtain inflorescence fresh weight (g/plant). To determine inflorescence dry weight (i.e., yield), ~25 g samples of fresh inflorescence material from the first three plants in each treatment were weighed, dried at 70C for 72 hours, and then re-weighed to obtain dry weight (DW). Yield was computed on a per-plant basis as the total inflorescence FW × (sample DW/sample FW). Cured 'whole-bud' cannabis inflorescence sold commercially normally contains 10 to 15 % water. Therefore, the marketable yield can be calculated from inflorescence DW by factoring in the appropriate water content.
The cannabinoid analysis was conducted using ultra performance liquid chromatography (UPLC) separation. The composite sub-sample of dried cannabis was milled to a fine powder; from which1.0g was extracted with an Acetontrile/H2O mixture with sonication and agitation for 20 minutes at ambient temperature. A 1.5 mL aliquot was diluted and filtered into a HPLC vial and analysed as per 7020006509EN (Layton and Aubin, 2019).

Statistical design and analysis
RStudio software (RStudio Team, 2020) was used for data analysis. Normality and homoscedasticity of the data was assessed, and the data met these assumptions. The RStudio package 'rsm' (Lenth, 2009) was used to analyze inflorescence yield and to generate threedimensional and contour plots to represent the response surface. To improve precision of yield estimates, the average yield of the five replicates in each treatment was used. Two sets of three surface and contour plots were created, each while holding one of the nutrient concentrations fixed at its centre point. These surface and contour plots, along with canonical analysis, were then used to determine the optimal rate of all three factors. Correlation analysis of yield and vegetative parameters was performed using the RStudio software package 'ggplot2' (Wickham, 2016). To determine if there were differences in inflorescence cannabinoid content attributable to treatment, data from cannabinoid analysis was tested with a one-way ANOVA followed by Tukey's HSD post-hoc test.

Inflorescence yield response
Cannabis inflorescence yield responded to increasing N and P supply but did not respond to K within the tested range (Figures 2 and 3). Based on the surface response model, the highest average yield of 144 g/plant was achieved with N and P concentrations of 194 mg L -1 and 59 mg L -1 , respectively. Visual analysis of contour graphs (with a 5 gram resolution) show that yield responded to N best in the range of 160-230 mg L -1 , and P in the range of 40-80 mg L -1 (Figure  2).

Figure 2. Three-dimensional response surfaces for inflorescence yield (g/plant) at a range of nutrient solution N, P, and K concentrations (mg L -1 ) of Cannabis sativa grown in deep water culture (P ≤ 0.05, R 2 = 0.57). (A) Surface plot of K vs. P at N = 180 mg L -1 (B) Surface plot of K vs. N at P = 60 mg L -1 (C) Surface plot of P vs. N at K = 200 mg L -1 .
A B Preprints (www.preprints.org) | NOT PEER-REVIEWED | Posted: 24 August 2021 doi:10.20944/preprints202108.0460.v1

Cannabinoid content
There were no nutrient treatment effects on the inflorescence cannabinoid content. The average cannabinoid contents are listed in Table 4. In addition to those cannabinoids listed, the following were below the detection limits (i.e., < 0.5 mg/g): CBC, CBD, CBDA, CBN,  8 THC.
B C

Relationships between inflorescence yield and vegetative growth attributes
No nutrient deficiency or toxicity symptoms were observed on any plants. Inflorescence yield was linearly and positively correlated with the measured vegetative growth attributes. Inflorescence yield had significant correlations with aboveground plant fresh weight (Figure 4), plant growth index (Figure 5), and root dry weight (Figure 6).

Discussion
The goal of this study was to determine the optimal concentration of N, P, and K in the nutrient solution for the flowering stage of soilless cannabis production using RSM. The optimal concentrations of nutrient solution N and P was predicted to be approximately 194 mg L -1 N, 59 mg L -1 P, respectively. Based on analysis of the response surface model, it was found N and P were the most important factors in predicting inflorescence yield. Inflorescence yield decreased markedly outside of the range of 160-240 mg L -1 N, and 40-80 mg L -1 P. These findings suggest that drug-type cannabis responds well to nitrogen and phosphorus during the flowering stage. Inflorescence yield did not respond to nutrient solution K concentration within the tested range, indicating the K currently supplied (300-400 mg L -1 ) by some commercial cultivators are likely too high.
Inflorescence yield had a strong positive correlation with a number of vegetative growth attributes. The strong correlation between inflorescence yield and plant growth index indicates that larger plant size can result in higher inflorescence yield. Nutrient supply, especially N, can determine cannabis plant size as N is an essential component of plant chlorophyll and ribulose-1,5-bisphosphate carboxylase-oxygenase (Rubisco). Low levels of N can reduce plant photosynthetic capacity and limit plant growth (Saloner and Bernstein, 2020). For flowering drug-type cannabis in soilless culture, supply of 30 and 80 mg L -1 N restricted whole plant and inflorescence growth, but plants performed optimally with supply of 160-320 mg L -1 N (Saloner and Bernstein, 2021). The optimal N supply (194 mg L -1 ) found in our study is within their range, despite the two studies using two different growing methods and plants with different genetic backgrounds. For drug-type cannabis during the flowering stage in an organic-based soilless production system, the optimal N supply was slightly higher (212 to 261 mg L -1 ; Caplan et al., 2017a) than the optimal level found in the present study. A possible explanation for the higher optimal N supply in the organic fertilizer study is that N from organic-based fertilizers may not always be readily available, as the release of N from organic fertilizers depends on the speed and extent of the mineralization process (Hartz et al., 2010;Dion et al., 2020). Though it is unclear what source of organic nitrogen was used in their study, factoring in organic N availability of around 60% would put our findings in line with those by Caplan et al. (2017a). Along with aboveground growth, root growth also contributes to overall plant size. We found that inflorescence yield had a strong positive correlation with root dry weight, supporting our conclusion that larger plants produce higher yields. The context of where plants spend their energy is important. For industrial hemp, increasing N supply increased plant growth, but this growth was partitioned more toward stem material rather than valuable inflorescence material (Campiglia et al., 2017). Further investigations of cannabis response to nitrogen should consider product quality, and the distribution of biomass to various plant organs to maximize inflorescence growth and quality.
While modeling of cannabis inflorescence yield response to N, P, and K with surface analysis accounts for interaction between nutrients, the surface response model demonstrated that K, within the tested range of 60-340 mg L -1 , had no effect on inflorescence yield. This lack of response may suggest that 60 mg L -1 K is not low enough to cause nutrient deficiency, and 340 mg L -1 K is not high enough to cause toxicity. Moreover, cannabis responses to K may be cultivar specific. Plants of one cannabis cultivar supplied with 240 mg L -1 K had 25% reduced fresh shoot and root biomass by compared to those fed 175 mg L -1 , while plants of another cultivar had up to 40% increased shoot and root biomass (Saloner et al., 2019). Plant height, number of nodes on the main stem, and stem diameter of these two cultivars remained similar, so this difference in biomass was caused by one cultivar becoming "bushier" than the other under high K supply. These difference in response to K supply may be due to differences in plant tissue (e.g., main stem vs side branch) sensitivity to K. Plant phenological stage (i.e., vegetative or flowering stage) may also be a factor in cannabis response to K supply. In a previous study of flowering aquaponic cannabis response to K, inflorescence yield increased when plants were provided with K up to 150 mg L -1 (Yep and Zheng, 2020). Cultivar and plant phenological stage should be considered in future studies looking at cannabis response to nutrients, especially K.
Many commercial cannabis cultivation operations currently use fertilizer formulations that contain very high levels of P (more than 200 mg L -1 P in some cases). This practice is based on anecdotal evidence that P enhances inflorescence production. These concentrations are much higher than the optimal rate of 60 mg L -1 P found in our study, and at the higher range could cause reduction of both plant growth and inflorescence yield. In addition to reducing plant growth and yield, excessive supply of nutrients is a potential source of environmental pollution. Though, cannabis does appear to have the ability to store and mobilize certain amount of P when required. When provided with P higher than 30 mg L -1 in the vegetative stage, cannabis sequestered excess P in root tissue to prevent excess accumulation in the shoots (Shiponi and Bernstein, 2021). A greater understanding of what cannabis P requirements are, and whether there is any truth to the practice of supplying high concentrations of P, should be a priority for making cannabis production more sustainable. However, based on our data it appears that the levels of P found in many cannabis specific commercial fertilizers are far higher than needed and could lead to negative environmental impacts.
While the cannabinoid concentrations in the floral tissues in our study did not respond to nutrient solution NPK concentrations, other studies indicate that plant mineral nutrition can affect production of secondary metabolites in cannabis (Caplan et al., 2017a;Saloner and Bernstein, 2021). There appears to be an inverse relationship between cannabis yield and potency, with cannabinoid concentrations decreasing as plant inflorescence yield increases. Inflorescence from plants supplied with 160 mg L -1 N had approximately 30% and 20% lower concentrations of THCA and CBDA than plants supplied with 30 mg L -1 N (Saloner and Bernstein, 2021). But while nutrient stress and deficiency may enhance inflorescence cannabinoid content, this method is not ideal for optimizing overall plant productivity as plants supplied with 160 mg L -1 N yielded twice that of those supplied with 30 mg L -1 N. Cannabis grown in two organic growing media with different organic fertilizer rates (i.e., 57, 113, 170, 226 and 283 mg L -1 N) had negative linear relationships between the concentrations of inflorescence THCA and CBGA and the fertilizer application rate for some of the treatment combinations (i.e., growing media and fertilizer rate) (Caplan et al., 2017a). However, for majority of the treatment combinations, fertilizer rates from 57-226 mg L -1 N did not have any effects on THCA or CBGA concentrations; and the cannabinoid concentrations only dropped when the fertilizer rate increased to the highest level of 283 mg L -1 N. The context of yield is again important when analyzing differences in cannabinoid content as THCA concentrations dropped by ~20% in the highest fertilizer rate, but inflorescence yield almost doubled vs. lowest fertilizer rate. As noted by Bernstein et al. (2019), an understanding of how nutrient supply influences cannabinoid concentrations would be an important step toward controlling and standardizing the cannabinoid contents of medical cannabis. Cannabinoid concentrations are also important to recreational consumers, who rank THC and CBD concentrations among the most important factors when making purchasing decisions (Zhu et al., 2020). Given that cannabinoids are the compounds that make cannabis so uniquely valuable, more work needs to be done to investigate the effect of mineral nutrition on cannabis yield, and the relationship between yield and potency.
The use of central-composite design allows experimenters to account for potential interactions between the different nutrients. This is important as nutrient interactions have been shown to affect plant nutrient uptake (Fageria, 2001;Rietra et al., 2017). A recent study found that high K supply decreased concentrations of Ca and Mg in cannabis leaf tissue, indicating antagonistic relationships between these positively charged ions (Saloner et al., 2019). An understanding of how combinations of nutrients at different concentrations affect crop growth and yield is important for the development of recommendations for the commercial cannabis industry. Had the same number of nutrients and nutrient levels as were included in this study been investigated with a traditional full-factorial design, many more nutrient solution treatment groups would have been required, compared to the number of treatment groups used in this study. The difference in number of treatment groups needed can be more pronounced as more factors (i.e., Ca, Mg) are included. Considering the high cost of cannabis and growing space in controlled environments, the response surface approach allowed us to complete this study where another experimental design may have been prohibitive.
No matter the experimental design used, an inherent problem in nutrient solution experiments is that nutrients cannot be added individually but must be added as a compound containing both anions and cations. Further, the ionic balance constraint requires the sum of the charges of cations and anions in solution to be equal (De Rijck and Schrevens, 1999b). The implication for formulating experimental treatment solutions is that it is practically impossible to change the level of one nutrient while keeping concentrations of all other nutrients the same. In this study, we focused on N, P, and K concentrations while attempting to keep all other nutrients at reasonable levels using commonly available horticultural fertilizer compounds. For example, potassium nitrate and calcium nitrate usually contribute the bulk of nitrogen, potassium, and calcium in horticultural nutrient solutions (Resh, 2012). Formulating a high N, low K nutrient solution with these fertilizers results in higher levels of Ca than other nutrient solution treatments. Likewise, a low N, high K nutrient solution necessitates an additional source of K such as KCl, which would increase solution Cl concentration. Higher concentrations of nutrients such as Ca and Cl bring the potential for nutrient interactions which may affect experimental results. The lack of response to K in the range of 60-340 mg L -1 observed in our trial may be partially due to competition for uptake from Ca. Regarding experimental Cl levels, hydroponic cannabis has been shown to tolerate rates of 180 mg L -1 Cl with no impact on yield or potency (Yep et al., 2020a) so it is unlikely Cl levels limited plant growth in this trial. Though less than ideal in an experimental setting, there is no perfect solution for the problem of keeping all nutrient concentrations the same when formulating treatment solutions.
While this trial determined the theoretical optimum levels of N and P for the DWC growing method, these levels may not be definitive for all production methods. Our trial was conducted in solution culture with weekly nutrient solution changes, and the EC and pH dynamics of our DWC units are likely different than other growing methods, meaning that plant nutrient availability and overall salinity of the nutrient solution would also likely be considerably different. Many commercial cannabis operations utilize substrate-based soilless cultivation systems, such as coir in containers, that may offer more nutrient and pH buffering capacity (Zheng, 2020). Having said that, our trial does represent or closely resemble some common soilless production practices, such as growing cannabis in rockwool, in the current cannabis production industry (Zheng, 2021). Another limitation of our study was that we only used a single cannabis cultivar. Similar experiments should be performed on different cultivars, with disparate growth habits and cannabinoid compositions to investigate how individual cultivars may respond to NPK treatment levels. Additionally, this study only looked at inflorescence yield and cannabinoid composition and did not evaluate the impact of NPK on inflorescence terpene content or organoleptic properties.
Drug-type cannabis is still a relatively new crop in the legal setting, especially for large-scale commercial production, and many aspects of its cultivation are relatively unknown. We found that response surface methodology was a suitable experimental approach for investigation of cannabis responses to NPK, and that modeling of yield response to these nutrients aided us in achieving our experimental objective. Based on the results of this study, we recommend providing plants with a nutrient solution containing N and P at approximately 194 and 59 mg L -1 , respectively, to achieve maximal inflorescence yield. Future studies should investigate the inflorescence yield and vegetative growth response of genetically diverse cultivars to macronutrients and include more quality parameters to ensure that plant yields do not compromise product quality. Improving our understanding of cannabis responses to mineral nutrients is an essential step towards the effective and sustainable cultivation of this high-value horticultural crop.