Real-time estimations of blood glucose concentrations from sweat measurements using the local density random walk model

TL;DR

This paper introduces a faster model for estimating blood glucose from sweat, enabling real-time monitoring without blood draws.

Contribution

The Local Density Random Walk model enables efficient real-time glucose estimation from sweat, replacing slower biophysical simulations.

Findings

The LDRW model achieves similar accuracy to the original biophysical model with a root mean square difference of 0.04 mmol/L.

The LDRW model reduces computational time to 2.6 seconds per data point, 0.7% of the original method's time.

The LDRW model has a better balance of fit and complexity, as indicated by a lower corrected Akaike Information Criterion value.

Abstract

Sweat provides a non-invasive alternative to blood draws, enabling glucose-concentration monitoring for both healthy individuals and diabetic patients. In our previous work, we demonstrated a strategy that accurately estimates blood glucose concentrations from sweat measurements. However, this method involves time-consuming simulations using a biophysical model, limiting its application to offline use. The goal of this study is to propose an approach that increases computational efficiency, thereby facilitating real-time estimation of blood glucose concentrations using sweat-sensing technology. To this end, we propose replacing the original biophysical model with the Local Density Random Walk (LDRW) model. This is justified because both models describe the pharmacokinetics of glucose transport through a convective-diffusion process. The performance of the LDRW model and the original biophysical model are compared in terms of estimation accuracy, computational efficiency, and model complexity, using seven datasets from the literature. The estimation of blood glucose concentrations using the LDRW model closely approximates that of the original model, with a root mean square difference of just 0.04 mmol/L between the two models' estimates. Remarkably, the LDRW model significantly reduces the average computational time to 2.6 s per data point, representing only 0.7% of the time required by the original method. Furthermore, the LDRW model demonstrates a smaller corrected Akaike Information Criterion value than the original method, indicating an improved balance between goodness of fit and model complexity. The proposed novel approach paves the way for the clinical adoption of sweat-sensing technology for non-invasive, real-time monitoring of diabetes.

Figures

Schematic illustration of the biophysical model for glucose transport

Impact of model parameters on the LDRW model response. (a) Effect of varying $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A$$\end{document}$ ( $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A=1, 2, 4$$\end{document}$ ) with $\docu

Flowchart of double-loop optimization strategy for inverse estimation of blood glucose concentrations based on sweat measurements. The optimization objective of loop 1 is to estimate the blood glucose concentration, and the optimization objective of loop 2 is to fine-tune the parameters of the LDRW model. The loops operate in an alternating manner, with one cycle involving loop 1 and the subsequent one involving loop 2

Comparison of experimental blood glucose concentrations with estimates from the biophysical and LDRW models using double-loop optimization. Experimental data come from Lee et al. [[40](#CR40)]

First sensitivity analysis: coefficient of variation (CV) of model parameters for (a) LDRW and (b) biophysical modelsFig. 7Second sensitivity analysis: coefficients of variation (CV) for LDRW model parameters in response to variations in biophysical model parameters

Full text

Introduction

Diabetes is a metabolic disorder characterized by elevated blood glucose concentrations due to insufficient production or ineffective use of insulin [1]. In recent years, this disease has gained increasing attention due to its prevalence across populations. In 2021, approximately 537 million people were living with diabetes, representing 10.5% of the overall population worldwide [2]. The mortality rate associated with diabetes and its complications was as high as 6.7 million deaths in 2021 [2]. Diabetes-related vascular complications such as cardiovascular disease, diabetic kidney disease, and diabetic retinopathy are major factors contributing to illness and death [3–5]. These complications not only reduce the overall quality of life and lead to higher mortality rates but also result in substantial economic burdens [4, 6], reaching 966 billion U.S. dollars globally in 2021 [2].

To prevent and mitigate diabetic complications, early diagnosis and timely medical intervention are pivotal [1, 7, 8]. A current challenge in diabetes management involves customizing medication doses and timing to meet the individual patient needs [9]. Specifically, medication doses should be adjusted in response to an individual’s current blood glucose concentrations. Real-time monitoring of blood glucose is essential, as it allows precise and timely adjustments of medication doses. However, common methods for monitoring blood glucose concentrations are typically invasive. The fingerstick test is the most prevalent method and it requires pricking the fingertip to obtain a small blood sample [10]. Another sampling concept employs microneedles for transdermal sampling of the interstitial fluid (ISF) [11]. These invasive techniques may cause discomfort and potential harm, and are prone to occlusions and challenges in obtaining adequate sample volume, limiting the frequency of sampling and potentially impacting patient compliance [12, 13]. The limitations of existing methods for blood glucose monitoring highlight the need for real-time, non-invasive alternative methods.

To enable non-invasive real-time monitoring of blood glucose concentrations, sweat-sensing technology emerges as a valuable option, offering an alternative to traditional blood and ISF sampling. However, research studying the relationship between blood and sweat glucose concentrations is still very limited. The study by Moyer et al. achieved only a moderate determination coefficient ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${R}^{2}$$\end{document} ) of 0.69 from a total of 115 experimental samples, utilizing linear regression analysis [14]. Klous and colleagues reported a correlation coefficient of 0.73 from a total of 48 data points using linear regression analysis [15]. Our previous research proposed a biophysical model to characterize glucose transport across various compartments as a convective-diffusion process, simulating the pharmacokinetic transport of glucose from blood capillaries to sweat excretion. Building on this model, we introduced a novel double-loop optimization strategy. This approach effectively solves the inverse problem, allowing for the estimation of blood glucose concentrations based on measured sweat glucose concentrations [16]. This double-loop strategy consists of two intertwined optimization loops, both employing an optimizer to reduce the error between measured and estimated sweat glucose concentrations. The first loop is dedicated to optimizing blood glucose concentrations, which are subsequently used as inputs for the biophysical model. The second loop aims at refining the parameters of the biophysical model in a personalized way. This approach achieved a high Pearson correlation coefficient ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R$$\end{document} ) of 0.98 using 108 experimental data points for estimating blood glucose from sweat glucose measurements.

