Next Article in Journal
Investigation of Ocean Sub-Surface Processes in Tropical Cyclone Phailin Using a Coupled Modeling Framework: Sensitivity to Ocean Conditions
Previous Article in Journal
Clinical Observations Associated with Phenobarbital Serum Monitoring to Manage Epilepsy in a California Sea Lion with Domoic Acid Toxicosis
Previous Article in Special Issue
Impact of Egg Exposure to UV Filter-Spiked Sediment on the Survival, Hatching Success, Cardiac Frequency, and Metabolic Scope of Zebrafish Embryos
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Analysis and Modeling of Sunscreen Ingredients’ Behavior in an Aquatic Environment

1
Green Engineering & Resources Research Group (GER), Departamento de Química e Ingeniería de Procesos y Recursos, ETSIIT, Universidad de Cantabria, Avda. de los Castros 46, 39005 Santander, Cantabria, Spain
2
Departamento de Química Analítica, Facultad de Ciencias del Mar y Ambientales, Universidad de Cádiz, Campus Universitario Río San Pedro, 11519 Puerto Real, Spain
3
Departamento de Ecología y Gestión Costera, Instituto de Ciencias Marinas de Andalucía (CSIC), Campus Universitario Río San Pedro, 11519 Puerto Real, Spain
*
Author to whom correspondence should be addressed.
Oceans 2022, 3(3), 340-363; https://doi.org/10.3390/oceans3030024
Submission received: 11 May 2022 / Revised: 14 July 2022 / Accepted: 28 July 2022 / Published: 2 August 2022
(This article belongs to the Special Issue Skincare Chemicals and Marine Life)

Abstract

:
Sunscreens have become a product based on increasingly complex formulations that include, among many ingredients, a mixture of UV filters to provide optimal sun ultraviolet radiation protection. A significant group of scientific works deals with the impact of UV filters in aquatic media. However, the knowledge of the mechanism and kinetics of the compound’s direct release, fate, and its transformation and interaction with living organisms is necessary to assess its environmental occurrence and behavior and to predict potential and real impacts on the aquatic environment. This review outlines the existing analysis and modeling of the release and behavior of sunscreen’s ingredients in the marine environment, including aquatic organisms. The physical-chemical properties, photodegradation, and release kinetics of particles and chemicals into the water are studied by hydrodynamic and kinetic models. Direct photolysis of chemicals is modeled as pseudo-first-order kinetics, while the indirect pathway by the reaction of sunscreen with reactive oxygen species is described as second-order kinetics. The interaction of UV filters with marine biota is studied mainly by toxicokinetic models, which predict their bio-accumulation in the organisms’ tissues. These models consider the chemicals’ uptake and excretion, as well as their transfer between different internal animal organs, as a first-order kinetic process. The studies analyzed in the present work represent a driver of change for the beauty and personal care industry, in order to seek new ecological alternatives through the application of R&D tactics.

Graphical Abstract

1. Introduction

The growing awareness of beauty and personal care has increased in recent years, reaching historical growth in the products used in this area manufactured by the organic fine chemical sector. Despite the fact that theglobal cosmetic market decreased in 2020 by eight percent, in comparison to the previous year, as a consequence of the COVID-19 pandemic, skincare was the leading category of cosmetic products accounting for about 42% of this global market. Skincare is one of the most profitable products, as its revenue is projected to generate roughly USD 177 billion in 2025. Skincare encompass a large number of products, such as exfoliants, mineral-based facial treatments, serums, oils, or sunscreens. Solar protection is one of the key drivers behind the fast-growing suncare global market, which generated a revenue of USD 9.3 billion in 2020. The global market for sunscreen ingredients is forecasted to grow from USD 691.6 million in 2021 to USD 787.0 million by 2026, at a compound annual growth rate of 2.6% for the 2021–2026 period [1,2].
Sunscreens are among the many emerging pollutants that enter the sea and cause adverse ecological effects [3,4]. In recent years, the environmental effects of some ingredients used in the preparation of sunscreens once they are released into the aquatic environment has been discussed by the scientific community [5,6]. The increase in tourism associated with the sun and coast has increased the presence of these products in the marine environment and, consequently, the concern about the potential dangers to marine ecology, which are mainly caused by the UV chemical filters, both organic and inorganic, present in the sunscreens [5,7,8]. Following the verification of these adverse effects on marine life and ecosystems, specific ingredients used in these products have been banned in some countries and regions, mostly coastal and intensively touristic areas (e.g., Hawaii, Palau, and Thailand).
The analysis and modeling of the release, behavior, and fate of ingredients contained in sunscreen products when they come into contact with seawater allow us to predict their potential negative impacts on the marine environment. The development of a general modeling method, which allows us to model these complex (from an environmental point of view) matrices based on chemical cocktails, requires the minimization and management of the effects caused by the use of sunscreens in aquatic ecosystems. In the present work, we first present a brief perspective on sunscreen products. Secondly, we critically review our current knowledge of models studied in the literature regarding the behavior and transformation of these environmental harmful substances and their associated ingredients in aquatic media. This review may represent an advance in the current state of knowledge, given that it adds a different point of view on sunscreen release in aquatic ecosystems, investigating the chemical behaviour of sunscreen compounds depending on the experiments performed (in-lab, in-field, with organisms, without organisms). The study also reports the mechanisms and kinetic models used, as well as the requirement for the more precise and realistic modeling of skincare ingredients and, in general, for all the emerging contaminants.

1.1. Sunscreens: Function, Market and Legislation

Sunscreens are cosmetic products belonging to the skincare products group. Their main function is to protect individuals, in both short and long terms, from solar UV radiation that the human skin receives. These products, primarily designed to be used in high-UV-radiation-exposure situations, have evolved into multifunctional products capable of preventing skin aging due to their moisturizing effect [9].
The UV filter market dominates the sector, as UV filters are an ingredient used in cosmetic products as well as being used exclusively in the formulation of sunscreen. Higher-UV-protection products are demanded by both sectors, being a blend of UV filters the basic ingredients that provide this protection. From a global market volume of 27,000 tons of UV filters in 2016, an annual growth rate of 4% is expected [10].
The complex formulations of sunscreens, their use, and the possible release into the environment of the ingredients that make up their photodegradation and aging properties make it necessary to develop and apply some type of legislation or regulation to both the manufacturing and use of these products. In recent years, legislation proposals have been developed for these components, but not all of them have been enacted due to the lack of adequate research that demonstrates their harmful effects on the environment. Several governments (e.g., the states of Hawaii and Florida (Key West), US Virgin Islands, Palau, and Bonaire) have approved legislations banning the sale of sunscreens containing some UV organic filters, such as oxybenzone and octinoxate, due to their reported effects on aqueous environments and potential for bleaching coral reefs. Additionally, a recent study [11] conducted in the Mexican Caribbean highlighted the need to create more restrictive regulations on sunscreens. Table 1 includes IUPAC chemical name, chemical structure, regulation, and main properties of the most commonly used UV filters, including 10 organic and 2 inorganic compounds.

1.2. Design and Manufacturing of Sunscreen

There is a large variety of possible UV filters to be used in the formulation of a sunscreen. These are divided into two groups: organic (synthetic) and inorganic (minerals). Although the formulation of each product is different, all of them combine several UV filters, whether organic or inorganic. This is because the maximum concentrations allowed are lower than those that a protector would require to provide effective protection against UV radiation. In addition, by using different types of UV filters, a broader spectrum of protection is achieved [12]. The most common UVA organic filters are benzophenones 1, 2, 3, and 4, which are aromatic ketones only effective against short-wave UVA radiation. Benzophenones 4, 5, and 9 have sulfur and sodium atoms in their formulation. Ecamsule, avobenzone, and meradimate (with nitrogen atoms included in their compositions) are other commonly used filters. An important group of UVB filters are cinnamates (found in 90% of sunscreens). Octyl methoxycinnamate (OCM) is widely used, but presents an volatility to sun exposure. p-aminobenzoic acids (PABAs) and derivatives, salicylates, camphor derivatives, octocrylene, ensulizol, and urocanic acid are other groups of important UVB filters.
Inorganic filters used in sunscreen formulations are titanium dioxide and zinc oxide and, in turn, can be included as nanoparticles. The photocatalytic properties of nanoparticles cause their photopassivation in order to prevent any type of oxidative stress on the skin or the product itself from losin its protective characteristics [13,14]. In general, this photopassivation is achieved by the precipitation of an inert mineral layer of aluminum hydroxide or silica on the surface of the nanoparticles [13]. Another way to increase the photostability of TiO2 nanoparticles (TiO2-NPs) is by “doping” the particle with certain metals, such as vanadium, manganese, iron, chromium, or cobalt [15].
Sunscreen is a product classified within the organic fine chemical and pharmaceutical industry where manufacturers produce a wide range of chemical substances, which are typically of a high added value. The most important task in their manufacturing process by recipe-driven batch processes in multipurpose plants is to achieve compliance with all the desired functional characteristics and fulfil the requirements for the approval of the product by different administrations.
The different combinations of the formulated ingredients focus on obtaining a specific sun-protection factor, but also on meeting the needs of a specific group of consumers; in this sense, cream, lotion, gel, stick, or aerosol media allow different forms of application for different people. All of this gives rise to a complex matrix where, in addition to UV-filter combinations, aqueous and organic dispersing agents, emulsifiers, antioxidants, humectants, preservatives, stabilizers, and antimicrobial agents are included in the formulations.
When sunscreens are used, part of the many ingredients used in the sunscreen’s formulation are released into the surrounding environment, especially into the marine environment [8]. This release depends on the sunscreen ingredients’ properties, the receiving medium characteristics, and a series of physical–chemical interactions between both of them [16,17,18,19,20]. All these processes can cause the ingredients to remain unchanged in the medium, change their speciation, or even transform into other compounds at variable rates [21,22,23,24,25,26]. Once in the environment, the original or transformed ingredients can interact directly or indirectly with the environment and marine life, causing a series of adverse effects [27,28,29,30,31,32]. Different mechanisms and mathematical models are proposed in the present study to analyze and predict the compound’s release, its transformation, and interaction with the environment and living organisms, e.g. [33,34,35,36,37,38].

2. Analysis and Modeling of the Behavior of Sunscreen’s Ingredients in an Aquatic Environment in the Absence of Organisms

The behavior, transformation, aging, and modeling of sunscreen’s ingredients in the aquatic environment in the absence of organisms were identified through a systematic literature review. Advanced oxidation processes (AOPs), that is, processes that were not spontaneous in nature, were not considered. Relevant keywords were combined to form a broad search conducted through the ScienceDirect and Scopus databases. The search was limited to the papers published in the year 2000 and onwards. This process identified 25 published articles that fulfiled the objective of this section. Figure 1 classifies the publications founded according to the scale and ingredients studied. Table 2 summarizes the modeling approach used in the different studies reviewed in Section 2.

2.1. Modeling of the Sunscreen Ingredients’ Behavior in Field Scenarios

