Data to Decision in a Dynamic Ocean: Robust Species Distribution Models and Spatial Decision Frameworks

Human use of the oceans is increasingly in conflict with conservation of endangered species. Methods for managing the spatial and temporal placement of industries such as military, fishing, transportation and offshore energy, have historically been post hoc; i.e. the time and place of human activity is often already determined before assessment of environmental impacts. In this dissertation, I build robust species distribution models in two case study areas, US Atlantic (Best et al. 2012) and British Columbia (Best et al. 2015), predicting presence and abundance respectively, from scientific surveys. These models are then applied to novel decision frameworks for preemptively suggesting optimal placement of human activities in space and time to minimize ecological impacts: siting for offshore wind energy development, and routing ships to minimize risk of striking whales. Both decision frameworks relate the tradeoff between conservation risk and industry profit with synchronized variable and map views as online spatial decision support systems.

For siting offshore wind energy development (OWED) in the U.S. Atlantic (chapter 4), bird density maps are combined across species with weights of OWED sensitivity to collision and displacement and 10 km2 sites are compared against OWED profitability based on average annual wind speed at 90m hub heights and distance to transmission grid. A spatial decision support system enables toggling between the map and tradeoff plot views by site. A selected site can be inspected for sensitivity to a cetaceans throughout the year, so as to capture months of the year which minimize episodic impacts of pre-operational activities such as seismic airgun surveying and pile driving.

Routing ships to avoid whale strikes (chapter 5) can be similarly viewed as a tradeoff, but is a different problem spatially. A cumulative cost surface is generated from density surface maps and conservation status of cetaceans, before applying as a resistance surface to calculate least-cost routes between start and end locations, i.e. ports and entrance locations to study areas. Varying a multiplier to the cost surface enables calculation of multiple routes with different costs to conservation of cetaceans versus cost to transportation industry, measured as distance. Similar to the siting chapter, a spatial decisions support system enables toggling between the map and tradeoff plot view of proposed routes. The user can also input arbitrary start and end locations to calculate the tradeoff on the fly.

Essential to the input of these decision frameworks are distributions of the species. The two preceding chapters comprise species distribution models from two case study areas, U.S. Atlantic (chapter 2) and British Columbia (chapter 3), predicting presence and density, respectively. Although density is preferred to estimate potential biological removal, per Marine Mammal Protection Act requirements in the U.S., all the necessary parameters, especially distance and angle of observation, are less readily available across publicly mined datasets.

In the case of predicting cetacean presence in the U.S. Atlantic (chapter 2), I extracted datasets from the online OBIS-SEAMAP geo-database, and integrated scientific surveys conducted by ship (n=36) and aircraft (n=16), weighting a Generalized Additive Model by minutes surveyed within space-time grid cells to harmonize effort between the two survey platforms. For each of 16 cetacean species guilds, I predicted the probability of occurrence from static environmental variables (water depth, distance to shore, distance to continental shelf break) and time-varying conditions (monthly sea-surface temperature). To generate maps of presence vs. absence, Receiver Operator Characteristic (ROC) curves were used to define the optimal threshold that minimizes false positive and false negative error rates. I integrated model outputs, including tables (species in guilds, input surveys) and plots (fit of environmental variables, ROC curve), into an online spatial decision support system, allowing for easy navigation of models by taxon, region, season, and data provider.

For predicting cetacean density within the inner waters of British Columbia (chapter 3), I calculated density from systematic, line-transect marine mammal surveys over multiple years and seasons (summer 2004, 2005, 2008, and spring/autumn 2007) conducted by Raincoast Conservation Foundation. Abundance estimates were calculated using two different methods: Conventional Distance Sampling (CDS) and Density Surface Modelling (DSM). CDS generates a single density estimate for each stratum, whereas DSM explicitly models spatial variation and offers potential for greater precision by incorporating environmental predictors. Although DSM yields a more relevant product for the purposes of marine spatial planning, CDS has proven to be useful in cases where there are fewer observations available for seasonal and inter-annual comparison, particularly for the scarcely observed elephant seal. Abundance estimates are provided on a stratum-specific basis. Steller sea lions and harbour seals are further differentiated by ‘hauled out’ and ‘in water’. This analysis updates previous estimates (Williams & Thomas 2007) by including additional years of effort, providing greater spatial precision with the DSM method over CDS, novel reporting for spring and autumn seasons (rather than summer alone), and providing new abundance estimates for Steller sea lion and northern elephant seal. In addition to providing a baseline of marine mammal abundance and distribution, against which future changes can be compared, this information offers the opportunity to assess the risks posed to marine mammals by existing and emerging threats, such as fisheries bycatch, ship strikes, and increased oil spill and ocean noise issues associated with increases of container ship and oil tanker traffic in British Columbia’s continental shelf waters.

