Significant rainfall decreases and variations of the atmospheric circulation in the Mediterranean (1950–2000)

This paper aims at identifying the link between significant rainfall decreases (1950–2000) in the Mediterranean basin and the atmospheric circulation at the 500 hPa level. The months and seasons of the subregions with significant rainfall decrease during this period have been identified previously (Norrant and Douguédroit, Theor Appl Climatol 83(1–4):89–106, 2006): October in the Mediterranean Iberia, March in the Atlantic Iberia, January and winter in Greece, and winter in the Near East. Canonical Correlation Analyses based on the monthly and daily data records from 62 rainfall stations and 138 grid points at the 500 hPa level over a Euro-Atlantic window were first calculated to define the TeleConnection Patterns explaining significant regional rainfall decreases. Then, 500 hPa level weather types (ZWTs) of the rainy days with important or little rainfall associated with each Teleconnection Pattern were identified in each subregion. Rainfall-causing disturbances from the Atlantic reach Iberia directly; some of them are regenerated if they reach the Mediterranean. Other disturbances are generated locally near Greece and the Near East (Meteorological Office in Weather in the Mediterranean I: general meteorology, Her Majesty’s Stationery Office, London, 1962). The relationship between significant rainfall decreases and the corresponding 500 hPa level appears to be a nonlinear phenomenon. In all of the studied subregions, a break during the 1970s separates two subperiods differing significantly from each other. Rainfall decrease is due to the higher frequency of important rainfall ZWTs over low rainfall ZWTs, during the first period, which the opposite is true during the second period. Such an inversion could be partially linked with the prevailing North Atlantic Oscillation-positive phase during the last quarter of the twentieth century.     

How changes in diet and trade patterns have shaped the N cycle at the national scale: Spain (1961–2009)

During the last five decades (1961–2009), Spain has experienced a considerable expansion in the nutrient cycle of its agricultural sector and, in particular, a threefold increase in anthropogenic reactive nitrogen inputs, from 536 Gg N year^?1 in 1961–1965 to 1673 Gg N year^?1 in 2005–2009. Import of feed (soybean, cereals, and cakes) from America and Europe to supply a growing livestock population constitutes the largest share of this increase, along with intensification of synthetic fertilizer use. While in the early 1960s, Spain was nearly self-sufficient in terms of food and feed supply, the net import of agricultural products presently equals domestic crop production, when expressed in terms of nitrogen content (ca. 650 Gg N year^?1). The most important driver of this shift appears to be the rapid change in domestic consumption patterns, which evolved from a typical Mediterranean diet to an animal-protein-rich diet similar to the North European and American diets. Besides livestock production mostly for national consumption, the Spanish agricultural system has specialized in vegetal products with low N content such as olive oil, wine, vegetables, and citrus fruit, which are for the most part exported. The nitrogen load exported outside the Spanish borders by rivers is very low (6.5 % of the total net N input). As a result of the high import and low export of reactive nitrogen, the Spanish mainland is suffering from considerable pollution by local emissions of reactive nitrogen forms to air and water.     

Three decades (1983–2010) of contaminant trends in East Greenland polar bears (Ursus maritimus). Part 1: Legacy organochlorine contaminants

Legacy organochlorine contaminants were determined in adipose tissues from 294 polar bears (Ursus maritimus) sampled in East Greenland in 23 of the 28 years between 1983 and 2010. Of 19 major legacy contaminants and congeners (ΣPCB, 4 PCB congeners (CB153, 180, 170/190), ΣDDT, p,p′-DDE, p,p′ -DDD and p,p′-DDT, α- and β-hexachlorocyclohexane (HCH), HCB, octachlorostyrene, dieldrin, oxychlordane, cis- and trans-chlordane, cis- and trans-nonachlor, heptachlor epoxide and BB-153), 18 showed statistically significant average yearly declines of − 4.4% (range: − 2.0 to − 10.8%/year) among subadult polar bears (i.e. females < 5 years, males < 6 years). For example, the ∑ PCB concentrations declined 2.7 fold from 22 730 ng/g lw (95% C.I.: 12 470–32 990 ng/g lw) in 1983–1986 to 8473 ng/g lw (95% C.I.: 6369–9776 ng/g lw) in 2006–2010. Similar but fewer statistically significant trends were found for adult females and adult males likely due to smaller sample size and years. Despite declines as a result of international regulations, relatively high levels of these historic pollutants persist in East Greenland polar bear tissues.     