A large number of studies focus on the concentration and fate of sunscreen ingredients in different environment scenarios that are commonly used for recreational activities, such as rivers, lakes, estuaries, and beaches [8]. Substance or material flow analysis (SFA or MFA); particle flow analysis (PFA); and stochastic Lagrangian and hydrodynamic models, which are sometimes conducted together, have been previously used to model the behavior of sunscreen UV filters in the field [11,34,37,39,40,41,42,43].
Material flow analysis models describe the accumulation of the studied compound in different compartments, both productive and environmental, considering each compartment as a mixer (MFA models) or considering the probabilistic behaviors (PMFA models) [35,36,37,38,39,40,41,43]. Hydrodynamic models (HMs) describe the transport, diffusion, and transformation of the studied compounds in a natural aquatic environment, generally in marine environments, such as a gulf, beach, bay, or an estuary [19,36,40].
Recently, Zheng and Nowack [39] integrated the size distribution of the particle flows into the MFA model by a size-specific, dynamic, probabilistic, MFA (ss-DPMFA) model that computes size information about the TiO2-NPs flows released into the environment. The results obtained from this study could be used as inputs for environmental multimedia fate models.
Lindo-Atichati et al. [40] incorporated the fate and transport of 4-MBC, OC, and BP-4 UV filters into a stochastic Lagrangian model that was coupled with a high-resolution hydrodynamic model; this study generated the expected trajectories of water parcels that transport chemicals between pollution sources (wastewater-treatment plants) and mussel rafts that were exposed to pollution. The predicted concentrations of organic UV filters in the coastline and at the outfalls of urban wastewater-treatment plants were provided.
Arvidsson et al. [41], using a global approach, developed a particle-flow-analysis (PFA) methodology to determine the emissions of TiO2-NPs from several industrial products, such as paint, sunscreen, and self-cleaning cement. Their results indicate that the current, highest emissions of TiO2-NPs originate from the use of sunscreen.
The approach of Johnson et al. [19] was a combination of measurements and modeling of the fate of titanium in the influent and effluent in an active sludge plant. Their study calculated river-water concentrations using hydrological models estimating the magnitude and variability of flows across a catchment in combination with a range of water-quality models, including a catchment scale water-quality model. The range of concentrations obtained included other predicted concentration ranges using probabilistic material flow analysis to model the concentrations of nanoparticles in specific environments, including sediment [34,35,36,37,42,43]. In addition to these works, the review conducted by Yuan et al. [6] highlighted the lack of information on the environmental fate of inorganic UV filters in the real environment.

2.2. Modeling of the Sunscreen Ingredients’ Behavior in Lab-Scale Studies

Before analysing the laboratory-scale experiments, we highlighted the works of Tovar-Sánchez et al. [3] and Sánchez-Quiles and Tovar-Sánchez [8], which combined laboratory and field experiments. These works demonstrated that transformations of sunscreen nanoparticles in seawater and under solar radiation produce compounds and reactive oxygen species (ROS) as singlet oxygen (1O2), hydroxyl radicals (OH), and superoxide radicals (O2). They suggested that sunscreens in coastal waters may have deleterious effects on the coastal ecosystem, either by inhibiting the growth of some species of marine phytoplankton or by adding essential micronutrients (N, P, and Si compounds) that may stimulate the growth of others.
Laboratoryexperiments aim to propose reaction mechanisms and kinetic models to describe the physico-chemical transformations of sunscreen ingredients under controlled conditions. The results of such studies can be inputs to the model applied to real environmental scenarios. In these studies, aqueous solutions of the target chemical were prepared by dissolving individual sunscreen components directly in ultrapure water. Additions of selected concentrations of ions (nitrate, bicarbonate, or chloride) and organic matter are common to examine how their presence may influence the transformations of the studied chemicals [44,45,46,47,48]. In addition to ultrapure water and synthetic seawater, other water matrices, such as freshwater, estuarine, and seawater samples from bathing areas, as well as water from swimming pools, were used as bulk solutions in laboratory-scale experiments [24,31,48,49,50]. Laboratory experiments were conducted by submitting the liquid mixture to magnetic stirring under simulated sunlight mimicking daily natural cycles using irradiation from laboratory UV lamps. Water evaporated under irradiation was supplemented regularly in some experiments, while controlled temperature representative of seawater was maintained by external cooling in others [5,16,51].
Three sets of laboratory studies can be distinguished: (i) studies that focused on the photodegradation kinetics and mechanism of sunscreen UV organic filters; (ii) works studying the release, fate, transformation, and aging in water of sunscreen ingredients and their degradation products; and (iii) works that highlighted the mobilization and photoproduction of inorganic nutrients (i.e., SiO2, P-PO43−, and N-NO3) from sunscreen products. These works are analyzed in the following sections in the present work with Figure 2 showing the modeled transformations related to the release and transformation of the sunscreen ingredients in the aquatic environment in laboratory-scale studies.

2.2.1. Photodegradation of Sunscreen UV Organic Filters under Simulated Natural Conditions

Several works [24,25,32,44,45,47,48,50,52] modeled and determined the photodegradation kinetics in laboratory-simulated environmental conditions of organic UV filters as benzophenone-derived compounds (benzophenone-3 (BP-3), benzophenone-4 (BP-4), 4-hydroxybenzophenone (H-BP), 2-hydroxy-4-methoxybenzophenone (HM-BP), 2,4-dihydroxybenzophenone (DH-BP), and 2,2’-dihydroxy-4-methoxybenzophenone (DHM-BP)), and the acid agents p-aminobenzoic acid (PABA), octyl-dimethyl-p-aminobenzoic acid (ODPABA), and 2-phenylbenzimidazole-5-sulfonic acid (PBSA).
Photodegradation processes may occur in water by direct and/or indirect photolysis. The direct photolysis pathway implies self-decomposition after the organic UV filters absorb UV radiation, following a pseudo-first-order kinetics; this model allows a manageable solution for the kinetic system leading ro an apparent or observed kinetic constant. The indirect photolysis pathway involves different photosensitizer agents, such as dissolved organic matter (DOM), reactive oxygen species (ROS), or other inorganic substances (e.g., nitrates, bicarbonates, NaCl, or TiO2), usually modeled as a second-order reaction rate.
The transformation-rate constant is a strong function of significant environmental parameters, such as pH [45,52], as well as dissolved organic matter, nitrate, and chloride concentrations in water [44,48,53], which slow down the transformation rates of organic UV filters [48]. The above-mentioned works performed experiments with synthetic solutions of sunscreen ingredients, and therefore did not take into account the chemical release process from commercial sunscreens into an aqueous medium. However, the works of Sendra et al. [31] and Tovar-Sánchez et al. [54] studied the photolysis of organic UV filters, starting from solutions of commercial sunscreens in seawater. Although they did not propose reaction models, both works showed the evolution of H2O2 generation as a consequence of the degradation of both TiO2-NPs and organic filters.
In addition to modeling the UV filters’ transformations under simulated natural conditions in thelaboratory, there are works that modeled the degradation of UV filters through the named advanced oxidation processes (AOPs). In these treatment processes, more extreme conditions were used in order to achieve maximum degradation yields. However, the kinetic modeling proposed in these cases can be useful to compare with natural processes [21]. AOP processes for the removal of sunscreen UV filters in water established pseudo-first-order-rate degradation, obtaining an apparent kinetic constant, and second-order-rate constants of the reaction between UV filter and reactive species [49,55,56,57,58,59,60,61,62].

2.2.2. Degradation and Aging of TiO2-NPs Used in Sunscreen

Sunscreens consist of a complex oil-in-water or water-in-oil emulsion in which UV filters are incorporated. Titanium dioxide nanoparticles (TiO2-NPs) are usually functionalized with combinations of surface coatings, such as inorganic oxides (SiO2, Al2O3), poly-dimethylsyloxane (PDMS), or stearic acid, to satisfy the practical application as an inorganic UV filter being a complex matrix.
Several authors studied the physicochemical processes related to TiO2-NPs degradation and aging in aquatic media during a contact time of up to seven days using aqueous solutions of TiO2-NPs [16,20,53,63] or aqueous solutions of commercial sunscreen containing TiO2-based nanocomposites [51]. The results obtained show that large parts of TiO2-NPs are released from the sunscreen with fast degradation, losing their coating and, therefore, their protective function. The kinetics and thermodynamics of the global release and aging process occurs under a set of processes, with complex and interrelated mechanisms, which depend on the type of TiO2-NPs and the properties of the aquatic environment. TiO2-NP transformations as dissolution, adsorption–desorption, aggregation/agglomeration, sedimentation, sulfidation, mineralization, redox reactions (e.g., photo-oxidation and photoreduction), and interaction with macromolecules can affect its fate, toxicity/bioavailability, and persistence in the environment [7,17,23]. In these processes, the characteristics of TiO2-NPs and the environmental conditions play a crucial role. Properties, such as the crystalline structure, morphology, size distribution, surface charge, type of coating materials and stabilizing agents (SiO2, Al(OH)3, PDMS, and/or stearic acid), must be taken into account, together with the characteristics and conditions of the aqueous medium simulating natural-environment conditions (freshwater, seawater, estuarine water), such as ionic strength; pH; temperature; light; and the presence of ingredients, such as natural organic matter, macromolecules, radicals, and surfactants.
Several works related to the aging and transformation of TiO2-NPs [16,20,51,53,63] proposed the agglomeration of nanoparticles during the aging process in stable colloidal residue that could include macroscopic aggregates, agglomerates, and submicronic fractions; these colloidal residues interact with the organic phase of the sunscreen, and the variety of ingredients in the sunscreen can be released into the aqueous environement. Metals, elements associated with Ti or Zn nanoparticles, as well as their coatings and metal impurities from raw materials [64,65,66] can partially dissolve in water and partially remain as nanoparticles associated with organomineral agglomerates and aggregates.
Furthermore, Auffan et al. [63] and Slomberg et al. [16] showed that after 48 h of TiO2-NPs aging in both fresh- and seawater, the protective SiO2 layer was almost completely degraded with the high release of Si; however, Al(OH)3 was more stable with a low release of Al into the water. A review by Peng et al. [67] discussed the environmental behavior of metal-based NPs with an in-depth analysis of the mechanisms and kinetics of aggregation, dissolution, and ROS generation; the review analyzed proposals, such as the Arrhenius-based kinetic model for NP dissolution, interaction force boundary layer theory approaches for describing adsorption, and different diffusion-limited and reaction-limited aggregation models, to predict the aggregation kinetics of NPs in solutions.

2.2.3. Inorganic Nutrients of Sunscreen