While the double-loop optimization strategy is currently the most accurate method for estimating blood glucose concentrations via sweat-sensing technology, the substantial computational time requirements hamper its application for real-time monitoring. This strategy involves time-intensive simulations using a biophysical model implemented in COMSOL Multiphysics, requiring complete simulations for each iteration within both optimization loops. With multiple iterations typically necessary, the computational times generally range from 5 to 10 min per data point. Consequently, this duration far exceeds the thresholds for real-time estimation, necessitating more efficient computing. To address this challenge, our objective is to reduce the computational time of this strategy without compromising its accuracy. A simplified mathematical model, offering fewer parameters and reduced complexity, could be advantageous in improving computational efficiency through model reduction, providing a more feasible alternative to the original, complex biophysical model.

The biophysical model we employ to describe the convective-diffusion process of glucose transport is inherently complex. However, the Local Density Random Walk model (LDRW) can provide a simpler alternative for modeling convection–diffusion transport processes in pharmacokinetics. This model, originally developed to interpret and fit Indicator-Dilution Curves (IDC) in thermo- or dye-dilution studies, has been adapted for contrast-enhanced imaging applications [17–19]. The LDRW model operates by analyzing the probability density function of an indicator detected at a specific distance downstream following its upstream injection into an infinitely long tube. This transport is modeled as a convective-diffusion process, which describes the spatial and temporal evolution of the indicator concentration, solved under given boundary conditions. The LDRW model accommodates multiple passages of the indicator [17], providing an accurate representation of biomarker transport processes and making it well-suited for our study. Unlike the biophysical model, which uses 18 parameters and 10 equations as characterized in [16], the LDRW model simplifies the description of the biomarker transport process, e.g. for indicator molecules or particles, using a single analytical function with 3 parameters only. This model reduction, achieved through the LDRW model, is expected to accelerate the model identification by parameter estimation and enhance its reliability.

This study aims at developing a method for real-time estimation of blood glucose concentrations based on sweat glucose measurements. Building upon our previously introduced double-loop optimization strategy in [16], this study proposes replacing the original biophysical model with the LDRW model to enhance computational efficiency. To assess the effectiveness of our approach, we compare the results from this study with those of our previous work, in terms of estimation accuracy as well as computational time. Additionally, we conducted two separate sensitivity analyses: the first to evaluate the sensitivity of the glucose concentration estimates to the model parameters, and the second to explore the interlink between the parameters of the LDRW and the original biophysical model, thus investigating the physiological relevance of the LDRW model parameters.

Methods

Biophysical model and its impulse response

In our previous work, we proposed a compartmental biophysical model using COMSOL Multiphysics to simulate the transport of glucose from blood capillaries to sweat along a sweat gland [16], as illustrated in Fig. 1.Fig. 1. Schematic illustration of the biophysical model for glucose transport

In this section, we summarize the model working principle and we introduce the 18 model parameters used. The full model details are provided elsewhere [16]. The glucose transport is described considering three main steps. Firstly, glucose in the blood capillary compartment is transported to the interstitial fluid (ISF) compartment. This glucose flow, denoted as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${J}_{source}$$\end{document} in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$mol {s}^{-1}$$\end{document} , is driven by the difference in glucose concentrations between the blood capillary \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${C}_{p}$$\end{document} and the ISF compartment \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${C}_{ISF}$$\end{document} , both in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$mol {m}^{-3}$$\end{document} , given as

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${J}_{source}={k}_{DE}({C}_{p}-{C}_{ISF}){V}_{p},$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${V}_{p}$$\end{document} is the blood capillary compartment volume related to a single sweat gland in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${m}^{3}$$\end{document} [20], and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${k}_{DE}$$\end{document} represents the glucose clearance constant in the dermis in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${s}^{-1}$$\end{document} [21].

At the same time, water moves from the blood capillary compartment to the ISF compartment due to hydrostatic pressure gradients between these two compartments. The water flow rate, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Q}_{water}$$\end{document} , is governed by the Starling equation as

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Q}_{water}={L}_{p,c}{A}_{c}\left({P}_{c}-{P}_{ISF}\right),$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${P}_{c}$$\end{document} represents the hydrostatic pressure in the capillaries in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$mmHg$$\end{document} [22], \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${P}_{ISF}$$\end{document} represents the hydrostatic pressure in the ISF in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$mmHg$$\end{document} [22], \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${L}_{p,c}$$\end{document} is the hydraulic conductivity in the capillaries in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m{s}^{-1}{mmHg}^{-1}$$\end{document} [23], and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${A}_{c}$$\end{document} is the blood capillary compartment area related to a single sweat gland in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${m}^{2}$$\end{document} [20].

The water flow velocity in the ISF compartment, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${u}_{ISF}$$\end{document} ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m{s}^{-1}$$\end{document} ), is obtained by dividing the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Q}_{water}$$\end{document} by the ISF compartment's cross-sectional area \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${A}_{ISF}$$\end{document} ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${m}^{2}$$\end{document} ),

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${u}_{ISF}=\frac{{Q}_{water}}{{A}_{ISF}}.$$\end{document}

