Visual Abstract
Abstract
In oncologic PET, the SUV and standardized uptake ratio (SUR) of a viable tumor generally increase during the postinjection period. In contrast, the net influx rate (Ki), which is derived from dynamic PET data, should remain relatively constant. Uptake-time-corrected SUV (cSUV) and SUR (cSUR) have been proposed as uptake-time-independent, static alternatives to Ki. Our primary aim was to quantify the intrascan repeatability of Ki, SUV, cSUV, SUR, and cSUR among malignant lesions on PET/CT. An exploratory aim was to assess the ability of cSUR to estimate Ki. Methods: This prospective, single-center study enrolled adults undergoing standard-of-care oncologic PET/CT. SUV and Ki images were reconstructed from dynamic PET data obtained before (∼35–50 min after injection) and after (∼75–90 min after injection) standard-of-care imaging. Tumors were manually segmented. Quantitative metrics were extracted. cSUVs and cSURs were calculated for a 60-min postinjection reference uptake time. The magnitude of the intrascan test–retest percent change (test–retest |%Δ|) was calculated. Coefficients of determination (R2) and intraclass correlation coefficients (ICC) were also computed. Differences between metrics were assessed via the Wilcoxon signed-rank test (α, 0.05). Results: This study enrolled 78 subjects; 41 subjects (mean age, 63.8 y; 24 men) with 116 lesions were analyzed. For both tracers, SUVmax and maximum SUR (SURmax) had large early-to-late increases (i.e., poor intrascan repeatability). Among [18F]FDG-avid lesions (n = 93), there were no differences in intrascan repeatability (median test–retest |%Δ|; ICC) between the maximum Ki (Ki,max) (13%; 0.97) and either the maximum cSUV (cSUVmax) (12%, P = 0.90; 0.96) or the maximum cSUR (cSURmax) (13%, P = 0.67; 0.94). For DOTATATE-avid lesions (n = 23), there were no differences in intrascan repeatability between the Ki,max (11%; 0.98) and either the cSUVmax (13%, P = 0.41; 0.98) or the cSURmax (11%, P = 0.08; 0.94). The SUVmax, cSUVmax, SURmax, and cSURmax were all strongly correlated with the Ki,max for both [18F]FDG (R2, 0.81–0.92) and DOTATATE (R2, 0.88–0.96), but the cSURmax provided the best agreement with the Ki,max across early-to-late time points for [18F]FDG (ICC, 0.69–0.75) and DOTATATE (ICC, 0.90–0.91). Conclusion: Ki,max, cSUVmax, and cSURmax had low uptake time dependence compared with SUVmax and SURmax. The Ki,max can be predicted from cSURmax.
In oncologic imaging, SUV changes between scans are critical for treatment response assessment (1). However, SUV depends on uptake times, as many tumors accumulate tracer continuously (2,3). The logistic demands of busy clinical PET services often preclude precise scan-to-scan reproduction of uptake times, reducing the reliability of SUV as an oncologic biomarker (4). The Patlak model, which attempts to ameliorate this shortcoming, assumes that circulating tracer is trapped irreversibly, allowing tracer uptake to be quantified via the net influx rate (Ki) (5,6). Several clinically used PET tracers, including [18F]FDG, approximate this behavior, permitting Patlak modeling of dynamic PET data. Once steady-state conditions are achieved between the blood and tissue compartments, the Ki should remain relatively constant, whereas the SUV is expected to increase with time. Furthermore, Ki-based metrics are promising prognostic biomarkers for several cancer types, occasionally outperforming SUV-based metrics (7,8).
However, Ki derivation requires direct measurement or estimation of arterial input functions (AIFs) and dynamic acquisitions to generate tissue time–activity curves (5). The required modifications to PET protocols may increase imaging time or introduce motion-related quantitative errors (9). The uptake-time-corrected SUV (cSUV), which involves retrospectively modifying an observed SUV on the basis of actual versus targeted uptake times, is an alternative means of addressing the uptake time dependence of SUVs (10). Furthermore, the uptake-time-corrected tumor-to-blood standardized uptake ratio (cSUR) may allow for Ki estimation, without the need for AIFs or dynamic imaging (11).
To our knowledge, no prior studies have assessed the relative temporal stabilities of SUV, SUR, Ki, cSUV, and cSUR in a broad oncologic PET population. Thus, our primary aim was to quantify the intrascan repeatability of these metrics and thereby determine which approach provides the most time-independent assessment of tracer avidity on [18F]FDG PET and DOTATATE PET. An exploratory aim was to determine the ability of cSUR to estimate Ki.
MATERIALS AND METHODS
Study Design
This prospective, Institutional Review Board–approved, Health Insurance Portability and Accountability Act–compliant study (NCT04283552) enrolled 78 subjects from a pool of consecutive patients scheduled to undergo standard-of-care (SOC) oncologic PET/CT for various indications, using [18F]FDG, [68Ga]Ga-DOTATATE, [64Cu]Cu-DOTATATE, or [18F]piflufolastat (note that [68Ga]Ga-DOTATATE and [64Cu]Cu-DOTATATE are hereafter collectively called DOTATATE, as these scans were analyzed together). These tracers have been reported to satisfy the Patlak model’s assumptions (12–14). The 2 [18F]piflufolastat studies were excluded because of insufficient cases for tracer-specific analysis. All imaging occurred at a tertiary-care center between June 2020 and October 2022. Inclusion criteria included being at least 18 y of age, having the ability to provide written informed consent, and self-reporting the ability to tolerate approximately 90 min of near-motionless supine positioning. Study imaging was performed before and after SOC PET/CT using the same tracer dose.
Imaging Protocol
The study imaging protocol (details are available in the supplemental materials; available at http://jnm.snmjournals.org) is summarized in Figure 1. All patients were imaged on a single Biograph Vision 600 PET/CT scanner (Siemens Healthineers) equipped with commercially available software for direct reconstruction of multiparametric PET images (FlowMotion Multiparametric PET Suite; Siemens Healthineers).
PET/CT acquisition protocol. Brackets show timing of data used for early and late image reconstructions.
PET Image Reconstruction
Using automated scanner tools, volumes of interest were placed in the descending thoracic aorta on a 6-min dynamic chest acquisition and the subsequent 10 whole-body (WB) passes (15). Per default scanner software settings, AIFs were generated from measured blood activity concentrations via exponential ([18F]FDG) or linear piecewise (DOTATATE) curve fitting. After all WB PET passes were reviewed dynamically for large bulk motion events, early and late SUV and Ki images were reconstructed per manufacturer-recommended parameters (Supplemental Table 1). Each reconstruction used data from three 5-min WB passes with targeted acquisition times of 35–50 min (early) and 75–90 min (late) after injection. Note that the scanner software requires at least 3 WB passes for Patlak analysis. The 3 latest pre-SOC WB passes were selected for the early images, ensuring adequate time to achieve steady-state conditions. Importantly, subjects left the scanner to void immediately before SOC imaging per our standard clinical protocol, precluding automated scanner measurement of post-SOC blood tracer concentrations due to different patient positioning. Consequently, the AIF for the post-SOC Ki reconstructions was automatically derived by the scanner software from extrapolation of the pre-SOC AIF (i.e., no incorporation of measured post-SOC blood tracer concentrations). SUV was based on actual body weight with units of grams per milliliter. Ki had units of milliliter per minute per 100 mL. In contrast to [18F]FDG, intravascular DOTATATE does not enter red blood cells, requiring correction of measured Ki values for the subjects’ hematocrit levels (16):
Quantitative Analysis
Tracer-avid lesions deemed to represent sites of viable malignancy on the SOC PET/CT interpretation were selected by one author. In cases of numerous lesions, the largest or most tracer-avid lesions were selected (5 per subject maximum). Each lesion was manually segmented in MIM version 7.1.5 (MIM Software) on 4 PET image sets (Ki-early, SUV-early, Ki-late, SUV-late) to generate volumes of interest, using coregistered CT images for guidance. Maximum and peak values were extracted. Additionally, a cylindric volume of interest (1-cm diameter, 6-cm length) was placed in the descending thoracic aorta (avoiding vessel walls) to extract a mean value for SUR calculation: SUVmax and SUVpeak were used to calculate maximum SUR (SURmax) and peak SUR (SURpeak), respectively; the SUVmean of blood was used in both cases.
Uptake Time Correction
Actual uptake time ranges were extracted for each image set, with the mid point defining the effective uptake time (e.g., 44.5 min for 37–52 min after injection). cSUV and cSUR were calculated as follows (10,11):
SUV and SUR are measured values, T0 is the actual uptake time, and Tc is the correction time reference. Tc was set to 60 min, reflecting a commonly targeted uptake time in [18F]FDG and DOTATATE protocols (1). Figure 2 shows the SUV and SUR correction procedure for a representative case. For [18F]FDG, the value of the parameter b of 0.313 was based on a prior study (10). For DOTATATE, we empirically derived a b of 0.63 by determining the value (averaged across all subjects) that best reproduced the observed late SUVmax from the observed early SUVmax. The early and late values were then corrected to 60 min with the cSUV equation.
Uptake time correction procedure for SUV and SUR. Axial [18F]FDG PET and fused [18F]FDG PET/CT images are shown for Ki (top) and SUV (bottom) reconstructions at early and late time points. Ki,max, SUVmax, and SURmax of [18F]FDG-avid mediastinal lymph node (arrows) are shown. Test–retest %Δ values were 1.5%, 49.8%, and 78.6% for Ki,max, SUVmax, and SURmax, respectively, indicating much better intrascan repeatability for Ki,max. Procedure for correcting SUVmax and SURmax to 60 min after injection is also shown. Test–retest %Δ values were −1.4% for cSUVmax and +7.0% for cSURmax, similar to Ki,max results.
Manual Patlak Analysis
To explore apparent temporal variations in Ki, we selected 6 [18F]FDG, 1 [68Ga]Ga-DOTATATE, and 2 [64Cu]Cu-DOTATATE cases with at least 1 lesion exhibiting a large (>20% or <20%) test–retest percent change (%Δ) in maximum Ki (Ki,max) for further analysis. For all 9 cases, extrapolated AIF curve fits were compared with manually measured blood activity concentrations on the WB passes. Areas under the time–activity curve were compared for extrapolated AIF curve fits versus manual measurements via trapezoidal integration. For 4 cases, full manual Patlak analysis was performed for selected lesions and reference organs (supplemental materials) (17–20).
Statistical Analysis
Statistical analysis was conducted in Prism 9 (GraphPad) and Excel 2016 (Microsoft) by one author with statistician guidance. Participant and scan characteristics were summarized descriptively. Because of the anticipated pharmacokinetic differences, the [18F]FDG cases were analyzed separately from the DOTATATE cases. The 2-tailed Wilcoxon signed-rank test was used for pairwise comparisons of quantitative metrics. Intrascan test–retest changes were computed:
To facilitate comparisons across metrics of different magnitudes, intrascan test–retest %Δ was also computed:
Results were displayed via Bland–Altman plots and box-and-whisker plots (21). The mean (μ) and SD (σ) of the test–retest Δ and test–retest %Δ distributions were determined for each metric. The 95% limits of repeatability were defined as follows:
Given that near-zero %Δ values could be due to averaging of large negative and positive changes, absolute test–retest %Δ (test–retest |%Δ|) values were also computed. Intraclass correlation coefficients (ICC) and coefficients of determination (R2) were also used to assess test–retest repeatability and to quantify the accuracy of Ki prediction by other metrics. A P value of less than 0.05 defined statistical significance. More detailed statistical methods are available in the supplemental materials (22).
RESULTS
Study Cohort
Of the 78 study subjects, 41 subjects with 116 lesions (93 on [18F]FDG; 23 on DOTATATE) were included in the analysis (Fig. 3). This study cohort was 58.5% men (24/41) with a mean age of 63.8 y. Additional patient and scan characteristics are captured in Supplemental Table 2.
Study flowchart.
Intrascan Repeatability of Tumor Uptake Metrics
Test–retest repeatability results are summarized in Table 1 ([18F]FDG) and Supplemental Table 3 (DOTATATE). Scatterplots of late versus early metric values are shown in Supplemental Figures 1 and 2 ([18F]FDG) and Supplemental Figures 3 and 4 (DOTATATE). Bland–Altman plots for Ki, cSUV, and cSUR are displayed in Figure 4 ([18F]FDG) and Supplemental Figure 5 (DOTATATE). Table 2 ([18F]FDG) and Supplemental Table 4 (DOTATATE) summarize test–retest |%Δ| values for all metrics. Box-and-whisker plots of test–retest %Δ and test–retest |%Δ| distributions are presented in Figure 5 ([18F]FDG) and Supplemental Figure 6 (DOTATATE). The remainder of this section focuses on the [18F]FDG results, except as noted (detailed DOTATATE results are provided in the supplemental materials).
Intrascan Repeatability of SUV, cSUV, SUR, cSUR, and Ki Metrics Among [18F]FDG-Avid Lesions
Intrascan test–retest repeatability for [18F]FDG-avid lesions. Bland–Altman plots are shown for Ki,max (A), cSUVmax (B), cSURmax (C), Ki,peak (D), cSUVpeak (E), and cSURpeak (F). Horizontal orange lines represent test–retest %Δ means. Horizontal dashed and solid black lines represent ±1σ and ±2σ for test–retest %Δ distributions, respectively. Each open circle represents [18F]FDG-avid lesion. 95% limits of repeatability (LOR) are computed as −2σ, +2σ.
Test–retest Absolute %Δ for SUV, cSUV, SUR, cSUR, and Ki Metrics Among [18F]FDG-Avid Lesions
Test–retest %Δ and |%Δ| distributions for [18F]FDG-avid lesions. Box-and-whisker plots show test–retest %Δ (A and B) and |%Δ| (C and D) distributions for maximum (A and C) and peak (B and D) values of Ki, SUV, cSUV, SUR, and cSUR. All P values are based on comparison to Ki,max or Ki,peak. Table 2 provides descriptive statistics.
As expected, SUV and SUR metrics showed large, statistically significant early-to-late increases. In contrast, cSUV and cSUR metrics were similar at early and late time points, though with some small but significant early-to-late changes. For the maximum cSUV (cSUVmax), the early and late values were statistically equivalent (median, 8.0 vs. 7.3; P = 0.17). Surprisingly, the Ki metrics exhibited significant early-to-late increases (median, 1.8 vs. 2.3; P < 0.001). The early and late values of each metric were strongly correlated (R2, 0.90–0.96). However, the ICCs showed substantially better agreement between early and late values for Ki, cSUV, and cSUR metrics (range, 0.91–0.97) than for SUV and SUR metrics (range, 0.27–0.75). For DOTATATE, the results were similar, except that the Ki metrics exhibited significant (or nearly significant) early-to-late decreases.
In the Bland–Altman analysis, cSUVmax and maximum cSUR (cSURmax) showed the least bias between early and late values, with mean test–retest %Δ values of −6% and 7%, respectively, compared with 11% for Ki,max. In contrast, the mean test–retest %Δ values for SUVmax and SURmax were 47% and 81%, respectively, indicating large early-to-late increases. Regarding the magnitude of deviation from perfect repeatability (i.e., test–retest |%Δ| = 0), the test–retest |%Δ| of Ki,max (median, 13%) was similar to those of cSUVmax (median, 12%; P = 0.90) and cSURmax (median, 13%; P = 0.67) but significantly less than those of SUVmax (median, 48%; P < 0.001) and SURmax (median, 81%; P < 0.001). The test–retest |%Δ| of the peak Ki (Ki,peak) (median, 15%) was significantly lower than that of all other relevant metrics except for the peak cSUR (cSURpeak) (median, 13%; P = 0.36). For DOTATATE, the results were similar to those of [18F]FDG for the Ki,max analysis, though the median test–retest |%Δ| of the Ki,peak was similar to that of SUVpeak (rather than cSURpeak).
Prediction of Ki by cSUR
Supplemental Figures 7 and 8 ([18F]FDG) and Supplemental Figures 9 and 10 (DOTATATE) show correlation results for Ki,max versus SUVmax, SURmax, cSUVmax, and cSURmax. For [18F]FDG, the maximum values of all 4 metrics strongly predicted the corresponding early (R2, 0.81–0.92) and late (R2, 0.81–0.91) Ki,max values. However, agreement with Ki,max by ICCs was the highest for SURmax (0.94) and cSURmax (0.75) at the early time point and for cSURmax (0.69) at the late time point. Similarly, for DOTATATE, the maximum values of all 4 metrics strongly predicted the corresponding early (R2, 0.88–0.96) and late (R2, 0.88–0.93) Ki,max values. However, in contrast to [18F]FDG, agreement with Ki,max by ICCs was high for all 4 metrics at the early time point (range, 0.85–0.93) but only for cSUVmax (0.78) and cSURmax (0.90) at the late time point.
Manual Patlak Analysis
Supplemental Figures 11 and 12 show Ki,max test–retest %Δ values for each subject’s lesions for [18F]FDG and DOTATATE, respectively. Supplemental Tables 5 and 6 present manual Patlak analyses for several [18F]FDG and DOTATATE subjects, respectively. Supplemental Figures 13 and 14 capture AIF and tissue-response curves and manual Patlak plots for representative [18F]FDG and DOTATATE cases, respectively. Supplemental Figure 15 illustrates the effects of motion and image noise on Ki and SUV. For [18F]FDG, the AIF curve-fit extrapolation mildly underestimated late blood activity concentrations, contributing to higher late Ki values. Furthermore, motion of small lesions across WB passes contributed to Ki errors that were ameliorated by manual frame-by-frame segmentations. For DOTATATE, the AIF curve-fit extrapolation moderately overestimated late blood activity concentrations, resulting in lower late Ki values; additionally, DOTATATE binding appeared to be reversible at late time points for some cases. More details are provided in the supplemental materials.
DISCUSSION
In this study, we quantified the intrascan repeatability of Ki, SUV, cSUV, SUR, and cSUR metrics for [18F]FDG-avid and DOTATATE-avid lesions on WB PET in a general oncology population. For both tracers, SUVmax and SURmax showed large early-to-late increases (i.e., poor intrascan repeatability). For [18F]FDG, there were no significant differences in intrascan repeatability between Ki,max (median test–retest |%Δ|, 13%; ICC, 0.97) and either cSUVmax (median test–retest |%Δ|, 12%; P = 0.90; ICC, 0.96) or cSURmax (median test–retest |%Δ|, 13%; P = 0.67; ICC, 0.94). The intrascan repeatability of the Ki,peak was better than that of the peak cSUV (cSUVpeak) but similar to that of cSURpeak. For DOTATATE, there were no significant differences in intrascan repeatability between Ki,max (median test–retest |%Δ|, 11%; ICC, 0.98) and either cSUVmax (median test–retest |%Δ|, 13%; P = 0.41; ICC, 0.98) or cSURmax (median test–retest |%Δ|, 11%; P = 0.08; ICC, 0.94). Again, intrascan repeatability of the Ki,peak was better than that of cSUVpeak or cSURpeak.
Early-to-late increases are a well-known limitation of SUVs and SURs for tumor response assessments (2,3,23). As such, the Quantitative Imaging Biomarkers Alliance recommends that uptake times for baseline and follow-up scans be approximately 60 min with a no more than 10 min difference between scans (1). However, differences greater than 10 min are not uncommon. Methods to correct SUV and SUR for uptake time (i.e., cSUV, cSUR) have been published (10,23). For example, a study reported that correcting SUVs and SURs from 20 min to 55 min after injection reduced differences with actual values at 55 min from −30% to 2% for SUV and from −52% to −3% for SUR (10). This study, which used data from 9 male patients with colorectal liver metastases, proposed the simple SUV and SUR correction equations used in our work.
We verified that cSUV, using the published time parameter b of 0.313, is a relatively time-independent marker of tumoral [18F]FDG avidity. For [18F]FDG, our mean test–retest %Δ values of −6% and 7% for cSUVmax and cSURmax, respectively, are slightly greater in magnitude than the values cited above, possibly because of our longer early-to-late intervals or heterogeneous patient cohort. We empirically derived a b value of 0.63 for DOTATATE and found that cSUV is also a relatively time-independent marker of tumoral DOTATATE avidity, with mean test–retest %Δ values of 2% and −7% for cSUVmax and cSURmax, respectively. Compared with cSUVmax and cSURmax, cSUVpeak and cSURpeak showed worse intrascan repeatability, with sizeable negative test–retest %Δ values for both tracers. The reason for this somewhat surprising finding is unclear, as peak measurements (because of their larger sampling volumes and lower potential for noise-related errors) are generally considered more repeatable than maximum measurements (24).
In terms of test–retest |%Δ| and ICC, the intrascan repeatability was similar across Ki,max, cSUVmax, and cSURmax for both tracers. However, we observed small but statistically significant early-to-late increases and decreases in Ki,max for [18F]FDG and DOTATATE cases, respectively. In contrast, cSUVmax and cSURmax showed no significant early-to-late changes for either tracer, with the exception of a small significant increase in cSURmax for DOTATATE. For both tracers, the observed early-to-late Ki changes were partially attributable to inaccurate AIF curve-fit extrapolations, the need for which arose from incorporating SOC imaging into our study design. A protocol using nonextrapolated image-derived AIFs or population-based AIFs might reduce these apparent temporal changes in Ki. Several cases suggested late reversibility of DOTATATE binding, also contributing to the observed early-to-late Ki decreases. Overall, cSUVmax and cSURmax provided intrascan repeatability similar to that of Ki,max, without dynamic imaging or AIF estimation.
Ki images may still be worth their inherent complexities, as Ki metrics appear useful for guiding treatment decisions and predicting oncologic outcomes (7,8,25,26). One study showed that Ki correlated with SUR (R2, 0.96) much more strongly than with SUV (R2, 0.37), with all metrics measured at 50–60 min after injection (11). In contrast, we found that SUVmax, cSUVmax, SURmax, and cSURmax all strongly correlated with Ki,max for both [18F]FDG (R2, 0.81–0.92) and DOTATATE (R2, 0.88–0.96), though cSURmax had the best agreement with Ki,max across early and late time points for [18F]FDG (ICC, 0.69–0.75) and DOTATATE (ICC, 0.90–0.91). Our findings indicate that Ki can be predicted from cSUR and that cSURmax exhibits a nearly 1:1 proportionality to Ki,max. To this point, cSUR and Ki appear to predict postchemoradiation lung cancer outcomes better than does SUV (27). That said, Patlak images may provide higher lesion conspicuity and fewer false positives than with SUV images (28,29).
Our study has limitations, including its single-center, single-scanner design. The results should be corroborated at other centers on other scanners. Our patient cohort was heterogeneous; the relatively small sample size precluded subgroup analysis by cancer type or imaging indication. The b parameter of 0.63 for DOTATATE was derived empirically (rather than from AIF curve fitting) and needs to be validated in other cohorts. Again, the AIF curve-fit extrapolations created late Ki errors. A more thorough investigation of potential causes of the observed temporal variability in Ki is still warranted. Finally, our study excluded subjects who anticipated difficulty with a 90-min imaging period, potentially enriching our cohort for patients capable of remaining relatively motionless; as such, Ki images may be more degraded by motion in an unselected oncologic population.
CONCLUSION
Ki,max, cSUVmax, and cSURmax exhibit comparably high intrascan repeatability in a general oncologic population undergoing PET with [18F]FDG or DOTATATE, with significantly less uptake time dependence compared with SUVmax and SURmax. cSURmax can predict Ki,max without dynamic acquisitions.
DISCLOSURE
This work was supported by a research grant from Siemens Healthineers to Washington University, including salary support for Tyler Fraum. Richard Wahl has received consulting income from Siemens Healthineers. All participants were imaged on a Siemens PET/CT scanner. Saeed Ashrafinia and Anne Smith are Siemens employees. These authors participated in the initial study design, provided occasional technical support, and critically reviewed the manuscript. However, all data collection, analysis, and manuscript preparation were performed by Washington University authors. No other potential conflict of interest relevant to this article was reported.
KEY POINTS
QUESTION: In oncologic PET, which quantitative metric of tracer avidity is least dependent on uptake times?
PERTINENT FINDINGS: In this prospective, cross-sectional study of general oncology subjects undergoing PET with [18F]FDG or DOTATATE, Ki,max, cSUVmax, and cSURmax exhibited comparably high intrascan repeatability, whereas SUVmax and SURmax substantially increased at later uptake times.
IMPLICATIONS FOR PATIENT CARE: Tracer-avidity metrics with high intrascan repeatability should be used to avoid quantitative interpretive errors related to different uptake times across scans.
Footnotes
Published online Aug. 14, 2024.
- © 2024 by the Society of Nuclear Medicine and Molecular Imaging.
Immediate Open Access: Creative Commons Attribution 4.0 International License (CC BY) allows users to share and adapt with attribution, excluding materials credited to previous publications. License: https://creativecommons.org/licenses/by/4.0/. Details: https://jnm.snmjournals.org/page/permissions.
REFERENCES
- Received for publication July 28, 2023.
- Accepted for publication July 10, 2024.