Organic UV-filter degradation in aqueous media in the presence of TiO2-NPs and light, as well as the degradation of nanoparticle surface coatings, can generate nutrients as nitrogen compounds (NO2, NO3, NH4+), phosphorous compounds (PO43−), and SiO2, which can have a considerable environmental impact [3,5,68].
The photocatalytic degradation of nitrogen-containing compounds from sunscreen can result in nitrate, nitrite, ammonia, and nitrogen production, whose relative concentrations and kinetic processes depend mainly on (i) the nature of nitrogen-containing organic compounds (i.e., the functional group in which nitrogen atoms are presented, the initial oxidation state of nitrogen, the nitrogen-existing type in compounds, and the number and position of nitrogen atoms); (ii) the TiO2-NPs properties as particle-size distribution and/or crystalline structure; and (iii) the reaction conditions, such as the light intensity and duration, concentration of organic compounds, TiO2-NPs concentration, and pH [69,70].
The study conducted by Low et al. [69] presented the formation of ammonia (NH3) with a basic pH or ammonium ion (NH4+) in acidic conditions; both species transformed into nitrite ions (NO2) and, later at much lower rates and following a first-order reaction, into nitrate ions (NO3), with chloride and sulphate ions being the inhibitors of the reaction. Jing et al. [70] classified different scenarios according to the types of nitrogen groups in organic compounds of sunscreen: amino groups were predominantly mineralized into NH4+, while the nitro group was mainly converted into NO3. For the heterocyclic structures, the nitrogen atom was transformed into either one or both of the NH4+ and NO3 species. N2 was mainly generated from the photodegradation of the –N=N– double-bond moieties. In relation to the kinetic modeling of the nitrogen compounds’ transformation, Jing et al. [70] suggested zero-order kinetics at high concentrations (5 × 10−3 mol L−1) because diffusion was the dominant process; in contrast, at lower concentrations, the degradation efficiency was proportional to the substance concentration in accordance with an apparent first-order kinetics.
Moreover, Ji et al. [50] indicated that the sulfonic moiety from 2-phenylbenzimidazole-5-sulfonic acid (PBSA) was primarily converted to SO42−, while nitrogen atoms were predominantly released as NH4+ and, to a lesser extent, as NO3. Along the same lines, Abdelraheem et al. [61] showed that the mineralization of hetero N-atoms from PBSA and 1H-benzimidazole-2-sulfonic acid (BSA) were released in solution as NH4+. Ammonia could be oxidized into NO3 by hydroxyl radicals generated from visible light/TiO2 present in the aqueous media or resulting from the autocatalyzed photo-oxidation process of NH4+ by UV/TiO2.
In addition to the nitrogen compounds, organic sunscreen ingredients with phosphorous, as well as phosphorous impurities, from TiO2-NPs, can release species of phosphorous as inorganic nutrients into the aquatic environment with a probably significant negative environmental impact. Liu et al. [68] studied the phosphorus-release kinetics from five types of TiO2-NPs under the influence of varying solution chemistries. Similar to the previous literature, the authors reported a rapid, initial decline (2 h) in the amount of adsorbed phosphate on the surface, followed by an extended slow desorption rate, and after phosphorous dissolution, re-adsorption on the particulate surface was reported. This interaction between aqueous PO43− and nanoparticle surfaces (ZnO-NPs) was studied by Lv et al. [71], who showed the dissolved Zn2+ formed nanoparticles of zinc (ZnO-NPs) complexed with aqueous PO43−, resulting in the formation of amorphous zinc phosphate hydrate precipitate. In both the studies previously mentioned, the introduction of phosphorus into the aquatic environment increased the concentration of nutrients, in addition to being bioavailable, generating an impact on the environment.
Similar to nitrogen and phosphorus, silicon compounds can release species containing this element, constituting a new nutrient element in the medium. Liu et al. [72] investigated the release kinetics of Si impurity frequently introduced during NP synthesis, and the resulting effect on TiO2-NPs transformation in aqueous solutions. The Si-release kinetics follow the parabolic kinetics mode, which is similar to the diffusion-controlled dissolution of minerals. Slomberg et al. [16] evaluated the transformation of the UV TiO2-NPs filter with a surface SiO2 coating, showing a high level of SiO2-layer dissolution after 48 h in fresh and marine waters.
On the other hand, Rodríguez-Romero et al. [5] studied the combined release of metals and nutrients from a commercial sunscreen. A pseudo-first-order kinetics model was proposed to describe the metal-release process, considering that metals are released from the organic phase of the sunscreen into the aquatic environment and can then be adsorbed on the stable colloid formed by the aging of the initial sunscreen. The photodegradation of organic compounds from sunscreen with phosphorus and silica, releases phosphates and silica into the seawater, which is described as first-order kinetics.
Table 2. Models and the modeling approach used in the different reviewed studies described in Section 2.
Table 2. Models and the modeling approach used in the different reviewed studies described in Section 2.
ModelModeling Approach and Model ExpressionsChemicalRef.
Sunscreen ingredients’ modeling in field scenarios
Substance flow analysis (SFA), material flow analysis (MFA), probabilistic and dynamic probabilistic MFA (ss-DPMFA). Particle-flow analysis (PFA).
  • Mass-balance multi-compartment models.
  • Particle number is used as a flow and stock metric instead of mass.
  • Probabilistic material flow analysis to model nanoparticle concentration in compartments.
TiO2-NPs[19,35,36,37,38,39,40,41,43]
Hydrodynamic models (HMs).Hydrodynamic numerical models describing and estimating the velocity field of three-dimensional currents and turbulent vertical mixing. Combination with water-quality models.TiO2-NPs, 4-MBC, OC, BP4[19,40]
Kinetic model of laboratory-scale photodegradation of sunscreen UV organic filters
  • Direct photolysis: pseudo-first order.
  • Indirect photolysis: second order.
  • Langmuir–Hinshelwood model.
r i , j = k i , j   C i
r i , j = k i , j   C i 2
r i , j = k i , j   K   C i 1 + K   C i
BP3, BP4, H-BP, HM-BP, DH-BP, DHM-BP, PABA, OD-PAB, PBSA, 4-MBC[24,25,32,44,45,47,48,50,52]
Degradation and aging of TiO2-NPs used in sunscreen
Several kinetic models for NP dissolution coupled with adsorption and aggregation.Pseudo-first- and second-dissolution kinetic orders. Diffusion or reaction limited agregation models.TiO2-NPs[16,20,51,63,67]
Inorganic nutrients from sunscreen
Pseudo-first-order kinetic model.Zero-order kinetics at high nitrogen concentrations. r i , j = k i , j   C i
r i , j = k i , j
Organic UV filters with N2 and P compounds. P, SiO2impurities/coatings.[5,50,68,69,70,71,72]

3. Analysis and Modeling of the Behavior of Sunscreen Ingredients in Contact with Aquatic Organisms

The mechanism and mathematical models proposed to analyze and predict the compound’s release, transformation, and interaction with living aquatic organisms were identified through a systematic literature review. Figure 3 presents the number of published articles that fulfil the objective of this section, collected according to the type of substance (UV filters or metals), differentiating between field (in green) and laboratory studies (in blue). It also presents the number of publications according to the aquatic organism studied. The most studied organisms were bivalves, which are good biomonitors of coastal pollution because they are filter feeders and exposed to high volumes of water for feeding and respiration [73,74]. Furthermore, their commercial and ecological importance in coastal environments make bivalves an excellent potential target species for assessing the effects of marine pollution [75].
The research conducted on metal accumulation is more extensive than for UV filters (Figure 3). High concentrations of metals in the marine environment near industrialized areas have led to an interest in their bioaccumulation and possible influence on the food chain, as well as the availability of tools for determining their effects. However, existing studies with UV filters have focused more on determining the toxic effects on aquatic organisms [67,76,77,78,79,80].
Due to the increasing use of sunscreens, their chemical components have reached the aquatic environment by direct sources due to various leisure and recreational activities (e.g., swimming, diving and wave-riding, and/or indirect sources through effluents from wastewater-treatment plants [81,82,83]). Once in the aquatic environment, these compounds can bioaccumulate in living organisms. Furthermore, this bioaccumulation can lead to biomagnifications where species higher in the food chain may be exposed to all the chemicals that lower-order species accumulate [84].
Organic UV filters are expected to accumulate in aquatic biota due to their lipophilic character and solubility in water. Many organic UV filters have high lipophilicity, which makes them susceptible to accumulate in the organisms of the biota, being stored somewhat faster than they are metabolized or excreted, thus primarily fixing them in soil and sediments. These matrices are considered a sink for organic pollutants in aquatic ecosystems [27].
Additionally, inorganic filters cause environmental toxicity and enzymatic changes in aquatic organisms. Several studies identified and evaluated the toxicity and ecotoxicological effects (different biological responses) of TiO2 and ZnO nanoparticles on different aquatic species [28,30,76,77,85]. Both of them can bioaccumulate and cause oxidative stress in clams [31,77].
In general, the resaerch collected from the literature study the harmful effects on living organisms of UV filters (only one filter or mixtures of several filters prepared in the laboratory). Only a few studies on the effects of sunscreens on marine organisms have been conducted using commercial sunscreen products, e.g., [29,79,86,87].
In environmental studies, the modeling itself is used for quantification, synthesis, and prediction (or hypothesis testing). Bio-toxicokinetic models, based on the description of the bioconcentration rate of chemical substances in living organisms, play an important role in driving the field and provide a greater understanding of the bioaccumulation and toxicity of chemical compounds in aquatic organisms. Therefore, there have been many attempts to develop bio-toxicokinetic models to predict the bioaccumulation, biotransformation, and toxicity of metals in aquatic organisms, the studies dealing with organic compounds being scarce, such as organic UV filters. This is probably because these organic compounds in organisms can be digested and transformed into different ones, which in some cases can be more toxic to the organism than the substance ingested. These transformations make modeling difficult, especially in organisms higher up the food chain.
Only two studies have presented a biokinetic model to describe the bioconcentration of organic UV filters. Gomez et al. [88] proposed a model to describe the bioconcentration of pharmaceutical and personal-care products (organic UV filters 2-ethyl-hexyl-4-trimethoxycinnamate (EHMC) and octocrylene (OC)) and benzodiazepines (diazepam and tetrazepam) in the marine mussel Mytilus galloprovincialis. Vidal-Liñán et al. [78] modeled the bioaccumulation kinetics of UV organic filters (4-methylbenzylidene camphor (4 MBC), benzophenone-3 (BP3), benzophenone-4 (BP4), OC, and octyl dimethyl p-aminobenzoate (ODPABA)) also using M.galloprovinciallis, suggesting that lipids are the final destination of the accumulated chemical substances.
However, sunscreens also contain trace metals and inorganic UV filters (TiO2-NPs/TiO2, ZnO-NPs/ZnO). The stability of inorganic filters in the aquatic medium is the subject of numerous studies, e.g., [72,89,90,91]. ZnO nanoparticles in the aquatic environment are transformed into their ionic form (Zn2+) [92,93]. Both ZnO-NPs and Zn2+ have been suggested to induce relevant toxicity effects, highlighting the release of Zn2+ ions from ZnO-NPs as important contributors to toxic effects [22,94,95]. On the other hand, Fan et al. [96] modeled the bioconcentration of six different types of TiO2-NPs in freshwater cladoceran Daphnia magna, observing significant effects on the exposure concentration and property of the nanoparticles, concluding that the studied nanoparticles could be considered as highly bioaccumulative.
Metals from sunscreen, once released into the aquatic environment, can bioconcentrate in organisms following the same processes as dissolved metals in the marine environment and, therefore, their bioconcentration can be described with models similar to those used in the bioaccumulation studies of dissolved metals obtained from sources other than sunscreen or just from synthetic solutions.

Biokinetic and Toxicokinetic Models