After glucose enters the ISF compartment, it is subject to diffusion and convection processes. During these processes, a fraction of the glucose flow, denoted as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${J}_{source}$$\end{document} , is absorbed by cells at a rate \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${r}_{uptake}$$\end{document} ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$mol{m}^{-3}{s}^{-1}$$\end{document} ) [24], given as

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{\partial {C}_{ISF}}{\partial t}=\frac{{J}_{source}}{{V}_{ISF}}-{r}_{uptake}+{D}_{ISF}\frac{{\partial }^{2}{C}_{ISF}}{\partial {y}_{ISF}^{2}}-{u}_{ISF}\frac{\partial {C}_{ISF}}{\partial {y}_{ISF}},$$\end{document}

with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${V}_{ISF}$$\end{document} the ISF compartment's effective volume in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${m}^{3}$$\end{document} [20], \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D}_{ISF}$$\end{document} is the glucose diffusion coefficient in the ISF in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${m}^{2}{s}^{-1}$$\end{document} [25], \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${y}_{ISF}$$\end{document} is the distance covered by glucose from its entry point into the ISF compartment in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m$$\end{document} . Notice that the glucose diffusion coefficient, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D}_{ISF}$$\end{document} , is considered spatially constant over the ISF in Eq. (2).

Due to the concentration gradient existing between the ISF compartment and the sweat gland compartment, glucose in the ISF subsequently transitions into the sweat gland. Diffusion primarily governs the passage of glucose across the sweat gland wall, which can be characterized using Fick's first law as.

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${J}_{ISF-sg}={D}_{sg,wall}\frac{\partial \left({C}_{ISF}-{C}_{sg}\right)}{\partial {h}_{sg}},$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${J}_{ISF-sg}$$\end{document} is the glucose flux entering the sweat gland in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$mol {s}^{-1}$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D}_{sg,wall}$$\end{document} represents the glucose diffusion coefficient in the sweat gland wall ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${m}^{2}{s}^{-1}$$\end{document} ) [26], \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${h}_{sg}$$\end{document} is the sweat gland wall thickness in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m$$\end{document} [27], and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${C}_{sg}$$\end{document} is the concentration of glucose in the sweat gland in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${mol m}^{-3}$$\end{document} .

Water also flows into the sweat gland compartment at a flow rate \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Q}_{water,sg}$$\end{document} in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${m}^{3}{s}^{-1}$$\end{document} , and subsequently travels along the sweat gland toward the skin surface. This process is governed by a pressure difference \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta P$$\end{document} between the ISF and the skin surface through the sweat gland in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$mmHg$$\end{document} [28] and can be given using Darcy's law as

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Q}_{water}=\frac{\Delta P}{R},$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R$$\end{document} represents the sweat gland's hydraulic resistance in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Pa\bullet s {m}^{-3}$$\end{document} [27], which can be given as

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R=\frac{128\eta L}{\pi {d}^{4}},$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta$$\end{document} is the water viscosity in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Pa\bullet s$$\end{document} [29], \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L$$\end{document} is the sweat gland length in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m$$\end{document} [30], and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d$$\end{document} is the sweat gland's luminal diameter in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m$$\end{document} [31].

The water flow velocity in the sweat gland, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${u}_{sg}$$\end{document} ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m{s}^{-1}$$\end{document} ), can be given as

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${u}_{sg}=\frac{{Q}_{water,sg}{u}_{sweat,n}}{{A}_{sg}},$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Q}_{water,sg}$$\end{document} is the ratio of the water flow rate, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${A}_{sg}$$\end{document} is the sweat gland's luminal area, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${u}_{sweat,n}$$\end{document} is the experimental sweat velocity normalized by the passive sweat velocity, which is \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3\times {10}^{-4} m{ s}^{-1}$$\end{document} [32]. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${u}_{sweat,n}$$\end{document} servers as a correction factor to adjust \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${u}_{sg}$$\end{document} for differences between passive sweating and experimental conditions, such as sauna, high-intensity exercise, and stimulated sweating.

Lastly, the glucose is transported through the sweat gland and ultimately reaches the human skin surface. This process can also be described as a convective-diffusion process, similar to Eq. (4), with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${C}_{sg,dil}$$\end{document} denoting the glucose concentration in the sweat gland compartment in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$mol {m}^{-3}$$\end{document} ,

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{\partial {C}_{sg,dil}}{\partial t}={D}_{sw}\frac{{\partial }^{2}{C}_{sg,dil}}{\partial {y}_{sg}^{2}}-{u}_{sg}\frac{\partial {C}_{sg,dil}}{\partial {y}_{sg}},$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D}_{sw}$$\end{document} is the spatially-invariant glucose diffusion coefficient in water in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${m}^{2}{s}^{-1}$$\end{document} [26].

The dilution effect of water influxes on glucose concentration at different sweat rates is given as

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${C}_{sg,dil}=\frac{{C}_{sg}}{1+{K}_{w/g}{u}_{sg,n}},$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{w/g}$$\end{document} represents the ratio between the volumetric flow rate of water and glucose [33], and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${u}_{sg,n}$$\end{document} is the normalized water velocity in the sweat gland ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${u}_{sg}$$\end{document} ) relative to the passive sweat velocity.