Starting with marine animal observations at specific coordinates and times, I combine these data with environmental data, often satellite derived, to produce seascape predictions generalizable in space and time. These habitat-based models enable prediction of encounter rates and, in the case of density surface models, abundance that can then be applied to management scenarios. Specific human activities, OWED and shipping, are then compared within a tradeoff decision support framework, enabling interchangeable map and tradeoff plot views. These products make complex processes transparent for gaming conservation, industry and stakeholders towards optimal marine spatial management, fundamental to the tenets of marine spatial planning, ecosystem-based management and dynamic ocean management.

Passive Acoustics: A Multifaceted Tool for Marine Mammal Conservation

Marine mammals exploit the efficiency of sound propagation in the marine environment for essential activities like communication and navigation. For this reason, passive acoustics has particularly high potential for marine mammal studies, especially those aimed at population management and conservation. Despite the rapid realization of this potential through a growing number of studies, much crucial information remains unknown or poorly understood. This research attempts to address two key knowledge gaps, using the well-studied bottlenose dolphin (Tursiops truncatus) as a model species, and underwater acoustic recordings collected on four fixed autonomous sensors deployed at multiple locations in Sarasota Bay, Florida, between September 2012 and August 2013. Underwater noise can hinder dolphin communication. The ability of these animals to overcome this obstacle was examined using recorded noise and dolphin whistles. I found that bottlenose dolphins are able to compensate for increased noise in their environment using a wide range of strategies employed in a singular fashion or in various combinations, depending on the frequency content of the noise, noise source, and time of day. These strategies include modifying whistle frequency characteristics, increasing whistle duration, and increasing whistle redundancy. Recordings were also used to evaluate the performance of six recently developed passive acoustic abundance estimation methods, by comparing their results to the true abundance of animals, obtained via a census conducted within the same area and time period. The methods employed were broadly divided into two categories – those involving direct counts of animals, and those involving counts of cues (signature whistles). The animal-based methods were traditional capture-recapture, spatially explicit capture-recapture (SECR), and an approach that blends the “snapshot” method and mark-recapture distance sampling, referred to here as (SMRDS). The cue-based methods were conventional distance sampling (CDS), an acoustic modeling approach involving the use of the passive sonar equation, and SECR. In the latter approach, detection probability was modelled as a function of sound transmission loss, rather than the Euclidean distance typically used. Of these methods, while SMRDS produced the most accurate estimate, SECR demonstrated the greatest potential for broad applicability to other species and locations, with minimal to no auxiliary data, such as distance from sound source to detector(s), which is often difficult to obtain. This was especially true when this method was compared to traditional capture-recapture results, which greatly underestimated abundance, despite attempts to account for major unmodelled heterogeneity. Furthermore, the incorporation of non-Euclidean distance significantly improved model accuracy. The acoustic modelling approach performed similarly to CDS, but both methods also strongly underestimated abundance. In particular, CDS proved to be inefficient. This approach requires at least 3 sensors for localization at a single point. It was also difficult to obtain accurate distances, and the sample size was greatly reduced by the failure to detect some whistles on all three recorders. As a result, this approach is not recommended for marine mammal abundance estimation when few recorders are available, or in high sound attenuation environments with relatively low sample sizes. It is hoped that these results lead to more informed management decisions, and therefore, more effective species conservation.

Bridging the Scale Gap from Leaf to Canopy in Biosphere-Atmosphere Gas and Particle Exchanges

Terrestrial ecosystems, occupying more than 25% of the Earth's surface, can serve as

`biological valves' in regulating the anthropogenic emissions of atmospheric aerosol

particles and greenhouse gases (GHGs) as responses to their surrounding environments.

While the signicance of quantifying the exchange rates of GHGs and atmospheric

aerosol particles between the terrestrial biosphere and the atmosphere is

hardly questioned in many scientic elds, the progress in improving model predictability,

data interpretation or the combination of the two remains impeded by

the lack of precise framework elucidating their dynamic transport processes over a

wide range of spatiotemporal scales. The diculty in developing prognostic modeling

tools to quantify the source or sink strength of these atmospheric substances

can be further magnied by the fact that the climate system is also sensitive to the

feedback from terrestrial ecosystems forming the so-called `feedback cycle'. Hence,

