[go: up one dir, main page]

 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (2,554)

Search Parameters:
Keywords = data validation and calibration

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
15 pages, 1463 KiB  
Article
Integration of FTIR Spectroscopy, Volatile Compound Profiling, and Chemometric Techniques for Advanced Geographical and Varietal Analysis of Moroccan Eucalyptus Essential Oils
by Aimen El Orche, Abdennacer El Mrabet, Amal Ait Haj Said, Soumaya Mousannif, Omar Elhamdaoui, Siddique Akber Ansari, Hamad M. Alkahtani, Shoeb Anwar Ansari, Ibrahim Sbai El Otmani and Mustapha Bouatia
Sensors 2024, 24(22), 7337; https://doi.org/10.3390/s24227337 (registering DOI) - 17 Nov 2024
Viewed by 256
Abstract
Eucalyptus essential oil is widely valued for its therapeutic properties and extensive commercial applications, with its chemical composition significantly influenced by species variety, geographical origin, and environmental conditions. This study aims to develop a reliable method for identifying the geographical origin and variety [...] Read more.
Eucalyptus essential oil is widely valued for its therapeutic properties and extensive commercial applications, with its chemical composition significantly influenced by species variety, geographical origin, and environmental conditions. This study aims to develop a reliable method for identifying the geographical origin and variety of eucalyptus oil samples through the application of advanced analytical techniques combined with chemometric methods. Essential oils from Eucalyptus globulus and Eucalyptus camaldulensis were analyzed using Gas Chromatography–Flame Ionization Detection (GC–FID) and Fourier Transform Infrared (FTIR) Spectroscopy. Chemometric analyses, including Orthogonal Partial Least Squares-Discriminant Analysis (O2PLS-DA) and Hierarchical Cluster Analysis (HCA), were utilized to classify the oils based on their volatile compound profiles. Notably, O2PLS-DA was applied directly to the raw FTIR data without additional spectral processing, showcasing its robustness in handling unprocessed data. For geographical origin determination, the GC–FID model achieved a Correct Classification Rate (CCR) of 100%, with 100% specificity and 100% sensitivity for both calibration and validation sets. FTIR spectroscopy achieved a CCR of 100%, specificity of 100%, and sensitivity of 100% for the calibration set, while the validation set yielded a CCR of 95.83%, specificity of 99.02%, and sensitivity of 94.44%. In contrast, the analysis based on species variety demonstrated 100% accuracy across all metrics CCR, specificity, and sensitivity—for both calibration and validation using both techniques. These findings underscore the effectiveness of volatile and infrared spectroscopy profiling for quality control and authentication, providing robust tools for ensuring the consistency and reliability of eucalyptus essential oils in various industrial and therapeutic applications. Full article
Show Figures

Figure 1

Figure 1
<p>The geographical positions of the samples.</p>
Full article ">Figure 2
<p>Key Constituents of Eucalyptus Leaf Essential Oils.</p>
Full article ">Figure 3
<p>Comprehensive O2PLS-DA and HCA Analysis of Essential Oils: (<b>a</b>) Varietal Differentiation Score Plot, (<b>b</b>) Geographical Origin Score Plot, (<b>c</b>) Key Volatile Compound Loading Plot, and (<b>d</b>) HCA Dendrogram for Regional Classification (KE: Kenitra, TA: Taounate, TI: Tiflet, CA: Casablanca, KH: Khemisset, ME: Marrakech).</p>
Full article ">Figure 4
<p>FTIR Spectral Analysis of Essential Oils: Comparison of Spectral Profiles for Camaldulensis and Globulus Varieties (R1: Kenitra, R2: Taounate, R3: Tiflet for Camaldulensis; and R4: Casablanca, R5: Khemisset, R6: Marrakech for Globulus).</p>
Full article ">Figure 5
<p>Classification and Clustering of Eucalyptus Essential Oils: O2PLS-DA and HCA Analysis. (<b>a</b>): Classification of Eucalyptus Essential Oils by Variety: LV1 vs. LV2 Plot, (<b>b</b>): Geographical Origin Classification of Eucalyptus Essential Oils: LV1 vs. LV2 Plot, (<b>c</b>): Geographical Origin Classification of Eucalyptus Essential Oils: LV1 vs. LV3 Plot, and (<b>d</b>): Hierarchical Cluster Analysis (HCA) of Eucalyptus Essential Oils (KE: Kenitra, TA: Taounate, TI: Tiflet, CA: Casablanca, KH: Khemisset, ME: Marrakech).</p>
Full article ">
16 pages, 2161 KiB  
Article
Validation of EuroSCORE II, ACEF Score, CHA2DS2-VASc, and CHA2DS2-VA in Patients Undergoing Left Main Coronary Artery Angioplasty: Analysis from All-Comers BIA-LM Registry
by Emil Julian Dąbrowski, Paweł Kralisz, Konrad Nowak, Kamil Gugała, Przemysław Prokopczuk, Grzegorz Mężyński, Michał Święczkowski, Sławomir Dobrzycki and Marcin Kożuch
J. Clin. Med. 2024, 13(22), 6907; https://doi.org/10.3390/jcm13226907 (registering DOI) - 16 Nov 2024
Viewed by 201
Abstract
Background: Simple surgical and clinical risk scores are useful in mortality prediction. Aims: The study’s aim was to validate three scores in real-world registry of percutaneous coronary intervention (PCI) for the left main coronary artery (LMCA). Methods: All data were [...] Read more.
Background: Simple surgical and clinical risk scores are useful in mortality prediction. Aims: The study’s aim was to validate three scores in real-world registry of percutaneous coronary intervention (PCI) for the left main coronary artery (LMCA). Methods: All data were obtained from the BIA-LM Registry. Discrimination and calibration of EuroSCORE II, ACEF, CHA2DS2-VASc, and CHA2DS2-VA were assessed with receiver operating characteristic (ROC) curves analysis and Hosmer–Lemeshow (HL) test. Results: The final cohort included 851 patients, median age was 71, and 156 patients had history of previous coronary artery bypass grafting (CABG). Median EuroSCORE II, ACEF, CHA2DS2-VASc, and CHA2DS2-VA were 3.1% (IQR 5.4%), 1.56 (IQR 0.9), 4 (IQR 2), and 4 (IQR 2), respectively. In the short- (30 days) and long-term (mean 4.1 years), there were 27 and 318 deaths. In short-term, EuroSCORE II showed the best discrimination in the overall population and subgroup with unprotected LMCA [area under the curve (AUC) 0.804, 95% CI 0.717–0.890 and AUC 0.826, 95% CI 0.737–0.913, respectively, p < 0.001 for comparisons with other models), with the best cut-off value at 7.1%. In long-term observation, EuroSCORE II and ACEF showed good predictive value (overall population: AUC 0.716, 95% CI 0.680–0.750 and AUC 0.725, 95% CI 0.690–760, respectively). In short- and long-term observation, EuroSCORE II and ACEF showed poor calibration (HL test p < 0.05) as compared to CHA2DS2-VASc (HL test p = 0.40 and 0.18). Conclusions: EuroSCORE II showed good mortality prediction in short-term observation; however, its predicted risk should be interpreted with caution due to poor calibration. ACEF and EuroSCORE II may be useful in long-term mortality prediction. Full article
Show Figures

Figure 1

Figure 1
<p>Flow chart of study design. CABG, coronary artery bypass grafting; LMCA, left main coronary artery disease; PCI, percutaneous coronary intervention.</p>
Full article ">Figure 2
<p>The receiver operating characteristic curve of EuroSCORE II, CHA<sub>2</sub>DS<sub>2</sub>-VASc, and ACEF score as a predictor of death. Abbreviations: AUC, area under curve; LMCA, left main coronary artery; ROC, receiver operating characteristic.</p>
Full article ">Figure 3
<p>Comparison of CHA<sub>2</sub>DS<sub>2</sub>-VASc with CHA<sub>2</sub>DS<sub>2</sub>-VA receiver operating characteristic curves. Abbreviations: AUC, area under curve; LMCA, left main coronary artery; ROC, receiver operating characteristic.</p>
Full article ">Figure 4
<p>Calibration plots of the analyzed models in the overall population undergoing left main coronary artery PCI in short- and long-term observation.</p>
Full article ">
18 pages, 3865 KiB  
Article
Rapid Authentication of Intact Stingless Bee Honey (SBH) by Portable LED-Based Fluorescence Spectroscopy and Chemometrics
by Diding Suhandy, Dimas Firmanda Al Riza, Meinilwita Yulia, Kusumiyati Kusumiyati, Mareli Telaumbanua and Hirotaka Naito
Foods 2024, 13(22), 3648; https://doi.org/10.3390/foods13223648 (registering DOI) - 16 Nov 2024
Viewed by 433
Abstract
Indonesian stingless bee honey (SBH) of Geniotrigona thoracica is popular and traded at an expensive price. Brown rice syrup (RS) is frequently used as a cheap adulterant for an economically motivated adulteration (EMA) in SBH. In this study, authentic Indonesian Geniotrigona thoracica SBH [...] Read more.
Indonesian stingless bee honey (SBH) of Geniotrigona thoracica is popular and traded at an expensive price. Brown rice syrup (RS) is frequently used as a cheap adulterant for an economically motivated adulteration (EMA) in SBH. In this study, authentic Indonesian Geniotrigona thoracica SBH of Acacia mangium (n = 100), adulterated SBH (n = 120), fake SBH (n = 100), and RS (n = 200) were prepared. In short, 2 mL of each sample was dropped directly into an innovative sample holder without any sample preparation including no dilution. Fluorescence intensity was acquired using a fluorescence spectrometer. This portable instrument is equipped with a 365 nm LED lamp as the fixed excitation source. Principal component analysis (PCA) was calculated for the smoothed spectral data. The results showed that the authentic SBH and non-SBH (adulterated SBH, fake SBH, and RS) samples could be well separated using the smoothed spectral data. The cumulative percentage variance of the first two PCs, 98.4749% and 98.4425%, was obtained for calibration and validation, respectively. The highest prediction accuracy was 99.5% and was obtained using principal component analysis–linear discriminant analysis (PCA-LDA). The best partial least square (PLS) calibration was obtained using the combined interval with R2cal = 0.898 and R2val = 0.874 for calibration and validation, respectively. In the prediction, the developed model could predict the adulteration level in the adulterated honey samples with an acceptable ratio of prediction to deviation (RPD) = 2.282, and range error ratio (RER) = 6.612. Full article
(This article belongs to the Section Food Analytical Methods)
Show Figures

