Patent Publication Number: US-2022229036-A1

Title: System and method for determining aqueous nitrate concentration in solution containing dissolved organic carbon

Description:
CROSS-REFERENCE TO RELATED APPLICATIONS 
     This application is a National Stage application of International Patent Application No. PCT/IL2020/050645, filed on Jun. 11, 2020, which claims priority to U.S. Patent Application No. 62/860,273, filed on Jun. 12, 2019 each of which is hereby incorporated by reference in its entirety. 
    
    
     TECHNICAL FIELD 
     The invention relates in general to systems and methods for optimizing agricultural crops&#39; yields while reducing water contamination due to excess application of fertilizers. More specifically, the invention relates to a system and method for determining a level of nitrate in a water sample, thereby to determine accurate amount of fertilization necessary. 
     BACKGROUND 
     Contamination of rivers, lakes, fresh water, drinking water, groundwater and soil&#39;s porewater by nitrate is a global problem. The term “nitrate” is briefly referred to herein also as “N”. It is universally recognized that nitrate contamination of drinking water is a threat to human health. It is also reported that human health suffers from adverse effects, even cancer, due to continual exposure to nitrate above a certain level. A significant rate of water pollution results from an excess of fertilization by farmers, due to a lack of real-time and accurate information with respect to the availability of nutrient in the soil. Therefore, an excess of fertilization results not only in a waste of resources, but also with a pollution of the groundwater, particularly by nitrate. 
     During the second half of the 20th century, clear trends of rising nitrate concentration in groundwater have been observed in aquifers all over the globe. The World Health Organization (WHO) had determined that nitrate levels in drinking water should not exceed 50 ppm. When exceeding this level of concentration, nitrate is harmful to infants, where it may cause “blue-babies syndrome” (methemoglobinemia) and can lead to severe illness and even death. Unfortunately, nitrate contamination is the most dominant factor responsible to severe degradation of groundwater and surface resources. On a global scale, eutrophication and hypoxia of streams, rivers and lakes, is mostly attributed to subsurface return flow from nitrate contaminated groundwater, leaking from phreatic aquifers underlying agricultural fields. Moreover, the impact of nitrate contaminated groundwater is not limited to terrestrial water resources, and has a great impact on marine ecosystems as well. Eutrophication and hypoxia on a large scale has been found in the Gulf of Mexico and the Black Sea, as well as severe impact on the Great Barrier Reef, Australia. Overall, nitrate contamination had led to more groundwater disqualification and water well shutdowns than any other contaminant, worldwide. While nitrate is considered the most common non-point source pollutant in groundwater, numerous studies have linked the increase of nitrate concentration in groundwater to excess use of fertilizers in agriculture. As a result, a global regulatory takes place by environmental protection and water authorities to reduce excessive application of agricultural fertilizers. For example, the European Union had established the Nitrates Directive, and the US Environmental Protection Agency (EPA) regards nitrate contamination in groundwater as an event requiring immediate action. 
     At present, fertilizer application in agriculture relies primarily on farmer&#39;s experience, expert&#39;s recommendation, and sporadic soil testing. None of these techniques provide information that is in line with the time scale of N-fertilizers mobilization (nitrate&#39;s solution movement rates through the soils sediments), consumption and transformation dynamics in the soil. 
     Presently, the monitoring of chemical parameters in soils is performed in water samples, that may be obtained, for example, by a suction cup which is mounted in the soil, or by extracting soil samples. Water samples collected by this mechanism is typically transferred to a laboratory for further chemical analysis, or is analyzed on-site by means of an analytical kit. However, nitrate in the soil is highly soluble, mobile, and is obviously consumed by the crops. Moreover, nitrate concentration in the soil may fluctuate in time-scales of hours to days, as a response to different irrigation schemes, precipitation, fertilization, root uptake and different plant growth phases. As such, the monitoring of nitrate concentration by conventional tools does not meet the required time resolution for optimizing fertilization schemes, while preventing groundwater pollution due to excess application of fertilizers. Moreover, current techniques typically require handling of the sampling schemes by a devoted research team, not by farmers. 
     Ultraviolet (UV) absorption spectroscopy is one of the most common methods for nitrate analysis in aqueous solution. Nitrate in aqueous solution absorbs light at two main wavelength regions: (a) a high-absorbance band in between 200-240 nm; and (b) a low absorbance-band in between 280-320 nm. While the absorbance in both bands is known to have direct correlation to nitrate concentration, the absorbance at 220 nm is two orders of magnitude higher when compared to the absorbance at 300 nm. 
     Tuly et al., (2009), “In Situ Monitoring of Soil Solution Nitrate: Proof of Concept”, https://www.researchgate.net/publication/231523625, suggests a technique for continuous monitoring of nitrate concentrations in soil solution. An absorbance spectroscopy is applied on a sample within a stainless-steel porous cup, which is installed in the soil. The porous cup is filled with deionized water, which is placed in a reservoir of potassium nitrate solution. Once the solution inside of the cup achieves chemical equilibrium by diffusion between the porous cup and the surrounding medium, the absorption spectrum of the solution is measured by means of a UV dip probe. The proposed setup of Tuly et al, (2009) uses a UV light source and a dip probe that are connected to spectrophotometer via optical fibers, to continuously determine the nitrate concentration within the cup. However, this technique has a limited applicability, for two main reasons: (a) the obtaining of chemical equilibrium between the porous cup and the surrounding medium, especially in unsaturated sediment with a limited water storage, is rather slow, resulting in a time lag between the actual variation of the nitrate concentration in the soil to its actual measurement. Therefore, rapid concentration variations, as expected following intensive irrigation or fertilization events may not be recorded; and (b) the presence of natural soil Dissolved Organic Carbon (DOC) limit the accuracy of UV absorption by means of spectroscopy analysis, since both nitrate and soil DOC absorb UV light in overlapping wavelengths ranges. 
     WO 2018/104939, Yeshno et al., 2019 suggests a nitrate monitoring technique which is based on a continuous spectral analysis of soil porewater in an optical flow cell. The optical flow cell is connected to a porous interface which obtains a continuous flux of soil porewater. The absorption spectrum of the soil porewater is continuously recorded and analyzed to determine in real-time the nitrate concentration. The analysis involves a scan of the absorption spectrum of the soil porewater to identify an optimal wave length where DOC interference to nitrate measurement is minimal. The system for carrying out this technique requires for its operation a wide-band UV spectrophotometer and accordingly a wide-range, deuterium UV lamp. However, in spite of the capability of the system to continuously measure nitrate concentration in-situ, the system&#39;s large-dimensions, along with its high energy consumption, is too bulky and expensive for practical and commercial applications. 
     In brief, the prior art optical system of WO 2018/104939 for determining a level of nitrate in a cultivated soil basically requires the following components:
         a. One or more sampling cells, each enabling the extraction of a porewater sample from a specific region of the soil;   b. Each sampling cell is connected to an optical flow cell for continuous or frequent spectral analysis of the soil porewater;   c. A UV light source for applying a light beam in wavelengths between 200-240 nm through the water sample in the flow cell;   d. A photodetector or a spectrophotometer for accumulating light fromthe light beam, following its passage through the water sample;   e. A processing unit for determining the level of light-absorption (based on Beer-Lamb ert equation), and
           f. A processing unit to estimate nitrate concentrations based on a previously prepared empiric calibration equation. The calibration equation 1 s obtained by an algorithm which locates an optimal wavelength, where a minimal interference from DOC is found;   
           g. WO 2018/104939 also shows how the system can perform analyses on a plurality of regions, while each time a single region is selected for analyses.       

     It is an object of the present invention to provide a system for determining in real time and in-situ a level of nitrate in soil, which is simpler in structure, compact in size, and of lower cost compared to similar prior art systems. 
     It is a particular object of the invention to eliminate the effects of DOC existing in a porewater sample on the absorption spectrum, thus to enable the use of absorption spectroscopy technique to estimate aqueous nitrate concentration in the presence of DOC in the solution. 
     Other objects and advantages of the invention will become clear as the description proceeds. 
     SUMMARY 
     The invention relates to a system for determining a level of nitrate in a water sample, comprising: (a) an optical flow cell which is at least partially transparent and which is configured to contain a sample of water; (b) a first illuminator for illuminating the sample within the cell by light in a first wavelength, and a first photodetector for collecting the first-wavelength illumination, following the light passage through the sample; (c) a second illuminator for illuminating the sample within the cell by light in a second, fluorescence-excite wavelength, and a second photodetector for collecting illumination in a third, fluorescence-emission wavelength from the sample; and (d) an analysis unit for: (d.1) determining an overall effect of nitrate+DOC (hereinafter, the terms “nitrate+DOC” and “nitrate and DOC” are used interchangeably, and their essence is the same) within the sample on the absorbance of the sample, said effect of nitrate+DOC being proportional to a rate of absorbance of light due to said illumination by said first illuminator, said absorbance being determined from a difference between a level of illumination by said first illuminator and a level of collected illumination by said first photodetector; (d.2) determining a concentration of DOC within the sample, said DOC concentration being proportional to an intensity of said fluorescence emission from the sample due to said illumination by said second illuminator, and as collected by said second photodetector; and (d.3) subtracting said effect of DOC from said effect of nitrate+DOC on absorbance, thereby to determine the concentration of nitrate within the sample. 
     In an embodiment of the invention, the system further comprising a first look-up table, for converting the difference as measured in step (d.1) to a nitrate and DOC concentration level. 
     In an embodiment of the invention, the system further comprising a second look-up table for converting said fluorescence emission as measured in step (d.2) to a DOC concentration level. 
     In an embodiment of the invention, the system further comprising a third look-up table, for calibrating the subtraction result based on a specific type of DOC known to be in the specific tested sample, wherein said type of DOC reflects a specific chemical DOC composition. 
     In an embodiment of the invention, a mathematical equation is used to convert absorbance and/or fluorescence measurements to concentration levels. 
     In an embodiment of the invention, the system further comprising one or more additional illuminators, and one or more additional photodetectors, in order to measure absorbance and/or fluorescence emission in additional wavelengths, thereby to determine a specific type of DOC within the sample. 
     In an embodiment of the invention, the analysis unit comprises a mathematical model to extract the value of nitrate based on the measurements of absorption and fluorescence, wherein the mathematical model comprising: 
     
       
         
           
             
               
                 
                   
                     C 
                     
                       NO 
                       3 
                       
                         - 
                         1 
                       
                     
                   
                   = 
                   
                     
                       
                         
                           
                             C 
                             
                               D 
                               ⁢ 
                               O 
                               ⁢ 
                               C 
                             
                           
                           ⁡ 
                           
                             ( 
                             
                               
                                 
                                   
                                     A 
                                     
                                       N 
                                       ⁢ 
                                       
                                         O 
                                         3 
                                         
                                           - 
                                           1 
                                         
                                       
                                     
                                   
                                   ⁡ 
                                   
                                     ( 
                                     
                                       λ 
                                       1 
                                     
                                     ) 
                                   
                                 
                                 + 
                                 
                                   
                                     A 
                                     
                                       D 
                                       ⁢ 
                                       O 
                                       ⁢ 
                                       C 
                                     
                                   
                                   ⁡ 
                                   
                                     ( 
                                     
                                       λ 
                                       1 
                                     
                                     ) 
                                   
                                 
                               
                               
                                 
                                   
                                     A 
                                     
                                       N 
                                       ⁢ 
                                       
                                         O 
                                         3 
                                         
                                           - 
                                           1 
                                         
                                       
                                     
                                   
                                   ⁡ 
                                   
                                     ( 
                                     
                                       λ 
                                       2 
                                     
                                     ) 
                                   
                                 
                                 + 
                                 
                                   
                                     A 
                                     
                                       D 
                                       ⁢ 
                                       O 
                                       ⁢ 
                                       C 
                                     
                                   
                                   ⁡ 
                                   
                                     ( 
                                     
                                       λ 
                                       2 
                                     
                                     ) 
                                   
                                 
                               
                             
                             ) 
                           
                         
                         ⁢ 
                         x 
                         ⁢ 
                         
                           
                             ɛ 
                             
                               D 
                               ⁢ 
                               O 
                               ⁢ 
                               C 
                             
                           
                           ⁡ 
                           
                             ( 
                             
                               λ 
                               2 
                             
                             ) 
                           
                         
                       
                       - 
                       
                         
                           ɛ 
                           
                             D 
                             ⁢ 
                             O 
                             ⁢ 
                             C 
                           
                         
                         ⁡ 
                         
                           ( 
                           
                             λ 
                             1 
                           
                           ) 
                         
                       
                     
                     
                       
                         
                           ɛ 
                           
                             NO 
                             3 
                             
                               - 
                               1 
                             
                           
                         
                         ⁡ 
                         
                           ( 
                           
                             λ 
                             1 
                           
                           ) 
                         
                       
                       - 
                       
                         
                           ( 
                           
                             
                               
                                 
                                   A 
                                   
                                     N 
                                     ⁢ 
                                     
                                       O 
                                       3 
                                       
                                         - 
                                         1 
                                       
                                     
                                   
                                 
                                 ⁡ 
                                 
                                   ( 
                                   
                                     λ 
                                     1 
                                   
                                   ) 
                                 
                               
                               + 
                               
                                 
                                   A 
                                   
                                     D 
                                     ⁢ 
                                     O 
                                     ⁢ 
                                     C 
                                   
                                 
                                 ⁡ 
                                 
                                   ( 
                                   
                                     λ 
                                     1 
                                   
                                   ) 
                                 
                               
                             
                             
                               
                                 
                                   A 
                                   
                                     N 
                                     ⁢ 
                                     
                                       O 
                                       3 
                                       
                                         - 
                                         1 
                                       
                                     
                                   
                                 
                                 ⁡ 
                                 
                                   ( 
                                   
                                     λ 
                                     2 
                                   
                                   ) 
                                 
                               
                               + 
                               
                                 
                                   A 
                                   
                                     D 
                                     ⁢ 
                                     O 
                                     ⁢ 
                                     C 
                                   
                                 
                                 ⁡ 
                                 
                                   ( 
                                   
                                     λ 
                                     2 
                                   
                                   ) 
                                 
                               
                             
                           
                           ) 
                         
                         ⁢ 
                         x 
                         ⁢ 
                         
                           
                             ɛ 
                             
                               NO 
                               3 
                               
                                 - 
                                 1 
                               
                             
                           
                           ⁡ 
                           
                             ( 
                             
                               λ 
                               2 
                             
                             ) 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   3 
                   ) 
                 
               
             
           
         
       
     
     Where A is the measured absorbance at a given wavelength (λ 1,2  (nm)), 
     εNO 3   −I  DOC (λ 1,2  (nm)) is a molar attenuation coefficient for either the nitrate or the DOC (L mol −1  cm −1 ) at a given wavelength (λ 1,2  (nm)), C DOC  is the DOC concentration (mol L −1 ) as obtained from said second photodetector, C NO3  is a nitrate concentration (mol L −1 ), and l is an optical pathlength (cm). 
     In an embodiment of the invention, the analysis unit applies a machine learning technique comprising: (a) generating a plurality of absorption and fluorescence measurements for different values of nitrate concentration and various DOC types and respective concentrations; (b) selecting and adapting one or more deep learning networks; (c) training at least one of the selected deep learning networks; and (d) using the trained network to calculate the nitrate concentration based on absorption and fluorescence measurements. 
     In an embodiment of the invention, the first wavelength is selected from the bands of 200-250 nm and 280-320 nm. 
     In an embodiment of the invention, the second excite wavelength is within a band of 225 nm-600 nm. 
     In an embodiment of the invention, the third, fluorescence emission wavelength is within a band of 250 nm-700 nm. 
     In an embodiment of the invention, the system further comprising a first filter for assuring that radiation only in the first wavelength arrives the first photodetector. 
     In an embodiment of the invention, the system further comprising a second filter for assuring that radiation only in the third wavelength arrives the second photodetector. 
     In an embodiment of the invention, the water sample is taken from a soil or from a water reservoir. 
     In an embodiment of the invention, the water sample is collected from a cultivated soil by a porous interface, and is provided in a low flow-rate through the optical flow cell. 
     The invention further relates to a method for determining a concentration rate of nitrate in a water sample, comprising: (a) providing the sample; (b) illuminating the sample in a first wavelength, and determining the combined effect of nitrate+DOC within the sample on the absorbance, said concentration of nitrate+DOC being proportional to a rate of absorbance of light due to said illumination in said first wavelength, said absorbance being determined from a difference between a level of illumination in said first wavelength before passing the sample and a level of collected illumination in said first wavelength following passage through the water sample; (c) illuminating the water sample in a second, exciting wavelength, and determining an effect of DOC within the sample, said effect of DOC being proportional to an intensity of fluorescence emission from the sample in a third wavelength due to said illumination of the sample in said second wavelength; and (d) deducting the effect of DOC from the overall effect of nitrate+DOC on the absorbance, as determined, thereby to obtain the concentration of nitrate in the sample. 
     In an embodiment of the invention, the method further comprising use of a first look-up table for converting the absorbance to a nitrate+DOC concentration level. 
     In an embodiment of the invention, the method further comprising use of a second look-up table for converting said fluorescence emission to a DOC concentration level. 
     In an embodiment of the invention, the method further comprising use of a third look-up table, for calibrating the deduction result based on a specific type of DOC known to be in the specific tested sample, wherein said type of DOC reflects a specific chemical DOC composition. 
     In an embodiment of the invention, the method further comprising use of a mathematical equation to convert absorbance and/or fluorescence measurements to concentration levels. 
     In an embodiment of the invention, the method further comprising measuring absorbance and/or fluorescence emission in additional wavelengths, thereby to determine a specific type of DOC within the sample. 
     In an embodiment of the invention, the method comprising a mathematical model to extract the value of nitrate based on the measurements of absorption and fluorescence, wherein the mathematical model comprising: 
     
       
         
           
             
               
                 
                   
                     C 
                     
                       NO 
                       3 
                       
                         - 
                         1 
                       
                     
                   
                   = 
                   
                     
                       
                         
                           
                             C 
                             
                               D 
                               ⁢ 
                               O 
                               ⁢ 
                               C 
                             
                           
                           ⁡ 
                           
                             ( 
                             
                               
                                 
                                   
                                     A 
                                     
                                       N 
                                       ⁢ 
                                       
                                         O 
                                         3 
                                         
                                           - 
                                           1 
                                         
                                       
                                     
                                   
                                   ⁡ 
                                   
                                     ( 
                                     
                                       λ 
                                       1 
                                     
                                     ) 
                                   
                                 
                                 + 
                                 
                                   
                                     A 
                                     
                                       D 
                                       ⁢ 
                                       O 
                                       ⁢ 
                                       C 
                                     
                                   
                                   ⁡ 
                                   
                                     ( 
                                     
                                       λ 
                                       1 
                                     
                                     ) 
                                   
                                 
                               
                               
                                 
                                   
                                     A 
                                     
                                       N 
                                       ⁢ 
                                       
                                         O 
                                         3 
                                         
                                           - 
                                           1 
                                         
                                       
                                     
                                   
                                   ⁡ 
                                   
                                     ( 
                                     
                                       λ 
                                       2 
                                     
                                     ) 
                                   
                                 
                                 + 
                                 
                                   
                                     A 
                                     
                                       D 
                                       ⁢ 
                                       O 
                                       ⁢ 
                                       C 
                                     
                                   
                                   ⁡ 
                                   
                                     ( 
                                     
                                       λ 
                                       2 
                                     
                                     ) 
                                   
                                 
                               
                             
                             ) 
                           
                         
                         ⁢ 
                         x 
                         ⁢ 
                         
                           
                             ɛ 
                             
                               D 
                               ⁢ 
                               O 
                               ⁢ 
                               C 
                             
                           
                           ⁡ 
                           
                             ( 
                             
                               λ 
                               2 
                             
                             ) 
                           
                         
                       
                       - 
                       
                         
                           ɛ 
                           
                             D 
                             ⁢ 
                             O 
                             ⁢ 
                             C 
                           
                         
                         ⁡ 
                         
                           ( 
                           
                             λ 
                             1 
                           
                           ) 
                         
                       
                     
                     
                       
                         
                           ɛ 
                           
                             NO 
                             3 
                             
                               - 
                               1 
                             
                           
                         
                         ⁡ 
                         
                           ( 
                           
                             λ 
                             1 
                           
                           ) 
                         
                       
                       - 
                       
                         
                           ( 
                           
                             
                               
                                 
                                   A 
                                   
                                     N 
                                     ⁢ 
                                     
                                       O 
                                       3 
                                       
                                         - 
                                         1 
                                       
                                     
                                   
                                 
                                 ⁡ 
                                 
                                   ( 
                                   
                                     λ 
                                     1 
                                   
                                   ) 
                                 
                               
                               + 
                               
                                 
                                   A 
                                   
                                     D 
                                     ⁢ 
                                     O 
                                     ⁢ 
                                     C 
                                   
                                 
                                 ⁡ 
                                 
                                   ( 
                                   
                                     λ 
                                     1 
                                   
                                   ) 
                                 
                               
                             
                             
                               
                                 
                                   A 
                                   
                                     N 
                                     ⁢ 
                                     
                                       O 
                                       3 
                                       
                                         - 
                                         1 
                                       
                                     
                                   
                                 
                                 ⁡ 
                                 
                                   ( 
                                   
                                     λ 
                                     2 
                                   
                                   ) 
                                 
                               
                               + 
                               
                                 
                                   A 
                                   
                                     D 
                                     ⁢ 
                                     O 
                                     ⁢ 
                                     C 
                                   
                                 
                                 ⁡ 
                                 
                                   ( 
                                   
                                     λ 
                                     2 
                                   
                                   ) 
                                 
                               
                             
                           
                           ) 
                         
                         ⁢ 
                         x 
                         ⁢ 
                         
                           
                             ɛ 
                             
                               NO 
                               3 
                               
                                 - 
                                 1 
                               
                             
                           
                           ⁡ 
                           
                             ( 
                             
                               λ 
                               2 
                             
                             ) 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   3 
                   ) 
                 
               
             
           
         
       
     
     Where A is the measured absorbance at a given wavelength (λ 1,2  (nm)), 
     εNO 3   −I  DOC (λ 1,2  (nm)) is a molar attenuation coefficient for either the nitrate or the DOC (L mol −1  cm −1 ) at a given wavelength (λ 1,2  (nm)), C DOC  is the DOC concentration (mol L −1 ) as obtained from said second photodetector, C NO3  is a nitrate concentration (mol L −1 ), and l is an optical pathlength (cm). 
     In an embodiment of the invention, the first wavelength is selected from the bands of 200-250 nm, or 280-320 nm. 
     In an embodiment of the invention, the second, exciting wavelength is in the order of 225-400 nm. 
     In an embodiment of the invention, the third, fluorescence emission wavelength is in the order of 250 nm-500 nm. 
     In an embodiment of the invention, the water sample is taken from a soil or from a water reservoir. 
     In an embodiment of the invention, the method further applies a machine learning technique, comprising the steps of: (a) generating a plurality of absorption and fluorescence measurements for different values of nitrate concentration and different types of DOC, and their respective DOC concentrations; (b) selecting one or more deep learning networks; (c) training the network and selecting a one with a best performance; and (d) calculating the nitrate concentration based on said absorption and fluorescence measurements. 
    
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
       In the drawings: 
         FIG. 1  generally illustrates in a flow diagram form a method for determining a concentration of nitrate in a water sample, according to an embodiment of the invention; 
         FIG. 2  schematically illustrates a structure of a system for determining a concentration of nitrate in soil porewater, according to an embodiment of the invention; 
         FIG. 3  illustrates a method for determining nitrate concentration in soil solution (or aqueous solution), according to an embodiment of the invention; 
         FIG. 4  exemplifies a complex pattern in the relationship between nitrate concentration, DOC concentration and a UV absorption for soil water as obtained from conventional greenhouse soil water samples; 
         FIG. 5  shows a 3D projection of the experimental results; 
         FIG. 6  shows a calibration curve for a DOC concentration; 
         FIG. 7  shows observed (known) nitrate concentration vs. predicted nitrate concentration by the new analytical concept of the invention, for soil water samples from 6 agricultural sites. 
         FIG. 8  shows observed vs. predicted nitrate concentrations, for the hummus soil mixture water extract, when the calibration equation was obtained for absorption at 235 nm. 
     
    
    
     DETAILED DESCRIPTION 
     The invention provides an optical system for determining a level of nitrate in a cultivated soil, which overcomes drawbacks of similar prior art systems. In brief:
         a. The system can use low-cost and low-energy components, such as LED type light sources and semiconductor type photodetectors;   b. While prior art systems are based on absorption measurement only, the system of the invention measures: (i) absorption resulting from illumination by a first light source to measure a combined effect of both nitrate and DOC; and (ii) a florescence emission resulting from illumination by a second light source to measure the concentration of DOC, and cancel its effect on the absorption measurement by use of a signal processing algorithm. Such a dual-measurement technique enables separation of the effect of DOC from the combined effect of nitrate+DOC on the measurement, thereby the technique enables determination of the actual concentration of nitrate in the water sample. This structure and analysis eliminate the necessity to utilize a cumbersome system with a relatively high power and high-cost light-source, and the necessity to utilize a spectrometer in order to determine an optimal operational frequency in which the effect of DOC on the measurement is minimal;   c. The system applies algorithms to calculate the nitrate concentration;   d. The system applies one or more databases for calibrating the measurement. Optionally, several databases are used to accommodate for different types of soils (therefore different compositions of materials and optical properties within the DOC);   e. A second embodiment of the system of the invention applies a mathematical model of the DOC concentration as a function of fluorescence power and the light absorption as a function of nitrate+DOC concentration and solves an equation of two variables;   f. A third embodiment of the system of the invention applies many measurements to train a deep learning network, which is later used to calculate the nitrate concentration based on two measurements of absorption and fluorescence measurements.       

       FIG. 1  generally illustrates in a flow diagram form a method  100  for determining a concentration of nitrate in a water sample (for example, which is taken from cultivated soil, waste water treatment plant effluent, hydroponic/aquaponic agriculture systems, rivers, lakes, and other agriculture-environmental applications. Hereinafter, for a sake of simplicity, the description will refer to a sample which is taken from soil) according to an embodiment of the invention. In step  112 , a water sample is obtained from the soil, for example, by use of a porous interface. In one embodiment, the sample flows in a low flow-rate through an optical flow cell. In step  114 , the water sample is illuminated, and the combined effect of the nitrate+DOC (referred to in the drawings as “N+DOC”) on the absorbance is determined. In step  116 , and following application of an excitation illumination on the sample in the optical flow cell, the concentration of DOC-only in the sample is determined by measuring the intensity of fluorescence emission from the sample. Finally, in step  118 , the nitrate concentration in the sample is determined, by subtracting the effect of the DOC (as obtained in step  116 ) on the absorbance, from the combined absorbance at N+DOC concentration as obtained in step  114 . The term “subtracting”, which appears throughout this application does not necessarily refers to a pure mathematical operation, but rather to an operation which offsets the effect of DOC from the combination of N+DOC. 
       FIG. 2  schematically illustrates a structure of system  200  for determining a concentration of nitrate in soil porewater, according to an embodiment of the invention. A small dead-volume porous interface  202  is placed in the soil to obtain a continuous low flux stream of soil porewater solution. The soil porewater flows through an optical flow cell  206  via small diameter tubing (e.g., inner diameter mm)  204 . The sample extraction from the soil is driven, for example, by applying low pressure (vacuum) on the porous interface. The sample can be later discharged or accumulated for further analysis or system calibration at sample accumulation chamber  232 . The casing of flow cell  206  is at least partially transparent, to allow passage of light beams therethrough. A first light source  212 , preferably of LED type, illuminates the cell in a first UV wavelength, for example in a proximity of 300 nm. The light beam of the first light source  212  passes through the optical flow cell  206  in which the soil porewater flows, while some of the light-beam&#39;s energy is absorbed by the water constituents (which in turn contains nitrate and DOC). The remaining energy from the light beam of the first light source  212  is accumulated by photodetector  222 , forming an absorbance signal  252 . A processing unit  240  is used to calculate the absorbance signal which reflects the difference between the illumination intensity by the first light source  212 , and the light intensity which is accumulated by photodetector  222  (Beer Lambert equation). A second light source  214 , preferably of a LED type, illuminates the cell  206  in a second wavelength, for example, in an excitation wavelength in proximity of 350 nm. The light beam from the second light source  214  excites the DOC within the water sample, causing a DOC fluorescence emission at a secondary wavelength of, for example, 451 nm. The fluorescence emission results substantially only from the DOC, is proportional to the DOC concentration, and is independent from the nitrate concentration within the sample. The fluorescence emission resulting from the second light beam is accumulated by second photodetector  224 , forming a fluorescence signal  254 . The fluorescence signal  254  is conveyed to the processing unit  240  where a predetermined calibration equation is used to estimate the DOC concentration in the sample. As previously mentioned, the fluorescence emission intensity is in proportion to the DOC concentration in the solution. However, the intensity of the fluorescence emission, in itself, is also proportional to the intensity of the excitation illumination by the second light source  214 . Yet, the intensity of the excitation illumination is a parameter which is controlled by the operator of the system. Preferably, filters  242  and  244  are located in front of optical detectors  222  and  224 , respectively, to ensure passage of only the wavelengths of interest towards the respective detectors. For example, a filter  244  allowing only light at 451 nm is located in front of the fluorescence detector  224 , ensuring that the 350 nm light from the excitation beam will not be received and saturate the detector  224  or mask the fluorescence reading. In Addition, a filter for allowing only the passage of, for example, 300 nm is located in front of the absorbance detector  222 , to ensure that the measurement will only be carried at the region where nitrate has a maximum absorbance peak. 
       FIG. 3  illustrates a method  300  for determining nitrate concentration in soil (or in another medium containing nitrate, as discussed above), according to an embodiment of the invention. In step  312 , the soil porewater is illuminated in a 1 st  wavelength by an illumination beam, typically either in the 200-240 nm band, or preferably in the 280-320 nm where a LED type illuminator can operate. In step  322 , the total absorbance which is the outcome of the presence of both N+DOC in the solution is calculated in reference to a predetermined blank solution, for example, double distilled water (DDW). In step  340  the calculated absorbance is substituted within a predetermined polynomial equation. This equation is further substituted with values from a 1 st  look-up-table  332  that provides a specific N/DOC concentration level for each specific combined absorption of N+DOC in the solution. Look-up-table  332  is mainly obtained through the use of experiments as will be later elaborated under the example section. Substantially simultaneously, the water sample is illuminated  314  by light in a 2nd wavelength, for example, 350 nm or any other wavelength which is suitable to excite a fluorescence illumination of DOC from the sample. In step  324 , a photodetector, which is tuned to detect in the wavelength of the fluorescence emission, for example 451 nm, detects the intensity of the fluorescence emission from the sample. Fluorescence from excitation/emission at 350/451 nm results only from the existence of DOC in the sample. Moreover, it has been found that, for a given type of soil, the intensity of the fluorescence emission can be found in a linear or polynomial correlation with the concentration rate of DOC in the sample. In step  344 , and given the intensity of the fluorescence emission as determined in step  324 , the concentration of DOC in the sample is determined. This is also done by use of a calibration 2nd look-up-table  364 . The calibration 2nd look-up-table  364 , upon submission of a specific emission intensity, provides a specific concentration rate of DOC. The calibration look-up-table  364 , is also prepared beforehand, typically by use of experiments. Preferably, look-up-table  364  selectively refers to the specific soil-type which is actually used. It has been found that the chemical composition of the DOC somewhat changes from one type of soil to another. Therefore, it is preferable to prepare a plurality of calibration look-up-tables  364 , each for a specific type of soil, and to actually use the look-up-table (or data therefrom) which most suits the type of the examined soil. The DOC concentration as determined in step  344  is then substituted within the predetermined polynomial calibration equation 340. Finally, the polynomial equation uses the given DOC concentration to deduct its attribution to the combined absorption of N+DOC, as given by step  322 . In other words, the polynomial equation returns the N concentration as a function of the DOC concentration (or its effect) in the solution, and the overall absorbance caused by N+DOC. 
     Soils from five different agricultural fields of the coastal plain of Israel, were collected and analyzed for this study, including: organic and conventional greenhouses for vegetable crops, open field which is used for rotating mixed crop, and a citrus&#39;s orchard. These sites were chosen to represent a spectrum of typical agricultural practices on different soils. In addition, commercial hummus soil mixture from “Dovrat” commercial compost was also examined to represent a potential impact of compost application in agriculture on the soil water DOC. 
     Soil solution samples were obtained from a soil and DDW mixture. The mixtures were left to stand for 24 hours to allow the solution to achieve a chemical equilibrium with the soil natural DOC. The soil phase and liquid phase in each sample was separated by a standard laboratory centrifuge, and the suspended solids were removed by 0.22 μm membrane filters. The samples were then diluted to obtain a series of replicas reflecting different DOC concentrations. Each replica of DOC concentration was spiked with a specific volume of 10,000 ppm standard potassium nitrate solution, to obtain between 4 to 6 different nitrate concentrations per each level of DOC. As a result, a matrix composed of 25 to 30 samples of different combinations of DOC and nitrate concentrations, ranging from zero to about a 1000 ppm nitrate, and between zero to about 100 ppm DOC was created from the soil extracts at each agricultural site. 
     The initial values of DOC and Total Nitrogen (TN) in each sample were estimated by an Analytic Jena multi N/C 2100s TOC/TN analyzer, while the nitrate concentration in each sample was determined by Dionex ICS 5000 Ion chromatograph. The absorption of each of the samples at 300 nm was determined using TECAN Spark 10M multimode microplate reader spectrophotometer. The light absorbance was defined by the Lambert-Beer equation: 
     
       
         
           
             
               
                 
                   Absorbance 
                   = 
                   
                     
                       - 
                       
                         log 
                         10 
                       
                     
                     ⁢ 
                     
                       I 
                       
                         I 
                         0 
                       
                     
                   
                 
               
               
                 
                   ( 
                   1 
                   ) 
                 
               
             
           
         
       
     
     where I reflects the light intensity after passing through the examined solution and Io is the light intensity of the source, or the light intensity of the source after passing through DDW as a reference. 
     A fluorescence spectroscopy technique was applied to measure the DOC concentration in the examined solution. The fact that the DOC fluorescence spectroscopy is not affected by the presence of nitrate in the solution makes it easy to separate the unique effect of absorption by nitrate and DOC. Fluorescence measurements were performed by TECAN Infinite M200 spectrophotometer with excitation (EX)/Emission (EM) at 350/451 nm. The results of the chemical and spectral analyses were used to obtain a series of matrix databases containing nitrate and DOC concentrations and absorbance at 300 nm (one per sampling site). In theory, the application of a UV absorption technique on aqueous nitrate solutions should be resulted in a clear linear correlation between the absorption rate and nitrate concentration.  FIG. 4  exemplifies said complex pattern m the relationship between nitrate concentration, DOC concentration and the UV absorption for soil water obtained from the conventional greenhouse field station. In fact, since both nitrate and DOC absorb light in different intensities at 300 nm, their contribution to the overall absorbance at 300 nm can be described as their cumulative absorbance superposition. As such, a 2D model (as presented in  FIG. 4 ) is not so adequate to characterize the relation between N/DOC concentration and the UV absorption.  FIG. 4  shows nitrate concentration vs. absorption at 300 nm of matrix of water samples obtained from the conventional greenhouse. Each data point represents soil water sample with a known concentration of DOC (indicated by notations next to each point (ppm)) and with known nitrate concentration. It should be noted that the X-axis was converted to a log scale merely to improve the visualization. However, when projecting the data on a 3D domain, it can be seen that the relationship between the three variables creates a plan which can be quantified mathematically ( FIG. 5 ). The applying of a multivariate regression model (e.g., polynomial equation in MATLAB) has shown how the nitrate concentration in a sample can be estimated as a function of the DOC and absorbance at 300 nm. The following is a general polynomial equation for nitrate estimation: 
     
       
         
           
             
               
                 
                   
                     Nitrate 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     
                       ( 
                       
                         DOC 
                         , 
                         Abs 
                       
                       ) 
                     
                   
                   = 
                   
                     
                       P 
                       00 
                     
                     + 
                     
                       
                         P 
                         01 
                       
                       ⁢ 
                       xAbs 
                     
                     + 
                     
                       
                         P 
                         10 
                       
                       ⁢ 
                       xDOC 
                     
                   
                 
               
               
                 
                   ( 
                   2 
                   ) 
                 
               
             
           
         
       
     
     Where: Nitrate is the nitrate concentration (ppm), DOC indicates the DOC concentration (ppm), Abs indicates the absorbance as measured at 300 nm (arbitrary units), and P 00 , P 10 , and P 01  are the coefficients as obtained from the regressions. 
     Fluorescence emission at 451 nm from samples with known DOC concentrations were used to develop a calibration curve for each site, as shown in  FIG. 6 . An analysis for the relation between DOC concentration and fluorescence intensity in the various water samples from different fields, respectively, has shown that even though some curves follow a same trend, thus having identical calibration curves (citrus orchard and open field), each group of samples has its own unique calibration equation due to the variation in the DOC chemical and optical characteristics between the various sites. The composition of organic matter that compose the DOC is affected by the site-specific characteristics such as: type of agricultural crop, biological activity in the soil, the type of applied compost, and the type of material used to fertilize. Therefore, calibration curves for DOC concentrations have to be site-specific in order to be properly used in fluorescence spectroscopy. Moreover, since the chemical variability in DOC affects the absorbance spectrum as well, the previously described polynomial equation for estimating nitrate concentration should preferably be site-specific. Nonetheless, it was concluded that the impact of DOC on the absorption spectrum, as resulted from its chemical composition, remains relatively constant over time and is a site-specific feature. This shows that once an initial calibration equation is obtained, this equation can be used for that particular site repeatedly in the long term, making a real-time continuous measurement of nitrate concentration in soil feasible. 
     Measurements of DOC concentration, as achieved from fluorescence emission at 451 nm, along with the total absorbance at 300 nm, enabled the estimation of nitrate concentration in a series of solutions obtained from the previously mentioned field sites ( FIG. 7 ). The concept model successfully estimated nitrate concentrations in all cases, in spite of the presence of DOC in the water samples. The high correlation between the observed and the predicted nitrate measurements (general R 2  &gt;0.96 and average RMSE=66.4 ppm nitrate) shows that the technique is applicable for a continuous, real-time measurement of nitrate concentration in soil as well as in other water sources containing DOC. The method satisfies practical demands for real-time continuous monitoring of nitrate concentration in soil, where the nitrate concentration range may vary from tens to thousands ppm during fertilization and irrigation cycles. Accordingly, LED-based components operating at 280-320 nm, and corresponding sensors can be used for the online, in-situ nitrate monitoring system in agricultural soils. Nevertheless, the technique of the invention can gain even higher resolution and higher accuracy in the future, when LED UV light sources ( 212  in  FIG. 2 ) at a wavelength of 200-250 nm will be developed and will be commercially available. 
     In the embodiment above, a calibration with respect to the concentration of nitrate uses a prior knowledge of the type of DOC at the soil where the specific nitrate-concentration determination is actually performed (i.e., what is the additional effect of the specific mixture of the DOC at that area on the absorption). Based on this prior knowledge, a calibration is performed. 
     In still another embodiment, the identification of the type of DOC at the area may be obtained automatically. In one alternative, the system measures the nitrate concentration in various DOC types by use of a two-phase procedure, as follows:
         a. A first phase in which three optical measurements of the solution are performed. The optical measurement may include either two absorption measurements and one fluorescence measurement, or one absorption measurement and two fluorescence measurement. Other combinations of absorbance and fluorescence measurements at many wavelengths may be used to get more accurate determination of the DOC type. Each optical measurement is performed at a different wave length;   b. A second phase in which, based on the results of the above first-phase measurements and prior-knowledge from data-base look-up tables that were prepared for different DOCs in advance, the system identifies which specific calibration equation or calibration curve to use.       

     In an embodiment of the invention, the analysis unit comprises a mathematical model to extract the value of nitrate based on the measurements of absorption, based on the Beer Lambert equation: 
     
       
         
           
             
               A 
               ⁡ 
               
                 ( 
                 λ 
                 ) 
               
             
             = 
             
               
                 ɛ 
                 ⁡ 
                 
                   ( 
                   λ 
                   ) 
                 
               
               ⁢ 
               x 
               ⁢ 
               C 
               ⁢ 
               x 
               ⁢ 
               l 
             
           
         
       
     
     Where A is the measured absorbance at a given wavelength (λ(nm)), ε(λ) is the molar attenuation coefficient (L mol −1  cm −1 ) at a given wavelength (λ(nm)), C is the examined Ion concentration (mol L −1 ) and l is an optical pathlength (cm). 
     Nitrate concentration can be then obtained by a set of two equations with two unknowns as following: 
     
       
         
           
             
               
                 
                   
                     A 
                     ⁢ 
                     1 
                     ⁢ 
                     
                       ( 
                       
                         λ 
                         1 
                       
                       ) 
                     
                   
                   = 
                   
                     
                       
                         
                           A 
                           
                             N 
                             ⁢ 
                             
                               O 
                               3 
                               
                                 - 
                                 1 
                               
                             
                           
                         
                         ⁡ 
                         
                           ( 
                           
                             λ 
                             1 
                           
                           ) 
                         
                       
                       + 
                       
                         
                           A 
                           
                             D 
                             ⁢ 
                             O 
                             ⁢ 
                             C 
                           
                         
                         ⁡ 
                         
                           ( 
                           
                             λ 
                             1 
                           
                           ) 
                         
                       
                     
                     = 
                     
                       
                         
                           
                             ɛ 
                             
                               N 
                               ⁢ 
                               
                                 O 
                                 3 
                                 
                                   - 
                                   1 
                                 
                               
                             
                           
                           ⁡ 
                           
                             ( 
                             
                               λ 
                               1 
                             
                             ) 
                           
                         
                         ⁢ 
                         
                           x 
                           ⁡ 
                           
                             ( 
                             
                               C 
                               
                                 NO 
                                 3 
                                 
                                   - 
                                   1 
                                 
                               
                             
                             ) 
                           
                         
                         ⁢ 
                         xl 
                       
                       + 
                       
                         
                           
                             ɛ 
                             
                               D 
                               ⁢ 
                               O 
                               ⁢ 
                               C 
                             
                           
                           ⁡ 
                           
                             ( 
                             
                               λ 
                               1 
                             
                             ) 
                           
                         
                         ⁢ 
                         
                           x 
                           ⁡ 
                           
                             ( 
                             
                               C 
                               
                                 D 
                                 ⁢ 
                                 O 
                                 ⁢ 
                                 C 
                               
                             
                             ) 
                           
                         
                         ⁢ 
                         xl 
                       
                     
                   
                 
               
               
                 
                   ( 
                   1 
                   ) 
                 
               
             
             
               
                 
                   
                     A2 
                     ⁡ 
                     
                       ( 
                       
                         λ 
                         2 
                       
                       ) 
                     
                   
                   = 
                   
                     
                       
                         
                           A 
                           
                             N 
                             ⁢ 
                             
                               O 
                               3 
                               
                                 - 
                                 1 
                               
                             
                           
                         
                         ⁡ 
                         
                           ( 
                           
                             λ 
                             2 
                           
                           ) 
                         
                       
                       + 
                       
                         
                           A 
                           
                             D 
                             ⁢ 
                             O 
                             ⁢ 
                             C 
                           
                         
                         ⁡ 
                         
                           ( 
                           
                             λ 
                             2 
                           
                           ) 
                         
                       
                     
                     = 
                     
                       
                         
                           
                             ɛ 
                             
                               N 
                               ⁢ 
                               
                                 O 
                                 3 
                                 
                                   - 
                                   1 
                                 
                               
                             
                           
                           ⁡ 
                           
                             ( 
                             
                               λ 
                               2 
                             
                             ) 
                           
                         
                         ⁢ 
                         
                           x 
                           ⁡ 
                           
                             ( 
                             
                               C 
                               
                                 N 
                                 ⁢ 
                                 
                                   O 
                                   3 
                                   
                                     - 
                                     1 
                                   
                                 
                               
                             
                             ) 
                           
                         
                         ⁢ 
                         x 
                         ⁢ 
                         l 
                       
                       + 
                       
                         
                           
                             ɛ 
                             
                               D 
                               ⁢ 
                               O 
                               ⁢ 
                               C 
                             
                           
                           ⁡ 
                           
                             ( 
                             
                               λ 
                               2 
                             
                             ) 
                           
                         
                         ⁢ 
                         
                           x 
                           ⁡ 
                           
                             ( 
                             
                               C 
                               
                                 D 
                                 ⁢ 
                                 O 
                                 ⁢ 
                                 C 
                               
                             
                             ) 
                           
                         
                         ⁢ 
                         xl 
                       
                     
                   
                 
               
               
                 
                   ( 
                   2 
                   ) 
                 
               
             
           
         
       
     
     From consolidating equation (1)+(2) we can extract the nitrate concentration (equation (3)): 
     
       
         
           
             
               
                 
                   
                     C 
                     
                       NO 
                       3 
                       
                         - 
                         1 
                       
                     
                   
                   = 
                   
                     
                       
                         
                           
                             C 
                             
                               D 
                               ⁢ 
                               O 
                               ⁢ 
                               C 
                             
                           
                           ⁡ 
                           
                             ( 
                             
                               
                                 
                                   
                                     A 
                                     
                                       N 
                                       ⁢ 
                                       
                                         O 
                                         3 
                                         
                                           - 
                                           1 
                                         
                                       
                                     
                                   
                                   ⁡ 
                                   
                                     ( 
                                     
                                       λ 
                                       1 
                                     
                                     ) 
                                   
                                 
                                 + 
                                 
                                   
                                     A 
                                     
                                       D 
                                       ⁢ 
                                       O 
                                       ⁢ 
                                       C 
                                     
                                   
                                   ⁡ 
                                   
                                     ( 
                                     
                                       λ 
                                       1 
                                     
                                     ) 
                                   
                                 
                               
                               
                                 
                                   
                                     A 
                                     
                                       N 
                                       ⁢ 
                                       
                                         O 
                                         3 
                                         
                                           - 
                                           1 
                                         
                                       
                                     
                                   
                                   ⁡ 
                                   
                                     ( 
                                     
                                       λ 
                                       2 
                                     
                                     ) 
                                   
                                 
                                 + 
                                 
                                   
                                     A 
                                     
                                       D 
                                       ⁢ 
                                       O 
                                       ⁢ 
                                       C 
                                     
                                   
                                   ⁡ 
                                   
                                     ( 
                                     
                                       λ 
                                       2 
                                     
                                     ) 
                                   
                                 
                               
                             
                             ) 
                           
                         
                         ⁢ 
                         x 
                         ⁢ 
                         
                           
                             ɛ 
                             
                               D 
                               ⁢ 
                               O 
                               ⁢ 
                               C 
                             
                           
                           ⁡ 
                           
                             ( 
                             
                               λ 
                               2 
                             
                             ) 
                           
                         
                       
                       - 
                       
                         
                           ɛ 
                           
                             D 
                             ⁢ 
                             O 
                             ⁢ 
                             C 
                           
                         
                         ⁡ 
                         
                           ( 
                           
                             λ 
                             1 
                           
                           ) 
                         
                       
                     
                     
                       
                         
                           ɛ 
                           
                             NO 
                             3 
                             
                               - 
                               1 
                             
                           
                         
                         ⁡ 
                         
                           ( 
                           
                             λ 
                             1 
                           
                           ) 
                         
                       
                       - 
                       
                         
                           ( 
                           
                             
                               
                                 
                                   A 
                                   
                                     N 
                                     ⁢ 
                                     
                                       O 
                                       3 
                                       
                                         - 
                                         1 
                                       
                                     
                                   
                                 
                                 ⁡ 
                                 
                                   ( 
                                   
                                     λ 
                                     1 
                                   
                                   ) 
                                 
                               
                               + 
                               
                                 
                                   A 
                                   
                                     D 
                                     ⁢ 
                                     O 
                                     ⁢ 
                                     C 
                                   
                                 
                                 ⁡ 
                                 
                                   ( 
                                   
                                     λ 
                                     1 
                                   
                                   ) 
                                 
                               
                             
                             
                               
                                 
                                   A 
                                   
                                     N 
                                     ⁢ 
                                     
                                       O 
                                       3 
                                       
                                         - 
                                         1 
                                       
                                     
                                   
                                 
                                 ⁡ 
                                 
                                   ( 
                                   
                                     λ 
                                     2 
                                   
                                   ) 
                                 
                               
                               + 
                               
                                 
                                   A 
                                   
                                     D 
                                     ⁢ 
                                     O 
                                     ⁢ 
                                     C 
                                   
                                 
                                 ⁡ 
                                 
                                   ( 
                                   
                                     λ 
                                     2 
                                   
                                   ) 
                                 
                               
                             
                           
                           ) 
                         
                         ⁢ 
                         x 
                         ⁢ 
                         
                           
                             ɛ 
                             
                               NO 
                               3 
                               
                                 - 
                                 1 
                               
                             
                           
                           ⁡ 
                           
                             ( 
                             
                               λ 
                               2 
                             
                             ) 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   3 
                   ) 
                 
               
             
           
         
       
     
     Where A is the measured absorbance at a given wavelength (λ 1,2  (nm)), εNO 3   −1 , Doc (λ 1,2  (nm)) is a molar attenuation coefficient for either the nitrate or the DOC (L mol −1  cm −1 ) at a given wavelength (λ 1,2  (nm)), C DOC  is the DOC concentration (mol L −1 ) as obtained from said second photodetector, C NO3  is a nitrate concentration (mol L −1 ), and/is an optical pathlength (cm). 
     In still another embodiment, a machine learning technique may be used. The machine learning may include the following phases:
         a. Generation of many absorption and fluorescence measurements for different values of nitrate concentration and various DOC types and respective concentrations;   b. Selection and adaption of one or more deep learning networks from known results;   c. Training of the network and selection of a one promising a best performance;   d. Use of the trained network to calculate the nitrate concentration based on absorption and fluorescence measurements in a manner as described above.       

       FIG. 8  shows the observed vs. predicted nitrate concentration, for hummus soil mixture water extract, as obtained by the technique of the invention. In this specific experiment the predicted data was obtained for absorbance measurements at 235 nm. 
     While some embodiments of the invention have been described by way of illustration, it will be apparent that the invention can be carried into practice with many modifications, variations and adaptations, and with the use of numerous equivalents or alternative solutions that are within the scope of persons skilled in the art, without departing from the spirit of the invention or exceeding the scope of the claims.