Dielectrical Properties of Living Epidermis and Dermis in the Frequency Range from 1 kHz to 1 MHz.

We determine the in-vivo dielectric properties-resistivity and relative permittivity-of living epidermis and dermis of human skin soaked with a physiological saline solution for one minute between 1 kHz and 1 MHz. This is done by fitting approximate analytical solutions of a mechanistic model for the transport of charges in these layers to a training set comprising impedance measurements at two depth settings on stripped skin on the volar forearm of 24 young subjects. Here, the depth settings are obtained by varying the voltage at a second inject on the electrical-impedance-spectroscopy probe. The model and the dielectric properties are validated with a test set for a third depth setting with overall good agreement. In addition, the means and standard deviations of the thicknesses of living epidermis and dermis are estimated from a literature review as 61±7 μm and 1.0±0.2 mm respectively. Furthermore, extensions to resolve the skin layers in more detail are suggested.

Introduction

The human skin comprises several layers—stratum corneum (SC), living epidermis (LE), dermis (DE) and hypodermis— with different dielectric properties that are functions of not only the frequency of the applied voltage but also other factors such as age [1] and health condition [2] during electrical impedance spectroscopy (EIS) measurements.

The dielectric properties of human skin have been studied experimentally for different locations and frequencies. A database treating the skin as a one-layer entity for both wet and dry skin has been compiled for in-vivo EIS measurements in the frequency range 100 Hz to 20 GHz [3]. A higher resolution of the inherent layers was achieved through tape-stripping and in-vivo EIS, which established the dielectric properties of stratum corneum from 1 kHz to 1 MHz [4, 5, 6] and viable skin with living epidermis and dermis lumped together from 1 kHz to 1 MHz [5,6]. In addition, in-vitro EIS measurements of epidermis and devitalized dermis in a nutrient liquid saturation was carried out from 1 Hz to 100 kHz [7,8].

In view of the lack of in-vivo dielectric properties in the living epidermis and dermis between 1 kHz and 1 MHz, we apply and extend our old methodology of fitting a mechanistic mathematical model to experimental EIS measurements [6,9]. In essence, the methodology was applied to determine the dielectric properties of viable skin (living epidermis and dermis combined into one layer) and stratum corneum from two separate measurement series of stripped and intact skin respectively. The advantage of the mechanistic model is that it captures the local transport of charges throughout the skin layers and the EIS probe (see Fig. 1a-b); its main disadvantage, however, is the significant computational overhead since a complex-valued partial differential equation (PDE) needs to be solved numerically in several domains. We have recently managed to secure approximate analytical solutions [10,11] that reduce the multidomain PDE to a set comprising improper integrals and closed-form equations. This reduction is key when we here fit the dielectric properties of both living epidermis and dermis at the same time, because it becomes significantly more computationally difficult to fit four parameters at once—as compared to fitting only two for viable skin [9]. Besides introducing the approximate analytical solutions, we also tap into more information from the EIS measurements of stripped skin: the different charge pathways underneath the four-electrode, two-point-measurement EIS probe for the depth settings, which were obtained by varying the voltage at the second inject (see Fig. 1). A total of five depth settings are employed, out of which we choose depth settings 1 and 5 as the training set for parameter fitting and depth setting 3 as the test set for validation.

Fig.1

Schematic of (a) the four-electrode EIS probe with a sense (1), guard (2), second inject (3) and first inject (4); (b) the electrodes of the probe and the stripped skin comprising living epidermis and dermis; and (c) the reduced model consisting of dermis with the living epidermis and electrodes reduced to boundary conditions.

The layout of the paper is as follows: First, we summarize the experiments from our earlier study [9] and estimate the thicknesses of living epidermis and dermis from a literature survey. We then work through the main steps in the derivation towards the approximate analytical solution, which combines the idea of model reduction [10] and the generic, multi-electrode and -layer approximate solution [11]. After an outline of the parameter fitting, we discuss the experimental measurements at the three depth settings as well as the fitted dielectric parameters. Finally, we suggest extensions that seek to resolve the skin layers in more detail and end with conclusions.

Materials and methods

Experiments

The experimental EIS measurements are from a previous study, which was carried out after ethics approval and informed consent [9].

