Environmental characterization of the testing system
Before the in vivo evaluation, the physical parameters of the simulated seaside chamber were calibrated to ensure high environmental fidelity. The UV simulation unit utilized a 6 W lamp emitting a broad spectrum ranging from 300 to 420 nm with a peak wavelength at 340 nm. This specific spectral distribution ensures a dominant UVA emission accompanied by a biologically relevant portion of UVB, which effectively mirrors the typical ultraviolet spectrum of natural midday sunlight reaching the Earth's surface17. The wave maker operated at a flow rate of 10 L/min, which generates a mild fluid shear stress comparable to gentle coastal wave washout. Additionally, the airflow unit provided a constant wind speed of 35 CFM, translating to a wind velocity of approximately 3 to 4 m/s over the test area. This correlates with a gentle breeze (Beaufort scale 3) typical of seaside meteorological conditions18. These combined parameters reliably replicate the cyclic mechanical and environmental stressors experienced during typical seaside activities.
Sunscreen retention assessment
First, the simulated seaside apparatus was assembled. As shown in Figure 1A, the system primarily consists of a blower fan, a UV irradiation system, a simulated seawater tank, and a control system, which are respectively used for wind-drying, UV irradiation, seawater immersion, and centralized control. After placing the mice in the apparatus and starting the system, the conditions of immersion and emersion phases on the shore were successfully simulated (Figure 1B). Based on this setup, we first investigated sunscreen retention under these conditions. Initially, commercially available Sunscreens A, B, and C, labeled with a fluorescent marker, were applied to the dorsal skin of the mice. Subsequently, these mice were subjected to an alternating cycle in the simulation apparatus, consisting of 30 min of immersion (simulating swimming) and 30 min of emersion phase (simulating sunbathing on the shore). After 2 h, following a complete cycle of wave and wind, the retained fluorescence on the skin was detected using in vivo imaging. As shown in Figure 2A–B, compared to pre-simulation, all sunscreen groups exhibited a decrease in fluorescence intensity, indicating that the wash-off effect of seawater waves causes a significant amount of sunscreen to be removed. Among them, the quadrant treated with Sunscreen C retained a relatively high amount of fluorescence, Sunscreen B showed a moderate retention rate, and Sunscreen A exhibited the lowest retention rate (Figure 2C). This indicates that Sunscreens C and B possess relatively stronger retention capabilities. These data confirm that our apparatus can quantitatively differentiate the durability of different formulations under complex environmental stress.
Photoprotective efficacy
To further investigate the UV-protective effects of these sunscreens in a simulated seaside environment, macroscopic observations of the dorsal skin of mice were conducted. As shown in Figures 3A and 3B, the unprotected control exhibited UV-induced damage, characterized by pronounced erythema, redness, and sunburn. In contrast, the regions treated with sunscreens showed significantly alleviated sunburn symptoms. Furthermore, histological analysis using H&E staining was performed and was consistent with the aforementioned macroscopic observations (Figure 4A–B). When the unprotected control was exposed to the seaside simulator, significant thickening of the epidermis and marked infiltration of inflammatory cells (Figure 4C) in the dermis were observed, demonstrating that the simulated seaside testing system effectively replicates skin sunburn in a practical seaside environment. The application of sunscreens significantly mitigated these pathological changes. Notably, the epidermis thickness and inflammatory cell infiltration of the skin tissue treated with Sunscreen C were remarkably similar to those of the UV-shielded internal control subjected to the same environmental stress. This indicates that Sunscreen C provides the optimal protective effect in the simulated marine environment, likely attributable to its high retention rate and inherent photoprotective properties.

Figure 1: Overview of the simulated seaside testing system. (A) Schematic diagram of the apparatus. (B) Mice inside the testing chamber. The customized test chamber integrates a UV irradiation device, a water tank with a wave generator, and a blower fan, effectively simulating a practical seaside environment. No copyright concerns are associated with this figure as it is obtained by the authors. Please click here to view a larger version of this figure.

Figure 2: Sunscreen retention fluorescence analysis. (A) Representative fluorescence images of the mouse dorsal skin before and after the 120-minute experiment. (B) Quantitative analysis of fluorescence retention. (C) Relative fluorescence intensity analysis. Data are presented as mean ± SD (n = 5). Statistical analysis was performed using one-way ANOVA followed by Tukey’s post-hoc test. *p < 0.05, **p < 0.01. Please click here to view a larger version of this figure.

Figure 3: Macroscopic appearance of mouse skin after environmental exposure. (A) Representative macroscopic images of the mouse dorsal skin after a 120-minute exposure to the simulated seaside environment. (B) Locally enlarged views showing the occurrence of erythema and peeling in the unprotected control area compared to the sunscreen-protected areas and the UV-shielded internal control (Normal). Please click here to view a larger version of this figure.

Figure 4: Histological evaluation of photoprotective efficacy. (A) Representative H&E-stained sections of the mouse dorsal skin from different treatment areas following environmental exposure. Indicator lines identify the boundaries of the measured epidermis. Scale bar = 250 µm. (B) Quantitative analysis of epidermal thickness in each designated area. (C) Quantitative assessment of relative inflammatory cell counts (%). The UV-shielded internal control (Normal) was set as the baseline (100%). Data are presented as mean ± SD (n = 5). Statistical analysis was performed using one-way ANOVA followed by Tukey’s post-hoc test. *p < 0.05, **p < 0.01, ***p < 0.001. Please click here to view a larger version of this figure.