Three decades (1983–2010) of contaminant trends in East Greenland polar bears (Ursus maritimus). Part 2: Brominated flame retardants

Brominated flame retardants were determined in adipose tissues from 294 polar bears (Ursus maritimus) sampled in East Greenland in 23 of the 28 years between 1983 and 2010. Significant linear increases were found for sum polybrominated diphenyl ether (ΣPBDE), BDE100, BDE153, and hexabromocyclododecane (HBCD). Average increases of 5.0% per year (range: 2.9–7.6%/year) were found for the subadult polar bears. BDE47 and BDE99 concentrations did not show a significant linear trend over time, but rather a significant non-linear trend peaking between 2000 and 2004. The average ΣPBDE concentrations increased 2.3 fold from 25.0 ng/g lw (95% C.I.: 15.3–34.7 ng/g lw) in 1983–1986 to 58.5 ng/g lw (95% C.I.: 43.6–73.4 ng/g lw) in 2006–2010. Similar but fewer statistically significant trends were found for adult females and adult males likely due to smaller sample size and years. Analyses of δ¹⁵N and δ¹³C stable isotopes in hair revealed no clear linear temporal trends in trophic level or carbon source, respectively, and non-linear trends differed among sex and age groups. These increasing concentrations of organobromine contaminants contribute to complex organohalogen mixture, already causing health effects to the East Greenland polar bears.     

Time trends of perfluorinated alkyl acids in serum from Danish pregnant women 2008–2013

We aimed to estimate the levels and time trends of perfluorinated alkyl acids (PFAAs) in serum of 1533 Danish pregnant nulliparous women between 2008 and 2013. The selection criterion of only including nulliparous women was chosen to avoid confounding from parity. The serum samples were analyzed for sixteen PFAAs using solid phase extraction and liquid chromatography tandem mass spectrometry (LC-MS/MS). We investigated the time trends for seven PFAAs, which were detected in more than 50% of the samples: perfluorohexane sulfonate (PFHxS), perfluoroheptane sulfonate (PFHpS), perfluorooctane sulfonate (PFOS), perfluorooctanoic acid (PFOA), perfluorononanoic acid (PFNA), perfluorodecanoic acid (PFDA), and perfluoroundecanoic acid (PFUnA). We found that the serum levels of all seven PFAAs decreased during the period from 2008 to 2013; on average PFHxS decreased with 7.0% per year, PFHpS with 14.8%, PFOS with 9.3%, PFOA with 9.1%, PFNA with 6.2%, PFDA with 6.3%, and PFUnA with 7.1% per year. Adjustment for maternal age, body mass index (BMI), educational level and gestational age at blood sampling did not change the time trends much. To our knowledge, we are the first to report decreasing trends of PFNA, PFDA and PFUnA since year 2000, thereby indicating that the phase-out of these compounds are beginning to show an effect on human serum levels.     

Spatial identification by satellite imagery of the crop–fallow rotation cycle in northern Laos

In the mountainous regions of northern Laos, shifting cultivation, or slash-and-burn agriculture, is widely practiced. However, the crop–fallow rotation cycle is becoming shorter owing to forest conservation policies and population pressure, causing loss of productivity that deleteriously affects farmers’ livelihoods in the region. To investigate regional land use conditions, we have developed a method of identifying the crop–fallow rotation cycle from Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper+ (ETM+) data. We assessed the impact of the identified cycle on plant production measured by Normalized Difference Vegetation Index (NDVI). The study site was an area in Luang Prabang Province. Using eight TM and ETM+ images acquired annually from 1995 to 2003, except for 1998, when cloud-free data were not collected, we classified land use in each year as crop or fallow by the presence of vegetation in the late dry season. Conformity with fallow age determined by field investigation was 69.1%. The cultivation frequency from 1995 to 2002 showed that 77,000 ha (17.3% of the study site) had not been used for cropping during the period, but 41,000 ha (9.2%) had been used every year. Of the study site, 129,000 ha (29.1%) was cultivated one or two times, 83,000 ha (18.7%) was three or four times, and 54,000 ha (12.2%) was five or six times. The NDVI of crops in November did not provide sufficient evidence to prove the assumption that a longer fallow period would result in better crop yields. Instead, the regeneration of fallow vegetation was evidenced by the higher NDVI values after longer fallow. More than 8 years would be needed to reach the same NDVI as forest. From the produced maps indicating fallow age and cultivation frequency, we found that areas with high potential for regeneration decreased as cultivation frequency increased. Areas near rivers were intensively used, and fallow length was accordingly short. Low-potential areas were found in the western basin of the Mekong River. This spatial information can be used to detect areas where biomass productivity is at high risk of deteriorating. Readers should send their comments on this paper to: BhaskarNath@aol.com within 3 months of publication of this issue.     