the emergent need is to reduce uncertainties when assessing this complex and dynamic

feedback cycle that is necessary to support the decisions of mitigation and

adaptation policies associated with human activities (e.g., anthropogenic emission

controls and land use managements) under current and future climate regimes.

With the goal to improve the predictions for the biosphere-atmosphere exchange

of biologically active gases and atmospheric aerosol particles, the main focus of this

dissertation is on revising and up-scaling the biotic and abiotic transport processes

from leaf to canopy scales. The validity of previous modeling studies in determining

iv

the exchange rate of gases and particles is evaluated with detailed descriptions of their

limitations. Mechanistic-based modeling approaches along with empirical studies

across dierent scales are employed to rene the mathematical descriptions of surface

conductance responsible for gas and particle exchanges as commonly adopted by all

operational models. Specically, how variation in horizontal leaf area density within

the vegetated medium, leaf size and leaf microroughness impact the aerodynamic attributes

and thereby the ultrane particle collection eciency at the leaf/branch scale

is explored using wind tunnel experiments with interpretations by a porous media

model and a scaling analysis. A multi-layered and size-resolved second-order closure

model combined with particle

uxes and concentration measurements within and

above a forest is used to explore the particle transport processes within the canopy

sub-layer and the partitioning of particle deposition onto canopy medium and forest

oor. For gases, a modeling framework accounting for the leaf-level boundary layer

eects on the stomatal pathway for gas exchange is proposed and combined with sap

ux measurements in a wind tunnel to assess how leaf-level transpiration varies with

increasing wind speed. How exogenous environmental conditions and endogenous

soil-root-stem-leaf hydraulic and eco-physiological properties impact the above- and

below-ground water dynamics in the soil-plant system and shape plant responses

to droughts is assessed by a porous media model that accommodates the transient

water

ow within the plant vascular system and is coupled with the aforementioned

leaf-level gas exchange model and soil-root interaction model. It should be noted

that tackling all aspects of potential issues causing uncertainties in forecasting the

feedback cycle between terrestrial ecosystem and the climate is unrealistic in a single

dissertation but further research questions and opportunities based on the foundation

derived from this dissertation are also brie

y discussed.

(Table 1) Soil description of pits dug in coastal, inland and higland Antarctic Dry Valleys

In the past few years, it has become increasingly apparent that perchlorate (ClO4-) is present on all continents, except the polar regions where it had not yet been assessed, and that it may have a significant natural source. Here, we report on the discovery of perchlorate in soil and ice from several Antarctic Dry Valleys (ADVs) where concentrations reach up to 1100/µg/kg. In the driest ADV, perchlorate correlates with atmospherically deposited nitrate. Far from anthropogenic activity, ADV perchlorate provides unambiguous evidence that natural perchlorate is ubiquitous on Earth. The discovery has significant implications for the origin of perchlorate, its global biogeochemical interactions, and possible interactions with the polar ice sheets. The results support the hypotheses that perchlorate is produced globally and continuously in the Earth's atmosphere, that it typically accumulates in hyperarid areas, and that it does not build up in oceans or other wet environments most likely because of microbial reduction on a global scale.

Evaluation of the Genotoxic and Physiological Effects of Decabromodiphenyl Ether (BDE-209) and Dechlorane Plus (DP) Flame Retardants in Marine Mussels (Mytilus galloprovincialis)

Dechlorane Plus (DP) is a proposed alternative to the legacy flame retardant decabromodiphenyl ether (BDE-209), a major component of Deca-BDE formulations. In contrast to BDE-209, toxicity data for DP are scarce and often focused on mice. Validated dietary in vivo exposure of the marine bivalve (Mytilus galloprovincialis) to both flame retardants did not induce effects at the physiological level (algal clearance rate), but induced DNA damage, as determined by the comet assay, at all concentrations tested. Micronuclei formation was induced by both DP and BDE-209 at the highest exposure concentrations (100 and 200 mu g/L, respectively, at 18% above controls). DP caused effects similar to those by BDE-209 but at lower exposure concentrations (5.6, 56, and 100 mu g/L for DP and 56, 100, and 200 mu g/L for BDE-209). Moreover, bioaccumulation of DP was shown to be concentration dependent, in contrast to BDE-209. The results described suggest that DP poses a greater genotoxic potential than BDE-209

