Low-Cost Nitric Oxide Sensors: Assessment of Temperature and Humidity Effects

High equipment cost is a significant entry barrier to research for small organizations in developing solutions to air pollution problems. Low-cost electrochemical sensors show sensitivity at parts-per-billion by volume mixing ratios but are subject to variation due to changing environmental conditions, in particular temperature. In this study, we demonstrate a low-cost Internet of Things (IoT)-based sensor system for nitric oxide analysis. The sensor system used a four-electrode electrochemical sensor exposed to a series of isothermal/isohume conditions. When deployed under these conditions, stable baseline responses were achieved, in contrast to ambient air conditions where temperature and humidity conditions may be variable. The interrelationship between working and auxiliary electrodes was linear within an environmental envelope of 20–40 °C and 30–80% relative humidity, with correlation coefficients from 0.9980 to 0.9999 when measured under isothermal/isohume conditions. These data enabled the determination of surface functions that describe the working to auxiliary electrode offsets and calibration curve gradients and intercepts. The linear and reproducible nature of individual calibration curves for stepwise nitric oxide (NO) additions under isothermal/isohume environments suggests the suitability of these sensors for applications aside from their role in air quality monitoring. Such applications would include nitric oxide kinetic studies for atmospheric applications or measurement of the potential biocatalytic activity of nitric oxide consuming enzymes in biocatalytic coatings, both of which currently employ high-capital-cost chemiluminescence detectors.