Changes in the spatial patterns of human appropriation of net primary production (HANPP) in Europe 1990–2006

Understanding patterns, dynamics, and drivers of land use is crucial for improving our ability to cope with sustainability challenges. The human appropriation of net primary production (HANPP) framework provides a set of integrated socio-ecological indicators that quantify how land use alters energy flows in ecosystems via land conversions and biomass harvest. Thus, HANPP enables researchers to systematically and consistently assess the outcome of changes in land cover and land-use intensity across spatio-temporal scales. Yet, fine-scale HANPP assessments are so far missing, an information important to address site-specific ecological implications of land use. Here, we provide such an assessment for Europe at a 1-km scale for the years 1990, 2000, and 2006. The assessment was based on a consistent land-use/biomass flow dataset derived from statistical data, remote sensing maps, and a dynamic global vegetation model. We find that HANPP in Europe amounted to ~43 % of potential productivity, well above the global average of ~25 %, with little variation in the European average since 1990. HANPP was highest in Central Europe and lower in Northern and Southern Europe. At the regional level, distinct changes in land-use intensity were observed, most importantly the decline of cropland areas and yields following the breakdown of socialism in Eastern Europe and the subsequent recovery after 2000, or strong dynamics related to storm events that resulted in massive salvage loggings. In sum, however, these local dynamics cancelled each other out at the aggregate level. We conclude that this finding warrants further research into aspects of the scale-dependency of dynamics and stability of land use.     

Cancer mortality in the meat and delicatessen departments of supermarkets (1950–2006)

Meat cutters and meat wrappers in the meat department of supermarkets are exposed to oncogenic viruses present in raw meat from cattle, pigs, sheep, and poultry, and their products (unpasteurized milk and raw eggs). Up to the mid 1970s, meat wrappers were also exposed to carcinogens present in fumes emitted from the machine used to wrap meat. Because of this we studied cancer mortality in a cohort of 10,701 workers in the meat and delicatessen departments of supermarkets, and we report here the findings after the third follow-up. Standardized mortality ratios (SMR) were estimated in the cohort as a whole and in race/sex subgroups, using the US population for comparison. Study subjects were followed up from January 1950 to December 2006.Significantly increased SMRs of 1.3 (95% CI, 1.2–1.5), and 2.7 (95% CI, 1.2–5.3) were recorded for cancers of the lung, and tonsils/oropharynx, respectively, in the entire cohort, affecting nearly all race/sex subgroups. SMRs of 4.6 (95% CI, 1.0–13.6) for cancer of the floor of the mouth, and 2.8 (95% CI, 1.3–5.3) for cancer of the gall bladder and biliary tract were recorded only in White male meatcutters. Significantly decreased SMRs were observed for a few cancers.It is not known if the observed excess of cancers is a result of occupational exposures. However, substantial evidence points to fumes from the wrapping machine as a possible candidate for explaining the excess in female meat wrappers. Nested case–control studies that can examine risks from occupational exposures in greater detail, and adequately control for confounding factors are now needed, to permit specifically investigate the role of the oncogenic viruses, fumes and non-occupational risk factors in the occurrence of these cancers. The findings are important, not only occupationally but also because the general population may also experience these exposures, albeit to a lesser degree.     

Evolution and frequency (1970–2007) of combined temperature–precipitation modes in the Spanish mountains and sensitivity of snow cover

