Optimisation of Nitrogen, Phosphorus, and Potassium for Soilless Production of Cannabis sativa in the Flowering Stage Using Response Surface Analysis

Following legalisation, 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 optimisation.


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 to research on this crop. However, change in social attitudes towards 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 maximise crop yields, quality, and profits while minimising environmental impacts.
Proper supply of mineral nutrients is essential for efficient and sustainable cultivation of any crop. Among the most important 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 often rely on nutrient recipes developed by fertiliser 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., 2017a,b) 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 synchronise nutrient supply and demand to maximise 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 fertilisers is approximately 160 mg L −1 (Saloner and Bernstein, 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 fertilisers 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., 2017a,b). 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 high P 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 high concentrations of P needs to be evaluated.
While there are no published studies examining the 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). Although there is a lack of recommendations based on scientific research, some fertiliser 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 fertiliser 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., 2017a,b;Bernstein 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 optimising 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 by using fewer experimental units which conserves space, time, and resources. Nutrient solution optimisation has been approached by some researchers as a "mixture system" which is a type of multifactor optimisation similar to response surface analysis Schrevens, 1998, 1999a). However, the experimental design of a mixture system only optimises the nutrient composition of the solution but not the overall nutrient concentration as the design maintains a constant total nutrient supply in the solution. RSM allows the optimisation 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) C. 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 × 30.5 cm top outside diameter × 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 2-week-old cuttings (∼15 cm 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 × 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 × 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.
Plants were grown in the DWC systems vegetatively, under 18/6-h light/dark conditions, for 3 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 7 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 CO 2 supplementation.

Experimental Design and Treatments
A three-factor (i.e., N, P, and K), second order central rotatable composite design was used to model cannabis responses to TABLE 1 | Coded and un-coded factors for each treatment of the response surface analysis experiment investigating the effect of nitrogen (N), phosphorus (P), and potassium (K) rates on cannabis production in deep-water culture.

Coded factors
Un-coded factors (mg L −1 )  these mineral nutrients. Following a response surface design (Table 1), treatments combinations were defined by their concentrations of N (70, 120, 180, 250, and 290 mg L −1 ), P (20, 40, 60, 80, and 100 mg L −1 ), and K (60, 120, 200, 280, and 340 mg L −1 ) ( Table 2) achieved by using different amounts of varies straight fertilisers ( Table 3). The experimental unit (replicate) was one plant grown in an individual DWC unit. There were 15 different treatments, with at least 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  a N-NH 4 + N-NO 3 . b Radius adjustment factor for a three-factor design to make the design rotatable. 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 3-4 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 metre (BLU2300E Combo Metre, Bluelab Corporation, New Zealand). EC and pH of treatment feed solution and of the final drained solution are listed in Table 5.

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 × width 1 × width 2 )/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 h 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 h, 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 "wholebud" 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.

Statistical Design and Analysis
RStudio software (RStudio Team, 2020) was used for data analysis. Normality and homoscedasticity of the data were assessed, and the data met these assumptions. The RStudio package "rsm" (Lenth, 2009) was used to analyse inflorescence yield and to generate three-dimensional and contour plots to represent the response surface. To improve the 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, 3). Based on the surface response model, the estimated highest average yield of 144 g/plant would be achieved with N and P concentrations of 194 and 59 mg L −1 , respectively. Visual analysis of contour graphs (with a 5 g 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).

Cannabinoid Content
There were no nutrient treatment effects on the inflorescence cannabinoid content. The average cannabinoid contents are listed in Table 6. 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.

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, and 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-230 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,5bisphosphate 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 organicbased soilless production system, the optimal N supply was slightly higher (212-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 fertiliser study is that N from organic-based fertilisers may not always be readily available, as the release of N from organic fertilisers depends on the speed and extent of the mineralisation 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 towards 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 maximise inflorescence growth and quality.
While modelling 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 genotype specific. Plants of one cannabis genotype Royal Medic supplied with 240 mg L −1 K FIGURE 4 | Correlation between inflorescence yield (g/plant) and aboveground plant fresh weight (g/plant) of Cannabis sativa (r = 0.98, P < 0.001). Shaded area represents 95% confidence interval.
had 25% reduced fresh shoot and root biomass by compared to those fed with 175 mg L −1 , while plants of genotype Desert Queen 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 genotypes remained similar, so this difference in biomass was caused by one genotype becoming "bushier" than the other under high K supply. These differences in the 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). Genotype and plant phenological stage should be considered in future studies looking at cannabis response to nutrients, especially K.  Inflorescence yield (g) FIGURE 6 | Correlation between inflorescence yield and root dry weight of Cannabis sativa (r = 0.9, P < 0.001). Shaded area represents 95% confidence interval.
Many commercial cannabis cultivation operations currently use fertiliser 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 mobilise 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 existing data it appears that the levels of P found in many cannabis specific commercial fertilisers 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). However, while nutrient stress and deficiency may enhance inflorescence cannabinoid content, this method is not ideal for optimising 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 fertiliser 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 fertiliser application rate for some of the treatment combinations (i.e., growing media and fertiliser rate) (Caplan et al., 2017a). However, for the most of the treatment combinations, fertiliser rates from 57 to 226 mg L −1 N did not have any effects on THCA or CBGA concentrations; and the cannabinoid concentrations only dropped when the fertiliser rate increased to the highest level of 283 mg L −1 N. The context of yield is again important when analysing differences in cannabinoid content as THCA concentrations dropped by ∼20% in the highest fertiliser rate, but inflorescence yield almost doubled vs. lowest fertiliser rate. As noted by Bernstein et al. (2019), an understanding of how nutrient supply influences cannabinoid concentrations would be an important step towards controlling and standardising 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. Further work should also evaluate other compounds that are known to impact product quality.
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, yield, and quality 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 fertiliser compounds. For example, potassium nitrate and calcium nitrate usually contribute to the bulk of nitrogen, potassium, and calcium in horticultural nutrient solutions (Resh, 2012). Formulating a high N, low K nutrient solution with these fertilisers 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 or genotypes. 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 utilise 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). The treatments were applied only during the short-day period (i.e., flowering stage), and considering that plant nutrient requirement may vary at different development stages, the same experiment may also need to be conducted for the vegetative stage. 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 modelling 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.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.