Figure 1

Figure 1
<p>Visual information of authentic SBH, fake SBH, and brown rice syrup (RS).</p>
Full article ">Figure 2
<p>Visual information of adulterated SBH with six different adulteration levels.</p>
Full article ">Figure 3
<p>The front-face mode spectral acquisition system with portable LED-based fluorescence spectroscopy equipped with an innovative sample holder.</p>
Full article ">Figure 4
<p>Typical fluorescence spectral data in the range of 348.5–866.5 nm with a fixed excitation at 365 nm: (<b>a</b>) raw spectral data; (<b>b</b>) smoothed spectral data.</p>
Full article ">Figure 5
<p>The smoothed fluorescence spectral data of adulterated SBH with three adulteration levels (low, medium, and high adulteration) at a full spectrum of 348.5–866.5 nm.</p>
Full article ">Figure 6
<p>The result of PCA score plot calculation using a full spectrum of 348.5–866.5 nm (<b>a</b>) based on raw spectral data; (<b>b</b>) based on the smoothed spectral data.</p>
Full article ">Figure 7
<p>The calculated wavelength versus x-loadings for the first two PCs using a full spectrum of 348.5–866.5 nm.</p>
Full article ">Figure 8
<p>The calculated Hotelling’s T<sup>2</sup> versus Q-residual using a full spectrum of 348.5–866.5 nm.</p>
Full article ">Figure 9
<p>The result of classification model development with fewer selected variables: (<b>a</b>) the LDA method; (<b>b</b>) the PCA-LDA method.</p>
Full article ">Figure 10
<p>The scatter plot between actual and predicted adulteration levels in calibration and validation: (<b>a</b>) full-spectrum PLS model; (<b>b</b>) combined-interval PLS model.</p>
Full article ">Figure 11
<p>The scatter plot between actual and predicted adulteration levels.</p>
Full article ">
21 pages, 10593 KiB  
Article
Improved Phase Gradient Autofocus Method for Multi-Baseline Circular Synthetic Aperture Radar Three-Dimensional Imaging
by Shiliang Yi, Hongtu Xie, Yuanjie Zhang, Zhitao Wu, Mengfan Ge, Nannan Zhu, Zheng Lu and Pengcheng Qin
Remote Sens. 2024, 16(22), 4242; https://doi.org/10.3390/rs16224242 - 14 Nov 2024
Viewed by 238
Abstract
Multi-baseline circular synthetic aperture radar (MB CSAR) can be applied to obtain a three-dimensional (3D) image of the observation scene. However, the phase error caused by radar platform motion or atmospheric propagation delay restricts its 3D imaging capabilities. The phase error calibration of [...] Read more.
Multi-baseline circular synthetic aperture radar (MB CSAR) can be applied to obtain a three-dimensional (3D) image of the observation scene. However, the phase error caused by radar platform motion or atmospheric propagation delay restricts its 3D imaging capabilities. The phase error calibration of MB CSAR data is an essential step in the 3D imaging procedure due to the limited accuracy of positioning sensors. Phase gradient autofocus (PGA) is widely utilized to estimate the phase errors but is subject to shifts in the direction perpendicular to the line of sight and long iteration time in some sub-apertures. In this paper, an improved PGA method for MB CSAR 3D imaging is proposed, which can suppress the shifts and reduce computation time. This method is based on phase gradient estimation, but the prominent units are selected with an energy criterion. Then, weighted phase gradient estimation is presented to suppress the influence of prominent units with poor quality. Finally, a contrast criterion is adopted to reach faster convergence. The experimental results based on the measured MB CSAR data (Gotcha dataset) demonstrate the validity and feasibility of the proposed phase error calibration method for MB CSAR 3D imaging. Full article
(This article belongs to the Section Engineering Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Imaging geometry of MB CSAR system.</p>
Full article ">Figure 2
<p>The contrast-based convergence condition. (<b>a</b>) Convergence in a sub-aperture. (<b>b</b>) Iteration times of all sub-apertures.</p>
Full article ">Figure 3
<p>The flowchart of the proposed autofocus method.</p>
Full article ">Figure 4
<p>Radar tracks in the Gotcha dataset.</p>
Full article ">Figure 5
<p>Two-dimensional image of the whole scene in HH polarization, obtained at the first pass.</p>
Full article ">Figure 6
<p>Optical images of cars in chosen areas. (<b>a</b>) Ford Taurus; (<b>b</b>) Toyota Camry.</p>
Full article ">Figure 7
<p>Locations of points A, B, C and D in the 2D SAR image of the Ford Taurus.</p>
Full article ">Figure 8
<p>Imaging results of four selected points of the Ford Taurus area using different calibration methods. (<b>a</b>) Point A. (<b>b</b>) Point B. (<b>c</b>) Point C. (<b>d</b>) Point D.</p>
Full article ">Figure 9
<p>Projection of 3D images of the Ford Taurus. The left, middle and right columns are on the <math display="inline"><semantics> <mrow> <mi>x</mi> <mo>−</mo> <mi>y</mi> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>x</mi> <mo>−</mo> <mi>z</mi> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>y</mi> <mo>−</mo> <mi>z</mi> </mrow> </semantics></math> planes, respectively. (<b>a</b>) No calibration. (<b>b</b>) WLS method. (<b>c</b>) PGA method. (<b>d</b>) Proposed method.</p>
Full article ">Figure 10
<p>Three-dimensional images of Ford Taurus formed by different methods. (<b>a</b>) No calibration. (<b>b</b>) WLS method. (<b>c</b>) PGA method. (<b>d</b>) Proposed method.</p>
Full article ">Figure 11
<p>Locations of points E, F, G and H in the 2D SAR image of the Toyota Camry.</p>
Full article ">Figure 12
<p>Imaging results of four selected points of the Toyota Camry area using different calibration methods. (<b>a</b>) Point E. (<b>b</b>) Point F. (<b>c</b>) Point G. (<b>d</b>) Point H.</p>
Full article ">Figure 12 Cont.
<p>Imaging results of four selected points of the Toyota Camry area using different calibration methods. (<b>a</b>) Point E. (<b>b</b>) Point F. (<b>c</b>) Point G. (<b>d</b>) Point H.</p>
Full article ">Figure 13
<p>Three-dimensional images of the Toyota Camry formed by different methods. (<b>a</b>) No calibration. (<b>b</b>) WLS method. (<b>c</b>) PGA method. (<b>d</b>) Proposed method.</p>
Full article ">Figure 14
<p>The projection of 3D images of the Toyota Camry. The left, middle and right columns are on the <math display="inline"><semantics> <mrow> <mi>x</mi> <mo>−</mo> <mi>y</mi> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>x</mi> <mo>−</mo> <mi>z</mi> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>y</mi> <mo>−</mo> <mi>z</mi> </mrow> </semantics></math> planes, respectively. (<b>a</b>) No calibration. (<b>b</b>) WLS method. (<b>c</b>) PGA method. (<b>d</b>) Proposed method.</p>
Full article ">Figure 14 Cont.
<p>The projection of 3D images of the Toyota Camry. The left, middle and right columns are on the <math display="inline"><semantics> <mrow> <mi>x</mi> <mo>−</mo> <mi>y</mi> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>x</mi> <mo>−</mo> <mi>z</mi> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>y</mi> <mo>−</mo> <mi>z</mi> </mrow> </semantics></math> planes, respectively. (<b>a</b>) No calibration. (<b>b</b>) WLS method. (<b>c</b>) PGA method. (<b>d</b>) Proposed method.</p>
Full article ">
14 pages, 5245 KiB  
Article
Numerical Modeling of Instream Flow for Corbicula Habitat Preservation in Aquatic Ecosystem of Seomjin River Estuary, South Korea
by Chunggil Jung, Gayeong Lee and Jongyoon Park
Water 2024, 16(22), 3268; https://doi.org/10.3390/w16223268 - 14 Nov 2024
Viewed by 368
Abstract
South Korea’s River Act mandates the maintenance of instream flow to support river ecosystems. This regulation has evolved from early river management practices to more advanced, systematic approaches, including the Instream Flow Incremental Methodology (IFIM). Despite these advancements, river management in South Korea, [...] Read more.
South Korea’s River Act mandates the maintenance of instream flow to support river ecosystems. This regulation has evolved from early river management practices to more advanced, systematic approaches, including the Instream Flow Incremental Methodology (IFIM). Despite these advancements, river management in South Korea, particularly in the Seomjin River Basin, continues to face numerous challenges. In this study, a three-dimensional numerical model was developed to simulate the hydrodynamic and salinity conditions of the Seomjin River Estuary. This study proposes optimal instream flows to support critical habitats for the Corbicula bivalve, which has seen a significant decline due to salinity intrusion by environmental changes. Using the Environmental Fluid Dynamics Code (EFDC), the model simulates salinity and river discharge with calibration and validation by incorporating historical data. Subsequently, this study evaluates how river discharge affects salinity in four major Corbicula habitats (Dugok, Shinbi, Mokdo, and Hwamok). Finally, we determine the minimum flow (instream flow) needed to sustain Corbicula habitats. In short, this study found that the minimum flow rates (instream flow) required to meet target salinities varied significantly across these sites and under different tidal conditions. These findings highlight the necessity of adapting river flow management practices to preserve the ecological health for Corbicula in the Seomjin River Estuary. Furthermore, this study suggests integrating an additional water supply to be used with local water management plans by suggesting short-term and long-term alternatives in order to sustain adapting river minimum flow (instream flow). Full article
(This article belongs to the Special Issue Research on Watershed Ecology, Hydrology and Climate)
Show Figures

Figure 1

Figure 1
<p>Monitoring sites in focus area. Multi-purpose water supply at Seomjin river and Juam dams, instream flow notification station, and river water intake at Daap. (<b>a</b>) From map of South Korea, the study area has been zoomed out including Seomjin Rive Basin, (<b>b</b>) map of the whole Seomjin River Basin, (<b>c</b>) Daap water intake facility, (<b>d-1</b>) salinity monitoring site at Seomjin bridge, (<b>d-2</b>) salinity monitoring site at Seomjin River bridge.</p>
Full article ">Figure 2
<p>Process of model application for estimating instream flow.</p>
Full article ">Figure 3
<p>Hydrological data measured using current instream flow criterion: (<b>a</b>) daily Daap water intake, and (<b>b</b>) changes of riverbed cross-section at Gurye-gun (Songjeong-ri) station.</p>
Full article ">Figure 4
<p>Estimation of relationship between streamflow and salinity using Daap water intake and with streamflow at Gurye-gun (Songjeong-ri) station during (<b>a</b>) spring tide and (<b>b</b>) neap tides.</p>
Full article ">
29 pages, 5844 KiB  
Article
Early Modeling of the Upcoming Landsat Next Constellation for Soybean Yield Prediction Under Varying Levels of Water Availability
by Luís Guilherme Teixeira Crusiol, Marcos Rafael Nanni, Rubson Natal Ribeiro Sibaldelli, Liang Sun, Renato Herrig Furlanetto, Sergio Luiz Gonçalves, Norman Neumaier and José Renato Bouças Farias
Remote Sens. 2024, 16(22), 4184; https://doi.org/10.3390/rs16224184 - 9 Nov 2024
Viewed by 890
Abstract
The upcoming Landsat Next will provide more frequent land surface observations at higher spatial and spectral resolutions that will greatly benefit the agricultural sector. Early modeling of the upcoming Landsat Next products for soybean yield prediction is essential for long-term satellite monitoring strategies. [...] Read more.
The upcoming Landsat Next will provide more frequent land surface observations at higher spatial and spectral resolutions that will greatly benefit the agricultural sector. Early modeling of the upcoming Landsat Next products for soybean yield prediction is essential for long-term satellite monitoring strategies. In this context, this article evaluates the contribution of Landsat Next’s improved spectral resolution for soybean yield prediction under varying levels of water availability. Ground-based hyperspectral data collected over five cropping seasons at the Brazilian Agricultural Research Corporation were resampled to Landsat Next spectral resolution. The spectral dataset (n = 384) was divided into calibration and external validation datasets and investigated using three strategies for soybean yield prediction: (1) using the reflectance from each spectral band; (2) using existing and new vegetation indices developed based on three general equations: Normalized Difference Vegetation Index (NDVI-like), Band Ratio Vegetation Index (RVI-like), and Band Difference Vegetation Index (DVI-like), replacing the traditional spectral bands by all possible combinations between two bands for index calculation; and (3) using a partial least squares regression (PLSR) model composed of all Landsat Next spectral bands, in comparison to PLSR models using Landsat OLI and Sentienel-2 MSI bands. The results show the distribution of the new spectral bands over the most prominent changes in leaf reflectance due to water deficit, particularly in the visible and shortwave infrared spectrum. (1) Band 18 (centered at 1610 nm) had the highest correlation with yield (R2 = 0.34). (2) A new vegetation index, called Normalized Difference Shortwave Vegetation Index (NDSWVI), is proposed and calculated from bands 19 and 20 (centered at 2028 and 2108 nm). NDSWVI showed the best performance (R2 = 0.37) compared to traditional existing and new vegetation indices. (3) The PLSR model gave the best results (R2 = 0.65), outperforming the Landsat OLI and Sentinel-2 MSI sensors. The improved spectral resolution of Landsat Next is expected to contribute to improved crop monitoring, especially for soybean crops in Brazil, increasing the sustainability of the production systems and strengthening food security in Brazil and globally. Full article
Show Figures

Figure 1

Figure 1
<p>Spectral bandpasses for the sensors on all Landsat satellites [<a href="#B2-remotesensing-16-04184" class="html-bibr">2</a>].</p>
Full article ">Figure 2
<p>Location of Embrapa Soja in the context of Brazil, Paraná State, and the municipality of Londrina; experimental area overview; and description of the weather station and treatment plots: irrigated (IRR), non-irrigated (NIRR) and water deficit induced at vegetative (WDV) and reproductive (WDR) stages.</p>
Full article ">Figure 3
<p>Climatic water balance at 10-day periods in the WDV, WDR, NIRR, and IRR treatments in 2016/2017, 2017/2018, 2018/2019, 2022/2023, and 2023/2024 cropping seasons.</p>
Full article ">Figure 4
<p>Soil moisture content (%) at 0–20 cm and 20–40 cm depths in the 2016/2017, 2017/2018, 2018/2019, 2022/2023 and 2023/2024 cropping seasons at the transition from vegetative to reproductive stages (<b>a</b>) and at the R5 phenological stages towards the maturity stage (<b>b</b>). Means followed by the same letter among treatments within each depth and on each date do not differ by Tukey’s test (<span class="html-italic">p</span> ≤ 0.05).</p>
Full article ">Figure 5
<p>Soybean yield (kg ha<sup>−1</sup>) in the 2016/2017, 2017/2018, 2018/2019, 2022/2023 and 2023/2024 cropping seasons. Means followed by the same letter do not differ by the Tukey test (<span class="html-italic">p</span> ≤ 0.05).</p>
Full article ">Figure 6
<p>Spectral assessment in the field (<b>a</b>), detail of spectroradiometer (<b>b</b>), and the plant probe device (<b>c</b>). Photo by Décio de Assis—Embrapa Soja.</p>
Full article ">Figure 7
<p>Flowchart of the methodology adopted for spectral data processing and yield modeling.</p>
Full article ">Figure 8
<p>Yield values (<b>a</b>) and principal component analysis of the spectral response in the correspondent Landsat Next spectral bands (<b>b</b>) from samples collected in the 2016/2017, 2017/2018, 2018/2019, 2023/2023, and 2023/2024 cropping seasons pooled into the calibration (red squares) and external validation (green dots) datasets.</p>
Full article ">Figure 9
<p>Soybean spectral response across Vis-NIR-SWIR wavelengths in the IRR and WDR treatments (<b>a</b>) and the percentage of reflectance increasing from WDR treatment in relation to IRR with the delimitation of the correspondent Landsat Next spectral bands (color bars—(<b>b</b>)).</p>
Full article ">Figure 10
<p>Principal component analysis of the spectral response in the correspondent Landsat Next spectral bands of soybean crop under the evaluated water conditions in 2016/2017 (<b>a</b>), 2017/2018 (<b>b</b>), 2018/2019 (<b>c</b>), 2022/2023 (<b>d</b>) and 2023/2024 (<b>e</b>) cropping seasons.</p>
Full article ">Figure 11
<p>Correlation between soybean yield and Landsat Next spectral bands reflectance using the calibration dataset.</p>
Full article ">Figure 12
<p>Correlation between observed and predicted values of soybean yield through linear regression between yield and band 18 from Landsat Next.</p>
Full article ">Figure 13
<p>Coefficient of determination from the linear regression between soybean yield and all possible combinations for calculating two-band vegetation indices using Landsat Next spectral band reflectance from the training dataset (288 samples) under normalized difference (<b>a</b>), ratio (<b>b</b>), and difference (<b>c</b>) equations.</p>
Full article ">Figure 14
<p>Correlation between observed and predicted values of soybean yield through linear regression between yield and the outstanding NDVI (<b>a</b>), RVI (<b>b</b>), and DVI (<b>c</b>).</p>
Full article ">Figure 15
<p>Correlation between observed and predicted values of soybean yield through PLSR at the calibration and cross-validation (leave-one-out) stage (<b>a</b>) using 75% of data (288 samples—training dataset) and validated with the remaining 25% of the data (96 samples—testing dataset) at the external validation stage (<b>b</b>).</p>
Full article ">Figure 16
<p>Regression coefficients of PLSR for soybean grain yield prediction at R5 stage in a model developed using 75% of data (288 samples—calibration dataset).</p>
Full article ">Figure 17
<p>Correlation between observed and predicted values of soybean yield using band 18; NDVI calculated with bands 19 and 20, RVI calculated with bands 19 and 20; DVI calculated with bands 4 and 9; PLSR model using all Landsat Next bands; PLSR model using all Landsat OLI bands; and PLSR model using all Sentinel-2 MSI bands.</p>
Full article ">Figure 18
<p>Yield values from samples collected in the 2016–2017, 2017–2018, 2018–2019, 2023–2023 and 2023–2024 cropping seasons.</p>
Full article ">Figure 19
<p>Correlation between observed and predicted values of soybean yield through PLSR at the calibration and cross-validation (leave-one-out) steps using spectral response in the correspondent Landsat Next spectral bands in 2016–2017 (<b>a</b>), 2017–2018 (<b>b</b>), 2018–2019 (<b>c</b>), 2022–2023 (<b>d</b>) and 2023–2024 (<b>e</b>) cropping seasons.</p>
Full article ">Figure 20
<p>Statistics at the cross-validation step for soybean prediction using Landsat Next spectral bands under PLSR modeling for each soybean genotype within 2016/2017, 2017/2018, 2018/2019, 2022/2023, and 2023/2024 cropping seasons.</p>
Full article ">
25 pages, 7729 KiB  
Article
A Fast-Calibrated Computational Fluid Dynamic Model for Timber–Concrete Composite Ventilated Façades
by Sofia Pastori, Mohammed-Sadegh Salehi, Stefan Radl and Enrico Sergio Mazzucchelli
Buildings 2024, 14(11), 3567; https://doi.org/10.3390/buildings14113567 - 9 Nov 2024
Viewed by 489
Abstract
Timber–concrete composite (TCC) systems join the positive aspects of engineered wood products (good seismftaic behaviour, low thermal conductivity, environmental sustainability, good behaviour under fire if appropriately designed) with those of concrete (high thermal inertia, durability, excellent fire resistance). TCC facades are typically composed [...] Read more.
Timber–concrete composite (TCC) systems join the positive aspects of engineered wood products (good seismftaic behaviour, low thermal conductivity, environmental sustainability, good behaviour under fire if appropriately designed) with those of concrete (high thermal inertia, durability, excellent fire resistance). TCC facades are typically composed of an internal insulated timber-frame wall and an external concrete slab, separated by a ventilated air cavity. However, there is very limited knowledge concerning the performance of TCC facades, especially concerning their thermal behaviour. The present paper deals with the development and optimization of a 2D Computational Fluid Dynamic (CFD) model for the analysis of TCC ventilated façades’ thermal behaviour. The model is calibrated and validated against experimental data collected during the annual monitoring of a real TCC ventilated envelope in the north of Italy. Also, a new solver algorithm is developed to significantly speed up the simulation (i.e., 45 times faster simulation at an error below 3.5 °C compared to a typical CFD solver). The final model can be used for the time-efficient analysis (simulation time of approximately 23 min for a full day in real-time) and the optimization of the thermal performance of TCC ventilated facades, as well as other ventilated facades with external massive cladding. Our simulation strategy partially avoids the expensive and time-consuming construction of mock-ups, or the use of comparably slow (conventional) CFD solvers that are less suitable for optimization studies. Full article
(This article belongs to the Special Issue Thermal Fluid Flow and Heat Transfer in Buildings)
Show Figures

Figure 1

Figure 1
<p>Horizontal section of the TCC façade studied (units in cm). The layers are (1) a reinforced concrete slab; (2) ventilated air cavity; (3) OSB panel; (4) rockwool insulation (100 kg/m<sup>3</sup>); (5) timber-frame structure; and (6) an OSB panel. © Pastori S.</p>
Full article ">Figure 2
<p>Monitored building, and illustration of the TCC ventilated envelope area analyzed (red dashed box). © Pastori S.</p>
Full article ">Figure 3
<p>Naming of the regions used in the CFD model.</p>
Full article ">Figure 4
<p>Main boundary conditions used in the CFD model.</p>
Full article ">Figure 5
<p>Geometry of and mesh within the 2D CFD model (colors indicate different regions).</p>
Full article ">Figure 6
<p>Illustration of the new solver algorithm.</p>
Full article ">Figure 7
<p>Graphical representation of the values listed in <a href="#buildings-14-03567-t006" class="html-table">Table 6</a>.</p>
Full article ">Figure 8
<p>Prediction of (<b>a</b>) TO, (<b>b</b>) TC, and (<b>c</b>) TC_ext by model t2000 and t4500.</p>
Full article ">Figure 9
<p>Volume-averaged temperature of fluid.1, as well as the patch-averaged temperature of solid2.1 and solid.1.3 for case t0104 (<b>a</b>), and case t4500 (<b>b</b>).</p>
Full article ">Figure 10
<p>Comparison of experimental and simulation data at the surface of the OSB panel facing the cavity (“TO”, left) and the concrete slab facing the cavity (“TC”, right) for the base case (i.e., case t4500) and a case with high conductivity of the rockwool insulation. “T_exp,a” indicates the imposed outside air temperature, and “T_i” the indoor air temperature.</p>
Full article ">Figure A1
<p>Temporal evolution of the boundary conditions considered for summer sunny days: outdoor temperature (T_ext), indoor temperature (T_int), incident solar irradiation on the façade (Solar_irr), wind speed (Wind_speed), and air speed at the bottom opening of the ventilated cavity (Air_speed).</p>
Full article ">Figure A2
<p>Temporal evolution of the boundary conditions considered for summer cloudy days (the meaning of the colours is identical to that in <a href="#buildings-14-03567-f0A1" class="html-fig">Figure A1</a>).</p>
Full article ">Figure A3
<p>The algorithm for the initial period (<math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="bold-italic">τ</mi> </mrow> <mrow> <mi mathvariant="bold-italic">i</mi> <mi mathvariant="bold-italic">n</mi> <mi mathvariant="bold-italic">i</mi> <mi mathvariant="bold-italic">t</mi> <mi mathvariant="bold-italic">i</mi> <mi mathvariant="bold-italic">a</mi> <mi mathvariant="bold-italic">l</mi> </mrow> </msub> </mrow> </semantics></math>) in the new solver.</p>
Full article ">Figure A4
<p>An illustration of the steps in the solver algorithm defining the phases.</p>
Full article ">Figure A5
<p>Illustration of the “normal period” for the solver algorithm.</p>
Full article ">Figure A6
<p>Effect of a turbulence model in comparison to the base case results (i.e., case t4500) on the solid temperature distribution (data at 8:00 a.m. is shown; the figure has been scaled in the y-direction for better representation; in the rightmost panel the predicted turbulent viscosity is illustrated).</p>
Full article ">Figure A7
<p>Effect of using a 3D computational domain in comparison to the base case results (i.e., case t4500) on the solid and fluid temperature distribution (data at time <span class="html-italic">t</span> = 120 [min] is shown; the figure has been scaled in the y-direction for better representation).</p>
Full article ">Figure A8
<p>Effect of different TCC layer thicknesses on the temperatures facing the cavity (“base case” refers to case t4500).</p>
Full article ">
23 pages, 7255 KiB  
Article
Exploring the Relationship Between Very-High-Resolution Satellite Imagery Data and Fruit Count for Predicting Mango Yield at Multiple Scales
by Benjamin Adjah Torgbor, Priyakant Sinha, Muhammad Moshiur Rahman, Andrew Robson, James Brinkhoff and Luz Angelica Suarez
Remote Sens. 2024, 16(22), 4170; https://doi.org/10.3390/rs16224170 - 8 Nov 2024
Viewed by 658
Abstract
Tree- and block-level prediction of mango yield is important for farm operations, but current manual methods are inefficient. Previous research has identified the accuracies of mango yield forecasting using very-high-resolution (VHR) satellite imagery and an ’18-tree’ stratified sampling method. However, this approach still [...] Read more.
Tree- and block-level prediction of mango yield is important for farm operations, but current manual methods are inefficient. Previous research has identified the accuracies of mango yield forecasting using very-high-resolution (VHR) satellite imagery and an ’18-tree’ stratified sampling method. However, this approach still requires infield sampling to calibrate canopy reflectance and the derived block-level algorithms are unable to translate to other orchards due to the influences of abiotic and biotic conditions. To better appreciate these influences, individual tree yields and corresponding canopy reflectance properties were collected from 2015 to 2021 for 1958 individual mango trees from 55 orchard blocks across 14 farms located in three mango growing regions of Australia. A linear regression analysis of the block-level data revealed the non-existence of a universal relationship between the 24 vegetation indices (VIs) derived from VHR satellite data and fruit count per tree, an outcome likely due to the influence of location, season, management and cultivar. The tree-level fruit count predicted using a random forest (RF) model trained on all calibration data produced a percentage root mean squared error (PRMSE) of 26.5% and a mean absolute error (MAE) of 48 fruits/tree. The lowest PRMSEs produced from RF-based models developed from location, season and cultivar subsets at the individual tree level ranged from 19.3% to 32.6%. At the block level, the PRMSE for the combined model was 10.1% and the lowest values for the location, seasonal and cultivar subset models varied between 7.2% and 10.0% upon validation. Generally, the block-level predictions outperformed the individual tree-level models. Maps were produced to provide mango growers with a visual representation of yield variability across orchards. This enables better identification and management of the influence of abiotic and biotic constraints on production. Future research could investigate the causes of spatial yield variability in mango orchards. Full article
Show Figures

Figure 1

Figure 1
<p>Location of mango farms in the three mango growing regions of Australia.</p>
Full article ">Figure 2
<p>Flowchart showing the sequence of procedure steps used in this study to generate the results.</p>
Full article ">Figure 3
<p>Example of 18 tree locations on the classified NDVI map (<b>a</b>) and on the ESRI basemap image (<b>b</b>). The points with L, M and H prefixes represent the different tree vigour classes of low, medium and high, respectively.</p>
Full article ">Figure 4
<p>Summary of fruits counted (<b>a</b>) per farm and (<b>b</b>) heterogeneity of cultivar yield distribution from 2015 to 2021. The numerical values and black dots associated with each boxplot represent the number of trees of that particular cultivar and outliers, respectively.</p>
Full article ">Figure 5
<p>Correlation between fruit count and the 24 VIs using the entire datasets of 1958 datapoints. The green and red colour ramps show the strength and direction of the correlation being positive and negative, respectively.</p>
Full article ">Figure 6
<p>Distribution of slopes for CIRE_1 with average slope and standard deviation.</p>
Full article ">Figure 7
<p>Relationships identified between RENDVI and fruit count: (<b>a</b>) and (<b>b</b>) were positive for 2016 and 2017, (<b>c</b>) negative for 2020 and (<b>d</b>) non-existent for 2021.</p>
Full article ">Figure 8
<p>RF prediction of fruit count using all individual tree datasets (combined model). The different coloured points represent the sampled trees from the respective farms and regions. n = 390 represents the number of datapoints (20%) used for model validation.</p>
Full article ">Figure 9
<p>RF-based location (region) prediction of fruit count in the (<b>a</b>) Northern Territory (NT), (<b>b</b>) Northern Queensland (N–QLD) and (<b>c</b>) South East Queensland (SE–QLD). The different coloured points represent the sampled trees on a given farm in the respective regions.</p>
Full article ">Figure 10
<p>RF-based variable importance plots for models from (<b>a</b>) combined datasets, (<b>b</b>) Northern Territory (NT), (<b>c</b>) Northern Queensland (N–QLD) and (<b>d</b>) South East Queensland (SE–QLD) and the best (<b>e</b>) seasonal and (<b>f</b>) cultivar models.</p>
Full article ">Figure 11
<p>Comparison of total actual and predicted yield for the 51 validation points (blocks per season) obtained from 29 unique blocks with available actual harvest data from 2016 to 2021.</p>
Full article ">Figure 12
<p>An example of a tree-level yield variability map derived from the RF-based combined model (<b>right</b>). The RGB image of the mango orchard mapped is shown on the (<b>left</b>). The legend presents an industry-based categorization of yield variability ranging from low (0–55) to high (139–170) for this study.</p>
Full article ">
24 pages, 5699 KiB  
Article
Synthetic Wind Estimation for Small Fixed-Wing Drones
by Aman Sharma, Gabriel François Laupré, Pasquale Longobardi and Jan Skaloud
Atmosphere 2024, 15(11), 1339; https://doi.org/10.3390/atmos15111339 - 8 Nov 2024
Viewed by 387
Abstract
Wind estimation is crucial for studying the atmospheric boundary layer. Traditional methods such as weather balloons offer limited in situ capabilities; besides an Air Data System (ADS) combined with inertial measurements and satellite positioning is required to estimate the wind on fixed-wing drones. [...] Read more.
Wind estimation is crucial for studying the atmospheric boundary layer. Traditional methods such as weather balloons offer limited in situ capabilities; besides an Air Data System (ADS) combined with inertial measurements and satellite positioning is required to estimate the wind on fixed-wing drones. As pressure probes are an important constituent of an ADS, they are susceptible to malfunctioning or failure due to blockages, thus affecting the capability of wind sensing and possibly the safety of the drone. This paper presents a novel approach, using low-fidelity aerodynamic models of drones to estimate wind synthetically. In our work, the aerodynamic model parameters are derived from post-processed flight data, in contrast to existing approaches that use expensive wind tunnel calibration for identifying the same. In sum, our method integrates aerodynamic force and moment models into a Vehicle Dynamic Model (VDM)-based navigation filter to yield a synthetic wind estimate without relying on an airspeed sensor. We validate our approach using two geometrically distinct drones, each characterized by a unique aerodynamic model and different quality of inertial sensors, altogether tested across several flights. Experimental results demonstrate that the proposed cross-platform method provides a synthetic wind velocity estimate, thus offering a practical backup to traditional techniques. Full article
(This article belongs to the Section Atmospheric Techniques, Instruments, and Modeling)
Show Figures

Figure 1

Figure 1
<p>VDM-based navigation system. Image courtesy: Ref. [<a href="#B71-atmosphere-15-01339" class="html-bibr">71</a>].</p>
Full article ">Figure 2
<p>Aerodynamic calibration procedure. Image courtesy: Ref. [<a href="#B4-atmosphere-15-01339" class="html-bibr">4</a>].</p>
Full article ">Figure 3
<p>TP2 payload: (<b>left</b>) CAD model; (<b>right</b>) practical realization.</p>
Full article ">Figure 4
<p>Drones: (<b>top</b>) <span class="html-italic">TP2</span>; (<b>bottom</b>) <span class="html-italic">Concorde S</span>.</p>
Full article ">Figure 5
<p>Comparison of wind estimated by VDM and INS/GNSS/Pitot fusion for <span class="html-italic">TP-2</span>.</p>
Full article ">Figure 5 Cont.
<p>Comparison of wind estimated by VDM and INS/GNSS/Pitot fusion for <span class="html-italic">TP-2</span>.</p>
Full article ">Figure 6
<p>Wind residual between VDM and INS/GNSS/Pitot fusion for <span class="html-italic">TP-2</span>.</p>
Full article ">Figure 7
<p>Zoomed-in view of the residual error.</p>
Full article ">Figure 8
<p>Comparison of wind estimated by VDM and INS/GNSS/Pitot fusion for <span class="html-italic">Concorde S</span>.</p>
Full article ">Figure 9
<p>Wind residual between VDM and INS/GNSS/Pitot fusion for <span class="html-italic">Concorde S</span>.</p>
Full article ">Figure A1
<p>Wind estimated using incorrect VDM parameters for <span class="html-italic">TP2</span>—STIM13.</p>
Full article ">Figure A2
<p>eBeeX drone.</p>
Full article ">Figure A3
<p>Wind estimated by eBeex.</p>
Full article ">
15 pages, 2425 KiB  
Article
The Demographic and Clinical Characteristics, Prognostic Factors, and Survival Outcomes of Head and Neck Carcinosarcoma: A SEER Database Analysis
by Wanting Hou, Ouying Yan and Hong Zhu
Biomedicines 2024, 12(11), 2556; https://doi.org/10.3390/biomedicines12112556 - 8 Nov 2024
Viewed by 337
Abstract
Background: Head and neck carcinosarcoma (HNCS) is a rare and highly aggressive malignancy with limited research, resulting in an incomplete understanding of disease progression and a lack of reliable prognostic tools. This study aimed to retrospectively analyze the clinical characteristics and outcomes of [...] Read more.
Background: Head and neck carcinosarcoma (HNCS) is a rare and highly aggressive malignancy with limited research, resulting in an incomplete understanding of disease progression and a lack of reliable prognostic tools. This study aimed to retrospectively analyze the clinical characteristics and outcomes of HNCS patients using data from the Surveillance, Epidemiology, and End Results (SEER) database and to develop a nomogram to predict overall survival (OS) and cancer-specific survival (CSS). Methods: Patients diagnosed with HNCS from 1975 to 2020 were identified in the SEER database. Univariate and multivariate Cox regression analyses were conducted to identify independent prognostic indicators, with the optimal model selected using the minimal Akaike Information Criterion (AIC). The identified prognostic factors were incorporated into nomograms to predict OS and CSS. Model performance was assessed using the concordance index (C-index), area under the curve (AUC), calibration curves, and decision curve analysis (DCA). Survival curves were generated using Kaplan–Meier analysis and compared via the log-rank test. Results: A total of 152 HNCS patients were included, with 108 assigned to the training cohort and 44 to the validation cohort in a 7:3 ratio. Prognostic factors including age, primary tumor site, marital status, radiotherapy, chemotherapy, tumor size, pathological grade, and tumor stage were incorporated into the nomogram models. The models demonstrated strong predictive performance, with C-index values for OS and CSS of 0.757 and 0.779 in the training group, and 0.777 and 0.776 in the validation group, respectively. AUC values for predicting 3-, 5-, and 10-year OS were 0.662, 0.713, and 0.761, and for CSS the values were 0.726, 0.703, and 0.693. Kaplan–Meier analysis indicated significantly improved survival for patients with lower risk scores. The 3-, 5-, and 10-year OS rates for the entire cohort were 54.1%, 45.6%, and 35.1%, respectively, and the CSS rates were 62.9%, 57.5%, and 52.2%, respectively. Conclusions: This study provides validated nomograms for predicting OS and CSS in HNCS patients, offering a reliable tool to support clinical decision-making for this challenging malignancy. These nomograms enhance the ability to predict patient prognosis and personalize treatment strategies. Full article
(This article belongs to the Section Molecular and Translational Medicine)
Show Figures

Figure 1

Figure 1
<p>Survival outcomes for HNCS by cohort (<b>a</b>,<b>b</b>), tumor site (<b>c</b>), and disease stage (<b>d</b>).</p>
Full article ">Figure 2
<p>Nomograms constructed to predict 3-, 5-, and 10-year OS (<b>a</b>) and CSS (<b>b</b>) rates in patients with HNCS.</p>
Full article ">Figure 3
<p>Receiver operating characteristic (ROC) curves for evaluating the predictive ability of independent prognostic factors for 3-, 5-, and 10-year OS and CSS in HNCS patients in the training cohort: (<b>a</b>) 3-year OS; (<b>b</b>) 5-year OS; (<b>c</b>) 10-year OS; (<b>d</b>) 3-year CSS; (<b>e</b>) 5-year CSS; (<b>f</b>) 10-year CSS.</p>
Full article ">Figure 4
<p>Calibration curves for the nomogram in the training group: (<b>a</b>) 3-year OS, (<b>b</b>) 5-year OS, (<b>c</b>) 10-year OS, (<b>d</b>) 3-year CSS, (<b>e</b>) 5-year CSS, and (<b>f</b>) 10-year CSS.</p>
Full article ">Figure 5
<p>Decision curve analysis for the nomograms in the training group: (<b>a</b>) 3-year OS, (<b>b</b>) 5-year OS, (<b>c</b>) 10-year OS, (<b>d</b>) 3-year CSS, (<b>e</b>) 5-year CSS, and (<b>f</b>) 10-year CSS. The <span class="html-italic">x</span>-axis represents the threshold probability, and the <span class="html-italic">y</span>-axis represents the net benefit.</p>
Full article ">Figure 6
<p>Kaplan–Meier (K-M) curves for the nomogram in the training group: (<b>a</b>) OS and (<b>b</b>) CSS.</p>
Full article ">
21 pages, 18420 KiB  
Article
High-Resolution Mapping of Litter and Duff Fuel Loads Using Multispectral Data and Random Forest Modeling
by Álvaro Agustín Chávez-Durán, Miguel Olvera-Vargas, Inmaculada Aguado, Blanca Lorena Figueroa-Rangel, Ramón Trucíos-Caciano, Ernesto Alonso Rubio-Camacho, Jaqueline Xelhuantzi-Carmona and Mariano García
Fire 2024, 7(11), 408; https://doi.org/10.3390/fire7110408 - 7 Nov 2024
Viewed by 519
Abstract
Forest fuels are the core element of fire management; each fuel component plays an important role in fire behavior. Therefore, accurate determination of their characteristics and spatial distribution is crucial. This paper introduces a novel method for mapping the spatial distribution of litter [...] Read more.
Forest fuels are the core element of fire management; each fuel component plays an important role in fire behavior. Therefore, accurate determination of their characteristics and spatial distribution is crucial. This paper introduces a novel method for mapping the spatial distribution of litter and duff fuel loads using data collected by unmanned aerial vehicles. The approach leverages a very high-resolution multispectral data analysis within a machine learning framework to achieve precise and detailed results. A set of vegetation indices and texture metrics derived from the multispectral data, optimized by a “Variable Selection Using Random Forests” (VSURF) algorithm, were used to train random forest (RF) models, enabling the modeling of high-resolution maps of litter and duff fuel loads. A field campaign to measure fuel loads was conducted in the mixed forest of the natural protected area of “Sierra de Quila”, Jalisco, Mexico, to measure fuel loads and obtain field reference data for calibration and validation purposes. The results revealed moderate determination coefficients between observed and predicted fuel loads with R2 = 0.32, RMSE = 0.53 Mg/ha for litter and R2 = 0.38, RMSE = 13.14 Mg/ha for duff fuel loads, both with significant p-values of 0.018 and 0.015 for litter and duff fuel loads, respectively. Moreover, the relative root mean squared errors were 33.75% for litter and 27.71% for duff fuel loads, with a relative bias of less than 5% for litter and less than 20% for duff fuel loads. The spatial distribution of the litter and duff fuel loads was coherent with the structure of the vegetation, despite the high complexity of the study area. Our modeling approach allows us to estimate the continuous high-resolution spatial distribution of litter and duff fuel loads, aligned with their ecological context, which dictates their dynamics and spatial variability. The method achieved acceptable accuracy in monitoring litter and duff fuel loads, providing researchers and forest managers with timely data to expedite decision-making in fire and forest fuel management. Full article
Show Figures

Figure 1

Figure 1
<p>Image depicting the study area location and showing sampling plots. The projection coordinate system used is Universal Transverse Mercator Zone 13 North (UTM 13N).</p>
Full article ">Figure 2
<p>Flowchart for estimating the spatial distribution of litter and duff fuel loads using variable selection with random forests (VSURF).</p>
Full article ">Figure 3
<p>Sampling plots for field data. (<b>A</b>) One-hectare sampling plot divided into 400 m<sup>2</sup> (20 × 20 m) subplots. (<b>B</b>) A subplot showing nanoplots of 0.09 m<sup>2</sup> (0.30 × 0.30 m).</p>
Full article ">Figure 4
<p>Field data collection and laboratory process. (<b>A</b>) Litter and duff material collection in situ. (<b>B</b>) Sample classification in laboratory.</p>
Full article ">Figure 5
<p>Multispectral orthomosaics for P1, P2, and P3 from the Parrot Sequoia sensor. False-color composition: Green, red, and red edge bands.</p>
Full article ">Figure 6
<p>Predicted vs. observed litter (<b>left</b>) and duff (<b>right</b>) fuel loads. The blue solid line represents the trend line and the shadowed area the 95% confidence interval. The black solid represents the 1:1 line.</p>
Full article ">Figure 7
<p>Litter fuel load spatial distribution. P1, P2, and P3 are 1 ha sampling plots. (<b>A</b>) Litter fuel loads (Mg/ha); (<b>B</b>) Litter fuel load uncertainty (%).</p>
Full article ">Figure 8
<p>Duff fuel load spatial distribution. P1, P2, and P3 are 1 ha sampling plots. (<b>A</b>) Duff fuel loads (Mg/ha); (<b>B</b>) Duff fuel loads uncertainty (%).</p>
Full article ">Figure A1
<p>Vegetation indices used to estimate the spatial distribution of litter fuel loads. P1, P2, and P3 are 1 ha sampling plots; the black solid lines represent the boundaries of each sampling plot. (<b>A</b>) GNDVI; (<b>B</b>) NDVI.</p>
Full article ">Figure A2
<p>Texture bands used to estimate the spatial distribution of litter fuel loads. P1, P2, and P3 are 1 ha sampling plots; the black solid lines represent the boundaries of each sampling plot. (<b>A</b>) Green band homogeneity texture; (<b>B</b>) Red band homogeneity texture; (<b>C</b>) Green band variance texture.</p>
Full article ">Figure A3
<p>Vegetation indices used to estimate the spatial distribution of duff fuel loads. P1, P2, and P3 are 1 ha sampling plots; the black solid lines represent the boundaries of each sampling plot. (<b>A</b>) GRVI; (<b>B</b>) Datt4; (<b>C</b>) NDVIre.</p>
Full article ">Figure A4
<p>Texture bands used to estimate the spatial distribution of duff fuel loads. P1, P2, and P3 are 1 ha sampling plots; the black solid lines represent the boundaries of each sampling plot. (<b>A</b>) Red edge band homogeneity texture; (<b>B</b>) Green band second moment texture; (<b>C</b>) Red edge band contrast texture.</p>
Full article ">
23 pages, 5261 KiB  
Article
Autonomous Underwater Pipe Damage Detection Positioning and Pipe Line Tracking Experiment with Unmanned Underwater Vehicle
by Seda Karadeniz Kartal and Recep Fatih Cantekin
J. Mar. Sci. Eng. 2024, 12(11), 2002; https://doi.org/10.3390/jmse12112002 - 7 Nov 2024
Viewed by 484
Abstract
Underwater natural gas pipelines constitute critical infrastructure for energy transportation. Any damage or leakage in these pipelines poses serious security risks, directly threatening marine and lake ecosystems, and potentially causing operational issues and economic losses in the energy supply chain. However, current methods [...] Read more.
Underwater natural gas pipelines constitute critical infrastructure for energy transportation. Any damage or leakage in these pipelines poses serious security risks, directly threatening marine and lake ecosystems, and potentially causing operational issues and economic losses in the energy supply chain. However, current methods for detecting deterioration and regularly inspecting these submerged pipelines remain limited, as they rely heavily on divers, which is both costly and inefficient. Due to these challenges, the use of unmanned underwater vehicles (UUVs) becomes crucial in this field, offering a more effective and reliable solution for pipeline monitoring and maintenance. In this study, we conducted an underwater pipeline tracking and damage detection experiment using a remote-controlled unmanned underwater vehicle (UUV) with autonomous features. The primary objective of this research is to demonstrate that UUV systems provide a more cost-effective, efficient, and practical alternative to traditional, more expensive methods for inspecting submerged natural gas pipelines. The experimental method included vehicle (UUV) setup, pre-test calibration, pipeline tracking mechanism, 3D navigation control, damage detection, data processing, and analysis. During the tracking of the underwater pipeline, damages were identified, and their locations were determined. The navigation information of the underwater vehicle, including orientation in the x, y, and z axes (roll, pitch, yaw) from a gyroscope integrated with a magnetic compass, speed and position information in three axes from an accelerometer, and the distance to the water surface from a pressure sensor, was integrated into the vehicle. Pre-tests determined the necessary pulse width modulation values for the vehicle’s thrusters, enabling autonomous operation by providing these values as input to the thruster motors. In this study, 3D movement was achieved by activating the vehicle’s vertical thruster to maintain a specific depth and applying equal force to the right and left thrusters for forward movement, while differential force was used to induce deviation angles. In pool experiments, the unmanned underwater vehicle autonomously tracked the pipeline as intended, identifying damages on the pipeline using images captured by the vehicle’s camera. The images for damage assessment were processed using a convolutional neural network (CNN) algorithm, a deep learning method. The position of the damage relative to the vehicle was estimated from the pixel dimensions of the identified damage. The location of the damage relative to its starting point was obtained by combining these two positional pieces of information from the vehicle’s navigation system. The damages in the underwater pipeline were successfully detected using the CNN algorithm. The training accuracy and validation accuracy of the CNN algorithm in detecting underwater pipeline damages were 94.4% and 92.87%, respectively. The autonomous underwater vehicle also followed the designated underwater pipeline route with high precision. The experiments showed that the underwater vehicle followed the pipeline path with an error of 0.072 m on the x-axis and 0.037 m on the y-axis. Object recognition and the automation of the unmanned underwater vehicle were implemented in the Python environment. Full article
(This article belongs to the Special Issue Autonomous Marine Vehicle Operations—2nd Edition)
Show Figures

Figure 1

Figure 1
<p>Unmanned underwater vehicle used in the experiment.</p>
Full article ">Figure 2
<p>Damaged pipe used in the experiment (red-colored region is defined as “damage”).</p>
Full article ">Figure 3
<p>Damaged pipe used in the experiment (red-colored region is defined as “damage”).</p>
Full article ">Figure 4
<p>Pipeline damage detection experiment results.</p>
Full article ">Figure 5
<p>CNN training consistency and loss curves.</p>
Full article ">Figure 6
<p>Pipeline tracking experiment scenario schematic representation.</p>
Full article ">Figure 7
<p>Pipeline tracking pool experiment.</p>
Full article ">Figure 8
<p>Vehicle’s pixhawk and thrusters used in the experiment.</p>
Full article ">Figure 9
<p>Reference surge speed (blue line) and measured (red line) surge speed.</p>
Full article ">Figure 10
<p>Reference yaw angle (blue line) and measured yaw angle (red line).</p>
Full article ">Figure 11
<p>Reference depth value (blue line) and measured depth value (red line).</p>
Full article ">Figure 12
<p>Location 1 of pipe line damage detection.</p>
Full article ">Figure 13
<p>Location 2 of pipe line damage detection.</p>
Full article ">Figure 14
<p>Location 3 of pipe line damage detection.</p>
Full article ">Figure 15
<p>Reference path (blue line) and followed path (red line).</p>
Full article ">
25 pages, 23247 KiB  
Article
Infrared and Visible Camera Integration for Detection and Tracking of Small UAVs: Systematic Evaluation
by Ana Pereira, Stephen Warwick, Alexandra Moutinho and Afzal Suleman
Drones 2024, 8(11), 650; https://doi.org/10.3390/drones8110650 - 6 Nov 2024
Viewed by 484
Abstract
Given the recent proliferation of Unmanned Aerial Systems (UASs) and the consequent importance of counter-UASs, this project aims to perform the detection and tracking of small non-cooperative UASs using Electro-optical (EO) and Infrared (IR) sensors. Two data integration techniques, at the decision and [...] Read more.
Given the recent proliferation of Unmanned Aerial Systems (UASs) and the consequent importance of counter-UASs, this project aims to perform the detection and tracking of small non-cooperative UASs using Electro-optical (EO) and Infrared (IR) sensors. Two data integration techniques, at the decision and pixel levels, are compared with the use of each sensor independently to evaluate the system robustness in different operational conditions. The data are submitted to a YOLOv7 detector merged with a ByteTrack tracker. For training and validation, additional efforts are made towards creating datasets of spatially and temporally aligned EO and IR annotated Unmanned Aerial Vehicle (UAV) frames and videos. These consist of the acquisition of real data captured from a workstation on the ground, followed by image calibration, image alignment, the application of bias-removal techniques, and data augmentation methods to artificially create images. The performance of the detector across datasets shows an average precision of 88.4%, recall of 85.4%, and [email protected] of 88.5%. Tests conducted on the decision-level fusion architecture demonstrate notable gains in recall and precision, although at the expense of lower frame rates. Precision, recall, and frame rate are not improved by the pixel-level fusion design. Full article
(This article belongs to the Special Issue Intelligent Image Processing and Sensing for Drones 2nd Edition)
Show Figures

Figure 1

Figure 1
<p>Decision-level data fusion stages.</p>
Full article ">Figure 2
<p>Pixel-level data fusion stages.</p>
Full article ">Figure 3
<p>FusionGAN algorithm application: (<b>a</b>) Input Electro-optical (EO) image. (<b>b</b>) Input Infrared (IR) image. (<b>c</b>) Pixel Fused output image.</p>
Full article ">Figure 4
<p>Experimental setup with highlight on the sensors.</p>
Full article ">Figure 5
<p>Calibration procedure using the calibration board created: (<b>a</b>) EO image at close range. (<b>b</b>) EO image at far range. (<b>c</b>) IR image at close range. (<b>d</b>) IR image at far range.</p>
Full article ">Figure 5 Cont.
<p>Calibration procedure using the calibration board created: (<b>a</b>) EO image at close range. (<b>b</b>) EO image at far range. (<b>c</b>) IR image at close range. (<b>d</b>) IR image at far range.</p>
Full article ">Figure 6
<p>UAVs captured during flight experiments: (<b>a</b>) FE A—Mini-E. (<b>b</b>) FE A—DJI Mavic 2. (<b>c</b>) FE B—MIMIQ. (<b>d</b>) FE C—DJI Inspire 1. (<b>e</b>) FE D—Zeta FX-61 Phantom Wing. (<b>f</b>) FE D—DJI Mini 3 Pro.</p>
Full article ">Figure 7
<p>Schematics of flight paths: waypoints (blue) and workstation (red): (<b>a</b>) Flight Experiment A (Mini-E). (<b>b</b>) Flight Experiment A (DJI Mavic 2). (<b>c</b>) Flight Experiment C.</p>
Full article ">Figure 8
<p>Bias-removal algorithm application: (<b>a</b>) IR original image. (<b>b</b>) IR bias-corrected image. (<b>c</b>) Estimated bias.</p>
Full article ">Figure 9
<p>Artificial image pair creation algorithm: (<b>a</b>–<b>c</b>) EO images. (<b>d</b>–<b>f</b>) IR images.</p>
Full article ">Figure 10
<p>Dataset examples: (<b>a</b>–<b>c</b>) EO images. (<b>d</b>–<b>f</b>) IR images. (<b>g</b>–<b>i</b>) IR images with bias removed. (<b>j</b>–<b>l</b>) Pixel Fused images. (<b>m</b>–<b>o</b>) Pixel Fused images with bias removed.</p>
Full article ">Figure 11
<p>UAV recorded at twilight: (<b>a</b>) EO image. (<b>b</b>) IR image.</p>
Full article ">Figure 12
<p>Independent model detection and tracking on higher robustness target cases: (<b>a</b>) EO blurry UAV image. (<b>b</b>) IR blurry UAV image. (<b>c</b>) EO partially cut UAV image. (<b>d</b>) IR partially cut UAV image.</p>
Full article ">Figure 13
<p>Independent model detection and tracking on lower robustness target cases: (<b>a</b>) EO intra-class variation image. (<b>b</b>) IR intra-class variation image. (<b>c</b>) EO presence of birds image. (<b>d</b>) IR presence of birds image. (<b>e</b>) EO textured background image. (<b>f</b>) IR textured background image.</p>
Full article ">Figure 14
<p>Alignment failure on Pixel Fused images: (<b>a</b>) Vertical shift of input images to FusionGAN. (<b>b</b>) Significant vertical shift of input images leading to complete UAV overlap miss on Pixel Fused images.</p>
Full article ">Figure 15
<p>Data fusion detection and tracking on the intra-class variation target case: (<b>a</b>) EO-IR architecture. (<b>b</b>) IR-EO architecture. (<b>c</b>) Pixel-level fused architecture.</p>
Full article ">Figure 16
<p>Data fusion detection and tracking with the presence of birds target case: (<b>a</b>) EO-IR architecture. (<b>b</b>) IR-EO architecture.</p>
Full article ">Figure 17
<p>Data fusion detection and tracking on the textured background target case: (<b>a</b>) EO-IR architecture. (<b>b</b>) IR-EO architecture. (<b>c</b>) Pixel-level fused architecture.</p>
Full article ">
14 pages, 3706 KiB  
Article
The Characterization of Aquifer Parameters in Using Skimming Tubewells Through the Pumping Test Method: A Case Study of Tando Allahyar
by Xiufang Yang, Muhammad Uris Mirjat, Abdullah Baloch, Mashooque Ali Talpur, Shafi Muhammad Kori, Rajesh Kumar Soothar, Sher Ali Shaikh, Irshad Ali Mari and Farman Ali Chandio
Water 2024, 16(22), 3180; https://doi.org/10.3390/w16223180 - 6 Nov 2024
Viewed by 481
Abstract
Sindh is in the lower reaches of the Indus River; it is most vulnerable to a variety of upstream water development challenges. The aim of this research was to determine aquifer characteristics in the command area of Tando Allahyar-II distributary within the culmination [...] Read more.
Sindh is in the lower reaches of the Indus River; it is most vulnerable to a variety of upstream water development challenges. The aim of this research was to determine aquifer characteristics in the command area of Tando Allahyar-II distributary within the culmination of underground water potential. The hydraulic properties of the aquifer as well as the susceptibility of the formation to tedious extraction and saltwater upcoming were recognized. Three pumping tests were performed at head, middle, and tail reaches along the selected distributary. The drawdowns were measured at head reach (5.1667 h), at middle reach (6.0 h), and at tail reach (19.667 h) of the selected distributary by performing the pumping tests. Groundwater levels were lower at the tail reach compared to those at the head and middle reaches, likely due to a higher concentration of tubewells in the lower reach. The head and middle reaches showed higher groundwater levels, possibly due to constant head conditions promoting infiltration and recharge. The pumping test versus drawdown analysis revealed that the tubewells should be run with 7-h (on) and 4-h (off) operation. Further, the tubewells at all reaches (head, middle, and tail) should be closed for a minimum of 4 h between operations. This strategy would allow safe groundwater extraction, maintain water quality, and prevent water table depletion in the study area. The hydrodynamic and hydro-salinity behaviors were scrutinized in PWMIN 5.3 (version) by means of the MODFLOW mode. The results were estimated to compare the calibration and validation simulation outcomes using measured data. The model was successfully calibrated, and the root mean square (RMS) value of the head tubewell varied between 0.024 and 0.108, whereas it speckled between 0.0166 and 0.0349 for the middle tubewell and between 0.0659 and 0.0069 for the tail tubewell. The RMS values for hydrodynamic behavior for the head, middle, and tail reaches were less than 10%. These values represent a suitable match between the observed and simulated heads when a water table depletion of 1 to 2 m was observed due to extreme pumping. However, the average relative error values, for all validated procedures, were less than 10%. Full article
(This article belongs to the Special Issue Soil and Groundwater Quality and Resources Assessment)
Show Figures

Figure 1

Figure 1
<p>The red color line shows the command area map of Tando Allahyar-II distributary (blue line shows Naseer canal and light blue shows the Tando Allahyar-II distributary).</p>
Full article ">Figure 2
<p>Geological map of Tando Allahyar district.</p>
Full article ">Figure 3
<p>AQTESOLV results calculated by different testing methods at head reach.</p>
Full article ">Figure 4
<p>AQTESOLV results calculated by different testing methods at middle reach.</p>
Full article ">Figure 5
<p>Comparison of simulated and observed hydraulic heads at head reach tubewell.</p>
Full article ">Figure 6
<p>Comparison of simulated and observed hydraulic heads at middle reach tubewell.</p>
Full article ">Figure 7
<p>Comparison of simulated and observed hydraulic heads at tail reach tubewell.</p>
Full article ">
42 pages, 20876 KiB  
Article
Validating a Calibrated Model of a Groundwater Pump-And-Treat System Using Robust Multiple Regression
by Michael Rush, Leslie Gains-Germain, Lauren M. Foster and Tom Stockton
Water 2024, 16(22), 3178; https://doi.org/10.3390/w16223178 - 6 Nov 2024
Viewed by 517
Abstract
Validation of groundwater models is relatively challenging due to the need to reserve scarce water level data for calibration targets. In addition, traditional statistical validation metrics are unintuitive for non-technical audiences and do not directly identify model behaviors that require further refinement. We [...] Read more.
Validation of groundwater models is relatively challenging due to the need to reserve scarce water level data for calibration targets. In addition, traditional statistical validation metrics are unintuitive for non-technical audiences and do not directly identify model behaviors that require further refinement. We developed a novel model validation method that analyzes rate change events at pump-and-treat wells and statistically compares the water level responses at nearby monitoring wells between the data and model. The method takes advantage of events that occur alongside ambient pumping, unlike parameter estimation techniques that require independent drawdown or recovery events. The ability of the model to match well connections that are evident (or not evident) in the observations is characterized statistically, leading to four decision scenarios: model matches the observed connection (1) or lack thereof (2), model exhibits a connection that is not observed (3), or model over- or understates the observed connection (4). The method is applied to an FEHM-based groundwater flow and transport model that is shown to match 84.5% of the well connections analyzed. The method provides novel perspectives on the influence of calibration targets on the flow field and suggests that although the overall effect of drawdown targets was to improve the model, the choice of target well pairs creates flow pathways that may be inconsequential during normal operational conditions. The model adequately matches the flow over short spatial scales (<800 m) and over-represents the flow over greater distances (300–1200 m), suggesting the need for “null” drawdown targets in subsequent rounds of calibration. Full article
(This article belongs to the Section Hydrogeology)
Show Figures

Figure 1

Figure 1
<p>Example of CrEX-1 pumping event (<b>top</b>), corresponding CrPZ-2 monitoring well drawdown response in the data (<b>middle</b>) and in multiple model calibrations (<b>bottom</b>).</p>
Full article ">Figure 2
<p>Example of CrEX-1 pumping event (<b>top</b>), CrPZ-2 monitoring well water level recovery response in the data (<b>middle</b>) and multiple model calibrations (<b>bottom</b>).</p>
Full article ">Figure 3
<p>Flow chart describing well-pair classification scheme.</p>
Full article ">Figure 4
<p>Well-pair classifications based on <span class="html-italic">p</span>-value criteria; significance levels denoted as black lines and gray shading. Note that the effect of adjusting <math display="inline"><semantics> <mrow> <mi>α</mi> </mrow> </semantics></math> can be visually interpreted as shifts of the black line within the gray shading. Well pairs with data-model difference <span class="html-italic">p</span>-values greater than 0.01 are classified as matches (upper quadrants, 84.5%), with the model reproducing significant connections on the upper right (blue, 15.5%) and no connection on the upper left (green, 69%). Well pairs with data-model difference <span class="html-italic">p</span>-values less than 0.01 are classified as mismatched (lower quadrants, 15.5%), with the model overstating (purple, 6%) or understating (pink, 1.5%) a significant connection on the lower right and the model simulating connections that are not observed in the data on the lower left (orange, 7.5%).</p>
Full article ">Figure 5
<p>Well-pair classifications for CrEX-1 well pairs (<b>top</b>) and other wells (<b>bottom</b>); CrEX-2 and -3, CrIN-3, -4, and -5). CrEX-1 well pairs that were used to develop drawdown targets for the calibration are identified in cyan.</p>
Full article ">Figure 6
<p>Map-based representation of connection classifications for monitoring wells used in the development of drawdown targets during model calibration. Line colors summarize results across calibrations based on the most common outcome for the monitoring well to clarify the broad behavior of the model with respect to the data.</p>
Full article ">Figure 7
<p>Well-pair classifications based on <span class="html-italic">p</span>-value criteria against distance between monitoring and pumping wells. Instances of “Model Overly Connected” tend to occur at distances of 250–1750 m, whereas instances of “Model Matches Connection” occur at shorter distances (&lt;750 m).</p>
Full article ">Figure 8
<p>Map-based representation of response classifications for CrEX-1 well pairs. Similar to <a href="#water-16-03178-f006" class="html-fig">Figure 6</a>, line colors summarize results across calibrations based on the most common outcome for the monitoring well as a way to understand the broad behavior of the model with respect to the data. The general direction of flow is southeast (from the top left of the figure to the bottom right).</p>
Full article ">Figure A1
<p>Scatter plots and Huber regressions of pumping rate changes and monitoring well responses for CrEX-1 well pairs. Plot colors summarize results across calibrations based on the most common outcome for the monitoring well to clarify the broad behavior of the model with respect to the data. R-35a is excluded from the CrEX-1 analysis due to a low sample size (&lt;15 events). Each point on each panel represents <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math> for a single event for the given well pair. <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math> are unitless because they are normalized by their respective standard deviations.</p>
Full article ">Figure A2
<p>Scatter plots and Huber regressions of pumping rate changes and monitoring well responses for CrEX-2 well pairs. Plot colors summarize results across calibrations based on the most common outcome for the monitoring well to clarify the broad behavior of the model with respect to the data. Each point on each panel represents <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math> for a single event for the given well pair. <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math> are unitless because they are normalized by their respective standard deviations.</p>
Full article ">Figure A3
<p>Scatter plots and Huber regressions of pumping rate changes and monitoring well responses for CrEX-3 well pairs. Plot colors summarize results across calibrations based on the most common outcome for the monitoring well to clarify the broad behavior of the model with respect to the data. Each point on each panel represents <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math> for a single event for the given well pair. <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math> are unitless because they are normalized by their respective standard deviations.</p>
Full article ">Figure A4
<p>Scatter plots and Huber regressions of pumping rate changes and monitoring well responses for CrIN-3 well pairs. Plot colors summarize results across calibrations based on the most common outcome for the monitoring well to clarify the broad behavior of the model with respect to the data. Each point on each panel represents <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math> for a single event for the given well pair. <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math> are unitless because they are normalized by their respective standard deviations.</p>
Full article ">Figure A5
<p>Scatter plots and Huber regressions of pumping rate changes and monitoring well responses for CrIN-4 well pairs. Plot colors summarize results across calibrations based on the most common outcome for the monitoring well to clarify the broad behavior of the model with respect to the data. Each point on each panel represents <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math> for a single event for the given well pair. <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math> are unitless because they are normalized by their respective standard deviations.</p>
Full article ">Figure A6
<p>Scatter plots and Huber regressions of pumping rate changes and monitoring well responses for CrIN-5 well pairs. Plot colors summarize results across calibrations based on the most common outcome for the monitoring well to clarify the broad behavior of the model with respect to the data. Each point on each panel represents <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math> for a single event for the given well pair. <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math> are unitless because they are normalized by their respective standard deviations.</p>
Full article ">Figure A7
<p>Ratio of the standardized water level response (<math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math>) to the standardized pumping rate change (<math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math>) for CrEX-1 well pairs. Faint grey lines connect the same event to allow comparison of each individual event among the data and each of calibrations 1 through 4. The <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math>/<math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math> ratio is unitless because <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math> are normalized by their respective standard deviations.</p>
Full article ">Figure A8
<p>Ratio of the standardized water level response (<math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math>) to the standardized pumping rate change (<math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math>) for CrEX-2 well pairs. Faint grey lines connect the same event to allow comparison of each individual event among the data and each of calibrations 1 through 4. The <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math>/<math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math> ratio is unitless because <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math> are normalized by their respective standard deviations.</p>
Full article ">Figure A9
<p>Ratio of the standardized water level response (<math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math>) to the standardized pumping rate change (<math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math>) for CrEX-3 well pairs. Faint grey lines connect the same event to allow comparison of each individual event among the data and each of calibrations 1 through 4. The <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math>/<math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math> ratio is unitless because <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math> are normalized by their respective standard deviations.</p>
Full article ">Figure A10
<p>Ratio of the standardized water level response (<math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math>) to the standardized pumping rate change (<math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math>) for CrIN-3 well pairs. Faint grey lines connect the same event to allow comparison of each individual event among the data and each of calibrations 1 through 4. The <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math>/<math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math> ratio is unitless because <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math> are normalized by their respective standard deviations.</p>
Full article ">Figure A11
<p>Ratio of the standardized water level response (<math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math>) to the standardized pumping rate change (<math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math>) for CrIN-4 well pairs. Faint grey lines connect the same event to allow comparison of each individual event among the data and each of calibrations 1 through 4. The <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math>/<math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math> ratio is unitless because <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math> are normalized by their respective standard deviations.</p>
Full article ">Figure A12
<p>Ratio of the standardized water level response (<math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math>) to the standardized pumping rate change (<math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math>) for CrIN-5 well pairs. Faint grey lines connect the same event to allow comparison of each individual event among the data and each of calibrations 1 through 4. The <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math>/<math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math> ratio is unitless because <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>q</mi> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>∆</mo> <mi>h</mi> </mrow> </semantics></math> are normalized by their respective standard deviations.</p>
Full article ">Figure A13
<p>Histograms of residuals for Huber regressions of pumping rate changes and monitoring well responses for CrEX-1 well pairs. R-35a is excluded from the CrEX-1 analysis due to a low sample size (&lt;15 events). Residuals are unitless because the predicted values are water levels normalized by standard deviation.</p>
Full article ">Figure A14
<p>Histograms of residuals for Huber regressions of pumping rate changes and monitoring well responses for CrEX-2 well pairs. Residuals are unitless because the predicted values are water levels normalized by standard deviation.</p>
Full article ">Figure A15
<p>Histograms of residuals for Huber regressions of pumping rate changes and monitoring well responses for CrEX-3 well pairs. Residuals are unitless because the predicted values are water levels normalized by standard deviation.</p>
Full article ">Figure A16
<p>Histograms of residuals for Huber regressions of pumping rate changes and monitoring well responses for CrIN-3 well pairs. Residuals are unitless because the predicted values are water levels normalized by standard deviation.</p>
Full article ">Figure A17
<p>Histograms of residuals for Huber regressions of pumping rate changes and monitoring well responses for CrIN-4 well pairs. Residuals are unitless because the predicted values are water levels normalized by standard deviation.</p>
Full article ">Figure A18
<p>Histograms of residuals for Huber regressions of pumping rate changes and monitoring well responses for CrIN-5 well pairs. Residuals are unitless because the predicted values are water levels normalized by standard deviation.</p>
Full article ">Figure A19
<p>Residuals against fitted values for the CrEX-1 well pairs. R-35a is excluded from the CrEX-1 analysis due to a low sample size (&lt;15 events). Fitted values and residuals are unitless because the predicted values are water levels normalized by standard deviation.</p>
Full article ">Figure A20
<p>Residuals against fitted values for the CrEX-2 well pairs. Fitted values and residuals are unitless because the predicted values are water levels normalized by standard deviation.</p>
Full article ">Figure A21
<p>Residuals against fitted values for the CrEX-3 well pairs. Fitted values and residuals are unitless because the predicted values are water levels normalized by standard deviation.</p>
Full article ">Figure A22
<p>Residuals against fitted values for the CrIN-3 well pairs. Fitted values and residuals are unitless because the predicted values are water levels normalized by standard deviation.</p>
Full article ">Figure A23
<p>Residuals against fitted values for the CrIN-4 well pairs. Fitted values and residuals are unitless because the predicted values are water levels normalized by standard deviation.</p>
Full article ">Figure A24
<p>Residuals against fitted values for the CrIN-5 well pairs. Fitted values and residuals are unitless because the predicted values are water levels normalized by standard deviation.</p>
Full article ">Figure A25
<p>Map-based representation of connection classifications for CrEX-1 well pairs. Line colors summarize results across calibrations based on the most common outcome for the monitoring well to clarify the broad behavior of the model with respect to the data.</p>
Full article ">Figure A26
<p>Map-based representation of connection classifications for CrEX-2 well pairs. Line colors summarize results across calibrations based on the most common outcome for the monitoring well to clarify the broad behavior of the model with respect to the data.</p>
Full article ">Figure A27
<p>Map-based representation of connection classifications for CrEX-3 well pairs. Line colors summarize results across calibrations based on the most common outcome for the monitoring well to clarify the broad behavior of the model with respect to the data.</p>
Full article ">Figure A28
<p>Map-based representation of connection classifications for CrIN-3 well pairs. Line colors summarize results across calibrations based on the most common outcome for the monitoring well to clarify the broad behavior of the model with respect to the data.</p>
Full article ">Figure A29
<p>Map-based representation of connection classifications for CrIN-4 well pairs. Line colors summarize results across calibrations based on the most common outcome for the monitoring well to clarify the broad behavior of the model with respect to the data.</p>
Full article ">Figure A30
<p>Map-based representation of connection classifications for CrIN-5 well pairs. Line colors summarize results across calibrations based on the most common outcome for the monitoring well to clarify the broad behavior of the model with respect to the data.</p>
Full article ">
Back to TopTop