Snow cover in Spanish mountains is crucial for ensuring water availability in spring and summer months, for the success of winter tourism or for the maintenance of biodiversity in mountain ecosystems. A changing climate may affect the volume of snow cover even in high mountains, where weather conditions are usually favorable for snow accumulation. In this paper, we aim to investigate the evolution (1970–2007) of combined precipitation–temperature modes in the Spanish mountains, and the sensitivity of the snowpack to their occurrence. The climatic database “Spain02” and snow thickness data for Spanish mountains were used for this purpose. Results showed that the frequency of dry-warm and wet-warm days has increased over time in all mountain ranges, while the frequency of the “cold” modes has decreased. The thickness of the snowpack in the Pyrenees has also decreased and its evolution is negatively correlated with the frequency of dry-warm days, and positively correlated with the frequency of dry-cold and wet-cold days. This work constitutes the first approach that relates the evolution of climatic conditions favorable or unfavorable for snow accumulation and the evolution of the snowpack in Spanish mountains.     

Polybrominated diphenyl ethers (PBDEs) in marine mammals from Arctic and North Atlantic regions, 1986–2009

A selection of PBDE congeners was analyzed in pooled blubber samples of pilot whale (Globicephala melas), ringed seal (Phoca hispida), minke whale (Balaenoptera acutorostrata), fin whale (Balaenoptera physalus), harbor porpoise (Phocoena phocoena), hooded seal (Cystophora cristata) and Atlantic white-sided dolphin (Lagenorhynchus acutus), covering a time period of more than 20 years (1986–2009). The analytes were extracted and cleaned-up using open column extraction and multi-layer silica gel column chromatography, and the analysis was performed on a GC-MS system operating in the NCI mode. The highest PBDE levels were found in the toothed whale species pilot whale and white-sided dolphin, and the lowest levels in fin whales and ringed seals. One-sided analyses of variance (ANOVA) followed by Tukey comparisons of means were applied to test for differences between years and sampling areas. Due to inter-year sampling variability, only general comparisons of PBDE concentrations between different sampling areas could be made. Differences in PBDE concentrations between three sampling periods, from 1986 to 2007, were evaluated in samples of pilot whales, ringed seals, white-sided dolphins and hooded seals. The highest PBDE levels were found in samples from the late 1990s or beginning of 2000, possibly reflecting the increase in the global production of technical PBDE mixtures in the 1990s. The levels of BDE #153 and #154 increased relative to the total PBDE concentration in some of the species in recent years, which may indicate an increased relative exposure to higher brominated congeners. In order to assess the effect of measures taken in legally binding international agreements, it is important to continuously monitor POPs such as PBDEs in sub-Arctic and Arctic environments.     

Temporal trends of polychlorinated dibenzo-p-dioxins and dibenzofurans and dioxin-like polychlorinated biphenyls in mothers' milk from Sweden, 1972–2011

Temporal trends of polychlorinated dibenzo-p-dioxins (PCDDs), dibenzofurans (PCDFs) and dioxin-like polychlorinated biphenyls (DL-PCBs) in mothers' milk are still quite rare. Data are particularly scarce when it comes to concentrations from the last decade, 2000 and onwards. The aims of the present study were to assess temporal trends of PCDD, PCDF and DL-PCB in mothers' milk from Stockholm, 1972–2011 and to compare the results with previous analysis of some of the older samples.The samples were analyzed by high resolution GC/MS and results were statistically evaluated for the periods, 1972–2011 and 2002–2011. The rate of which ∑PCDDs, ∑DL-PCBs and the ∑TEQ are decreasing (on pg/g fat WHO-TEQ2005) is higher in the last decade compared to the 40 year period, 1972–2011. A similar trend is indicated, but not confirmed, for ∑TEQ of PCDFs, probably due to too many PCDF congeners below LOQ in the period 2002–2011. Concentrations of ∑PCDDs, ∑PCDFs, ∑DLPCBs and ∑TEQ, all expressed as pg/g fat on TEQ-WHO2005-basis, show a statistically significant decline over time, 5.8–6.8% per year, 1972–2011. The last ten years the annual declines for ∑PCDDs, ∑DL-PCBs and ∑TEQ are 9.2–11% and for ∑PCDF, 5.4%. Congener specific trend analysis, 2002–2001, of PCDDs and DL-PCBs showed the same pattern, while the PCDF congeners showed no such general trend. The results from the re-analysis showed good agreement with slightly lower ∑TEQ1998 pg/g fat concentrations in six out of seven samples and mean difference of 13% in ∑TEQ1998. The study shows that time series can be elongated from previous studies, as long as the sample population remains the same.     