For a more detailed description of the biophysical model, please refer to our previous work [16]. All model parameters are listed in Table 1.Table 1. Parameter of the biophysical modelParameterUnitValueRefCapillary hydrostatic pressure: P_c_mmHg30[22]Capillary hydraulic conductivity: L_p,c_m s^−1^ mmHg^−1^6.5 × 10^–10^[23]Clearance constant of glucose: k_DE_s^−1^6.51 × 10^–4^[21]Effective volume of capillary: V_p_m^3^3.02 × 10^–13^[20]Effective volume of ISF: V_ISF_m^3^6.0 × 10^–13^[20]Hydrostatic pressure in ISF: P_ISF_mmHg−3[22]Glucose diffusion coefficient for sweat gland wall: D_sg,wall_m^2^ s^−1^6.46 × 10^–10^[26]Glucose diffusion coefficient in ISF: D_ISF_m^2^ s^−1^2.64 × 10^–10^[25]Glucose diffusion coefficient in water: D_sw_m^2^ s^−1^6.7 × 10^–10^[26]Luminal diameter of sweat gland: dm5 × 10^–6^[31]ISF cross-sectional area: A_ISF_m^2^2.2 × 10^–8^[20]Pressure difference across sweat gland: ΔPmmHg10[28]Ratio of volumetric flow rate of water to glucose: K_w/g_-12[33]Sweat gland length: Lm4 × 10^–3^[30]Sweat gland wall thickness: h_sg_m5 × 10^–5^[27]Surface area of capillary: A_c_m^2^1.5 × 10^–8^[20]Uptake rate of glucose: r_uptake_mol m^−3^ s^−1^2.78 × 10^–2^[24]Water viscosity: ηPa·s1 × 10^–3^[29]

The mechanism described by the biophysical model can be approximated as linear and time-invariant (LTI) within a short time frame. This approximation allows us to characterize the system through its impulse response, which in LTI systems describes how the system's output changes over time in response to an impulse. This response not only captures the characteristics of the system but also enables prediction of system's response to any input signal. Specifically, the system's output can be determined by the convolution of its input with the impulse response [34].

To investigate the impulse response of the biophysical model, we input a Dirac delta function, denoted as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta (t)$$\end{document} , into the model [35, 36]. The response of the biophysical model to this input is termed its impulse response, which is represented by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h(t)$$\end{document} .

LDRW model and model identification

As mentioned in the previous section, our biophysical model employs 18 parameters to characterize glucose transport, primarily as a process dominated by convection and diffusion, as described in Eqs. (4) and (9). The convective-diffusion process also forms the basis of extensive research focused on modeling the time evolution of the concentration of an indicator, represented by its IDC, following its rapid injection (see Fig. 2(a)) [19]. In general, a diffusion-convection process is described as

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{\partial {c}_{i}\left(z,t\right)}{\partial t}=\frac{\partial }{\partial z}\left({D}_{i}\frac{\partial {c}_{i}\left(z,t\right)}{\partial z}\right)-v\frac{\partial {c}_{i}\left(z,t\right)}{\partial z},$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${c}_{i}(z,t)$$\end{document} is the indicator concentration in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$mol {m}^{-3}$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D}_{i}$$\end{document} is the diffusion coefficient of the indicator in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${m}^{2}{s}^{-1}$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v$$\end{document} is the fluid flow rate in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m {s}^{-1}$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$z$$\end{document} is the coordinate along the longitudinal dimension in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m$$\end{document} .Fig. 2(a) Schematic illustration of the LDRW experimental model. (b) Flowchart for Identifying the Initial Parameters of the LDRW Model Using the Impulse Response of the Biophysical Model

The LDRW model provides a solution to the diffusion-convection Eq. (11) under specific boundary conditions [37]. The fast injection hypothesis is represented by the initial condition as

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${c}_{i}\left(z,0\right)=\frac{m}{S}\delta \left(z\right),$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m$$\end{document} is the solute dose in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$mol$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S$$\end{document} is the cross-sectional area of the tube in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${m}^{2}$$\end{document} . The boundary condition ensuring mass conservation during transport is expressed as

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\int }_{0}^{\infty }{c}_{i}\left(z,t\right)dz=\frac{m}{S}.$$\end{document}

Given these conditions, the concentration–time profile of the LDRW model at the detection section is given as:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$C\left(t\right)=A{e}^{\lambda }\sqrt{\frac{\lambda }{2\pi \mu t}}{e}^{-\frac{\lambda }{2}\left(\frac{t}{\mu }+\frac{\mu }{t}\right)},$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A$$\end{document} in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$mol s {m}^{-3}$$\end{document} is the ratio of the injected dose of solute \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m$$\end{document} in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$mol$$\end{document} to the volumetric flow of the carrier fluid \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Q$$\end{document} in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${m}^{3}{s}^{-1}$$\end{document} ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Q=vS$$\end{document} ), \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu$$\end{document} is the mean transit time of the solute from the injection to the detection section in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$s$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda$$\end{document} is a dimensionless parameter that indicates the skewness of the curve and is equal to half of the Peclet number ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Pe/2$$\end{document} ), which represents the ratio of convection to diffusion of a solute. These three parameters define the LDRW model. The impact of variations in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu$$\end{document} , and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda$$\end{document} on the LDRW model response is illustrated in Fig. 3.Fig. 3. Impact of model parameters on the LDRW model response. (a) Effect of varying \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A$$\end{document} ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A=1, 2, 4$$\end{document} ) with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu =1 s$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda =1$$\end{document} . (b) Effect of varying \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu$$\end{document} ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu =1, 2, 4$$\end{document} ) with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A=1 mol {s}^{-1}{m}^{-3}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda =1$$\end{document} . (c) Effect of varying \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda$$\end{document} ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda =1, 2, 4$$\end{document} ) with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A=1 mol{ s}^{-1}{ m}^{-3}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu =1 s$$\end{document}