Introduction
Whilst it is now accepted that increased atmospheric pollutant concentrations are driving climate change, they also have a range of pronounced impacts on the environment and human health. Approximately 2.4 billion people reside in dwellings that use biomass fuel as the primary source of fuel for heating, cooking, or both. Estimates suggest that 3.2 million people die prematurely each year from illness directly attributable to the indoor air pollution generated from these biomass fuels [1]. A Taiwanese study in 2020 showed a high correlation between outpatient visits and increased concentrations of NO 2 , O 3 , particulate matter with a diameter ≤ 10 µm (PM 10 ), CO 2 , and NO [2].
Oxides of nitrogen originating from combustion processes in internal combustion engines are found in the photochemically generated smog that envelopes many modern cities, along with O 3 . Therefore, monitoring these compounds is important to inform appropriate health advice and assess the effectiveness of management strategies to reduce air pollution.
Many techniques are available for measuring NO concentrations, ranging from simple diffusive samplers [3][4][5][6] to absorption and emission spectroscopy [7][8][9], gas chromatography [10][11][12], and mass spectroscopy [13,14]. Chemiluminescence, an example of an emission spectroscopy technique, is considered the standard for the measurement of nitrogen oxides for air quality monitoring [6,[15][16][17][18][19]. Whilst capable of detecting nitrogen oxides at the parts per billion level (ppbV), drawbacks to the use of chemiluminescence detectors include the initial capital equipment cost (i.e., in the order of tens of thousands of US dollars per unit), ongoing maintenance and calibration costs, as well as the requirement for an ozone source for the oxidation reaction of NO to NO 2 .
The high capital cost of equipment is a significant entry barrier to research for small organizations in developing solutions to air pollution problems. This entry barrier and the steadily increasing concerns regarding the impact of air quality on human health have led to several commercial enterprises developing low-cost monitoring systems. Whilst these low-cost systems provide viable alternatives to high-cost government facilities for individuals to track their exposure, the accuracy and detection limits of these systems may not be suitable for all monitoring purposes [5].
Low-cost electrochemical sensors (LCESs) are electrochemical cells that, when exposed to the target gas, generate a current that is linearly proportional to the concentration of that gas. LCESs have been used for industrial safety applications since at least the 1970s [20]. LCES's used for industrial safety applications typically measure target gases at the parts per million (ppm) level, whilst air quality monitoring is typically measured in the parts per billion (ppb) range, some 2-3 orders of magnitude lower than the industrial applications. Following improvements in the signal-noise ratio, sensitivity, and conditioning circuitry, Mead et al. demonstrated that LCESs can provide accurate measurements in the parts per billion (ppb) range under certain conditions [21]. Meyerhoff et al. compared a commercially available LCES with a chemiluminescence detector to measure the nitrite concentration in food samples and found a linear regression curve close to unity (R = 0.999) [22].
The drawback with these improved sensors thus far has been the response, particularly of NO sensors, to variations in ambient environmental conditions [21,[23][24][25]. Significant effort has been made to develop methods to overcome these deficiencies, ranging from simple post-data corrections to machine learning algorithms [26][27][28].
Wei et al. [29] investigated the accuracy of LCESs under stable laboratory conditions and demonstrated the high potential of such sensors compared to conventional reference instruments. Wei et al. applied a linear correction methodology considering temperature and relative humidity, concluding that "careful data analysis and correction protocols are essential to guarantee good data quality" [29] (p. 73).
Wei et al. studied NO 2 , CO, NO and O x sensors, whereas Rogulski et al. [30] focused solely on Alphasense NO 2 sensors. finding that temperature affected the sensor accuracy more than relative humidity, but also that air temperatures > 30 • C could lead to absolute percentage errors > 150%. Thus, applying the correction methods recommended by Alphasense may not be appropriate under all environmental; conditions. The improved correction method proposed by Rogulski et al. [30] was based on a second-degree polynomial regression.
The main goal of our study was to assess the influences of temperature and humidity on an isolated low-cost NO electrochemical sensor and establish a polynomial surface function that combines the dependent variables of temperature and humidity for calibration of NO concentration. In contrast to previous studies, our focus was not on applying the derived functions to measure ambient air quality under changing environmental conditions. Rather, we focused on employing the sensors under isothermal/isohume conditions for applications in areas such as NO kinetic studies for atmospheric applications or measurement of potential biocatalytic activity of NO-consuming enzymes in biocatalytic coatings.

Materials and Methods
A Chemglass Life Sciences CG-1929-29 5-L double-walled reaction vessel fitted with a CG-1941 4-neck reaction vessel lid was used as the atmospheric chamber. The reaction vessel was equipped with a 3 4 -inch beaded inlet/outlet pipe for water circulation between the vessel walls. The four inlet necks allowed a cable to power the Raspberry Pi and for inlet of zero air and nitric oxide gas. When assembled, the vessel had a total volume of 7.87 L (Figure 1a). A custom ABS, 3D printed frame was constructed to support the LCES. From a cost perspective, the above system provides significant capital expenditure savings compared to chemiluminescence detectors (Table 2, Figure 1b).

Software Control System
User control of the system was achieved through a graphical user interface written in C#. To enable the compiled code to run on the Raspbian operating system, Mono, a software framework built as an open-source and free alternative to the. NET Framework, was installed on the Raspberry Pi.
The C# program calls a series of python scripts that handle the interrogation of the attached sensors. Open-source python libraries were imported into the scripts at runtime to provide access to system functionality such as the ASD1115 bus and GPIO pins on the Raspberry Pi4. Data were written to comma-separated value (csv) files which were subsequently imported into Microsoft Excel for analysis. The data files were initially written to the micro-SD card inserted into the Raspberry Pi before being uploaded to a cloud storage account. Data were visualized in real time using an Adafruit IO cloud-based dashboard.  Chamber temperature was controlled by a Thermoline TU3 water circulator/heaterequipped water reservoir. Instrument-grade air from BOC Gases (Table 1) was used as a source of zero air for purging the reaction chamber. Following the purge, humidity was first generated and then maintained using water/glycerol solutions according to ASTM D5036-2011. The use of water/glycerol solutions for generating and maintaining relative humidity conditions was preferred over gas bubblers due to the limitations in generating humidity's greater than approximately 50-55% with gas bubbler systems [29]. Using this approach has demonstrated temperature and relative humidity stability of ±0.1 • C and ±0.5% relative humidity, respectively.
Once the desired temperature and humidity conditions were established, 60 ppm NO gas in nitrogen was injected into the chamber using a Swagelok VAF-G2-01M-1-0 Variable Area Flowmeter at a flow rate of approximately 1.0 L/h. NO gas flow was continued until the desired concentration was achieved. The sensor system was constructed using an Alphasense Limited (Alphasense) NO NO-A4 mounted on an Alphasense analog front end (AFE). The AFE measured the electrode voltages from the working and auxiliary electrodes of the electrochemical sensor and provided an analog output of the values. The AFE was then connected to a South Coast Science (SCS) Alpha Pi Eng. Board (APEB) which conditioned the analog output of the Alphasense AFE and provided digitization of the electrode voltages through the use of ADS1115 analog to digital converters (ADCs). The voltage range of these converters is in the order of 2 volts (V) with 16-bit resolution, thereby providing a resolution of 3 × 10 −5 V.
The APEB was equipped with a 40-pin header compatible with microcomputers such as the Raspberry Pi 3 and 4. A Raspberry Pi 4, 4GB was attached to the above system. Raspbian "Buster", a Linux-based distribution, was used as the operating system, allowing a Microsoft Windows appearance to the microcomputer.
A Bosch BME280 sensor was attached to the 40-pin header to measure temperature and relative humidity. The BME280 is an integrated environmental sensor that has been specifically developed for applications such as home automation control, personalized weather stations, and sports fitness tools. The temperature sensor incorporated into the BME280 was calibrated against a National Institute of Standards and Technology (NIST)-calibrated thermometer. Calibration of relative humidity was undertaken using ASTM D5032, wherein a water/glycerol solution was prepared with the concentration required to achieve the desired relative humidity. This solution was added to Petri dishes and placed inside the chamber along with the LCES system. For a typical calibration, three Petri dishes each containing 50 ± 1 mL of solution were used. Following sealing, the chamber was purged with zero air at a 5 L/min flow rate for approximately 5 min. Relative humidity was then allowed to equilibrate and the final refractive index of the solution was used to determine the stabilized relative humidity value due to water evaporation or absorption into or from the airspace above the solution during the stabilization period.
A custom ABS, 3D printed frame was constructed to support the LCES. From a cost perspective, the above system provides significant capital expenditure savings compared to chemiluminescence detectors (Table 2, Figure 1b).

Software Control System
User control of the system was achieved through a graphical user interface written in C#. To enable the compiled code to run on the Raspbian operating system, Mono, a software framework built as an open-source and free alternative to the. NET Framework, was installed on the Raspberry Pi.
The C# program calls a series of python scripts that handle the interrogation of the attached sensors. Open-source python libraries were imported into the scripts at runtime to provide access to system functionality such as the ASD1115 bus and GPIO pins on the Raspberry Pi4. Data were written to comma-separated value (csv) files which were subsequently imported into Microsoft Excel for analysis. The data files were initially written to the micro-SD card inserted into the Raspberry Pi before being uploaded to a cloud storage account. Data were visualized in real time using an Adafruit IO cloud-based dashboard.

Control of NO Concentrations
NO at a mixing ratio of 60 ppm ±2% in nitrogen was obtained from CAC Gas and Instrumentation. A 0.5 L/min fixed flow rate regulator was attached to the NO cylinder. Gas was fed from this regulator to a calibrated Swagelok VAF-G2-01M-1-0 Variable Area Flowmeter which was set to deliver a 1.0 L/h flow rate. Gas exited the flowmeter and entered a 2-position, 3-way electrically controlled solenoid valve (refer Supplementary Materials) controlled from the microcomputer. The response time of the solenoid valve was approximately 30 ms.
At predetermined intervals, the solenoid valve was energized, and gas flow was diverted into the reaction chamber. Typically, at a 1 L/h flow rate of 60 ppmV ±2% NO gas, the solenoid valve required energizing for 9.444 s to deliver sufficient gas to increase the NO concentration in the chamber by 20 ppbV (refer Supplementary Materials for a typical electrode response to 20 ppbV stepwise nitric oxide addition and resultant calibration curve).
According to Alphasense, the NO electrochemical sensor has a response time (t90) of <25 s. To ensure that the sensor had stabilized following gas addition, a post stabilization time of 5 min was provided before the next gas addition.
The 30 ms response time of the solenoid valve, therefore, resulted in an error of ±0.064 ppbV under the conditions described above. This error is approximately 3.18% of the error associated with the stock gas concentration. In experiments with stepwise additions of NO gas, the total combined error associated with this addition method was ±0.464 ppbV per 20 ppbV step change. When a typical 100 ppbV, single-step gas addition was used, the total error was estimated to be ±2.064 ppbV, compared to ±2.320 ppbV for a multi-step addition.

Electrode Output Stability
The output of LCESs, particularly NO sensors, are susceptible to variations in ambient environmental conditions [28]. When the environment can be precisely controlled, such as under isothermal/isohume conditions in a sealed chamber, stable outputs could be obtained from the LCESs (Figure 2). Bosch BME280 sensor 7.00 Total 767.50

Control of NO Concentrations
NO at a mixing ratio of 60 ppm ±2% in nitrogen was obtained from CAC Gas and Instrumentation. A 0.5 L/min fixed flow rate regulator was attached to the NO cylinder. Gas was fed from this regulator to a calibrated Swagelok VAF-G2-01M-1-0 Variable Area Flowmeter which was set to deliver a 1.0 L/h flow rate. Gas exited the flowmeter and entered a 2-position, 3-way electrically controlled solenoid valve (refer Supplementary Materials) controlled from the microcomputer. The response time of the solenoid valve was approximately 30 ms.
At predetermined intervals, the solenoid valve was energized, and gas flow was diverted into the reaction chamber. Typically, at a 1 L/h flow rate of 60 ppmV ±2% NO gas, the solenoid valve required energizing for 9.444 s to deliver sufficient gas to increase the NO concentration in the chamber by 20 ppbV (refer Supplementary Materials for a typical electrode response to 20 ppbV stepwise nitric oxide addition and resultant calibration curve) According to Alphasense, the NO electrochemical sensor has a response time (t90) of <25 s. To ensure that the sensor had stabilized following gas addition, a post stabilization time of 5 min was provided before the next gas addition.
The 30 ms response time of the solenoid valve, therefore, resulted in an error of ±0.064 ppbV under the conditions described above. This error is approximately 3.18% of the error associated with the stock gas concentration. In experiments with stepwise additions of NO gas, the total combined error associated with this addition method was ±0.464 ppbV per 20 ppbV step change. When a typical 100 ppbV, single-step gas addition was used, the total error was estimated to be ±2.064 ppbV, compared to ±2.320 ppbV for a multi-step addition.

Electrode Output Stability
The output of LCESs, particularly NO sensors, are susceptible to variations in ambient environmental conditions [28]. When the environment can be precisely controlled, such as under isothermal/isohume conditions in a sealed chamber, stable outputs could be obtained from the LCESs (Figure 2).  In the latest generation four-electrode LCESs, such as those used during this study, the auxiliary electrode is isolated from the environment containing the target gas. Therefore, relative humidity fluctuations only affect the working electrode, leaving the responses from the auxiliary electrode unaltered. Figure 3 demonstrates a scenario in which the chamber was purged with zero air and allowed to equilibrate before again purging with zero air at the 13 electrode, which was exposed to the chamber environment, whereas the auxiliary electrode remained unaffected.
In the latest generation four-electrode LCESs, such as those used during this study, the auxiliary electrode is isolated from the environment containing the target gas. Therefore, relative humidity fluctuations only affect the working electrode, leaving the responses from the auxiliary electrode unaltered. Figure 3 demonstrates a scenario in which the chamber was purged with zero air and allowed to equilibrate before again purging with zero air at the 13:40 timestamp, to simulate a humidity transient. The result confirmed a transient response on the working electrode, which was exposed to the chamber environment, whereas the auxiliary electrode remained unaffected. Where temperature fluctuations occurred, both the working and auxiliary electrodes were affected, albeit to slightly different extents (Figure 4). From these data, the temperature dependence for the auxiliary and working electrode responses could be determined. Both the auxiliary and working electrodes displayed high correlation coefficients with temperature ( Figure 5). This in turn allowed the interrelationship between the auxiliary and working electrode responses resulting from minor temperature fluctuations to be plotted.  Where temperature fluctuations occurred, both the working and auxiliary electrodes were affected, albeit to slightly different extents (Figure 4). From these data, the temperature dependence for the auxiliary and working electrode responses could be determined. Both the auxiliary and working electrodes displayed high correlation coefficients with temperature ( Figure 5). This in turn allowed the interrelationship between the auxiliary and working electrode responses resulting from minor temperature fluctuations to be plotted.
In the latest generation four-electrode LCESs, such as those used during this study, the auxiliary electrode is isolated from the environment containing the target gas. Therefore, relative humidity fluctuations only affect the working electrode, leaving the responses from the auxiliary electrode unaltered. Figure 3 demonstrates a scenario in which the chamber was purged with zero air and allowed to equilibrate before again purging with zero air at the 13:40 timestamp, to simulate a humidity transient. The result confirmed a transient response on the working electrode, which was exposed to the chamber environment, whereas the auxiliary electrode remained unaffected. Where temperature fluctuations occurred, both the working and auxiliary electrodes were affected, albeit to slightly different extents (Figure 4). From these data, the temperature dependence for the auxiliary and working electrode responses could be determined. Both the auxiliary and working electrodes displayed high correlation coefficients with temperature ( Figure 5). This in turn allowed the interrelationship between the auxiliary and working electrode responses resulting from minor temperature fluctuations to be plotted.   Figure 6 shows a highly linear interrelationship between the auxiliary and working electrodes over a minor temperature range variation. In this situation, chamber environmental conditions varied by approximately ±1.2 • C, and relative humidity was maintained at ±0.5%. After the temperature and auxiliary electrode responses were measured, the expected zero value for the working electrode response could be calculated.  Figure 6 shows a highly linear interrelationship between the auxiliary and working electrodes over a minor temperature range variation. In this situation, chamber environmental conditions varied by approximately ±1.2 °C, and relative humidity was maintained at ±0.5%. After the temperature and auxiliary electrode responses were measured, the expected zero value for the working electrode response could be calculated. Using this calculated zero air value for the working electrode response, the target gas concentration under the controlled conditions could be determined from the difference between the expected zero air baseline response at the current sensor temperature and the actual measured output, based on a calibration curve that had previously been determined under the same temperature and relativity humidity conditions. By maintaining isohume conditions and varying the chamber temperature, we determined that this linear relationship extended across a temperature range from 20 to 40 °C. The calculated gradient was greater than unity, indicating that the working electrode output response was increasing at a more rapid rate than that of the auxiliary electrode. When data from varying humidity's were combined, it was also noted that whilst the linear relationship between the auxiliary and working electrode zero air responses was maintained, the   Figure 6 shows a highly linear interrelationship between the auxiliary and working electrodes over a minor temperature range variation. In this situation, chamber environmental conditions varied by approximately ±1.2 °C, and relative humidity was maintained at ±0.5%. After the temperature and auxiliary electrode responses were measured, the expected zero value for the working electrode response could be calculated. Using this calculated zero air value for the working electrode response, the target gas concentration under the controlled conditions could be determined from the difference between the expected zero air baseline response at the current sensor temperature and the actual measured output, based on a calibration curve that had previously been determined under the same temperature and relativity humidity conditions. By maintaining isohume conditions and varying the chamber temperature, we determined that this linear relationship extended across a temperature range from 20 to 40 °C. The calculated gradient was greater than unity, indicating that the working electrode output response was increasing at a more rapid rate than that of the auxiliary electrode. When data from varying humidity's were combined, it was also noted that whilst the linear relationship between the auxiliary and working electrode zero air responses was maintained, the Using this calculated zero air value for the working electrode response, the target gas concentration under the controlled conditions could be determined from the difference between the expected zero air baseline response at the current sensor temperature and the actual measured output, based on a calibration curve that had previously been determined under the same temperature and relativity humidity conditions. By maintaining isohume conditions and varying the chamber temperature, we determined that this linear relationship extended across a temperature range from 20 to 40 • C. The calculated gradient was greater than unity, indicating that the working electrode output response was increasing at a more rapid rate than that of the auxiliary electrode. When data from varying humidity's were combined, it was also noted that whilst the linear relationship between the auxiliary and working electrode zero air responses was maintained, the gradient and intercept of each data series changed considerably, suggesting a polynomial rather than linear surface function (Table 3; Figure 7). gradient and intercept of each data series changed considerably, suggesting a polynomial rather than linear surface function (Table 3; Figure 7). Table 3. Trendline equations and correlation coefficients for auxiliary and working electrode relationships across varying humidity's. The surface function described by Figure 7 was obtained from the MATLAB Curve Fitting App using a third-degree polynomial fit (Poly33) in temperature and relative humidity. The function was of the following form :   ,  00  10  01  20  11  02  30  21  12 03

Relative Humidity (%) Trendline Equation
(1) where x is temperature, and y is relative humidity. The coefficients for the equation are reported in Table 4. Table 4. Coefficients for Figure 8 (with 95% confidence bounds). The surface function described by Figure 7 was obtained from the MATLAB Curve Fitting App using a third-degree polynomial fit (Poly33) in temperature and relative humidity. The function was of the following form:

Value 95% Confidence Bounds
f (x, y) = p00 + p10x + p01y + p20x 2 + p11xy + p02y 2 + p30x 3 + p21x 2 y + p12xy 2 + p03y 3 (1) where x is temperature, and y is relative humidity. The coefficients for the equation are reported in Table 4. Table 4. Coefficients for Figure 8 (with 95% confidence bounds). The third-degree polynomial fit described by Equation (1) yielded an R 2 value of 0.9697, a sum of squares due to error (SSE) value of 5.003 × 10 −5 , and a root mean square error (RMSE) of 0.001582. Higher degree polynomials did not notably improve the R 2 , SSE, or RMSE values. The residuals plot derived from the polynomial fit (Figure 9 showed higher deviations with the surface fit with increasing temperature and relative humidity conditions, supporting observations by Rogulski et al. [30].

Effect of Varying Environmental Conditions on NO Calibration Curves
Individual calibration curves using the method described previously were obtained across the environmental envelope of 20-40 • C and 30-80% relative humidity. Highly linear correlation coefficients were found for all the stepwise experiments (Table 5). Several sets of data were obtained when calibration curves were prepared under isothermal/isohume conditions across an environmental envelope. Firstly, the electrode responses for both auxiliary and working electrodes could be determined for any NO concentration. Third-degree polynomial surface fits similar to those shown in Figure 7 could also be constructed from these data. More importantly, a third-degree polynomial surface fit of the form described by Equation (1), for the offset between the working and auxiliary electrode, could be constructed for these data for each measured environmental condition (Figure 7). This surface fit also displayed a high correlation, with an R 2 value of 0.963, an SSE value of 6.213 × 10 −5 , and RMSE of 0.001762. The residual plot (not shown) demonstrated similar deviations from the surface fit with increasing temperature and relative humidity.
Further surface fits could likewise be constructed for the gradients and intercepts of the calibration curves described above (Figure 8-Surface Fits, Figure 9-Residuals plot).

Reproducibility and Transferability of Surface Functions
Having defined surface functions for the working to auxiliary electrode offset, and gradients and intercepts for calibration graphs, a series of experiments were conducted to determine the reproducibility with the original sensor and the transferability to new NOA4 sensors. Four environmental conditions with three NO concentrations were employed, and the results were calculated using both the initial calibration curve for the specific environmental condition and the gradient and intercept surface functions.
The original sensor produced results within ±3% of the expected value using the calibration curve method suggesting that the initially constructed calibration curves had remained stable over a period of approximately three months (Figure 10a,b). Sensors 2 and 3, however, showed variation from expected values ranging from −40% to +115% of the known concentrations (Figure 10c-f). These variations typically increased under higher temperature and relative humidity. Higher levels of deviation from expected values resulted from the use of the surface functions. It is unclear if this was an artifact of the surface functions themselves, the small sample set of sensors, or a combination of both. Further work increasing the sample size of sensors and/or individual environmental conditions may assist in determining the cause of the deviation from expected NO values when using the surface functions. The results highlight that while the calibration for individual sensors is stable over monthly timescales, the calibration parameters are not transferable R 2 values of 0.9321 and 0.9741 were obtained for the gradient and intercept fits, respectively; SSE and RMSE values of 3.017 × 10 −4 and 44.85, respectively, were obtained for the gradient fit, and of 122.48 and 2.4747, respectively, for the intercept fit.

Reproducibility and Transferability of Surface Functions
Having defined surface functions for the working to auxiliary electrode offset, and gradients and intercepts for calibration graphs, a series of experiments were conducted to determine the reproducibility with the original sensor and the transferability to new NOA4 sensors. Four environmental conditions with three NO concentrations were employed, and the results were calculated using both the initial calibration curve for the specific environmental condition and the gradient and intercept surface functions.
The original sensor produced results within ±3% of the expected value using the calibration curve method suggesting that the initially constructed calibration curves had remained stable over a period of approximately three months (Figure 10a,b). Sensors 2 and 3, however, showed variation from expected values ranging from −40% to +115% of the known concentrations (Figure 10c-f). These variations typically increased under higher temperature and relative humidity. Higher levels of deviation from expected values resulted from the use of the surface functions. It is unclear if this was an artifact of the surface functions themselves, the small sample set of sensors, or a combination of both. Further work increasing the sample size of sensors and/or individual environmental conditions may assist in determining the cause of the deviation from expected NO values when using the surface functions. The results highlight that while the calibration for individual sensors is stable over monthly timescales, the calibration parameters are not transferable across sensors, and therefore the calibration procedure needs to be carried out for each sensor.

Conclusions
Electrochemical sensors, particularly NO sensors, are prone to variations in response due to variations in ambient environmental conditions. We demonstrated that under well controlled environmental conditions, i.e., isothermal/isohume conditions within a sealed environmental (reaction) chamber, stable responses can be achieved by NO electrochemical sensors. Furthermore, linear (although environmental-condition-dependent) calibration curves were obtained. Polynomial surface functions constructed from the individual linear calibration curves displayed good reproducibility under certain environmental conditions but resulted in increasing deviations under higher temperature and relative humidity conditions. Furthermore, we demonstrated that neither the linear nor the polynomial surface function calibration parameters could be transferred from the original sensors to alternate units.
Whilst the results reported in this paper are not immediately applicable to applications such as air quality monitoring without post-data corrections or use of machine learning algorithms, the linearity and reproducibility under isothermal/isohume conditions suggested that the sensors are suitable for applications where the experimental conditions are well controlled. Such applications would include NO kinetic studies for atmospheric applications, or the measurement of potential biocatalytic activity of NO -consuming enzymes in biocatalytic coatings, both of which currently employ high-capital-cost chemiluminescence detectors. We will focus on these areas in subsequent studies.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. These data are not publicly available due to data collection format and authentication requirements of the current cloud storage location.