Evaluation of the Genotoxic and Physiological Effects of Decabromodiphenyl Ether (BDE-209) and Dechlorane Plus (DP) Flame Retardants in Marine Mussels (Mytilus galloprovincialis)

Dechlorane Plus (DP) is a proposed alternative to the legacy flame retardant decabromodiphenyl ether (BDE-209), a major component of Deca-BDE formulations. In contrast to BDE-209, toxicity data for DP are scarce and often focused on mice. Validated dietary in vivo exposure of the marine bivalve (Mytilus galloprovincialis) to both flame retardants did not induce effects at the physiological level (algal clearance rate), but induced DNA damage, as determined by the comet assay, at all concentrations tested. Micronuclei formation was induced by both DP and BDE-209 at the highest exposure concentrations (100 and 200 mu g/L, respectively, at 18% above controls). DP caused effects similar to those by BDE-209 but at lower exposure concentrations (5.6, 56, and 100 mu g/L for DP and 56, 100, and 200 mu g/L for BDE-209). Moreover, bioaccumulation of DP was shown to be concentration dependent, in contrast to BDE-209. The results described suggest that DP poses a greater genotoxic potential than BDE-209

Review of the flood risk management system in Germany after the major flood in 2013

Widespread flooding in June 2013 caused damage costs of €6 to 8 billion in Germany, and awoke many memories of the floods in August 2002, which resulted in total damage of €11.6 billion and hence was the most expensive natural hazard event in Germany up to now. The event of 2002 does, however, also mark a reorientation toward an integrated flood risk management system in Germany. Therefore, the flood of 2013 offered the opportunity to review how the measures that politics, administration, and civil society have implemented since 2002 helped to cope with the flood and what still needs to be done to achieve effective and more integrated flood risk management. The review highlights considerable improvements on many levels, in particular (1) an increased consideration of flood hazards in spatial planning and urban development, (2) comprehensive property-level mitigation and preparedness measures, (3) more effective flood warnings and improved coordination of disaster response, and (4) a more targeted maintenance of flood defense systems. In 2013, this led to more effective flood management and to a reduction of damage. Nevertheless, important aspects remain unclear and need to be clarified. This particularly holds for balanced and coordinated strategies for reducing and overcoming the impacts of flooding in large catchments, cross-border and interdisciplinary cooperation, the role of the general public in the different phases of flood risk management, as well as a transparent risk transfer system. Recurring flood events reveal that flood risk management is a continuous task. Hence, risk drivers, such as climate change, land-use changes, economic developments, or demographic change and the resultant risks must be investigated at regular intervals, and risk reduction strategies and processes must be reassessed as well as adapted and implemented in a dialogue with all stakeholders.

Influence of environmental conditions on the age, hatch dates, and growth of juvenile Atlantic menhaden in the Choptank River, MD

Since 1993 Atlantic menhaden has experienced sustained low juvenile production (recruitment) in the Chesapeake Bay. Factors controlling growth, abundance, and mortality of larval and juvenile menhaden change throughout ontogeny such that larval growth rates could carry over to juvenile growth and survival. The effects of winter thermal conditions on the hatch dates and growth of larval and juvenile Atlantic menhaden in Atlantic shelf and Chesapeake Bay habitats were examined using otolith (ear-stone) increment analyses and growth models. For 2010-2013, truncated hatch-date distributions provided evidence for a winter recruitment bottleneck in Atlantic menhaden caused by cold temperatures. Hatch-dates of surviving juveniles were skewed towards warmer months for years characterized by colder temperatures. Reduced larval growth rates, influenced by reduced temperature and food availability, carried over to juvenile growth rates. A growing degree-day model performed well in simulating observed juvenile growth rates in the Choptank River tributary of Chesapeake Bay.

Assessing the influence of abiotic factors and leaf-level properties on the stability of growing-season canopy greenness in a deciduous forest