The initial parameter values for the LDRW model, denoted as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\theta }_{0}=\left\{{A}_{0}, {\lambda }_{0}, {\mu }_{0}\right\}$$\end{document} , were identified by fitting the impulse response, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h(t)$$\end{document} , derived from the biophysical model. This fitting was accomplished using a multiple linear regression (MLR) algorithm in the logarithmic domain, developed by our group [37], which has demonstrated optimal performance for LDRW fitting, as detailed in our prior works [17, 37]. The whole process is illustrated in Fig. 2(b).

Inverse estimation by double-loop optimization

Inverse estimation is a method used to infer the input of a system based on its output. In our study, this involves estimating blood glucose concentrations from sweat measurements. This approach is particularly important for enabling non-invasive monitoring of blood glucose concentrations, given the unclear relationship between sweat and blood glucose concentrations. The overall strategy for inverse estimation using the LDRW model is outlined in Fig. 4 and is explained in detail below:Fig. 4. Flowchart of double-loop optimization strategy for inverse estimation of blood glucose concentrations based on sweat measurements. The optimization objective of loop 1 is to estimate the blood glucose concentration, and the optimization objective of loop 2 is to fine-tune the parameters of the LDRW model. The loops operate in an alternating manner, with one cycle involving loop 1 and the subsequent one involving loop 2

Initially, the complex biophysical model is used only once to characterize the system via its impulse response \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h(t)$$\end{document} . \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h(t)$$\end{document} is obtained by applying \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta (t)$$\end{document} to the biophysical model, which is parameterized using values from the literature listed in Table 1. The initial parameters of the LDRW model, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\theta }_{0}=\left\{{A}_{0}, {\lambda }_{0}, {\mu }_{0}\right\}$$\end{document} , are then identified by fitting \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h(t)$$\end{document} with the LDRW model expression. Subsequently, the double-loop optimization strategy focuses on refining the three parameters of the more computationally efficient LDRW model, rather than the 18 parameters of the complex biophysical model. In the double-loop strategy, the first loop optimizes the model input, specifically the estimated blood glucose concentration, starting with a typical initial value \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${C}_{blood glu,0}$$\end{document} of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$5.5 mmol {L}^{-1}$$\end{document} , based on average values in healthy individuals [38]. The second loop focuses on fine-tuning the model parameters, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta$$\end{document} , to tailor the LDRW model in a personalized way. These loops are intertwined, with one iteration for input optimization followed by one iteration for parameter optimization. This iterative process refines the estimated blood glucose concentrations and the model parameter values. For each iteration \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i$$\end{document} , the error \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${e}_{i}$$\end{document} , defined as the mean square error between the estimated sweat glucose \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\widehat{C}}_{sweat glu,i}$$\end{document} and the experimental sweat glucose \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${C}_{sweat glu,i}$$\end{document} , is calculated. The optimization process continues until \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${e}_{i}$$\end{document} falls below 0.001. To mitigate the risk of convergence to local optima, both loops in the algorithm are executed only once per optimization iteration. For a more comprehensive explanation of the double-loop optimization strategy, readers can refer to our previous work [16].

In our implementation, MLR in the logarithmic domain provides a fast and robust initialization of the LDRW model parameters. This linear approach avoids the convergence issues faced by traditional nonlinear methods when dealing with noisy data [37]. Yet, applying a log transformation, the noise characteristics change too, making a least-squares optimization sub-optimal. Therefore, following this initialization step, the double-loop strategy employs a sparse nonlinear optimizer (SNOPT) [39] without log-transformation in both optimization loops. SNOPT’s sequential quadratic programming framework efficiently handles large-scale, sparse problems, making it particularly suitable for multi-parameter estimations. Specifically, it is ideal for iteratively refining the biophysical model's 18 parameters and the estimated blood glucose concentration \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\widehat{C}}_{blood glu}$$\end{document} . Additionally, its line search approach ensures that the algorithm dynamically adjusts the step length to converge to an optimal solution even with suboptimal initialization.

Data sources

The described strategy is tested on seven sets of experimental data from the literature [14, 15, 40–43], as detailed in Table 2. The experiments were conducted under different conditions, including meal intake, fasting, and glucose ingestion, to induce variations in blood glucose concentrations. The participants included 18 healthy individuals and 2 diabetic patients, resulting in a total of 108 paired blood and sweat glucose measurements. Such diversity in experimental conditions and subject profiles is critical for validating the robustness and generalizability of the proposed estimation strategy across different physiological states. Specifically, datasets corresponding to Experiments 1–6 (as listed in Table 2) are publicly available through the following references: [14, 40–43]. The original data and experimental procedures can be accessed in the respective publications. The dataset corresponding to Experiment 7 was acquired through collaboration with Vrije Universiteit Amsterdam, as detailed in [15]. This study was approved by the Medical Ethics Committee of Erasmus MC Rotterdam (MEC-2019–0202) and conducted in accordance with the revised Declaration of Helsinki (2013).Table 2. Experimental dataset of sweat and blood measurementsData SourcesParticipanthealth conditionSweat glucose(Mean ± Std)[μmol L^−1^]Blood glucose(Mean ± Std)[mmol L^−1^]MeasuringpointsnumberExp 1[40]1 healthy97.1 ± 13.65.6 ± 0.915Exp 2[41]2 healthy37.5 ± 22.411.4 ± 5.719Exp 3[14]1 diabetic178.3 ± 34.315.4 ± 2.16Exp 4[42]1 healthy19.3 ± 12.15.0 ± 0.17Exp 5[42]1 healthy42.4 ± 26.44.7 ± 0.65Exp 6[43]1 Healthy and 1 Diabetic69.4 ± 36.36.4 ± 1.28Exp 7[15]12 healthy12.8 ± 4.94.6 ± 0.648

Model evaluation and performance metrics