Bioaccumulation is an important framework in ecotoxicology research, serving as a link between the environment and organisms, being of crucial importance to assess the chemical status of water masses [97,98]. Any biological effect manifested by aquatic organisms due to a chemical is directly related to that dose at the site of action, and is the result of its bioaccumulation [97,99], which is modeled using biokinetic or toxicokinetic models.
Bio-toxicokinetic models apply the science and techniques of chemical kinetics to chemical behavior in biological systems [38]. The application of these models allows for the mathematical characterization of the processes of chemical absorption, distribution, and elimination. Bio-toxicokinetic models incorporate bioconcentration through expressions based on the material equilibrium between influx and efflux rates through differential equations that describe their evolution over time [100]. They are compartmental models that describe the movement of toxic substances between compartments. A compartment represents the quantity of a compound that behaves as if it exists in a homogeneously well-mixed container, and moves across the compartment boundary with a single absorption- or elimination-rate coefficient [101]. Biokinetic models consider the organism as a whole (a single compartment). However, toxicokinetic models consider different compartments within the organism, providing information on the bioconcentration in each compartment and a prediction of the most affected organs and/or tissues [102,103,104]. Both biokinetic and toxicokinetic models are modeled using mass balances, generally considering that the mass transfer between compartments follows first-order kinetics [33,38,105].
Figure 4 schematically presented the processes that can be considered in the approach of these models. All the possible flows of the biokinetic models are presented in Figure 4a. The uptake of the substance by the organism can be performed through the gills, when the substance is dissolved or in a colloidal form, or through the digestive system if it is associated with particles or food. Once in the organism, it can bioconcentrate and finally be excreted into the aquatic environment.
Figure 4b illustrates the possible flows of matter considering a toxicokinetic model. The difference present in the biokinetic model is that it considers the flows of matter within the organism, the bioconcentration of the substance in the different compartments, and the interrelation between them.
Considering first-order kinetics, the bioconcentration of a substance in an aquatic organism can be modeled as biokinetic according to the following general expression:
d C t dt = k u · C W + k u , col · C Col + k u · C aq ,   int + AE · IR food   C food +   AE · IR sed   C sed k e · C t  
C t is the bioaccumulated concentration of the substance in the organism, C W the concentration of the substance in water, C Col is the concentration of the substance in a colloidal form, C aq ,   int is the concentration of the substance in the pore water of the sediment, C food is the concentration of the substance present in food, and C sed is the concentration of the substance in the sediment particles. k u , k u , col , and k u represent the uptake-rate coefficients of the substance from water, colloids, and interstitial water, respectively. k e is the efflux-rate coefficient. AE · IR is the product of the efficiency of food assimilation in the digestive system ( AE ) produced by the rate of the specific ingestion of matter ( IR ) . The subscript “food” refers to food and the subscript “sed” to sediment particles. Equation (1) describes all the possible flows; however, depending on the study, this equation can be simplified. If the uptake from the sediments and colloids is not considered, Equation (1) is simplified into Equations (2) and (3), with food and food and sedimemt as the uptake inputs respectively:
d C t dt = k u · C W + AE · IR food   C food k e · C t
d C t dt = k u · C W + AE · IR food   C food + AE · IR sed   C sed k e · C t
When the substance is dissolved in water, its uptake by food can be neglected and this term in Equation (2) can be considered as negligible, obtaining the following expression presented in Equation (4):
d C t dt = k u · C W k e · C t
The formulations mentioned above assume no organism growth, which is a reasonable assumption during short laboratory experiments, but may be violated when the models are applied to predict accumulation in the field. When the size of the organism increases, “growth dilution” occurs, as the new tissue mass dilutes the toxicant concentration [101]. The efflux-rate coefficient, k e , of a growing organism overestimates the actual elimination. The growth-rate coefficient (g) must be incorporated into the model to correct the value of the efflux-rate coefficient. Considering the growth, Equation (5) subsequently describes the bioaccumulated concentration of the substance in the organism:
d C t dt = k u · C W + AE · IR food   C food k e + g · C t
Table 3 presents the different studies compiled from the bibliography that proposed the biokinetic models. Bivalves were the main organisms used in these studies because they are considered as excellent candidates that biomonitor coastal pollution [106], as they are filter feeders and exposed to high volumes of water to feed and breathe. In this process, the gills and other exchange surfaces are exposed to metals and other contaminants present in the environment, accumulating them in their soft tissues and shells [105,107].
The simplest model (Equation (4)) was proposed, mainly to describe the bioaccumulation kinetics of short laboratory experiments, since it neglects the effect of the growth of the organisms and does not consider the effect of food or the presence of sediment. The coefficients of the uptake and output rates were determined in order to obtain the bioconcentration factors of substances as organic UV filters [78,88]. In addition, it was also used as a rough estimate to serve as a starting point for more complex models [108].
Lengthy experiments conducted at field and laboratory levels need to incorporate the effect of the organism’s growth (Equation (5)). The laboratory studies analyzed the effects of different variables, such as the change in salinity [109] or changes in organism allometry [112]. Field studies allowed the validation of laboratory results or predicted their bioconcentration and effects on organisms in the aquatic environment. Some studies incorporated sediment particles as an uptake pathway (Table 3) with the aim of analyzing the main uptake pathway in different organisms [114,115].
The bioconcentration of a substance in different compartments within the organism can be modeled using toxicokinetic models (Figure 4b). In these models, the uptake of the substance can be condcuted in two ways. Dissolved substances are mainly captured by the gills and solid substances (food or particles) of the digestive canal. Once in the organism, the substance can be bioconcentrated in different interrelated compartments. Considering the first-order kinetics, Equations (6) and (7) describe all the possible flows:
d C gill dt = k u · C W + k 4   C org + k u · C aq ,   int k 3 · C gill
d C org dt = k 3 · C braq + AE · IR food   C food + AE · IR sed   C sed ± dC i dt k e · C i , e
where C gill and C org are the concentrations taken up by the gills and digestive canal, respectively; C i is the concentration of the substance in each compartment, i, of the organism; and C i , e is the concentration of the substance from which the substance is excreted. Kinetic constants k 3 and k 4 represent the uptake-rate coefficients from the gills and other compartments, respectively.
The term dC i dt depends on the compartment scheme (number of compartments and the relationship between them). It represents the bioaccumulation of the substance in the compartment and depends on the input and output flows considered. Each flux is usually assimilated to a first-order kinetics. For the model proposal, as many additional equations as the number of compartments must be indicated. The resolution of the system of differential equations is solved through matrix developments and, generally, through numerical methods implemented in specific software.
Equations 6 and 7 can be simplified if the uptake from sediments and feed is not considered, obtaining the following expressions presented in Equations (8) and (9):
d C gill dt = k u · C W + k 4   C org k 3 · C gill
d C org dt = k 3 · C braq ± dC i dt k e · C i , e
Sánchez-Marín et al. [105] proposed two different three- or four-compartment models to clarify the entry pathway (gills or digestive canal) and elimination of copper complexes with different contents of humic acid in mussels. Modeling demonstrated that the colloidal metal was available for uptake through the gut of marine mussels. Redeker et al. [116] proposed a five-compartment model (water, sediment, gut, central, and peripheral) to model the bioaccumulation of cadmium and zinc from water and sediment by the aquatic oligochaete, Tubifex tubifex. Simulations of different conditions showed that both dissolved and sediment-associated metal can be important sources of metal exposure for worms, and that the relative importance was highly dependent on the metal and exposure conditions, including the lability of metals in the sediment phase.
Chen and Liao [108] proposed a biokinetic model to understand the effect of dissolved copper on tilapia (Oreochromis mossambicus), a freshwater fish, when subjected to pulsed variations of the metal in the aquatic environment. It proposed a 10-compartment toxicokinetic model to predict the damage of the different organism’ organs.
Increasing numbers of studies have been designed to examine the effects of UV filters, mainly TiO2 and ZnO nanoparticles, on different aquatic organisms, ranging from bacteria to fish. These studies identified and evaluated the toxicity and ecotoxicological effects, through several biological responses, and bioaccumulation obtained mainly in laboratory-based experiments under different environmental conditions, e.g., [3,28,30,31,85,117]. However, as it has been highlighted in this section, a reduced number of studies proposed mathematical models to describe the interaction of sunscreen with living organisms, which is necessary to assess and predict impacts more accurately.

4. Conclusions

Sunscreen has become a product widely used by people all over the world, both in everyday life and associated with sun and beach tourism. Since its invention, it has evolved through different stages to achieve optimal protection against the entire UV spectrum with increasingly complex formulations. However, restrictions of the use of certain components, due to the impact produced by their release and transformation in an aquatic environment, is a driver for the design of environmentally friendly formulations.
The analysis and modeling of the release and behavior of sunscreen ingredients into the marine environment, including aquatic organisms, allowed us to predict the fate and levels of the chemicals under study, and therefore allowed us to predict the potential negative impacts on the aquatic environment. A large group of models studied the behavior of the particles and chemicals released into the water by hydrodynamic and kinetic models in the laboratory, while the knowledge of the fate of sunscreen components in the real environment is still limited. The laboratory-scale works analyzed and model: (i) the physical–chemical property behavior of sunscreens once they are incorporated into the aquatic environment, (ii) the photodegradation of sunscreen components, and (iii) the release kinetics of metals present in sunscreens and the formation of nutrients in the receiving aquatic environment. The direct photolysis of sunscreens is modeled as pseudo-first-order kinetics, while the indirect pathway by the reaction of sunscreen with reactive oxygen species is described as second-order kinetics.
Another significant group of works analyzed the effect of the sunscreen components on marine biota (toxicokinetic models), which caused the bioaccumulation of certain compounds in the organisms’ tissues. These works considered that the processes of uptake and excretion of chemicals to and from organisms, as well as their transfer between different organs, followed a first-order kinetic. The proposed models differ in the following ways: (i) the chemical uptake and (ii) the number of internal compartments of the organism considered, ranging from a single compartment to several compartments, such as the gills, digestive system, and internal organs. The kinetic parameters were obtained from the models for each species and the chemical allowed us to monitor the pollution levels in coastal areas, as well as the quality of the ecosystem.
The knowledge of the mechanisms of release and behaviour of the ingredients in sunscreens is essential to understand and predict their potential detrimental effects on marine ecosystems, as well as to identify their most environmentally harmful ingredients. The modeling approach was a useful prediction tool for these substances and can be of great value to environmental managers and stakeholders. However, the complexity (multitude of ingredients) and diversity (according to brand or format) of the sunscreen matrix, and the multitude of compartments and environmental variables that comprise the marine ecosystem, limited the application of these models. In this sense, future investigations should focus on the development of models, which consider the aforementioned factors, in order to obtain predictions of the potential impacts of these emerging products, and to inform future effective coastal-protection strategies and conservation policies based on sustainable tourism.

Author Contributions

All of the authors contributed to the study’s conceptualization and successive reviews. G.R.-G. and J.R.V.F.: writing—original draft preparation and supervision. A.R.-R. and A.T.-S.: writing—review and editing. A.R.-R.: funding acquisition and project administration. All authors have read and agreed to the published version of the manuscript.

Funding

This work is funded by the 2014–2020 ERDF Operational Programme and by the Department of Economy, Knowledge, Business and University of the Regional Government of Andalusia (project reference: FEDER-UCA18-106672).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