In short, the study comprised 26 healthy, non-smoking subjects aged 27±6 (mean ± standard deviation) with an equal distribution of men and women without any known skin diseases or allergies; on the day of the measurements, they were asked to abstain from applying moisturizers. Before taking the EIS measurements with an in-house impedance spectrometer, SciBase I, the subjects' skin were stripped eighty times with Scotch® MagicTM Cellulose Tape on the volar forearm and then soaked with a physiological saline solution (at 0.9% salt concentration by mass) for 1 minute. Among the 26 measurements, two outliers were identified and excluded from further analysis: Their measured impedances were 15 kΩ and 50 kΩ at 1 kHz with no signs of double dispersion.

The circular noninvasive EIS probe—two voltage injection electrodes, a current detector and a guard ring— allows for 2-point measurements between 1 kHz and 1 MHz at 31 frequencies logarithmically distributed for five different current penetration depth settings, which arise from varying the applied voltage at the second injection electrode from 5 to 50 mV. As mentioned in the Introduction, in our earlier studies [6,9], only the depth setting 5 was employed; here, we use three depth settings: 1, 3 and 5.

Skin Thickness

For the mathematical model, we need the thickness of living epidermis as well as dermis for young subjects on the volar forearm. We calculated these in terms of weighted means and pooled standard deviations from a literature review summarized in Table 1, assuming that all data sets are normal-distributed:

Tab. 1

Skin thickness (mean value ± standard deviation) at the volar forearm for epidermis, stratum corneum and dermis. Some measurement techniques are: Confocal Microscopy (CM), Confocal laser scanning microscopy (CLSM), Optical coherence tomography (OCT), Transepidermal water loss (TEWL), Confocal Raman spectrometer (CRS) and Ultrasound (US).

ED (μm)	Subj	Age	Measurement	Source
73.54±2.84	10	23 - 47	CM	[12]
89±9	15	19 - 24	CLSM	[13]
71.8±10	30	21 - 35	OCT	[14]
102±7.4	5	20 - 35	OCT	[15]
66.39±5.82	19	20 - 29	OCT	[16]
65.1±8.9	8	31 - 37	OCT	[17]

Here, is the weighted mean, X is the mean in study i, n is the number of subjects in study i, m is the total number of studies, σ is the standard deviation of study i, σ is the pooled standard deviation.

First, we estimated the thickness of epidermis (ED) as 75±9 μm and the thickness of stratum corneum as 14±4 μm. By subtracting the mean thickness of stratum corneum from the thickness of epidermis and pooling their standard deviations, we found an estimate for the living epidermis thickness as 61±7 μm. Finally, we quantified the thickness of dermis as 1.0±0.2 mm.

Note that we combined the measurements from both males and females since we found no significant differences between them.

Mathematical Model

We consider conservation of charge—ohmic and displacement currents—in the electrodes of the EIS probe and two layers—living epidermis and dermis—in the stripped skin as illustrated schematically in Fig. 1b and summarized in Appendix A. We further employ scaling arguments [10] to reduce the living epidermis and electrodes to boundary conditions, as illustrated in Fig. 1c, which yields a reduced model, outlined in Appendix B. The reduced model is then solved with a Hankel transform (see Appendix C), resulting in the following approximate analytical solution for the predicted impedance, Z, for a given depth setting l:

where

with

where A refers to the the i-th row and j-th column of the 4×4 matrix A, I is an integral from the Hankel transform, J1() refers to Bessel function of the first kind.

The parameters for the model are summarized in Table 2.

Tab. 2

Dimensions and material parameters.

Current detection width, w1	1 mm
Ceramic width, w2, w4, w6	0.15, 0.15, 1.9 mm
Guard width, w3	0.3 mm
Secondary inject width, w5	0.5 mm
Primary inject width, w7	0.5 mm
Living epidermis thickness, hLE	61 μm
Dermis thickness, hDE	1.0 mm
Inject voltage, V0	0.05 V
Depth setting, α	0.1, 0.58, 1
Electrical permittivity in vacuum,	8.85 x10-12 Fm-1
ɛ0 [28]

Parameter Fitting

We employed the genetic optimization algorithm in Matlab [29] together with the measured and predicted impedance from Eq. 2 to fit the relative permittivity and resistance of the living epidermis and dermis throughout the frequency interval of 1 kHz to 1 MHz. The optimization was carried out by minimizing the objective function, defined as