The evaluation procedure of both the biophysical model and the LDRW model was organized into four aspects: estimation performance, computational efficiency, model complexity, and parameter sensitivity. The detailed assessment methods are described below.

Estimation performance

Firstly, we computed several performance metrics to evaluate and compare the estimation results of the biophysical model and the LDRW model against the corresponding experimental blood glucose values. These metrics included the Root Mean Square Error (RMSE), Mean Absolute Error (MAE), Root Mean Square Percentage Error (RMSPE), the Pearson correlation coefficient ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R$$\end{document} ), and the coefficient of determination ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${R}^{2}$$\end{document} ) [44]. Secondly, to assess the agreement between the two models directly, we calculated the Root Mean Square Difference (RMSD) between their respective estimated blood glucose concentrations for each experimental data point. To evaluate whether the differences in estimation results between the two models were statistically significant, a paired two-tailed Wilcoxon signed-rank test was performed.

Computational efficiency

To evaluate the computational speed of the biophysical and LDRW models, the computational time for each model under identical testing conditions was recorded. Specifically, we documented the computational time per data point for the two models and calculated the ratio of the computational time required by the biophysical model to that of the LDRW model. A ratio exceeding 1 indicates that the biophysical model is more time-consuming, while a ratio less than 1 implies that it is faster. A paired two-tailed Wilcoxon signed-rank test was conducted on the computational time per data point for each model to assess whether the observed differences in computational efficiency was statistically significant.

Model complexity

To assess the goodness-of-fit and complexities of both models, the corrected Akaike Information Criterion (AICc) values were calculated. While the standard Akaike Information Criterion (AIC) is commonly used for model selection [45], the AICc is preferred in cases of limited sample sizes, as it minimizes the risk of overfitting [46]. In our study, we computed the AICc values by comparing the estimated blood glucose concentrations from the double-loop optimization strategy for both the LDRW and biophysical models against the experimental values. The AIC is formulated as

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$AIC=2k-2\text{ln}\left(L\right),$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k$$\end{document} is the number of estimated parameters in the model, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L$$\end{document} is the maximum value of the likelihood function for the model. Given the limited sample size \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n$$\end{document} , the AICc was employed, which is given as

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$AICc=AIC+\frac{2k\left(k+1\right)}{n-k-1}.$$\end{document}

In model selection, the model with lower AICc value is considered optimal, as it represents the better trade-off between goodness of fit and model complexity.

Sensitivity analysis

Two sensitivity analyses were conducted on the two models'parameters. The first analysis aimed to evaluate how variations in the parameters of both the LDRW and biophysical models influence the simulated sweat glucose concentrations, providing insights into the extent to which these parameters can be accurately estimated. For the LDRW model, we analyzed three parameters: the ratio of the solute to the volumetric flow ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A$$\end{document} ), the mean transit time ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu$$\end{document} ), and the skewness of the curve ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda$$\end{document} ). The biophysical model, being more complex, consisting of 18 parameters. For each parameter in both models, 100 parameter samples were generated using a random number generator that followed a Gaussian distribution. These samples were centered around values sourced from the literature as listed in Table 1, with a standard deviation of 10%. The two models'estimation results were then evaluated to calculate coefficients of variation (CV).

The second analysis aimed at evaluating the relationship between the biophysical and LDRW models, exploring how changes in the biophysical model influenced the derived LDRW model parameters. To clarify, this analysis focused on varying only those biophysical model parameters that exhibited a CV greater than 0.1% in the results of the first sensitivity analysis. Biophysical model parameters with a CV less than 0.1% were deemed to have limited impact on the model's estimation performance and were not considered for further analysis. For each selected biophysical model parameter, 100 samples were generated using a random number generator following a Gaussian distribution centered around referenced literature values listed in Table 1, with a standard deviation of 10%. We then determined the corresponding impulse response of the biophysical model and used the LDRW model to fit these impulse responses. This approach yielded a corresponding distribution of LDRW model parameters, enabling us to assess the interplay between the two models by analyzing the calculated CV.

Implementation details

The biophysical model simulations were conducted using COMSOL Multiphysics version 6.1. The LDRW model fitting and the inverse estimation by double-loop optimization were implemented in MATLAB R2022a. All simulations and analyses were performed on a personal workstation equipped with an Intel Core i7-12700 CPU, 16 GB RAM, and Windows 11 (64-bit) operating system. No GPU acceleration was used.

Results

Estimation performance

The root mean square difference quantifies the difference between the estimated blood glucose concentrations obtained from the biophysical model and those estimated by the LDRW model, amounting to 0.04 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$mmol {L}^{-1}$$\end{document} across all 108 experimental data points. Both models'estimated blood glucose concentrations showed a Pearson correlation coefficient \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R$$\end{document} of 0.98 with the experimental measurements across all data points. The RMSE, MAE, RMSPE, and coefficient of determination further quantify the degree of alignment between the model's estimates and the experimental data, as detailed in Table 3. Shown in Fig. 5 are the estimated blood glucose concentrations from one illustrative experiment, chosen from among the seven experimental datasets considered, obtained by applying the double-loop optimization strategy with both models. A paired two-tailed Wilcoxon signed-rank test was conducted on the estimation results obtained using each model. The result (p = 0.30) indicates that there is no statistically significant difference between the two models in terms of estimation results.Table 3. Comparison of estimation performance metrics for the LDRW and biophysical modelsMetricsLDRWBiophysical ModelRMSE (mmol L^−1^)0.850.85MAE (mmol L^−1^)0.450.44RMSPE (%)12.41%12.36%Pearson correlation coefficient R0.980.98Coefficient of determination R^2^0.960.96Fig. 5Comparison of experimental blood glucose concentrations with estimates from the biophysical and LDRW models using double-loop optimization. Experimental data come from Lee et al. [40]