A. Rodríguez-Romero was supported by the Spanish grant Juan de la Cierva Incorporación referenced as IJC2018-037545-I.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sabanoglu, T. Cosmetics Industry—Statistics & Facts. Available online: https://www.statista.com/topics/3137/cosmetics-industry/#dossierKeyfigures (accessed on 12 January 2022).
  2. Market Report. Global Sunscreen Ingredients Market Report 2021, Featuring Key Players Ashland Inc., BASF, Evonik Industries, l’Oreal SA and Royal DSM. 30 September 2021 08:43 ET | Source: Research and Markets. 2021. Available online: https://www.researchandmarkets.com/categories/personal-care#hmc (accessed on 1 November 2021).
  3. Tovar-Sánchez, A.; Sánchez-Quiles, D.; Basterretxea, G.; Benedé, J.L.; Chisvert, A.; Salvador, A.; Moreno-Garrido, I.; Blasco, J. Sunscreen products as emerging pollutants to coastal waters. PLoS ONE 2013, 8, e65451. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Tsui, M.M.P.; Leung, H.W.; Wai, T.-C.; Yamashita, N.; Taniyasu, S.; Liu, W.; Lam, P.K.S.; Murphy, M.B. Occurrence, distribution and ecological risk assessment of multiple classes of UV filters in surface waters from different countries. Water Res. 2014, 67, 55–65. [Google Scholar] [CrossRef]
  5. Rodríguez-Romero, A.; Ruiz-Gutiérrez, G.; Viguri, J.R.; Tovar-Sánchez, A. Sunscreens as a new source of metals and nutrients to coastal waters. Environ. Sci. Technol. 2019, 53, 10177–10187. [Google Scholar] [CrossRef] [PubMed]
  6. Yuan, S.; Huang, J.; Jiang, X.; Huang, Y.; Zhu, X.; Cai, Z. Environmental fate and toxicity of sunscreen-derived inorganic ultraviolet filters in aquatic environments: A Review. Nanomaterials 2022, 12, 699. [Google Scholar] [CrossRef] [PubMed]
  7. Labille, J.; Catalano, R.; Slomberg, D.; Motellier, S.; Pinsino, A.; Hennebert, P.; Santaella, C.; Bartolomei, V. Assessing sunscreen lifecycle to minimize environmental risk posed by nanoparticulate UV-filters–A Review for safer-by-design products. Front. Environ. Sci. 2020, 8, 1–25. [Google Scholar] [CrossRef]
  8. Sánchez-Quiles, D.; Tovar-Sánchez, A. Are sunscreens a new environmental risk associated with coastal tourism? Environ. Int. 2015, 83, 158–170. [Google Scholar] [CrossRef] [Green Version]
  9. FDA. U.S Food & Drug Administration. Protector Solar. 2019. Available online: https://www.fda.gov/drugs/understanding-over-counter-medicines/sunscreen-how-help-protect-your-skin-sun (accessed on 27 February 2021).
  10. Coronado, M. Ingredient Trends and Innovation in Sun Protection. In-Cosmetic Global. 2017. Available online: https://www.researchgate.net/publication/317545165_Ingredients_Trends_and_Innovation_in_Sun_Protection (accessed on 5 April 2017).
  11. Casas-Beltrán, D.A.; Febles-Moreno, K.; Hernandez-Yac, E.; Gallaher, C.M.; Alvarado-Flores, J.; Leal-Bautista, R.M.; Lenczewski, M. Impact of tourist behavior on the discharge of sunscreen contamination in aquatic parks, sinkholes, and beaches of the Mexican caribbean. Appl. Sci. 2021, 11, 6882. [Google Scholar] [CrossRef]
  12. Watkins, Y.S.D.; Sallach, J.B. Investigating the exposure and impact of chemical UV filters on coral reef ecosystems: Review and research gap prioritization. Integr. Environ. Assess. Manag. 2021, 17, 967–981. [Google Scholar] [CrossRef]
  13. McCall, M.J.; Gulson, B.; Andrews, D. Consumer use of sunscreens containing nanoparticles. In Nanotechnology Environmental Health and Safety: Risks, Regulation, and Management; Hull, M., Bowman, D., Eds.; Elsevier: Cambridge, MA, USA, 2018; pp. 389–423. [Google Scholar]
  14. Osmond-McLeod, M.J.; Oytam, Y.; Rowe, A.; Sobhanmanesh, F.; Greenoak, G.; Kirby, J.; McInnes, E.F.; McCall, M.J. Long-term exposure to commercially available sunscreens containing nano-particles of TiO2 and ZnO revealed no biological impact in a hairless mouse model. Part. Fibre Toxicol. 2015, 13, 44. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. De la Calle, I.; Menta, M.; Klein, M.; Séby, F. Screening of TiO2 and Au nanoparticles in cosmetics and determination of elemental impurities by multiple techniques (DLS, SP-ICP-MS, ICP- MS and ICP-OES). Talanta 2017, 171, 291–306. [Google Scholar] [CrossRef] [PubMed]
  16. Slomberg, D.L.; Catalano, R.; Bartolomei, V.; Labille, J. Release and fate of nanoparticulate TiO2 UV filters from sunscreen: Effects of particle coating and formulation type. Environ. Pollut. 2021, 271, 116263. [Google Scholar] [CrossRef] [PubMed]
  17. Donia, D.T.; Carbone, M. Fate of the nanoparticles in environmental cycles. Int. J. Environ. Sci. Technol. 2019, 16, 583–600. [Google Scholar] [CrossRef]
  18. Tsui, M.M.P.; Lam, J.C.W.; Ng, T.Y.; Ang, P.O.; Murphy, M.B.; Lam, P.K.S. Occurrence, distribution, and fate of organic UV filters in coral communities. Environ. Sci. Technol. 2017, 51, 4182–4190. [Google Scholar] [CrossRef] [PubMed]
  19. Johnson, A.C.; Bowes, M.J.; Crossley, A.; Jarvie, H.P.; Jurkschat, K.; Jürgens, M.D.; Lawlor, A.J.; Park, B.; Rowland, P.; Spurgeon, D.; et al. An assessment of the fate, behaviour and environmental risk associated with sunscreen TiO2 nanoparticles in UK field scenarios. Sci. Total Environ. 2011, 409, 2503–2510. [Google Scholar] [CrossRef]
  20. Labille, J.; Feng, J.; Botta, C.; Borschneck, D.; Sammut, M.; Cabie, M.; Auffan, M.; Rose, J.; Bottero, J.Y. Aging of TiO2 nanocomposites used in sunscreen. Dispersion and fate of the degradation products in aqueous environment. Environ. Pollut. 2010, 158, 3482–3489. [Google Scholar] [CrossRef]
  21. Wang, N.; He, L.; Sun, X.; Li, X.; Li, M. The transformation of Benzophenone-3 in natural waters and AOPs: The roles of reactive oxygen species and potential environmental risks of products. J. Hazard. Mater. 2022, 427, 127941. [Google Scholar] [CrossRef]
  22. Lead, J.R.; Batley, G.E.; Alvarez, P.J.J.; Croteau, M.-N.; Handy, R.D.; McLaughlin, M.J.; Judy, J.D.; Schirmerh, K. Nanomaterials in the environment: Behavior, fate, bioavailability, and effects—An updated review. Environ. Toxicol. Chem. 2018, 37, 2029–2063. [Google Scholar] [CrossRef] [PubMed]
  23. Amde, M.; Liu, J.-F.; Tan, Z.-Q.; Bekana, D. Transformation and bioavailability of metal oxide nanoparticles in aquatic and terrestrial environments. A review. Environ. Pollut. 2017, 230, 250–267. [Google Scholar] [CrossRef] [PubMed]
  24. Li, Y.; Qiao, X.; Zhou, C.; Zhang, Y.; Fu, Z.; Chen, J. Photochemical transformation of sunscreen agent benzophenone-3 and its metabolite in surface freshwater and seawater. Chemosphere 2016, 153, 494–499. [Google Scholar] [CrossRef]
  25. Vione, D.; Caringella, R.; De Laurentiis, E.; Pazzi, M.; Minero, C. Phototransformation of the sunlight filter benzophenone-3 (2-hydroxy-4-methoxybenzophenone) under conditions relevant to surface waters. Sci. Total Environ. 2013, 463–464, 243–251. [Google Scholar] [CrossRef] [Green Version]
  26. Chen, B.C.; Chen, W.Y.; Ju, Y.R.; Tsai, J.W.; Jou, L.J.; Singh, S.; Liao, C.M. Combining bioaccumulation and coping mechanism to enhance long-term site-specific risk assessment for zinc susceptibility of bivalves. Chemosphere 2011, 84, 707–715. [Google Scholar] [CrossRef] [PubMed]
  27. da Silva, A.C.P.; Santos, B.A.M.C.; Castro, H.C.; Rodrigues, C.R. Ethylhexyl methoxycinnamate and butyl methoxydibenzoylmethane: Toxicological effects on marine biota and human concerns. J. Appl. Toxicol. 2022, 42, 73–86. [Google Scholar] [CrossRef] [PubMed]
  28. Cunningham, B.; Torres-Duarte, C.; Cherr, G.; Adams, N. Effects of three zinc-containing sunscreens on development of purple sea urchin (Strongylocentrotus purpuratus) embryos. Aquat. Toxicol. 2020, 218, 105355. [Google Scholar] [CrossRef] [PubMed]
  29. He, T.; Tsui, M.M.P.; Tan, C.J.; Ma, C.Y.; Yiu, S.K.F.; Wang, L.H.; Chen, T.H.; Fan, T.Y.; Lam, P.K.S.; Murphy, M.B. Toxicological effects of two organic ultraviolet filters and a related commercial sunscreen product in adult corals. Environ. Pollut. 2019, 245, 462–471. [Google Scholar] [CrossRef] [PubMed]
  30. Naasz, S.; Altenburger, R.; Kühnel, D. Environmental mixtures of nanomaterials and chemicals: The Trojan-horse phenomenon and its relevance for ecotoxicity. Sci. Total Environ. 2018, 635, 1170–1181. [Google Scholar] [CrossRef] [PubMed]
  31. Sendra, M.; Sánchez-Quiles, D.; Blasco, J.; Moreno-Garrido, I.; Lubián, L.M.; Pérez-García, S.; Tovar-Sánchez, A. Effects of TiO2 nanoparticles and sunscreens on coastal marine microalgae: Ultraviolet radiation is key variable for toxicity assessment. Environ. Int. 2017, 98, 62–68. [Google Scholar] [CrossRef]
  32. Kotnik, K.; Kosjek, T.; Žegura, B.; Filipič, M.; Heath, E. Photolytic fate and genotoxicity of benzophenone-derived compounds and their photodegradation mixtures in the aqueous environment. Chemosphere 2016, 147, 114–123. [Google Scholar] [CrossRef]
  33. Grech, A.; Brochot, C.; Dorne, J.-L.; Quignot, N.; Bois, F.Y.; Beaudouin, R. Toxicokinetic models and related tools in environmental risk assessment of chemicals. Sci. Total Environ. 2017, 578, 1–15. [Google Scholar] [CrossRef] [PubMed]
  34. Gottschalk, F.; Sun, T.; Nowack, B. Environmental concentrations of engineered nanomaterials: Review of modeling and analytical studies. Environ. Pollut. 2013, 181, 287–300. [Google Scholar] [CrossRef]
  35. Gottschalk, F.; Scholz, R.W.; Nowack, B. Probabilistic material flow modeling for assessing the environmental exposure to compounds: Methodology and an application to engineered nano-TiO2 particles. Environ. Modell. Softw. 2009, 25, 320–332. [Google Scholar] [CrossRef]
  36. Gottschalk, F.; Sonderer, T.; Scholz, R.W.; Nowack, B. Modeled environmental concentrations of engineered nanomaterials (TiO2, ZnO, Ag, CNT, Fullerenes) for different regions. Environ. Sci. Technol. 2009, 43, 9216–9222. [Google Scholar] [CrossRef] [PubMed]
  37. Mueller, N.C.; Nowack, B. Exposure modeling of engineered nanoparticles in the environment. Environ. Sci. Technol. 2008, 42, 4447–4453. [Google Scholar] [CrossRef] [PubMed]
  38. Barron, M.G.; Stehly, G.R.; Hayton, W.L. Pharmacokinetic modeling in aquatic animals l. Models and concepts. Aquat. Toxicol. 1990, 18, 61–86. [Google Scholar] [CrossRef]
  39. Zheng, Y.; Nowack, B. Size-specific, dynamic, probabilistic material flow analysis of titanium dioxide releases into the environment. Environ. Sci. Technol. 2021, 55, 2392–2402. [Google Scholar] [CrossRef] [PubMed]
  40. Lindo-Atichati, D.; Montero, P.; Rodil, R.; Quintana, J.B.; Miró, M. Modeling dispersal of UV filters in estuaries. Environ. Sci. Technol. 2019, 53, 1353–1363. [Google Scholar] [CrossRef] [PubMed]
  41. Arvidsson, R.; Molander, S.; Sandén, B.A. Particle flow analysis: Exploring potential use phase emissions of titanium dioxide nanoparticles from sunscreen, paint, and cement. J. Ind. Ecol. 2012, 16, 343–351. [Google Scholar] [CrossRef]
  42. Musee, N. Simulated environmental risk estimation of engineered nanomaterials: A case of cosmetics in Johannesburg city. Hum. Exp. Toxicol. 2010, 30, 1181–1195. [Google Scholar] [CrossRef]
  43. Gottschalk, F.; Ort, C.; Scholz, R.W.; Nowack, B. Engineered nanomaterials in rivers—Exposure scenarios for Switzerland at high spatial and temporal resolution. Environ. Pollut. 2011, 159, 3439–3445. [Google Scholar] [CrossRef] [PubMed]
  44. De Laurentiis, E.; Minella, M.; Sarakha, M.; Marrese, A.; Minero, C.; Mailhot, G.; Brigante, M.; Vione, D. Photochemical processes involving the UV absorber benzophenone-4 (2-hydroxy-4-methoxybenzophenone-5-sulphonic acid) in aqueous solution: Reaction pathways and implications for surface waters. Water Res. 2013, 47, 5943–5953. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Zhou, L.; Ji, Y.; Zeng, C.; Zhang, Y.; Wang, Z.; Yang, X. Aquatic photodegradation of sunscreen agent p-aminobenzoic acid in the presence of dissolved organic matter. Water Res. 2013, 47, 153–162. [Google Scholar] [CrossRef]
  46. Thio, B.J.R.; Zhou, D.X.; Keller, A.A. Influence of natural organic matter on the aggregation and deposition of titanium dioxide nanoparticles. J. Hazard. Mater. 2011, 189, 556–563. [Google Scholar] [CrossRef] [PubMed]
  47. Mao, L.; Meng, C.; Zeng, C.; Ji, Y.; Yang, X.; Gao, S. The effect of nitrate, bicarbonate and natural organic matter on the degradation of sunscreen agent p-aminobenzoic acid by simulated solar irradiation. Sci. Total Environ. 2011, 409, 5376–5381. [Google Scholar] [CrossRef] [PubMed]
  48. Sakkas, V.A.; Giokas, D.L.; Lambropoulou, D.A.; Albanis, T.A. Aqueous photolysis of the sunscreen agent octyl-dimethyl-p-aminobenzoic acid: Formation of disinfection byproducts in chlorinated swimming pool water. J. Chromatogr. A. 2003, 1016, 211–222. [Google Scholar] [CrossRef]
  49. Celeiro, M.; Vignola Hackbarth, F.; Selene, M.A.; de Souza, G.U.; Llompart, M.; Vilar, V.J.P. Assessment of advanced oxidation processes for the degradation of three UV filters from swimming pool water. J. Photochem. Photobiol. A Chem. 2018, 351, 95–107. [Google Scholar] [CrossRef]
  50. Ji, Y.; Zhou, L.; Zhang, Y.; Ferronato, C.; Brigante, M.; Mailhot, G.; Yang, X.; Chovelon, J.M. Photochemical degradation of sunscreen agent 2-phenylbenzimidazole-5-sulfonic acid in different water matrices. Water Res. 2013, 47, 5865–5875. [Google Scholar] [CrossRef] [PubMed]
  51. Botta, C.; Labille, J.; Auffan, M.; Borschneck, D.; Miche, H.; Cabié, M.; Masion, A.; Rose, J.; Bottero, J.Y. TiO2-based nanoparticles released in water from commercialized sunscreens in a life-cycle perspective: Structures and quantities. Environ. Pollut. 2011, 159, 1543–1550. [Google Scholar] [CrossRef]
  52. Ji, Y.; Zhou, L.; Ferronato, C.; Salvador, A.; Yang, X.; Chovelon, J.M. Degradation of sunscreen agent 2-phenylbenzimidazole-5-sulfonic acid by TiO2 photocatalysis: Kinetics, photoproducts and comparison to structurally related compounds. Appl. Catal. B Environ. 2013, 140–141, 457–467. [Google Scholar] [CrossRef]
  53. Lu, G.Y.; Ke, C.H.; Zhu, A.; Wang, W.X. Oyster-based national mapping of trace metals pollution in the Chinese coastal waters. Environ. Pollut. 2017, 224, 658–669. [Google Scholar] [CrossRef]
  54. Sánchez-Quiles, D.; Tovar-Sánchez, A. Sunscreens as a source of hydrogen peroxide production in coastal waters. Environ. Sci. Technol. 2014, 48, 9037–9042. [Google Scholar] [CrossRef] [Green Version]
  55. Yin, K.; Li, T.; Zhang, T.; Zhang, Y.; Yang, C.; Luo, S. Degradation of organic filter 2-Phenylbenzidazole-5-Sulfonic acid by light-driven free chlorine process: Reactive species and mechanisms. Chem. Eng. J. 2022, 430, 132684. [Google Scholar] [CrossRef]
  56. Wang, Z.; Khan, M.A.; Al-Othman, A.A.S.; Abdullah, Z.; Sillanpää, M. Pilot-scale study on photodegradation of benzophenone-3 and benzophenone-8 ultraviolet filters enriched synthetic effluent. J. Water Process. Eng. 2021, 44, 102327. [Google Scholar] [CrossRef]
  57. Lee, Y.M.; Lee, G.; Zoh, K.D. Benzophenone-3 degradation via UV/H2O2 and UV/persulfate reactions. J. Hazard. Mater. 2021, 403, 123591. [Google Scholar] [CrossRef] [PubMed]
  58. Yin, J.-Y.; Wang, H.; Yu, K.-P.; Lee, J.; Lin, K.-Y.A. Degradation of sunscreen agent 2-phenylbenzimidazole-5-sulfonic acid using monopersulfate activated by MOF-derived cobalt sulfide nanoplates. J. Water Process. Eng. 2021, 44, 102282. [Google Scholar] [CrossRef]
  59. Moradi, N.; Amin, M.M.; Fatehizadeh, A.; Ghasemi, Z. Degradation of UV-filter Benzophenon-3 in aqueous solution using TiO2 coated on quartz tubes. J. Environ. Health Sci. Eng. 2018, 16, 213–228. [Google Scholar] [CrossRef] [PubMed]
  60. Xue, Y.; Dong, W.; Wang, X.; Bi, W.; Zhai, P.; Li, H.; Nie, M. Degradation of sunscreen agent p-aminobenzoic acid using a combination system of UV irradiation, persulphate and iron(II). Environ. Sci. Pollut. Res. 2016, 23, 4561–4568. [Google Scholar] [CrossRef]
  61. Abdelraheem, W.H.M.; He, X.; Duan, X.; Dionysiou, D.D. Degradation and mineralization of organic UV absorber compound 2-phenylbenzimidazole-5-sulfonic acid (PBSA) using UV-254nm/H2O2. J. Hazard. Mater. 2015, 282, 233–240. [Google Scholar] [CrossRef] [PubMed]
  62. Gong, P.; Yuan, H.; Zhai, P.; Xue, Y.; Li, H.; Dong, W.; Mailhot, G. Investigation on the degradation of benzophenone-3 by UV/H2O2 in aqueous solution. Chem. Eng. J. 2015, 277, 97–103. [Google Scholar] [CrossRef]
  63. Auffan, M.; Pedeutour, M.; Rose, J.; Masion, A.; Ziarelli, F.; Borschneck, D.; Chaneac, C.; Botta, C.; Chaurand, P.; Labille, J.; et al. Structural degradation at the surface of a TiO2-based nanomaterial used in cosmetics. Environ. Sci. Technol. 2010, 44, 2689–2694. [Google Scholar] [CrossRef] [PubMed]
  64. Zachariadis, G.A.; Sahanidou, E. Multi-element method for determination of trace elements in sunscreens by ICP-AES. J. Pharm. Biomed. Anal. 2009, 50, 342–348. [Google Scholar] [CrossRef] [PubMed]
  65. Zmozinski, A.V.; Pretto, T.; Borges, A.R.; Duarte, Á.T.; Vale, M.G.R. Determination of Pb and Cr in sunscreen samples by high-resolution continuum source graphite furnace atomic absorption spectrometry and direct analysis. Microchem. J. 2016, 128, 89–94. [Google Scholar] [CrossRef]
  66. Brun, N.R.; Lenz, M.; Wehrli, B.; Fent, K. Comparative effects of zinc oxide nanoparticles and dissolved zinc on zebrafish embryos and eleuthero-embryos: Importance of zinc ions. Sci. Total Environ. 2014, 476–477, 657–666. [Google Scholar] [CrossRef] [PubMed]
  67. Peng, C.; Zhang, W.; Gao, H.; Li, Y.; Tong, X.; Li, K.; Zhu, X.; Wang, Y.; Chen, Y. Behavior and potential impacts of metal-based engineered nanoparticles in aquatic environments. Nanomaterials 2017, 7, 21. [Google Scholar] [CrossRef]
  68. Liu, X.; Chen, G.; Erwin, J.G.; Adam, N.K.; Su, C. Release of phosphorous impurity from TiO2 anatase and rutile nanoparticles in aquatic environments and its implications. Water Res. 2013, 47, 6149–6156. [Google Scholar] [CrossRef] [PubMed]
  69. Low, G.K.C.; McEvoy, S.R.; Matthews, R.W. Formation of nitrate and ammonium ions in titanium dioxide mediated photocatalytic degradation of organic compounds containing nitrogen atoms. Environ. Sci. Technol. 1991, 25, 460–467. [Google Scholar] [CrossRef]
  70. Jing, J.; Liu, M.; Colvin, V.L.; Li, W.; Yu, W.W. Photocatalytic degradation of nitrogen-containing organic compounds over TiO2. J. Mol. Catal. A Chem. 2011, 351, 17–28. [Google Scholar] [CrossRef]
  71. Lv, J.; Zhang, S.; Luo, L.; Han, W.; Zhang, J.; Yang, K.; Christie, P. Dissolution and microstructural transformation of ZnO nanoparticles under the influence of phosphate. Environ. Sci. Technol. 2012, 46, 7215–7221. [Google Scholar] [CrossRef] [PubMed]
  72. Liu, X.; Chen, G.; Erwin, J.G.; Su, C. Silicon impurity release and surface transformation of TiO2 anatase and rutile nanoparticles in water environments. Environ. Pollut. 2014, 184, 570–578. [Google Scholar] [CrossRef] [PubMed]
  73. Casado-Martínez, M.C.; Blasco, J.; DelValls, T.A.; González-Castromil, M.A.; Riba, I. Interlaboratory assessment of marine bioassays to evaluate the environmental quality of coastal sediments in Spain. V. Whole sediment toxicity test using juveniles of the bivalve Ruditapes philippinarum. Cienc. Mar. 2006, 32, 159–166. [Google Scholar] [CrossRef] [Green Version]
  74. Zhou, Q.; Zhang, J.; Fu, J.; Shi, J.; Jiang, G. Biomonitoring: An appealing tool for assessment of metal pollution in the aquatic ecosystem. Anal. Chim. Acta 2008, 606, 135–150. [Google Scholar] [CrossRef]
  75. Rodríguez-Romero, A.; Jiménez-Tenorio, N.; Basallote, M.D.; de Orte, M.R.; Blasco, J.; Riba, I. Predicting the impacts of CO2 leakage from subseabed storage: Effects of metal accumulation and toxicity on the model benthic organism Ruditapes philippinarum. Environ. Sci. Technol. 2014, 48, 12292–12301. [Google Scholar] [CrossRef] [PubMed]
  76. Gutiérrez, M.F.; Ale, A.; Andrade, V.; Bacchetta, C.; Rossi, A.; Cazenave, J. Metallic, metal oxide, and metalloid nanoparticles toxic effects on freshwater microcrustaceans: An update and basis for the use of new test species. Water Environ. Res. 2021, 93, 2505–2526. [Google Scholar] [CrossRef]
  77. Wong, S.W.Y.; Zhou, G.J.; Leung, P.T.Y.; Han, J.; Lee, J.S.; Kwok, K.W.H.; Leung, K.M.Y. Sunscreens containing zinc oxide nanoparticles can trigger oxidative stress and toxicity to the marine copepod Tigriopus japonicus. Mar. Pollut. Bull. 2020, 154, 111078. [Google Scholar] [CrossRef]
  78. Vidal-Liñán, L.; Villaverde-de-Sáa, E.; Rodil, R.; Quintana, J.B.; Beiras, R. Bioaccumulation of UV filters in Mytilus galloprovincialis mussel. Chemosphere 2018, 190, 267–271. [Google Scholar] [CrossRef] [PubMed]
  79. Corinaldesi, C.; Damiani, E.; Marcellini, F.; Falugi, C.; Tiano, L.; Brugè, F.; Danovaro, R. Sunscreen products impair the early developmental stages of the sea urchin Paracentrotus lividus. Sci. Rep. 2017, 7, 7815. [Google Scholar] [CrossRef] [Green Version]
  80. Miller, R.J.; Bennett, S.; Keller, A.A.; Pease, S.; Lenihan, H.S. TiO2 nanoparticles are phototoxic to marine phytoplankton. PLoS ONE 2012, 7, 1. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Moeller, M.; Pawlowski, S.; Petersen-Thiery, M.; Miller, I.B.; Nietzer, S.; Heisel-Sure, Y.; Kellermann, M.Y.; Schupp, P.J. Challenges in current coral reef protection—Possible impacts of UV filters used in sunscreens, a critical review. Front. Mar. Sci. 2021, 8, 665548. [Google Scholar] [CrossRef]
  82. Benedé, J.L.; Chisvert, A.; Salvador, A.; Sánchez-Quiles, D.; Tovar-Sánchez, A. Determination of UV filters in both soluble and particulate fractions of seawaters by dispersive liquid–liquid microextraction followed by gas chromatography–mass spectrometry. Anal. Chim. Acta 2014, 812, 50–58. [Google Scholar] [CrossRef] [PubMed]
  83. Gago-Ferrero, P.; Badia-Fabregat, M.; Olivares, A.; Piña, B.; Blánquez, P.; Vicent, T.; Caminal, G.; Díaz-Cruz, M.S.; Barceló, D. Evaluation of fungal- and photo-degradation as potential treatments for the removal of sunscreens BP3 and BP1. Sci. Total Environ. 2012, 427–428, 355–363. [Google Scholar] [CrossRef] [PubMed]
  84. Huang, Y.; Law, J.C.F.; Lam, T.K.; Leung, K.S.Y. Risks of organic UV filters: A review of environmental and human health concern studies. Sci. Total Environ. 2021, 755, 142486. [Google Scholar] [CrossRef] [PubMed]
  85. Schneider, S.L.; Lim, H.W. A review of inorganic UV filters zinc oxide and titanium dioxide. Photodermatol. Photoimmunol. Photomed. 2019, 35, 442–446. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Danovaro, R.; Bongiorni, L.; Corinaldesi, C.; Giovannelli, D.; Damiani, E.; Astolfi, P.; Greci, L. Sunscreens cause coral bleaching by promoting viral infections. Environ. Health. Perspect. 2008, 116, 441–447. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Araújo, C.V.M.; Rodríguez-Romero, A.; Fernández, M.; Sparaventi, E.; Medina, M.M.; Tovar-Sánchez, A. Repellency and mortality effects of sunscreens on the shrimp Palaemon varians: Toxicity dependent on exposure method. Chemosphere 2020, 257, 127190. [Google Scholar] [CrossRef] [PubMed]
  88. Gomez, E.; Bachelot, M.; Boillot, C.; Munaron, D.; Chiron, S.; Casellas, C.; Fenet, H. Bioconcentration of two pharmaceuticals (benzodiazepines) and two personal care products (UV filters) in marine mussels (Mytilus galloprovincialis) under controlled laboratory conditions. Environ. Sci. Pollut. Res. 2012, 19, 2561–2569. [Google Scholar] [CrossRef]
  89. Fang, J.; Xu, M.J.; Wang, D.J.; Wen, B.; Han, J.Y. Modeling the transport of TiO2 nanoparticle aggregates in saturated and unsaturated granular media: Effects of ionic strength and pH. Water Res. 2013, 47, 1399–1408. [Google Scholar] [CrossRef] [PubMed]
  90. Keller, A.A.; Wang, H.; Zhou, D.; Lenihan, H.S.; Cherr, G.; Cardinale, B.J.; Miller, R.; Ji, Z. Stability and aggregation of metal oxide nanoparticles in natural aqueous matrices. Environ. Science Technol. 2010, 44, 1962–1967. [Google Scholar] [CrossRef] [PubMed]
  91. Zhang, Y.; Chen, Y.; Westerhoff, P.; Crittenden, J. Impact of natural organic matter and divalent cations on the stability of aqueous nanoparticles. Water Res. 2009, 43, 4249–4257. [Google Scholar] [CrossRef] [PubMed]
  92. de Miranda, L.L.R.; Harvey, K.E.; Ahmed, A.; Harvey, S.C. UV-filter pollution: Current concerns and future prospects. Environ. Monit. Assess. 2021, 193, 840–856. [Google Scholar] [CrossRef] [PubMed]
  93. Ginzburg, A.L.; Blackburn, R.S.; Santillan, C.; Truong, L.; Tanguay, R.L.; Hutchison, J.E. Zinc oxide-induced changes to sunscreen ingredient efficacy and toxicity under UV irradiation. Photochem. Photobiol. Sci. 2021, 20, 1273–1285. [Google Scholar] [CrossRef]
  94. Vieira Sanches, M.; Oliva, M.; De Marchi, L.; Cuccaro, A.; Puppi, D.; Chiellini, F.; Freitas, R.; Pretti, C. Ecotoxicological screening of UV-filters using a battery of marine bioassays. Environ. Pollut. 2021, 290, 118011. [Google Scholar] [CrossRef] [PubMed]
  95. Franklin, N.; Rogers, N.J.; Apte, S.; Batley, G.; Gadd, G.; Casey, P. Comparative toxicity of tanoparticulate ZnO, Bulk ZnO, and ZnCl2 to a freshwater microalga (Pseudokirchneriella subcapitata): The importance of particle solubility. Environ. Sci. Technol. 2007, 41, 8484–8490. [Google Scholar] [CrossRef]
  96. Fan, W.; Liu, L.; Peng, R.; Wang, W.X. High bioconcentration of titanium dioxide nanoparticles in Daphnia magna determined by kinetic approach. Sci. Total Environ. 2016, 569–570, 1224–1231. [Google Scholar] [CrossRef]
  97. Beiras, R. Bioaccumulation. Marine Pollution Sources, Fate and Effects of Pollutants in Coastal Ecosystems; Elsevier: Amsterdam, The Netherlands, 2018; Chapter 11; pp. 187–204. [Google Scholar] [CrossRef]
  98. Schäfer, S.; Buchmeier, G.; Claus, E.; Duester, L.; Heininger, P.; Körner, A.; Mayer, P.; Paschke, A.; Rauert, C.; Reifferscheid, G.; et al. Bioaccumulation in aquatic systems: Methodological approaches, monitoring and assessment. Environ. Sci. Eur. 2015, 27, 5. [Google Scholar] [CrossRef] [PubMed]
  99. Wang, W.X.; Tan, Q.G. Applications of dynamic models in predicting the bioaccumulation, transport and toxicity of trace metals in aquatic organisms. Environ. Pollut. 2019, 252, 1561–1573. [Google Scholar] [CrossRef] [PubMed]
  100. Curis, E.; Nicolis, I.; Bensaci, J.; Deschamps, P.; Bénazeth, S. Mathematical modeling in metal metabolism: Overview and perspectives. Biochimie 2009, 91, 1238–1254. [Google Scholar] [CrossRef] [PubMed]
  101. Landrum, P.F.; Lee II, H.; Lydy, M.J. Toxicokinetics in aquatic systems: Model comparisons and use in hazard assessment. Environ. Toxicol. Chem. 1992, 11, 1709–1725. [Google Scholar] [CrossRef]
  102. Gestin, O.; Lacoue-Labarthe, T.; Coquery, M.; Delorme, N.; Garnero, L.; Dherret, L.; Ciccia, T.; Geffard, O.; Lopes, C. One and multi-compartments toxico-kinetic modeling to understand metals’ organotropism and fate in Gammarus fossarum. Environ. Int. 2021, 156, 106625. [Google Scholar] [CrossRef]
  103. Gao, Y.; Zhang, Y.; Feng, J.; Zhu, L. Toxicokinetic-toxicodynamic modeling of cadmium and lead toxicity to larvae and adult zebrafish. Environ. Pollut. 2019, 251, 221–229. [Google Scholar] [CrossRef]
  104. Chen, W.Y. Toxicokinetic modeling challenges for aquatic nanotoxicology. Front. Mar. Sci. 2016, 2, 114. [Google Scholar] [CrossRef] [Green Version]
  105. Sánchez-Marín, P.; Aierbe, E.; Lorenzo, J.I.; Mubiana, V.K.; Beiras, R.; Blust, R. Dynamic modeling of copper bioaccumulation by Mytilus edulis in the presence of humic acid aggregates. Aquat. Toxicol. 2016, 178, 165–170. [Google Scholar] [CrossRef] [PubMed]
  106. Lu, G.; Zhu, A.; Fang, H.; Dong, Y.; Wang, W.X. Establishing baseline trace metals in marine bivalves in China and worldwide: Meta-analysis and modeling approach. Sci. Total Environ. 2019, 669, 746–753. [Google Scholar] [CrossRef]
  107. Rainbow, P.S.; Phillips, D.J.H. Cosmopolitan biomonitors of trace metals. Mar. Pollut. Bull. 1993, 26, 593–601. [Google Scholar] [CrossRef]
  108. Chen, W.Y.; Liao, C.M. Interpreting copper bioaccumulation dynamics in tilapia using systems-level explorations of pulsed acute/chronic exposures. Ecotoxicology 2014, 23, 1124–1136. [Google Scholar] [CrossRef] [PubMed]
  109. Yin, Q.; Wang, W.X. Multiple trace element accumulation in the mussel Septifer virgatus: Counteracting effects of salinity on uptake and elimination. Environ. Pollut. 2018, 242, 375–382. [Google Scholar] [CrossRef] [PubMed]
  110. Kalman, J.; Smith, B.D.; Bury, N.R.; Rainbow, P.S. Biodynamic modelling of the bioaccumulation of trace metals (Ag, As and Zn) by an in faunal estuarine invertebrate, the clam Scrobicularia plana. Aquat. Toxicol. 2014, 154, 121–130. [Google Scholar] [CrossRef] [PubMed]
  111. Bourgeault, A.; Gourlay-Francéa, C.; Priadic, C.; Ayraultc, S.; Tusseau-Vuillemin, M.H. Bioavailability of particulate metal to zebra mussels: Biodynamic modelling shows that assimilation efficiencies are site-specific. Environ. Pollut. 2011, 159, 3381–3389. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  112. Pan, K.; Wang, W.X. Allometry of cadmium and zinc concentrations and bioaccumulation in the scallop Chlamys nobilis. Mar. Ecol. Prog. Ser. 2008, 365, 115–126. [Google Scholar] [CrossRef]
  113. Wang, W.X.; Fisher, N.S.; Luoma, S.N. Kinetic determinations of trace element bioaccumulation in the mussel Mytilus edulis. Mar. Ecol. Prog. Ser. 1996, 140, 91–113. [Google Scholar] [CrossRef]
  114. O’Mara, K.; Adams, M.; Burford, M.A.; Fry, B.; Cresswell, T. Uptake and accumulation of cadmium, manganese and zinc by fisheries species: Trophic differences in sensitivity to environmental metal accumulation. Sci. Total Environ. 2019, 690, 867–877. [Google Scholar] [CrossRef] [PubMed]
  115. Lee, J.H.; Birch, G.F.; Cresswell, T.; Johansen, M.P.; Adams, M.S.; Simpson, S.L. Dietary ingestion of fine sediments and microalgae represent thedominant route of exposure and metal accumulation for Sydney rockoyster (Saccostrea glomerata): A biokinetic model for zinc. Aquat. Toxicol. 2015, 167, 46–54. [Google Scholar] [CrossRef] [PubMed]
  116. Redeker, E.S.; Bervoets, L.; Blu, R. Dynamic Model for the Accumulation of Cadmium and Zinc from Water and Sediment by the Aquatic Oligochaete, Tubifex tubifex. Environ. Sci. Technol. 2004, 38, 6193–6200. [Google Scholar] [CrossRef]
  117. Haynes, V.N.; Ward, J.E.; Russell, B.J.; Agrios, A.G. Photocatalytic effects of titanium dioxide nanoparticles on aquatic organisms—Current knowledge and suggestions for future research. Aquat Toxicol. 2017, 185, 138–148. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Classification of publications based on the components of sunscreens, grouping them into organic and inorganic filters, metals and nutrients, and split between field and laboratory studies. Oceans 03 00024 i013 Field studies; Oceans 03 00024 i014 laboratory studies.