Spatiotemporal variation in rainfall erosivity on the Chinese Loess Plateau during the period 1956–2008

Water resources and soil erosion are the most important environmental concerns on the Chinese Loess Plateau, where soil erosion and sediment yield are closely related to rainfall erosivity. Daily rainfall data from 60 meteorological stations were used to investigate the spatiotemporal variations in annual rainfall, annual erosive rainfall and annual rainfall erosivity on the Chinese Loess Plateau during the period 1956–2008. The annual rainfall, erosive rainfall and rainfall erosivity decreased over the past five decades, as determined by the Mann–Kendall test. A comparison of the annual averaged rainfall, erosive rainfall and rainfall erosivity from 1980 to 2008 with that from 1956 to 1979 revealed a remarkable spatial difference in the rainfall trends on the Loess Plateau. Regions of the plateau with major decreases in rainfall were primarily in the Hekouzhen-Longmen section of the middle Yellow River, especially in the Wuding River basin, the Fenhe basin and the northern-central Shanxi province, where the annual rainfall and erosive rainfall decreased by more than 10% and the annual rainfall erosivity decreased by more than 15%. The rainfall erosivity also decreased more than the annual rainfall. Because the annual rainfall has decreased significantly on the Chinese Loess Plateau over the past 50 years, it is important to better understand the ecological and hydrological processes affected by this climate change.     

Historical carbon fluxes in the expanding deforestation frontier of Southern Brazilian Amazonia (1985–2012)

In tropical areas, pioneer occupation fronts steer the rapid expansion of deforestation, contributing to carbon emissions. Up-to-date carbon emission estimates covering the long-term development of such frontiers depend on the availability of high spatial–temporal resolution data. In this paper, we provide a detailed assessment of carbon losses from deforestation and potential forest degradation from fragmentation for one expanding frontier in the Brazilian Amazon. We focused on one of the Amazonia’s hot-spots of forest loss, the BR-163 highway that connects the high productivity agricultural landscapes in Mato Grosso with the exporting harbors of the Amazon. We used multi-decadal (1984–2012) Landsat-based time series on forested and non-forested area in combination with a carbon book-keeping model. We show a 36% reduction in 1984s biomass carbon stocks, which led to the emission of 611.5 TgCO₂ between 1985 and 1998 (43.6 TgCO₂ year⁻¹) and 959.8 TgCO₂ over 1999–2012 (68.5 TgCO₂ year⁻¹). Overall, fragmentation-related carbon losses represented 1.88% of total emissions by 2012, with an increasing relevance since 2004. We compared the Brazilian Space Agency deforestation assessment (PRODES) with our data and found that small deforestation polygons not captured by PRODES had increasing importance on estimated deforestation carbon losses since 2000. The comparative analysis improved the understanding of data-source-related uncertainties on carbon estimates and indicated disagreement areas between datasets that could be subject of future research. Furthermore, spatially explicit, annual deforestation and emission estimates like the ones derived from this study are important for setting regional baselines for REDD+ or similar payment for ecosystem services frameworks.

     

Comparison of temporal and Spatial Characteristics of Ozone Pollution at Ground Level in the Eastern China

Monitoring data from ozone（O3） automatic stations in three typical cities with different climatic areas in the southern and northern parts of eastern China are used to analyze temporal and spatial characteristics of ozone pollution at ground level. The results show that ozone pollution level has distinct regional differences and the concentration in the suburbs is higher than that in the urban areas. The seasonal variation of ozone concentration in different climatic areas is greatly affected by the variation of precipitation. Ozone concentration in Shenyang and Beijing , in the temperate zone, has one perennial peak concentration, occurring in early summer, May or June. Ozone concentration in Guangzhou, in sub-tropical zone, has two peak values year round. The highest values occur in October and the secondary high value in June. The ozone season in the south is longer than that in the north. The annual average daily peak value of ozone concentrations in different climates usually occur around 3 pm. The diurnal variation range of ozone concentration declines with the increase of latitude. Ozone concentration does not elevate with the increase of traffic flow. Ozone concentration in Guangzhou has a distinct reverse relation to CO and NOx. This complicated non-linearity indicates that the equilibrium of ozone photochemical reaction has regional differences. Exceeding the rate of Beijing＇s lh ozone concentration is higher than that of Guangzhou, whereas the average 8h ozone level is lower than that of Guangzhou, indicating that areas in low latitude are more easily affected by moderate ozone concentrations and longer exposure. Thus, China should work out standards for 8h ozone concentration.     