Computational efficiency

To compare the computational efficiency of the double-loop optimization strategy when applied to the LDRW and biophysical models, we recorded the computational times required by each model, as shown in Table 4. The biophysical model requires approximately 6 min per data point, whereas the LDRW model requires only 2.6 s, nearly achieving real-time estimation. A paired two-tailed Wilcoxon signed-rank test was conducted on the computational time per data point for each model to assess whether the observed difference in computational efficiency was statistically significant. The result (p =  \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2.55\times {10}^{-16}$$\end{document} ) indicates that the difference is highly significant, confirming that the LDRW model is significantly faster than the biophysical model.Table 4. Comparison of computational time per data point: LDRW vs. biophysical modelsData SourcesComputational Timeper Data point forLDRW model (s)Computational Time Ratio(Biophysical model/LDRW model)Exp 1[40]2.3 ± 0.3151.0Exp 2[41]3.8 ± 0.9128.8Exp 3[14]3.1 ± 0.7140.0Exp 4[42]1.8 ± 0.1162.3Exp 5[42]1.9 ± 0.2173.9Exp 6[43]3.6 ± 0.5127.7Exp 7[15]2.0 ± 0.2149.7Average2.6 ± 0.8151.2 ± 17.0

Model complexity

The goodness of fit and model complexity of the biophysical model and LDRW model, which were each separately utilized in the double-loop algorithm, were evaluated using AICc. The AICc value for the LDRW model, at 277.3, was lower compared to the biophysical model, which had an AICc value of 313.4.

Sensitivity analyses

The first sensitivity analysis was conducted to assess the impact of variations in individual parameters on the accuracy of sweat glucose concentration estimations by the LDRW and biophysical models. In presenting the results in Fig. 6, it should be noted that for the biophysical model, only its parameters exhibiting a CV greater than 0.1% are included. The analysis result indicates that within the LDRW model parameters, the parameter with the largest influence is the ratio of the solute to the volumetric flow ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A$$\end{document} ), which showed a CV of 9.08%. In the biophysical model, the ratio of the volumetric flow rate of water to glucose ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{w/g}$$\end{document} ) results in the largest CV, equal to 18.19%. The second sensitivity analysis was designed to explore the relationship between variations in individual parameters of the biophysical model and the corresponding changes in the LDRW model parameters. The results are presented in Fig. 7, where a heatmap visualizes the CVs. Within the heatmap, darker shades indicate greater changes in the LDRW parameters resulting from specific variations in the biophysical model parameters.Fig. 6. First sensitivity analysis: coefficient of variation (CV) of model parameters for (a) LDRW and (b) biophysical modelsFig. 7Second sensitivity analysis: coefficients of variation (CV) for LDRW model parameters in response to variations in biophysical model parameters

Discussion

In this study, we propose substituting our previously proposed biophysical model with the LDRW model to enhance computational-efficiency in estimating blood glucose concentrations from sweat measurements, employing a double-loop optimization strategy.

Model performance

Estimation accuracy

The estimation results obtained from the double-loop optimization strategy using both the LDRW model and the biophysical model show minimal discrepancies, with a low root mean square difference of 0.04 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$mmol {L}^{-1}$$\end{document} . Across all 108 data points, the two model's performance in estimating blood concentrations of glucose from sweat is almost identical, as evidenced by an RMSE of 0.85 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$mmol {L}^{-1}$$\end{document} and a Pearson correlation coefficients ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R$$\end{document} ) of 0.98. These findings support the replacement of the original biophysical model with the LDRW model, which is specifically designed for diffusion-convection processes. These results demonstrate that the LDRW model retains the estimation accuracy of the biophysical model.

Computational efficiency

The use of the LDRW model significantly enhances the computational efficiency of the double-loop optimization. As reported in Table 4, the average computational time per data point for the LDRW model across seven experiments was 2.6 s, compared to 378.9 s for the biophysical model. This substantial reduction in computational time underscores the efficiency of employing the LDRW model within the double-loop optimization strategy. It is worth noting that variations in computational time per data point across different experiments are largely due to differences in experimental conditions, such as the extent of subject's blood glucose concentration fluctuations. For instance, the second experiment (EXP 2 in Table 2), which recorded the longest computational times per data point, exhibited the greatest fluctuations in blood glucose concentrations, with a standard deviation of 5.7 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$mmol {L}^{-1}$$\end{document} . Conversely, the forth experiment (EXP 4 in Table 2), with the shortest computational times, had the most stable blood glucose concentrations, indicated by a standard deviation of only 0.1 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$mmol {L}^{-1}$$\end{document} . These factors contribute to the observed differences in computational times between these experiments. It is also important to note that in typical glucose monitoring scenarios, glucose measurements are required at relatively low frequencies, typically ranging from 5 to 15 min [47, 48]. The computational speed achieved by the LDRW-based approach ensures timely estimations, making it well-suited for real-time or quasi-real-time glucose monitoring applications.

Model complexity

With an AICc value of 277.3, the LDRW model demonstrates a superior balance between model fit quality and model complexity, compared to the biophysical model's higher AICc of 313.4. The lower AICc value indicates that the model is better at explaining the characteristics of the data with fewer parameters. The LDRW model, which utilizes 3 parameters only, achieves performance nearly identical to that of the biophysical model, which uses 18 parameters. This notable reduction in the number of model parameters is crucial, as it not only minimizes the risk of overfitting but also enhances the robustness and reliability of model estimation [49].

Sensitivity analysis

Parameter sensitivity