for every subject and each measured frequency; i.e. a total of 744 times for our 24 subjects and 31 frequencies. Here, Re and Im denote the real and imaginary parts of the predicted impedance, Z, and measured counterpart, at depth settings l=1 (α=1) and l=5 (α=0.1). In addition to the training set based on depth settings 1 and 5, a test set was taken from the middle depth, l=3 (α=0.58) for validation of the fitted dielectric properties. The depth settings 1 and 5 were chosen for the training set as they are the furthest apart in terms of the current injected; and the depth setting 3 was selected because it lies in between the two depths for training.

From a computational point of view, we note that the improper integrals, I , only need to be solved numerically once before the fitting of the relative permittivity and resistivity, which is computationally efficient.

Ethical approval

The experimental EIS measurements were gathered under ethical approval and patient-informed consent.

Results and Discussion

Training set

The experimentally measured impedance—magnitude and phase in the frequency range from 1 kHz to 1 MHz—for the training set is shown in Fig. 2a,b. Here, the magnitude of the impedance is higher for depth setting 5, because less current is injected as compared to depth setting 1; similarly, the phase is less for depth setting 5. We also note that while the measured magnitude remains smooth and monotonically decreasing, the phase starts to fluctuate for frequencies around 100 kHz and higher.

Fig. 2

The mean (depth 1 ●; depth 5 ■] and standard deviations (error bars) for the (a) magnitude and (b) phase of the experimentally measured impedances. The fitted mean resistivity and relative permittivity and their standard deviations can be found in (c) and (d) respectively for the living epidermis (♦) and dermis (▼); the fitted equations, Eqs. 7-8, are shown as lines. In (e) and (f), the validation set for depth setting 3 is shown for the experiments (●) and model predictions (lines).

According to a Lilliefors test at a significance level of 0.05 for each frequency, the magnitudes of depth setting 1 are all normal-distributed, whereas 71% of the magnitudes are normal-distributed for depth setting 5; similarly, 100% and 90% of the phases are normal-distributed for depth setting 1 and 5 respectively. We therefore assume that the measurements are normal-distributed.

In order to test that the means of the magnitude and phase of the impedance are not the same for depth 1 and 5, a paired t-Test at each of the 31 measured frequencies gives P-values, P≲10⁻8.Even taking into account that we are carrying out multiple t-Tests and applying, e.g., a Bonferroni correction, we can state that the means are significantly different for our purposes. Here, we have employed a paired t-Test, because we paired two depths for each calibration of a given subject at each frequency. Likewise, a paired t-Test gives similarly small P-values (P≲10⁻8for ν≲105Hz) throughout the frequency range for the imaginary and real parts except around 1 MHz where the P-values for the imaginary parts become O(10⁻1-10⁻2).

Dielectric Properties

The means and standard deviations of the fitted resistivities and relative permittivities for the 24 subjects and 31 frequencies are depicted in Fig. 2c-d. The aforementioned measurement fluctuations only propagate to the fitted relative permittivity of the dermis above 100 kHz; the resistivities and relative permittivity of the living epidermis do not exhibit any significant fluctuations.

The mean dielectric properties for LE and DE can conveniently be written as exponential functions of fifth-order polynomials:

Here, and are the material coefficients for the resistivity, ρ, and relative permittivity, ε, of the skin layers, given in Table 3; is a dimensionless frequency and N is the uppercase greek letter of ν, introduced for notational convenience; σ (= 1/ ρ), is the conductivity of layer k. These polynomials give overall smooth, monotonically decreasing dielectric properties without the fluctuations.

Tab. 3

Coefficients for the conductivity and relative permittivity of the skin layers for 1 kHz to 1 MHz.

j	cjLE	djLE
0	1.6306×10-2	-2.3922×10-2
1	-3.3591×10-1	5.5919×10-1
2	2.7019×100	-5.2033×100
3	-1.0667×101	2.4090×101
4	2.0531×101	-5.6056×101
5	-1.3260×101	5.7328×101
j	cjDE	djDE
0	-1.0475×10-2	-2.5399×10-1
1	2.2902×10-1	5.5554×100
2	-1.9609×100	-4.7903×101
3	8.2439×100	2.0330×102
4	-1.7209×101	-4.2527×102
5	1.4687×101	3.5741×102