The Ecological–Genetic Pattern of Variation in the Populations of Narrow-Leaved Hawk's Beard Crepis tectorum L. (Asteraceae)

The effect of long-term exposure to a stress factor on the pattern of variation and intrapopulation differentiation with respect to phenotypic plasticity was studied in Crepis tectorum L. by means of family analysis of various traits in seedlings. Trends in the changes of variation patterns in populations and the possible consequences of the decrease in the proportions of genotypes determining a low sensitivity to environmental changes and in genotypic diversity with respect to phenotypic plasticity are discussed.     

Effect of Summer Monthly Temperatures on Light Tree Ring Formation in Three Larch Species (Larix) in the Northern Forest–Tundra of Siberia

Russian Journal of Ecology - Studies in the northern forest–tundra of Siberia have been performed to estimate temperature conditions in individual months of the growing season that have an...     

Trends in rainfall regime over Israel, 1975–2010, and their relationship to large-scale variability

Variations and trends in the rain regime of Israel are analyzed for 1975–2010, when persistent global warming has been observed. Negative trend is observed over the majority of Israel, statistically significant only in the super-arid region. The decrease is significant over the majority of Israel only in the spring, reflecting a shortening of the rainy season, >3 days/decade. The dry spells are becoming longer, significantly in most of the stations. The factors affecting these variations, synoptic systems, large-scale oscillations and global temperature, were studied for extended period, 1953–2010. A simple multiple stepwise regression model applied for the inter-annual rainfall variations indicates that the occurrence of Cyprus lows is the dominant factor and the Mediterranean oscillation index, MOI2, is also a significant factor. In order to reduce the inter-annual noise and reveal inter-decadal variations, the time-series of the rainfall and its potential predictors were smoothed by 11-year window, showing an increase toward the 1990s, followed by a decrease, at a higher rate, onward. Correspondingly, the aridity lines propagated southward till the mid-1990s and then withdrew back, at a larger rate. The large-scale oscillations and the global temperature explain 83 % of the variance on the inter-decadal time-scale, half of it explained by the global temperature alone. The findings of this study support the expected poleward expansion of the Hadley cell due to global warming.     

Comparison of measured multi-decadal rainfall variability with farmers’ perceptions of and responses to seasonal changes in western Uganda

Regional Environmental Change - Smallholder farmers in Sub-Saharan Africa (SSA) are not only dealing with decreased production from land degradation, but are also impacted heavily by climate...     

Polybrominated diphenyl ethers (PBDEs) in human samples of mother–newborn pairs in South China and their placental transfer characteristics

There are limited data concerning the placenta transfer characteristics and accumulation of polybrominated diphenyl ethers (PBDEs) in infants. However, PBDEs received increasing health concerns due to their endocrine disrupt and neurodevelopment toxicity effects. The present study assessed the accumulation of PBDEs in 30 paired placenta, breast milk, fetal cord blood, and neonatal urine samples collected from five major cities of the South China. The age of mothers ranged from 21 to 39 (mean 27.6 ± 4.56). The ∑ PBDE concentrations were 15.8 ± 9.88 ng g^− 1 lipid in placenta, 13.2 ± 7.64 ng g^− 1 lipid in breast milk, 16.5 ± 19.5 ng g^− 1 lipid in fetal cord blood, and 1.80 ± 1.99 ng ml^− 1 in neonatal urine. BDE-47 was the predominant congener in all types of human sample. Octa-BDEs such as BDE-196/-197 were detected highly in placenta and cord blood while moderately in breast milk and neonatal urine. Significant (p < 0.01) correlations were observed for both total and most individual PBDEs in cord blood–maternal placenta and breast milk–urine paired individual samples. The extent of placental transfer of higher brominated BDEs such as BDE-196/-197 was greater than that of BDE-47. The estimated daily intake (EDI) analysis for breast-fed infants revealed that newborns in these areas were exposed to relatively high levels of PBDEs via breast milk. Our study not only provided systematic fundamental data for PBDE distribution but also revealed the placenta transfer characteristics of PBDE congeners in South China.