The results from the first sensitivity analysis demonstrate that parameters from both the LDRW and biophysical models substantially influence the estimation of sweat glucose concentrations. Within the LDRW model, the parameter \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A$$\end{document} , representing the ratio of solute to volumetric flow, exhibited a CV of 9.08%, making it the second most influential parameter after the biophysical model's \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{w/g}$$\end{document} . In contrast, the parameter \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda$$\end{document} in the LDRW model, with the smallest CV at 0.30%, ranks higher in sensitivity compared to the fifth most sensitive parameter in the biophysical model, where a total of 18 parameters were assessed but only those with a CV greater than 0.1% are displayed in Fig. 6. This finding underscores a comparable level of sensitivity across the two models.

Parameter interpretability

The second sensitivity analysis of the relationship between the parameters of the biophysical and LDRW models supports the possibility of a physiological interpretation of LDRW model parameters. Among all the biophysical model parameters tested, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{w/g}$$\end{document} was found to have a most notable influence on the LDRW model parameter \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A$$\end{document} , as evidenced by the highest CV of 19.0% for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A$$\end{document} . This correlation indicates a strong dependency of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A$$\end{document} on variations in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{w/g}$$\end{document} . This is theoretically expected, given that \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{w/g}$$\end{document} , representing the ratio of the volumetric flow rate of water to the volumetric flow rate of glucose, is conceptually similar to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A$$\end{document} , which represents the ratio of solute (glucose) to volumetric water flow.

Furthermore, changes in the biophysical model's luminal diameter of the sweat gland ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d$$\end{document} ) predominantly affect the LDRW model parameter \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu$$\end{document} . Following \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d$$\end{document} , the biophysical model's parameters \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta P$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta$$\end{document} , and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L$$\end{document} exert the next most substantial influence on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu$$\end{document} . This is expected, as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu$$\end{document} represents the transit time of glucose, which primarily depends on the velocity of water flow in the sweat gland ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${u}_{sg}$$\end{document} ). According to Eq. (8), \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${u}_{sg}$$\end{document} is calculated from the ratio of the water flow rate inside the sweat gland ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Q}_{water,sg}$$\end{document} ) to the luminal area ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${A}_{sg}$$\end{document} ), which is directly related to its diameter \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d$$\end{document} . Additionally, as per Eq. (6), \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Q}_{water,sg}$$\end{document} is determined by the pressure difference across the sweat gland ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta P$$\end{document} ) and the hydraulic resistance, which correlates with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L$$\end{document} as outlined in Eq. (7).

Similarly, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d$$\end{document} was also found to considerably influence the LDRW model parameter \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda$$\end{document} , with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda$$\end{document} exhibiting the highest CV of 25.6%. Following \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta P$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta$$\end{document} , and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L$$\end{document} also have the next most substantial impact on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda$$\end{document} . This influence is linked to the definition of the parameter \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda$$\end{document} , which represents half of the Peclet number (Pe/2), quantifying the ratio of convective to diffusive transport [17]. The Peclet number is calculated as the product of the compartment length ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L$$\end{document} ) and water flow velocity ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${u}_{sg}$$\end{document} ) divided by the diffusion coefficient. This relationship explains why changes in diffusion coefficients, such as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D}_{sw}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D}_{sg,wall}$$\end{document} , notably impact the LDRW model parameter \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda$$\end{document} . In addition, in the biophysical model, the glucose uptake by cells in the ISF compartment is represented by the term \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${r}_{uptake}$$\end{document} , acting as a sink term in the diffusion-convection Eq. (4). In contrast, the LDRW model simplifies the glucose transport by omitting the sink term. This simplification is supported by our sensitivity analysis, with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${r}_{uptake}$$\end{document} showing a CV of less than 0.1%, which indicates its negligible influence on the model outputs.

In summary, the parameter values of the LDRW model closely correspond with those of the biophysical model and are amenable to physiological interpretation. This connection and interpretability are crucial for supporting the application of the LDRW model in clinical settings.

Limitations and future perspectives

The results presented in this paper are based on the limited datasets available in the literature, highlighting the need for future research to incorporate additional data to enhance statistical power. Moreover, the described strategy applies only to biomarkers that rely solely on diffusion-convection transport mechanisms, as exemplified by glucose. Therefore, this method is not directly applicable to biomarkers involving more complex transport mechanisms, such as active transport. This limitation is due to the inherent characteristics of the LDRW model, which is primarily designed to address processes dominated by diffusion and convection. Specific modifications would be necessary to adapt the approach to accommodate these alternative transport mechanisms.

Future research can build upon this work by extending the biophysical model to encompass the transport processes of other biomarkers, such as urea and lactate. Adapting the existing model to these biomarkers would enable the proposed strategy to monitor a broader range of diseases. Such an expansion could significantly enhance the applicability of sweat-sensing technology in disease monitoring, allowing for semi-continuous and non-invasive healthcare solutions.

Conclusion

In this study, we introduced a computationally efficient method for estimating blood glucose concentrations from sweat measurements, by replacing the previously published 18-parameter biophysical model with a three-parameter LDRW model. Our earlier study demonstrated that the biophysical model significantly outperformed existing state-of-the-art methods. The present results confirm that the LDRW model remains comparable accuracy (R = 0.98) while significantly reducing computational time to only 0.7% of that required by the original method, with an average of 2.6 s per data point. This acceleration brings the method well within the sampling intervals required for current continuous glucose monitoring workflows, enabling semi-continuous blood glucose tracking and proving timely insights for both patients and clinicians. The proposed approach paves the way for prolonged diabetes management through non-invasive sweat-sensing technology.