To the dielectric properties for living epidermis and dermis, we can add the polynomial counterparts from Birgersson [9] for the stratum corneum and adipose tissue for a complete four-layer description of the skin on the volar forearm from 1 kHz to 1 MHz.

Test set

The fitted functions for the dielectric properties were validated with the test set, as illustrated in Fig. 2e-f. We see that the test set displays the same behavior above around 100 kHz as the training sets for the measured phase. Most importantly, the model predictions agree well with the measured impedance at depth setting 3 throughout the frequency range except for the phase above 100 kHz; the phase does, however, stay within or almost within the standard deviations of the measurements.

According to a Lilliefors test at a significance level of 0.05 for each frequency, 100% and 93% of the magnitudes and phases of depth setting 3 are normal-distributed. We therefore take the measurements to be normal-distributed.

Taking paired t-Tests between depth 1 and 3 as well as depth 3 and 5 reveals P≲10⁻5throughout our frequency range, except around 1 MHz where the P-values for the phases become O(10⁻2). Similarly, the paired t-Test gives similarly small P-values (P≲10⁻4for ν≲105Hz) throughout the frequency range for the imaginary and real parts, except again around 1 MHz where the P-values for the imaginary parts become O(10⁻1-10⁻2). Overall, we can consider the test set to be sufficiently different from the training set.

Resistivity

As expected, the mean resistivity is higher for the living epidermis than dermis (Fig. 2c): At 1 kHz, the mean LE resistivity is 110 Ωm whilst the mean DE resistivity is 1.9 Ωm; and at 1 MHz, the mean LE resistivity is 11 Ωm whilst the mean DE resistivity is 1.1 Ωm. Furthermore, the resistivity of the dermis does not vary as much as that of living epidermis.

Why did we expect this outcome? It is well known that a lower resistivity occurs when there is a larger number of mobile ions under the influence of an electric field [30]. In our case, there is a higher water content in the dermis than living epidermis [31], which contributes to a higher number of mobile ions in the presence of the applied electric field.

As an additional point of comparison for the dermis, we note that Tavernier [7] reported constant transversal and planar resistivities of around 0.6 Ωm and 0.4 Ωm respectively for in-vitro measurements between 1 kHz and 100 kHz. Now, while these values are lower than ours, mainly because they measured devitalized, fully saturated dermis tissue in a nutrient liquid rather than the in-vivo counterpart, which is influenced by the circadian rhythm and the persons’ overall water content, their resistivities remain near-to constant as well.

Relative permittivity

The relative permittivity is a proportionality factor between the electric charge and the electric field, and reflects how much the bound charges can be displaced or polarized under the influence of an electric field [32]. From Fig. 2d, our measured in-vivo permittivity for the dermis decreases from around 106 to 102 in the frequency range from 1 kHz to 1 MHz, while that of the living epidermis is relatively constant on the order of 104. This frequency dependence of dielectric properties, known as β-dispersion, of the relative permittivity of dermis in the kHz frequency range occurs due to the polarization of cellular membranes, proteins and other organic macromolecules [33,34].

Returning to the measurements by Tavernier [7] for the dermis, we note that they found the relative permittivity to be around 106 (planar) and 4 × 105 (transverse) at 1 kHz, which is on the same order as ours; at 100 kHz, they found that the relative permittivity had dropped by around one order of magnitude, which is lower than our drop of around two orders of magnitude. Again, this is most likely due to the nature of their measurements.

Validation

We validate the model with the test case comprising the EIS measurement at a depth setting of 3, as illustrated in Fig. 2e-f. Overall, the predicted impedance based on the functions for the dielectric constants in the mathematical model agrees well with the measured counterpart in the frequency range for the magnitude and phase—although less so for the phase above 100 kHz. We attribute the increase in phase for depth settings 1 and 3 that appears after 100 kHz to noise in the measurements when the phase falls below around 10°, because we only expect the β-dispersion in this frequency range [34]. This is further supported by the fact that the mean phase for depth 5 does not suddenly increase at around 1 MHz and that the magnitudes of the impedance for the three depth settings are starting to even out around 1 MHz.