Maps depicting spatial pattern in the stability of summer greenness could advance understanding of how forest ecosystems will respond to global changes such as a longer growing season. Declining summer greenness, or “greendown”, is spectrally related to declining near-infrared reflectance and is observed in most remote sensing time series to begin shortly after peak greenness at the end of spring and extend until the beginning of leaf coloration in autumn,. Understanding spatial patterns in the strength of greendown has recently become possible with the advancement of Landsat phenology products, which show that greendown patterns vary at scales appropriate for linking these patterns to proposed environmental forcing factors. This study tested two non-mutually exclusive hypotheses for how leaf measurements and environmental factors correlate with greendown and decreasing NIR reflectance across sites. At the landscape scale, we used linear regression to test the effects of maximum greenness, elevation, slope, aspect, solar irradiance and canopy rugosity on greendown. Secondly, we used leaf chemical traits and reflectance observations to test the effect of nitrogen availability and intrinsic water use efficiency on leaf-level greendown, and landscape-level greendown measured from Landsat. The study was conducted using Quercus alba canopies across 21 sites of an eastern deciduous forest in North America between June and August 2014. Our linear model explained greendown variance with an R2=0.47 with maximum greenness as the greatest model effect. Subsequent models excluding one model effect revealed elevation and aspect were the two topographic factors that explained the greatest amount of greendown variance. Regression results also demonstrated important interactions between all three variables, with the greatest interaction showing that aspect had greater influence on greendown at sites with steeper slopes. Leaf-level reflectance was correlated with foliar δ13C (proxy for intrinsic water use efficiency), but foliar δ13C did not translate into correlations with landscape-level variation in greendown from Landsat. Therefore, we conclude that Landsat greendown is primarily indicative of landscape position, with a small effect of canopy structure, and no measureable effect of leaf reflectance. With this understanding of Landsat greendown we can better explain the effects of landscape factors on vegetation reflectance and perhaps on phenology, which would be very useful for studying phenology in the context of global climate change

Estimating Population Trends in American Woodcock (Scolopax Minor) Using Population Reconstruction Models

The American woodcock (Scolopax minor) population index in North America has declined 0.9% a year since 1968 prompting managers to identify priority information and management needs for the species (Sauer et al 2008). Managers identified a need for a population model that better informs on the status of American woodcock populations (Case et al. 2010). Population reconstruction techniques use long-term age-at-harvest data and harvest effort to estimate abundances with error estimates. Four new models were successfully developed using survey data (1999 to 2013). The optimal model estimates sex specific harvest probability for adult females at 0.148 (SE = 0.017) and all other age-sex cohorts at 0.082 (SE = 0.008) for the most current year 2013. The model estimated a yearly survival rate of 0.528 (SE = 0.008). Total abundance ranged from 5,206,000 woodcock in 2007 to 6,075,800 woodcock in 1999. This study represents the first population estimates of woodcock populations.

An Ecological Analysis of the Potential for Moss-Based Green Roof Design

Green roofs are a maturing application of best management practices for controlling urban stormwater runoff. The majority of green roofs are planted with drought resistant, higher plant species, such as the genus Sedum. However, other plant varieties, such as mosses, may be equally applicable. Residential roofs and natural terrestrial communities were sampled in both Maryland and Tennessee to determine moss community structure and species water composition. This served as a natural analog for potential green roof moss communities. During sampling, 21 species of moss were identified throughout the 37 total sites. The average percent moss cover and water composition across all roof sites was 40.7% and 38.6%, respectively and across all natural sites, 76.7% and 47.7%, respectively. Additional maximum water holding capacity procedures were completed on sedum and 19 of the 21 sampled moss species to assess their individual potential for stormwater absorption. Sedum species on average held 166% of their biomass in water, while moss species held 732%. The results of this study are used as a basis to propose moss species that will improve green roof performance.

BALD EAGLES (HALIAEETUS LEUCOCEPHALUS) AS INDICATORS OF GREAT LAKES ECOSYSTEM HEALTH