Figure 1. Classification of publications based on the components of sunscreens, grouping them into organic and inorganic filters, metals and nutrients, and split between field and laboratory studies. Oceans 03 00024 i013 Field studies; Oceans 03 00024 i014 laboratory studies.
Oceans 03 00024 g001
Figure 2. Superstructure of the sunscreen’s transformations in the environment. References in the figure refer to works that propose process modeling. The most representative references for the process in Figure 2 are: Metals Release and Adsoption/Desorption [5]; Organic UV filters Release [40]; Organic UV filters Direct photolysis [24,25,32,44,45,47,48,50]; Organic UV filters Photocatalysis [52]; Organic UV filters Indirect photolysis [24,44,45,47,50].
Figure 2. Superstructure of the sunscreen’s transformations in the environment. References in the figure refer to works that propose process modeling. The most representative references for the process in Figure 2 are: Metals Release and Adsoption/Desorption [5]; Organic UV filters Release [40]; Organic UV filters Direct photolysis [24,25,32,44,45,47,48,50]; Organic UV filters Photocatalysis [52]; Organic UV filters Indirect photolysis [24,44,45,47,50].
Oceans 03 00024 g002
Figure 3. Classification of publications based on the components of sunscreens (UV filters and metals), according to the aquatic organisms studied and splitting between field and laboratory studies. Oceans 03 00024 i013 Field studies; Oceans 03 00024 i014 laboratory studies.
Figure 3. Classification of publications based on the components of sunscreens (UV filters and metals), according to the aquatic organisms studied and splitting between field and laboratory studies. Oceans 03 00024 i013 Field studies; Oceans 03 00024 i014 laboratory studies.
Oceans 03 00024 g003
Figure 4. Processes of uptake, bioconcentration, and excretion of a substance by the aquatic organism. (a) Biokinetic model and (b) toxicokinetic model. Where C t ,   C W , C Col , C aq ,   int , C food ,     C sed , C gill , C org , and C i , e are the concentrations of bioaccumulated substances in water, in colloid form, in the pore water of the sediment, in the food, in the sediment particles, uptake by the gills, and by the digestive canal and excreted, respectively. k u , k u , col , k u ,     k 3 , and k 4 represent the uptake-rate coefficients of the substance from water, colloids, interstitial water, gills, and other compartments, respectively; k e is the efflux-rate coefficient; AE · IR is the product of the efficiency of food assimilation in the digestive system ( AE ) by the rate of specific ingestion of matter ( IR ) ; the subscript “food” refers to food and the subscript “sed” to sediment particles..
Figure 4. Processes of uptake, bioconcentration, and excretion of a substance by the aquatic organism. (a) Biokinetic model and (b) toxicokinetic model. Where C t ,   C W , C Col , C aq ,   int , C food ,     C sed , C gill , C org , and C i , e are the concentrations of bioaccumulated substances in water, in colloid form, in the pore water of the sediment, in the food, in the sediment particles, uptake by the gills, and by the digestive canal and excreted, respectively. k u , k u , col , k u ,     k 3 , and k 4 represent the uptake-rate coefficients of the substance from water, colloids, interstitial water, gills, and other compartments, respectively; k e is the efflux-rate coefficient; AE · IR is the product of the efficiency of food assimilation in the digestive system ( AE ) by the rate of specific ingestion of matter ( IR ) ; the subscript “food” refers to food and the subscript “sed” to sediment particles..
Oceans 03 00024 g004
Table 1. IUPAC chemical name, chemical structure, regulation, and main properties of the most common UV filters used in sunscreens. Red color represents banned chemicals.
Table 1. IUPAC chemical name, chemical structure, regulation, and main properties of the most common UV filters used in sunscreens. Red color represents banned chemicals.
UV FilterIUPAC Chemical Name (Formula)Chemical StructureRegulation(*) (**)Properties
Benzophenone -3
Oxybenzone
(BP3, BZ-3)
2-Hydroxy-4-methoxybenzophenone
(C14H12O3)
Oceans 03 00024 i001Regulated by maximum contents in the US, Canada, the EU, ASEAN, MERCOSUR, China, India, Japan, Korea, Australia, NZ, and South Africa.
Banned in Palau, Thailand marine natural parks, Aruba, Bonaire, the US Virgin Islands, Key West (Florida, US), and Hawaii (US).
CAS number: 131-57-7; molecular weight (MW): 228.24
Density: 1.201 g/cm3
Solubility in water: <0.1 g/100 mL at 20 °C
pKa = 7.56; log Kow = 3.52
Benzophenone-4;
Sulisobenzone
(BP4)
5-Benzoyl-4-hydroxy-2-methoxybenzene-1-sulfonic acid (C14H12O6S) Oceans 03 00024 i002Regulated by maximum contents in the US, Canada, the EU, ASEAN, MERCOSUR, China, India, Japan, Korea, Australia, NZ, and South Africa.CAS number: 4065-45-6; MW: 308.31
Solubility in water: 1 g per 4 mL
pKa1 = −0.70 (sulfonic acid); pKa2 = 7.56 (hydroxyl); log Kow = 0.37
Octinoxate; Uvinul MC80; Octyl methoxycinnamate (OMC)(RS)-2-Ethylhexyl (2E)-3-(4-methoxyphenyl) prop-2-enoate (C18H26O3) Oceans 03 00024 i003Regulated by maximum contents in the US, Canada, the EU, ASEAN, MERCOSUR, China, India, Japan, Korea, Australia, NZ, and South Africa.
Banned in Palau, Thailand marine natural parks, US Virgin Islands, Key West (Florida), and Hawaii (US).
CAS number: 5466-77-3; MW: 290.40
Density: 1.010 g/cm3
Insoluble in water
log Kow = 5.80
Octocrylene; Uvinul N-539 (OC)2-Ethylhexyl 2-cyano-3,3-diphenylprop-2-enoate (C24H27NO2) Oceans 03 00024 i004Regulated by maximum contents in the US, Canada, the EU, ASEAN, MERCOSUR, China, India, Japan, Korea, and South Africa.
Banned in Palau and US Virgin Islands.
CAS number: 6197-30-4; MW: 361.48
Density: 1.055 g/cm3
Insoluble in water
log Kow = 7.35
Enzacamene; 4-Methylbenzylidene camphor (4-MBC)(3E)-1,7,7-Trimethyl-3-[(4-methylphenyl) methylene]-2-norbornanone (C18H22O) Oceans 03 00024 i005Regulated by maximum contents in Canada, the EU, ASEAN, MERCOSUR, China, India, Japan, Korea, Australia, NZ, and South Africa.
Banned in Palau and Thailand marine natural parks.
CAS number: 36861-47-9; MW: 254.37
Density: 1.064 g/cm3
Insoluble in water
log Kow =5.47
Ácido p-aminobenzoico (PABA)Ácido 4-aminobenzoico (C7H7NO2) Oceans 03 00024 i006Regulated by maximum contents in the US, Canada, the EU, ASEAN, MERCOSUR, China, India, Japan, Korea, Australia, NZ, and South Africa.CAS number: 150-13-0; MW: 137.14
Density: 1.374 g/cm3
Solubility in water: 1 g/170 mL (25 °C)
pKa1 = 2.38; pKa2 = 4.85; log Kow = 0.83
Padimate O; Escalol 507; octyldimethyl PABA (OD-PABA)2-ethylhexyl 4-(dimethylamino) benzoate (C17H27NO2) Oceans 03 00024 i007Permitted UV filters and regulated by maximum contents in the US, Canada, the EU, ASEAN, MERCOSUR, China, India, Japan, Korea, Australia, NZ, and South Africa.CAS number: 21245-02-3; MW: 277.40
Density: 0.990 g/cm3
Solubility in water: 0.54 mg/L
pKa = 2.9; log Kow = 5.77
Ensulizole; phenyl
benzimidazole sulfonic acid (PBSA)
2-Phenyl-3H-benzimidazole-5-sulfonic acid (C13H10N2O3S) Oceans 03 00024 i008Regulated by maximum contents in the US, Canada, the EU, ASEAN, MERCOSUR, China, India, Japan, Korea, Australia, NZ, and South Africa.CAS Number: 27503-81-7
Molecular M.: 274.29
Soluble in water
2,4-Dihydroxy
Benzophenone;
Benzophenone-1 (DH-BP)
(2,4-dihydroxyphenyl)phenyl-methanone (C13H10O3) Oceans 03 00024 i009Regulated by maximum contents in Japan and South Africa.CAS number: 131-56-6; MW: 214.22
Density: 1.302 g/cm3
Insoluble in water; log Kow = 3.17
Dioxybenzone;
Benzophenone-8 (DHM-BP)
2,2′-Dihydroxy-4-methoxybenzophenone (C14H12O4) Oceans 03 00024 i010Regulated by maximum contents in the US, Canada, ASEAN, MERCOSUR, Korea, Australia, NZ, and South Africa.CAS number: 131-53-3; MW: 244.24
Density: 1.38 g/cm3
Insoluble in water; pKa = 7.11; log Kow = 4.31
TiO2 nanoparticleTitanium dioxide (TiO2) Oceans 03 00024 i011Regulated by maximum contents in the US, Canada, the EU, ASEAN, MERCOSUR, China, Japan, Korea, Australia, and NZ.CAS number: 13463-67-7
XRD: rutile, anatase
Al(OH)3 or SiO2 coating; TiO2 %: 79–89
ZnO nanoparticleZinc oxide (ZnO) Oceans 03 00024 i012Regulated by maximum contents in the US, Canada, ASEAN, MERCOSUR, China, Japan, Korea, and NZ.CAS number: 1314-13-2; XRD: wurtzite
SiO2-coated or coated with a silicone derivative. ZnO %: 79.1–81.5
(*) In the Mexican states of Quintana Roo and Yucatán, biodegradable sunscreens are permitted and the use of sunscreen in sinkholes is regulated in several municipalities. (**) ASEAN (Brunei, Cambodia, Indonesia, Laos, Malaysia, Myanmar, Philippines, Singapore, Thailand, and Vietnam); MERCOSUR (Argentina, Brazil, Paraguay, and Uruguay).
Table 3. Biokinetic models used in different studies.
Table 3. Biokinetic models used in different studies.
UV Filters/MetalsWater TypeField/LaboratoryReference
d C t d t = k u · C W k e · C t
Mussel4-MBC, BP3, BP4, OC, OD-PABASeawaterLaboratory[78]
FishCuFresh waterLaboratory[108]
MusselEHMCSeawaterLaboratory[88]
Clams/musselsZnSeawaterField[26]
d C t d t = k u · C W + A E · I R f o o d   C f o o d k e + g · C t
Oysters/mussels/clamsAg, Cd, Cr, Cu, Ni, Pb, Ti, ZnSeawaterField[106]
MusselAs, Cd, Cr, Cu, Ni, Se, Pb, ZnSeawaterLaboratory[109]
ClamsAg, As, ZnSeawaterLaboratory[110]
MusselCd, Cr, Cu, Ni, ZnFresh waterField[111]
ScallopCd, ZnSeawaterLaboratory[112]
MusselAg, Am, Cd, Co, Se, ZnSeawaterField[113]
d C t d t = k u · C W + A E · I R f o o d   C f o o d + A E · I R s e d   C s e d k e · C t
Clams/prawns/whitingCd, Mn, ZnSeawaterLaboratory[114]
OysterZnSeawaterLaboratory[115]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ruiz-Gutiérrez, G.; Rodríguez-Romero, A.; Tovar-Sánchez, A.; Viguri Fuente, J.R. Analysis and Modeling of Sunscreen Ingredients’ Behavior in an Aquatic Environment. Oceans 2022, 3, 340-363. https://doi.org/10.3390/oceans3030024

AMA Style

Ruiz-Gutiérrez G, Rodríguez-Romero A, Tovar-Sánchez A, Viguri Fuente JR. Analysis and Modeling of Sunscreen Ingredients’ Behavior in an Aquatic Environment. Oceans. 2022; 3(3):340-363. https://doi.org/10.3390/oceans3030024

Chicago/Turabian Style

Ruiz-Gutiérrez, Gema, Araceli Rodríguez-Romero, Antonio Tovar-Sánchez, and Javier R. Viguri Fuente. 2022. "Analysis and Modeling of Sunscreen Ingredients’ Behavior in an Aquatic Environment" Oceans 3, no. 3: 340-363. https://doi.org/10.3390/oceans3030024

Article Metrics

Back to TopTop