Referring back to the paired t-Tests to see whether the information we used for fitting based on depth settings 1 and 5 is sufficiently different from the training set, depth 3, we note that P≲10⁻4throughout most of the frequency range suggests that this is indeed the case. Overall, the validation should be sufficient for our purposes: viz., fitting of the dielectric properties of living epidermis and dermis. We could have split up the cohort into two subsets for the training and test sets instead of taking depth setting 3 for all subjects as the test set, but would not have been able to use all the available information for the fitting in doing so. Furthermore, depth setting 3 was selected because it gives a reasonably different magnitude, phase and pathway of the current in the stripped skin as compared to the other two depth settings—these are the physical phenomena that depend on the dielectric properties during EIS measurements of human skin.

Limitations

We have modeled a two-layer entity of stripped skin consisting of living epidermis and dermis by assuming homogenous and isotropic dielectric parameters. We have thus not captured any variations inside these two layers nor any anisotropic behavior. It might be possible to add at least the latter by also exploiting the experimental measurements at depth settings 2 and 4, which were not used in this study.

We did not model the subcutaneous layer, because we have shown through a current-distribution analysis [10] that the amount of current passing through the adipose tissue in the subcutaneous layer is negligible for our frequencies. Any current passing through the muscle tissue should therefore also be negligible.

We also did not consider additional factors such as the skin-electrode-interface impedance, hair follicles and sweat ducts since they are known to have a minor effect in the 1 kHz to 1 MHz frequency range with soaking [35, 36, 37].

From an experimental point-of-view, we note that there is scope to improve the measurements by fine-tuning the electrical-impedance spectrometer to be less sensitive to noise as we approach 1 MHz, which should result in smoother fitted dielectric properties even before fitting polynomials. Furthermore, one could increase the number of subjects as well as try to measure the stripping of stratum corneum.

Multi-scale extensions

The current mathematical model solves for a volume-averaged potential and total current that currently does not resolve the details of the cells or other structural features in the skin layers; therefore, the fitted dielectric parameters are also averages. In order to capture smaller-scale features, two extensions of the current work could be pursued.

The first extension would be to solve the three-dimensional transport of ohmic and displacement currents in a computational domain that comprises the relevant skin components—such as, e.g., cell membranes, cytoplasm and blood capillaries—numerically with the complex-valued PDE in a given skin layer whilst treating the other skin layers as averages to keep the computational cost reasonable and simulations tractable. We could then compare the simulated impedance with the experimental findings. This approach, however, would not lend itself easily to fitting of unknown dielectric parameters and their functional form with regards to the applied frequency due to the computational cost of geometrically resolving multiple length scales from the skin thickness down to the thickness of cell membranes, which brings us to the second extension.

In the second extension, we would apply the theory of volume-averaging [38] through a paper-and-pen approach to derive effective dielectric parameters for intrinsic potentials and superficial fluxes in the intra- and extracellular phases and the cell membranes that should result in a volume-averaged description that captures the leading order physics at multiple length scales yet is still computational cheap to solve. The key feature here is that multiple length scales will manifest themselves inside effective dielectric parameters and the interaction between the phases mathematically in the resulting PDE(s), and can therefore be solved in simple layers that we consider here. In fact, the set of equations we have solved in this paper represent one outcome from volume-averaging that lumps up all phases into one effective phase. Furthermore, this paper-and-pen analysis should allow us to write the dielectric properties not only as fifth-order polynomials but also as functions of parameters that describe the skin layer content, geometries and material properties: e.g., volume fractions of cells, dielectric properties of cell membranes, cytoplasm, and blood (depending on the skin layer). We have demonstrated both the suggested extensions in other engineering applications such as, e.g., organic [39,40] and perovskite [41] solar cells; and refer the interested reader to these for more ideas and a structured way to derive the volume-averaged equations for multi-scale systems that involve charge transport.

Conclusions

We have obtained new resistivity and permittivity values for the living epidermis and dermis of the human volar forearm, which are expressed as easy-to-use functions for frequencies between 1 kHz and 1 MHz; as well as weighted means and pooled standard deviations of epidermis, stratum corneum, living epidermis and dermis from a literature review. These thicknesses were employed in the mathematical model for conservation of charge that was fitted to the impedance measurements of stripped skin. In addition, we have discussed possible extensions to resolve more structural features of each skin layer.