Environmental indicators have been proposed as a means to assess ecological integrity, monitoring both chemical and biological stressors. In this study, we used nestling bald eagles as indicators to quantify direct or indirect tertiary-level contaminant exposure. The spatial and temporal trends of polychlorinated biphenyl (PCB) congeners were evaluated in nestling plasma from 1999–2014. Two hexa-chlorinated congeners, PCB-138 and 153, were detected with the highest frequency and greatest concentrations throughout Michigan. Less-chlorinated congeners such as PCB-52 and 66 however, comprised a greater percentage of total PCB concentrations in nestlings proximate to urbanized areas, such as along the shorelines of Lake Erie. Toxic equivalents were greatest in the samples collected from nestlings located on Lake Erie, followed by the other Great Lakes spatial regions. Nestling plasma samples were also used to measure concentrations of the most heavily-used group of flame retardants, brominated diphenyl ethers (BDEs), and three groups of alternative flame retardants, non-BDE Brominated Flame Retardants (NBFRS), Dechloranes, and organophosphate esters (OPs). BDE-47, 99 and 100 contributed the greatest to total BDE concentrations. Concentrations of structurally similar NBFRs found in this study and recent atmospheric studies indicate that they are largely used as replacements to previously used BDE mixtures. A variety of Dechloranes, or derivatives of Mirex and Dechlorane Plus, were measured. Although, measured at lesser concentrations, environmental behavior of these compounds may be similar to mirex and warrant future research in aquatic species. Concentrations of OPs in nestling plasma were two to three orders of magnitude greater than all other groups of flame retardants. In addition to chemical indicators, bald eagles have also been proposed as indicators to identify ecological stressors using population measures that are tied to the fitness of individuals and populations. Using mortality as a population vitality rate, vehicle collisions were found to be the main source of mortality with a greater incidence for females during white-tailed deer (Odocoileus virginianus) hunting months and spring snow-melt. Lead poisoning was the second greatest source of mortality, with sources likely due to unretrieved hunter-killed, white-tailed deer carcasses, and possibly exacerbated by density-dependent effects due to the growing population in Michigan.

Caffeine As An Indicator Of Estrogenic Activity In Source Water.

Caffeine has already been used as an indicator of anthropogenic impacts, especially the ones related to the disposal of sewage in water bodies. In this work, the presence of caffeine has been correlated with the estrogenic activity of water samples measured using the BLYES assay. After testing 96 surface water samples, it was concluded that caffeine can be used to prioritize samples to be tested for estrogenic activity in water quality programs evaluating emerging contaminants with endocrine disruptor activity.

Electrochemical Oxidation Of Landfill Leachate In A Flow Reactor: Optimization Using Response Surface Methodology.

Response surface methodology based on Box-Behnken (BBD) design was successfully applied to the optimization in the operating conditions of the electrochemical oxidation of sanitary landfill leachate aimed for making this method feasible for scale up. Landfill leachate was treated in continuous batch-recirculation system, where a dimensional stable anode (DSA(©)) coated with Ti/TiO2 and RuO2 film oxide were used. The effects of three variables, current density (milliampere per square centimeter), time of treatment (minutes), and supporting electrolyte dosage (moles per liter) upon the total organic carbon removal were evaluated. Optimized conditions were obtained for the highest desirability at 244.11 mA/cm(2), 41.78 min, and 0.07 mol/L of NaCl and 242.84 mA/cm(2), 37.07 min, and 0.07 mol/L of Na2SO4. Under the optimal conditions, 54.99 % of chemical oxygen demand (COD) and 71.07 ammonia nitrogen (NH3-N) removal was achieved with NaCl and 45.50 of COD and 62.13 NH3-N with Na2SO4. A new kinetic model predicted obtained from the relation between BBD and the kinetic model was suggested.

Changes in optical properties caused by UV-irradiation of aquatic humic substances from the amazon river basin: Seasonal variability evaluation

Aquatic humic substances (AHS) isolated from two characteristic seasons of the Negro river, winter and summer corresponding to floody and dry periods, were structurally characterized by (13)C nuclear magnetic ressonance. Subsequently, AHS aqueous solutions were irradiated with a polychromatic lamp (290-475 nm) and monitored by its total organic carbon (TOC) content, ultraviolet-visible (UV-vis) absorbance, fluorescence and Fourier transformed infrared spectroscopy (FTIR). As a result, a photobleaching upto 80% after irradiation of 48 h was observed. Conformational rearrangements and formation of low molecular complexity structures were formed during the irradiation, as deduced from the pH decrement and the fluorescence shifting to lower wavelengths. Additionally a significant mineralization with the formation Of CO(2), CO, and inorganic carbon compounds was registered, as assumed by TOC losses of up to 70%. The differences in photodegradation between samples expressed by photobleaching efficiency were enhanced in the summer sample and related to its elevated aromatic content. Aromatic structures are assumed to have high autosensitization capacity effects mediated by the free radical generation from quinone and phenolic moieties.