• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Soundscape Dynamics at Anuran Reproductive Sites in Pannonian Biogeographical Region: Effects of Road Noise on Vocal Activity

    2016-09-28 06:55:22AndrWEIPERTHEdSMITHSzilviaSIMIGLAMiklPUKYandYezhongTANG
    Asian Herpetological Research 2016年1期

    András WEIPERTH, Ed SMITH, Szilvia SIMIGLA, Miklòs PUKYand Yezhong TANG

    1MTA Centre for Ecological Research, Danube Research Institute, H-1113 Budapest, Karolina út 29, Hungary

    2Department of Psychology, University of Maryland, College Park, MD 20742, USA

    3Chengdu Institute of Biology, Chinese Academy of Sciences, # 9 of Section 4, South Renmin Road, Chengdu 610041,Sichuan, China

    ?

    Soundscape Dynamics at Anuran Reproductive Sites in Pannonian Biogeographical Region: Effects of Road Noise on Vocal Activity

    András WEIPERTH1, Ed SMITH2*, Szilvia SIMIGLA1, Miklòs PUKY1and Yezhong TANG3*

    1MTA Centre for Ecological Research, Danube Research Institute, H-1113 Budapest, Karolina út 29, Hungary

    2Department of Psychology, University of Maryland, College Park, MD 20742, USA

    3Chengdu Institute of Biology, Chinese Academy of Sciences, # 9 of Section 4, South Renmin Road, Chengdu 610041,Sichuan, China

    The emerging field of soundscape ecology views ecosystems in terms of biophony, geophony and anthrophony. Soundscape ecology considers the effects of sound on fauna, and this research focuses on anuran breeding lek soundscapes. The sensitivity of anuran breeding leks to acoustic disturbances makes breeding leks an important venue for a comparative soundscape study. We made long-term (> 24 h) sound recordings in three representative wetlands and short-term (< 30 min) recordings in ten sites in the Pannonian Biogeographical Region of Hungary and around the Hungary and Slovakia border. Long-term soundscapes of the fl oodplain stretch, where there is relatively minor anthrophonical disturbance, showed an obvious circadian change in sound intensities. The site with moderate sound contamination exhibited a disturbed pattern of circadian sound variation, while the site with heavy traffi c noise displayed an apparently random temporal soundscape. At different amphibian breeding sites during mating season, our short-term recordings were dominated by anuran calls, bird songs and wind noises, while insect calls and rain were present to a lesser degree. Our study indicates that vehicle traffi c noise is a severe imposition to the natural soundscape,and suggests that soundscape monitoring can provide a reliable and sensitive index of environmental change for both short-term and long-term periods.

    soundscape monitoring, anuran breeding site, biophony, geophony, anthrophony

    1. Introduction

    Soundscape ecology, the science of sound in the landscape, is an emerging field which encompasses the causes and consequences of biological (biophony),geophysical (geophony), and human-produced(anthrophony) sounds (Pijanowski et al., 2011) to understand coupled animal and human dynamics across different scales of distance and time. It is an integrative framework that aims to describe how climate, land transformation, biodiversity patterns, and humanactivities interact through time to form dynamic acoustic landscapes. Monitoring and studying soundscapes may illuminate physical mechanisms for ecological processes and identify courses of landscape change accurately,sensitively, and economically.

    Furthermore, the landscape has been reconceived as a dynamic system composed of matter, structured energy,information and meaning (Cosgrove, 2003; Farina, 2010),thus expanding upon the more classical, geographicalecological oriented perspective (Risser et al., 1984;Forman and Godron, 1986; Pickett and Cadenasso, 1995;Wu and Hobbs, 2002; Turner, 2005). In detail, sound produced in the landscape derives from various sources including human, weather, geophysical, and bioacoustic sources (Francis et al., 2011). Soundscape ecology overlaps with landscape ecology since some ecological processes occurring within landscapes can be tightly linked to and refl ected in patterns of soundcape (Formanand Godron, 1981; Urban et al., 1987; Turner, 1989;Turner et al., 2001; Farina, 2006).

    In bioacoustics, four animal taxa are well known for intense acoustic emissions: birds, most anurans, some insects, and a few mammals. Most of these animals produce intense sounds during their breeding seasons to attract potential mates and to repel rivals, and they are usually silent at other times. Those sounds are generally the main components of the local soundscape in areas such as ponds and leks (Runkle et al., 1994; Catchpole and Slater, 2003; Farina et al., 2011; Wang et al., 2012). The severe consequences of anthropogenic noise on wildlife have been shown recently over a diverse array of taxa (Barber et al., 2010). Population density of frogs is negatively related to road traffi c which could be attributed partially to the reducement in attractiveness of vocal display by traffic noise since anuran chorus behavior might be affected by man-made acoustic interference either directly through modulating call rates of the chorus participants or indirectly, through suppressing calling behavior of one set of species which in turn stimulated calling in other species (Sun and Narins, 2005;Cunnington and Fahrig, 2010). Nevertheless, ecological changes in response to noise at broad-scales have not yet been examined or tracked over time. In the present study,we compare soundscapes of three wetlands where frogs reproduce in differing levels of anthropogenic noise, and report acoustic analyses of soundscapes at numerous sites around the Pannonian Biogeographical Region.

    2. Materials and Methods

    2.1 Acoustic recording Study areas include fl oodplains,wetlands and ponds in the Pannonian Biogeographical Region in Hungary and Slovakia where the fire bellied toad (Bombina bombina) is the dominant species in most sampling sites from fl oodplains along the River Ipoly to water habitats near Lake Balaton during our recording time. The common spadefoot toad (Pelobates fuscus),European green toad (Bufotes viridis), European tree frog(Hyla arborea) and water frogs (Pelophylax esculentus)are common anuran species while some of them advertised vocally earlier or later than the period when we recorded. Long-term recordings were made at two sites along the River Ipoly (Hugyag and Hont), northeast Hungary and one site at Lake Balaton (Balatonederics) in order to investigate the temporal changes in soundscapes in the term of circadian period. Short-term recordings were made at seven places in total with some places containing two ponds (i.e. sites) (Budapest, Hont, Hugyag, Ipolydamázsd, Ipolyság/Hont, Ipolysz?g,Letkés) (Figure 1) with the purpose of surveying changes among different studied sites.

    Figure 1 Acoustic recording sites in the Pannonian Biogeographical Region, triangle: only long-term recording (Balatonederics), circle:only short-term recordings (Budapest, Ipolydamázsd, Ipolyság/ Hont, Ipolysz?g, Letkés), rhombus: long- and short-term recordings.

    Long-term acoustic recordings were made with a portable recorder (Sony, Japan) placed in a nearby tree and oriented to the center of the wetland. The recording volume was fi xed at the level of 25 and the soundscapes were recorded continuously for more than 24 hours. The frequency response of the Sony recorder was from 100 Hz to 16 000 Hz with the sample rate was set to 44 000 Hz. Fifteen-minute short-term recordings were made with the Marantz PMD670/U1B recorder (USA) connected to an AE3300 microphone (Audio-Technica, USA) at each of ten sites. The recording sensitivity, via the gain knob,was set specifically for each recording to optimize the signal noise ratio. The frequency response of the Marantz recorder was fl at to 20 000 Hz, +/- 0.25 dB. Gain knob settings, address, time, temperature, relative humidity and GPS information were recorded at each site.

    All recordings were completed in June 1-15, 2013 with each site recorded 1-2 times. Five anuran species' calls were recorded in the Pannonian Biogeographical Region:B. bombina, H. arborea, P. esculentus, P. ridibundus,and P. lessonae, the latter two species contributing only slightly to the soundscapes.

    2.2 Data analyses For measuring the intensity, creating the sonogram and analyzing the acoustic component,long-term recordings were fi rst segmented manually into sixty-minute sections. Then, the fi rst fi ve-minute of each one-hour segment was analyzed further. Because the soundscape was changing quite slowly, the five-minute segments represented the complete hour to a large extent. PRAAT (an open-source program released by Universityof Amsterdam) was used to measure relative sound intensities in dB and to create sonograms. For intensity measurement, the “To Intensity” function was used with“Down to Intensity Tier” and “Down to Table Of Real”operations. These relative intensity data were saved as txt files which in turn were loaded into MS Excel in order to calculate means for each segment. PRAAT was also used to create sonograms using the “Analyse Spectrum” function with a Hanning window length of 30 milliseconds. Sonograms were displayed conventionally --as two-dimensional figures (x-axis: time; y-axis:frequency) with warmer hues indicating frequencyspecifi c energy. Sounds were identifi ed by experimenters through visual inspection of the sonograms and listening to the recordings.

    For each site, SPL was derived from the short-term recordings. The constant sensitivity of the microphone and the variable sensitivity of the recorder were determined in a calibration step, and then the SPL at each site was determined in a measurement step. In the calibration step,a 1-KHz tone was played into the microphone at 71.3 dB SPL unweighted, measured by Bruel & Kajer 2250G Integrating Sound Level Meter. Reference recordings were made for a range of recorder gain knob settings(“gain settings”) of 3 through 7, in steps of 0.5, and digital RMS levels were computed for each reference recording. Gain settings were plotted against the RMS levels, as shown in Figure 2. An interpolating polynomial (made with Matlab) was fi tted to the data to allow computation of digital RMS levels corresponding to the reference 71.3 dB tone for any gain knob setting between 3.0 and 7.5(Rossing, 1990).

    Figure 2 Plots of recorded root mean square (RMS) in volts to knob settings of the recorder used to compute SPL of site soundscapes.

    In the measurement step, representative segments of the site recordings were chosen, and average RMS levels were computed for each segment. For each segment, the linear ratio of the segment RMS level to the RMS level of the reference tone recording was computed, taking into account the gain setting used at each recording site. The site-specifi c unweighted SPL levels were then computed by converting the linear ratios to dB levels, and adding these dB levels to the original SPL measured during the calibration step (71.3 dB). A-weighted SPL levels were computed in exactly the same way, except that the sitespecifi c segments were fi ltered with an A-weighting fi lter(in Matlab) before the corresponding digital RMS levels were computed.

    3. Results

    3.1 Different circadian patterns of soundscap Circadian soundscape variations were examined at three wetland sites: the Hugyag site, the Hont site, and the Balaton site. The Hugyag site (N 48°05'874"; E 19°26'533"; H 147 m) is located in a remote border area with almost no anthropogenic noise. Similarly, the Hont site (N 48°03'494"; E 18°58'264"; H 119 m) is located along the floodplain of River Ipoly, near the border between Hungary and Slovakia but, unlike at Hugyag, it is only 120 meters away from highway E77, so there is heavy traffic noise. The Balaton site (N 46°48'231"; E 17° 24'226"; H 110 m) is situated along Lake Balaton at Balatonederics, and has a medium level of traffi c noise.

    The soundscape in Hugyag varied daily in intensity from 50 to 80 dB. The peak sound intensity occurred precisely at 21:30, but no obvious intensity valley could be seen. Low intensities started at 4:30 and lasted to 13:30, with variations between 50.4 and 62.5 dB ( 3A). The relative sound intensities in Hont varied irregularly between 70 and 90 dB (Figure 3B). The fact that the minimal intensity in Hont was higher than that in Hugyag was attributed to the difference in high background of traffic noise. At Lake Balaton a 24-hour variation in sound intensities could be found even with masking of the traffi c noise (Figure 3C). Because of the moderate traffi c noise, the lowest sound intensity at Balaton was higher than that at Hugyag but lower than that at Hont.

    Figure 3 Circadian changes in sound intensity from three wetlands where anurans reproduce. A. Evident variation in pattern of sine wave in Hugyag; B. Irregular pattern of alteration in Hont; C. A pattern with regular wave disturbed by random noises in Lake Balaton.

    3.2 Temporal changes in biological components The dominant species of anuran communities in the Pannonian Biogeographical Region is Bombina bombina,a poisonous toad. At Hugyag this species produced advertisement calls nocturnally and diurnally, while circadian variations in intensity peaked around 20:00. The main call energy was concentrated at 470 Hz (Figure 4A). The species of Hyla arborea contributed largely to the soundscape from 21:00 to 1:00, displacing B. bombina as the loudest call, with dominant frequency around 2600 Hz (Figure 4B). At Lake Balaton Pelophylax esculentus and P. ridibundus called simultaneously and formed a chorus consisting of two frog species and traffi c cars. Interestingly, vocal activities of the two Pelophylax species were evoked frequently by the traffic noise(Figure 4C). Crickets and rain contributed some energy to the soundscape at the Balaton site. Soundscape at Hont was consisted mainly of bird songs which were masked largely with traffic noises while P. ridibundus produced calls around and after the midnight, contributing slightly to the overall soundscape.

    3.3 Varied structures of soundscape at different sites Soundscapes at the ten sites in the Pannonian Biogeographical Region show variation in mean intensity,temperature, relative humidity and major components(Table 1). For the anuran breeding sites, the principal bioacoustic sources were some anuran species, while birds, insects, and other frogs were minor sound sources(Figure 4D). Other acoustic sources were direct wind and wind in plants. We found no correlations among biological sounds and environmental elements such as temperature and relative humidity (for homogeneity test for variance and binary regression, p values > 0.05); this was probably due to differences in species composition at each site.

    4. Discussion

    In spatial dimension, a diverse array of sounds produced by mammals, birds, amphibians, and insects might be the main components of the soundscape in forests,grasslands and wetlands (Marler and Slabbekoorn, 2004),while the urban soundscapes are composed of sounds generated by vehicles, machines and other humanproduced sounds (Botteldooren et al., 2004; Raimbault and Dubois, 2005). Abiological parts of the soundscape include gushing rivers flowing over terrain, rain falling through canopies, and wind (Swanson et al., 1988). Time scales of the soundscape vary daily, seasonally and annually in habitats (Tang et al., 2001; Wang et al.,2012), reflecting circadian, reproductive periods, and habitat and/or climate changes, respectively. Soundscapes change dramatically as environments change, and animal vocalizations account for most of these changes. In addition, long-term ecological changes in landscape, e.g.,those accompanying desertization, global climate change,construction of transportation thoroughfares and other human activities, are reflected by soundscape changes too.

    It is well known that most anurans, if not all, compete for mate selection through vocalization (Kelley, 2004). Wetlands and/or ponds are sites for lekking by anurans,where acoustic communication, competition, mating, egg laying, and tadpole development occur. There is usually high acoustic background noise, since other animals make sounds around these sites. Circadian and seasonal changes in anuran vocalization can be expected in order to mitigate the interference effects of bio-noise.

    Many frog-eating waders forage at the sites where we set up our study from 9:00 to 19:00. Bombina species secrete poison which protects them from birds, thusallowing this toad to make advertisement calls day and night with a slight decrease during bird predation. H. arborea, P. esculentus and P. ridibundus produced calls most intensively after midnight when birds are at rest. Why these two Pelophylax species overlap their calling times, and the breadth of sounds that are observed to elicit calling from the species are both enigmas. At the sites, calling was often initiated by traffi c noise, and we could cause males to start calling by orally mimicking their calls (unpublished data). Many frog species vocalize in the form of chorus in which one male's calls generally stimulate other males to produce calls (Ryan et al.,1981; Fang et al., 2013). It is likely that male frogs are easily inducible behaviorally with less discrimination of stimulation structures while females are usually fastidious, and thus males can be evoked to call by many sounds including the traffi c noise.

    Table 1 Sound intensities and acoustic components of the soundscape recorded at different times from different sites within the Pannonian region.

    Figure 4 Waveforms and spectrograms of advertisement calls from three anuran species. A. Calls of Bombina bombina; B. Calls of Hyla arborea and B. bombina; C. Calls of Pelophylax esculentus and P. ridibundus; D. Sound recorded from a small pond to show acoustic components. To clearly depict these different vocalizations, sonogram frequencies are from 0 to 5k Hz for A, 0 to 10k Hz for B and C, and 0 to 8.5k Hz for D.

    Man-made noise from vehicles and machines might vary in a circadian rhythm, but not with a seasonal period. No site of nature reserves in the continental US is free from this man-made noise (Barber et al., 2011), and the same is most probably true for Europe. An international road, busy with trucks day and night, crosses the border between Hungary and Slovakia near the Hont site. A national road, occupied with relatively few cars at night,runs along Lake Balaton, while no roads exist near the Hugyag site. Soundscapes at these three sites exhibit different temporal patterns, mostly correlated with the traffi c noise. Our study indicates a large infl uence of noise contamination on the anuran bioacoustic components in the soundscape, which would mask the auditory signals of anurans (Bee and Swanson, 2007). A parallel situation exists for birds (Francis et al., 2011).

    It has been demonstrated that anthrophony infl uenced negatively reproductive successes of vocal speices through masking acoustic signals for sexual displays(Sun and Narins, 2005; Lengagne, 2008; Cunnington and Fahrig, 2010; Halfwerk et al., 2011). In addition,migratory birds showed a change in ability to gain body condition during migratory stopover when stayed along a“phantom road” in spite of the noise (Ware et al., 2015). Some vocal animals could, however, adjust their call parameters in avoidance of the anthrophonic noise by upshifting frequency in bird (Slabbekoorn and Peet, 2003)and frog (Parris et al., 2009), enhancing intensity, i.e. Lombard effect in bird (Brumm, 2004) and prolonging call duration in frog (Love and Bee, 2010). In contrast,Hyla males were incapable of adjusting their temporal or frequency call structures to increase efficiency of the vocal communication in the noise environment(Lengagne, 2008). The further study should be necessary to investigate the soundscape by measuring acoustic parameters as many as possible not just intensity.

    Acknowledgement This work was supported fi nancially by the Program of Exchange Visit between Chinese Academy of Sciences and Hungarian Academy of Sciences and the National Natural Science Foundation of China (NSFC 31272304 to TYZ).

    References

    Barber J. R., Crooks K. R , Fristrup K. M. 2010. The costs of chronic noise exposure for terrestrial organisms. Trends Ecol Evol, 25: 180-189

    Barber J. P., Burdett C. L., Reed S. E., Warner K. A.,F(xiàn)ormichella C., Crooks K. R., Theobald D. M., Fristrup K. M. 2011. Anthropogenic noise exposure in protected natural areas: estimating the scale of ecological consequences. Landscape Ecol, 26: 1281-1295

    Bee M. A., Swanson E. M. 2007. Auditory masking of anuran advertisement calls by road traffi c noise. Anim Behav, 74: 1765-1776

    Botteldooren D., Coensel B. D., Meur T. D. 2004. The temporal structure of the urban soundscape. J Sound Vib, 292: 105-123

    Brumm H. 2004. The impact of environmental noise on song amplitude in a territorial bird. J Anim Ecol, 73: 434-440

    Catchpole C. K., Slater P. J. B. 2003. Bird song: Biological themes and variations. 2nd Edition. Cambradge Press

    Cosgrove D. 2003. Landscape: ecology and semiosis. In: Palang H,F(xiàn)ry G (eds), Landscape interfaces: cultural heritage in changing landscapes. Kluwer, Dordrecht, 15-20

    Cunnington G. M., Fahrig L. 2010. Plasticity in the vocalizations of anurans in response to traffi c noise. Acta Oecol, 36(5): 463-470

    Fang G. Z., Jiang F., Yang P., Cui J. G., Brauth S. E., Tang Y. Z. 2013. Male vocal competition is dynamic and strongly affected by social contexts in music frogs. Anim Cogn, 17(2): 483-494

    Farina A. 2006. Principles and methods in landscape ecology. Springer, NY

    Farina A. 2010. Ecology, cognition and landscape. Springer,Dordrecht

    Farina A., Lattanzi E., Malavasi R., Pieretti N., Piccioli L. 2011. Avian soundscapes and cognitive landscapes: Theory, application and ecological perspectives. Landscape Ecol, 26: 1257-1267

    Forman R. T. T., Godron M. 1981. Patches and structural components for a landscape ecology. BioScience, 31: 733-740

    Forman R. T. T., Godron M. 1986. Landscape ecology. John Wiley, New York

    Francis C. D., Paritsis J., Ortega C. P., Cruz A. 2011. Landscape patterns of avian habitat use and nesting success resulting from chronic gas well compressor noise in NW New Mexico, USA. Landscape Ecol. doi:10.1007/s10980-011-9609-z

    Halfwerk W., Holleman L. J. M., Lessells C. M., Lessells M.,Slabbekoorn H. 2011. Negative impact of traffi c noise on avian reproductive success. Journal of Applied Ecology, 48(1): 210-219.

    Slabbekoorn H., Peet M. 2003. Ecology: Birds sing at a higher pitch in urban noise. Nature, 424(6946): 267

    Lengagne T. 2008. Traffi c noise affects communication behaviour in a breeding anuran, Hyla arborea. Biol Conserv, 141(8): 2023-2031

    Love, E. K., Bee, M. A. 2010. An experimental test of noisedependent voice amplitude regulation in Cope's grey treefrog,Hyla chrysoscelis. Anim behav, 80(3): 509-515

    Kelley D. B. 2004. Vocal communication in frogs. Curr Opin Neurobiol, 14: 751-757

    Marler P., Slabbekoorn H. 2004. Nature's music: the science of birdsong. Elsevier Academic Press, San Diego, USA

    Parris K. M., Velik-Lord M., North J. M. A. 2009. Frogs Call at a Higher Pitch in Traffi c Noise. Ecol Soc, 14(1): 124-124

    Pickett S. T. A , Cadenasso M. L. 1995. Landscape ecology: spatial heterogeneity in ecological systems. Science, 269: 331-334

    Pijanowski B. C., Farina A., Gage S. H., Dumyahn S. L., Krause B. L. 2011. What is soundscape ecology? An introduction and overview of an emerging new science. Landscape Ecol, 26:1213-1232

    Raimbault M., Dubois D. 2005. Urban soundscapes: experiences and knowledge. Cities, 22(5): 339-350

    Risser P. G., Karr J. R., Forman R. T. T. 1984. Landscape ecology: Directions and approaches. Illinois Natural History Survey Special Publication 2, Champaign

    Rossing T. D. 1990. The Science of Sound. 2nd Edition. Addison-Wesley Publishing Company, 92-93

    Runkle L. S., Wells K. D., Robb C. C., Lance S. L. 1994. Individual, nightly, and seasonal variation in calling behavior of the gray tree frog, Hyla versicolor: implications for energy expenditure. Behav Ecol, 5: 318-325

    Ryan M. J., Tuttle M. D., Taft L. K. 1981. The costs and benefi ts of frog chorusing behavior. Behav Ecol Sociobiol, 8(4): 273-278

    Swanson F. J., Kratz T. K., Caine N., Woodmansee R. G. 1988. Landform effects on ecosystem patterns and processes. BioScience, 38(2): 92-98

    Sun J. W., Narins P. M. 2005. Anthropogenic sounds differentially affect amphibian call rate. Biol Conserv, 121(3): 419-427

    Tang Y. Z., Zhuang L. Z., Wang Z. W. 2001. Advertisement Calls and Their Relation to Reproductive Cycles in Gekko gecko(Reptilia, Lacertilia). Copeia, 1: 248-253

    Turner M. G. 1989. Landscape ecology: the effect of pattern on process. Annu Rev Ecol Syst, 20: 171-197

    Turner M. G., Gardner R. H., O'Neill R. V. 2001. Landscape ecology in theory and practice: pattern and process. Springer Press, New York

    Turner M. G. 2005. Landscape ecology: what is the state of the science? Annu Rev Ecol Syst, 36: 319-344

    Urban D. L., O'Neill R. V., Shugart H. H. 1987. Landscape ecology. BioScience, 37: 119-127

    Wang, J. C., Cui J. G., Shi H. T., Brauth S. E., Tang Y. Z. 2012. Effects of Body Size and Environmental Factors on the Acoustic Structure and Temporal Rhythm of Calls in Rhacophorus dennysi. Asian Herpetol Res, 3: 205-212

    Ware H. E., Mcclure C. J. W., Carlisle J. D., Barber J. R. 2015. A phantom road experiment reveals traffi c noise is an invisible source of habitat degradation. P Natl Acad Sci USA, 39: 12105-12109

    Wu J., Hobbs R. 2002. Key issues and research priorities in landscape ecology: an idiosyncratic synthesis. Landscape Ecol,17: 355-365

    *Corresponding authors: Dr. Ed SMITH, from Department of Psychology, University of Maryland, USA, with his research focusing on bioacoustics; Prof. Yezhong TANG, from Chengdu Institute of Biology, Chinese Academy of Sciences, Chengdu, China, with his research focusing on neurobiology of amphibians and reptiles.

    E-mail: fastcared@yahoo.com (E. SMITH); tangyz@cib.ac.cn (Y. Z. TANG)

    22 September 2015 Accepted: 25 December 2015

    亚洲欧美精品综合一区二区三区| 噜噜噜噜噜久久久久久91| 禁无遮挡网站| 男女床上黄色一级片免费看| 国产真人三级小视频在线观看| av在线蜜桃| 精品久久久久久久人妻蜜臀av| 亚洲成a人片在线一区二区| 哪里可以看免费的av片| 91av网一区二区| 国内精品久久久久精免费| 别揉我奶头~嗯~啊~动态视频| 露出奶头的视频| 成人高潮视频无遮挡免费网站| 亚洲av成人av| 婷婷六月久久综合丁香| 1000部很黄的大片| 一级毛片高清免费大全| 三级毛片av免费| 99re在线观看精品视频| 不卡av一区二区三区| 国内精品久久久久精免费| 成人三级做爰电影| 亚洲欧美日韩卡通动漫| 亚洲欧美日韩无卡精品| 国产精品国产高清国产av| 在线观看舔阴道视频| 精品国内亚洲2022精品成人| 婷婷丁香在线五月| 最近在线观看免费完整版| 成人无遮挡网站| 免费人成视频x8x8入口观看| 色在线成人网| 少妇的逼水好多| 国产人伦9x9x在线观看| 精品久久久久久成人av| 精品免费久久久久久久清纯| 村上凉子中文字幕在线| 国产亚洲精品久久久久久毛片| 一本综合久久免费| 非洲黑人性xxxx精品又粗又长| а√天堂www在线а√下载| 日本 av在线| 国内精品美女久久久久久| 国产熟女xx| 在线播放国产精品三级| 日韩有码中文字幕| 日日干狠狠操夜夜爽| 久久久久久久久久黄片| 欧美日本亚洲视频在线播放| 在线a可以看的网站| 丁香六月欧美| 精品午夜福利视频在线观看一区| av女优亚洲男人天堂 | 午夜福利在线观看吧| 日本a在线网址| 日本a在线网址| 天天躁狠狠躁夜夜躁狠狠躁| 久久国产精品影院| 中文资源天堂在线| 搡老熟女国产l中国老女人| 韩国av一区二区三区四区| 亚洲第一电影网av| 精品无人区乱码1区二区| 真人一进一出gif抽搐免费| 午夜福利视频1000在线观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲成人精品中文字幕电影| 亚洲自偷自拍图片 自拍| 国产av麻豆久久久久久久| 伦理电影免费视频| 性色av乱码一区二区三区2| 中文字幕人成人乱码亚洲影| 亚洲国产精品999在线| 免费在线观看影片大全网站| 小说图片视频综合网站| 波多野结衣巨乳人妻| 日日夜夜操网爽| 久久久久国产一级毛片高清牌| 97超级碰碰碰精品色视频在线观看| 91久久精品国产一区二区成人 | 亚洲在线观看片| 国产精品永久免费网站| 精品99又大又爽又粗少妇毛片 | 午夜精品久久久久久毛片777| 国模一区二区三区四区视频 | 欧美激情在线99| 国产伦精品一区二区三区视频9 | 久久草成人影院| 亚洲无线在线观看| 国产野战对白在线观看| 精品电影一区二区在线| 伊人久久大香线蕉亚洲五| 91av网一区二区| 久久热在线av| 久久午夜综合久久蜜桃| 级片在线观看| 亚洲人成网站高清观看| 欧美另类亚洲清纯唯美| 男人和女人高潮做爰伦理| 宅男免费午夜| 美女免费视频网站| 男人舔女人下体高潮全视频| 丰满的人妻完整版| 国产三级在线视频| 夜夜看夜夜爽夜夜摸| 亚洲一区二区三区色噜噜| 每晚都被弄得嗷嗷叫到高潮| 欧美不卡视频在线免费观看| 日日夜夜操网爽| 午夜成年电影在线免费观看| 老司机午夜十八禁免费视频| 午夜精品在线福利| 免费人成视频x8x8入口观看| 欧美性猛交╳xxx乱大交人| 亚洲av免费在线观看| 一区福利在线观看| 国产野战对白在线观看| 国产精品综合久久久久久久免费| 久久香蕉精品热| 白带黄色成豆腐渣| 国产精品亚洲美女久久久| 日日摸夜夜添夜夜添小说| 天堂影院成人在线观看| 少妇的丰满在线观看| 精品无人区乱码1区二区| 午夜福利视频1000在线观看| 欧美成人一区二区免费高清观看 | 级片在线观看| 欧美黄色淫秽网站| 嫁个100分男人电影在线观看| 观看美女的网站| 最近最新中文字幕大全免费视频| 在线观看午夜福利视频| 最近在线观看免费完整版| 一级作爱视频免费观看| 嫩草影院入口| tocl精华| 久久久久性生活片| 欧美丝袜亚洲另类 | 国产欧美日韩一区二区三| 日韩欧美三级三区| 日本 av在线| 少妇人妻一区二区三区视频| 无遮挡黄片免费观看| 久久欧美精品欧美久久欧美| 国内揄拍国产精品人妻在线| 久久久久久久久久黄片| 1024香蕉在线观看| 色综合站精品国产| 成年女人永久免费观看视频| 天天添夜夜摸| 国产成人精品久久二区二区91| 欧美激情在线99| a级毛片在线看网站| 听说在线观看完整版免费高清| 九色成人免费人妻av| 舔av片在线| 日日干狠狠操夜夜爽| 亚洲精品中文字幕一二三四区| av在线天堂中文字幕| 久久精品91无色码中文字幕| 在线观看免费午夜福利视频| 亚洲人成电影免费在线| 亚洲成av人片免费观看| 亚洲五月婷婷丁香| 高清毛片免费观看视频网站| bbb黄色大片| 精品久久蜜臀av无| 淫妇啪啪啪对白视频| 日日夜夜操网爽| 国产亚洲精品一区二区www| 女同久久另类99精品国产91| 伊人久久大香线蕉亚洲五| 男女做爰动态图高潮gif福利片| 美女大奶头视频| 好男人电影高清在线观看| 国产高清视频在线播放一区| 少妇人妻一区二区三区视频| 宅男免费午夜| 伊人久久大香线蕉亚洲五| 18禁观看日本| www.自偷自拍.com| 色在线成人网| 久久人人精品亚洲av| 亚洲国产精品sss在线观看| 免费看a级黄色片| 久久久国产欧美日韩av| 美女大奶头视频| 成人特级黄色片久久久久久久| 黄色成人免费大全| 欧美日韩综合久久久久久 | 午夜免费男女啪啪视频观看| 国内精品一区二区在线观看| 一个人看视频在线观看www免费| 色网站视频免费| 亚洲第一区二区三区不卡| 自拍偷自拍亚洲精品老妇| 免费观看a级毛片全部| 高清av免费在线| 国产成人a∨麻豆精品| 国产黄片视频在线免费观看| 91久久精品国产一区二区成人| 最近手机中文字幕大全| 国产人妻一区二区三区在| 18禁动态无遮挡网站| 国产免费又黄又爽又色| 国产成人一区二区在线| 国产成人freesex在线| 天天一区二区日本电影三级| 大话2 男鬼变身卡| 26uuu在线亚洲综合色| 国产精品一区二区三区四区久久| 亚洲av成人av| 亚洲国产欧美人成| 久久久久久大精品| 久久午夜福利片| 老司机影院成人| 中文字幕精品亚洲无线码一区| 99热这里只有是精品在线观看| 直男gayav资源| 一级毛片aaaaaa免费看小| 国产高清不卡午夜福利| 精品人妻视频免费看| 国产精品久久久久久精品电影| 亚洲国产精品sss在线观看| 亚洲怡红院男人天堂| 一级毛片电影观看 | 少妇熟女欧美另类| 午夜精品在线福利| 99热全是精品| 国产av一区在线观看免费| 亚洲在线自拍视频| 在线播放无遮挡| 日本黄大片高清| 亚洲欧美成人综合另类久久久 | 乱码一卡2卡4卡精品| 国产在视频线在精品| 国产探花在线观看一区二区| 麻豆成人av视频| 国产免费一级a男人的天堂| 最新中文字幕久久久久| 久久久成人免费电影| 国产一区亚洲一区在线观看| 最近视频中文字幕2019在线8| 色播亚洲综合网| 99热6这里只有精品| 日韩视频在线欧美| 一个人看的www免费观看视频| 国产午夜福利久久久久久| 日日啪夜夜撸| www日本黄色视频网| 狂野欧美激情性xxxx在线观看| 亚洲最大成人手机在线| 嘟嘟电影网在线观看| 欧美成人精品欧美一级黄| 中文字幕精品亚洲无线码一区| 亚洲精品乱码久久久v下载方式| 秋霞伦理黄片| 国产精品人妻久久久久久| 久久鲁丝午夜福利片| h日本视频在线播放| 日韩av不卡免费在线播放| 亚洲精品色激情综合| 嫩草影院精品99| 久久99热这里只频精品6学生 | 亚洲最大成人手机在线| 美女高潮的动态| 高清在线视频一区二区三区 | 精品国产露脸久久av麻豆 | 亚洲人成网站在线播| 国产在视频线在精品| 国产伦一二天堂av在线观看| 1000部很黄的大片| 蜜臀久久99精品久久宅男| 日韩精品有码人妻一区| 久久久久久国产a免费观看| 丰满少妇做爰视频| 午夜精品国产一区二区电影 | av.在线天堂| 精华霜和精华液先用哪个| 九九热线精品视视频播放| 午夜免费激情av| 午夜福利成人在线免费观看| 亚洲人成网站高清观看| 国产精品99久久久久久久久| 免费av不卡在线播放| 国产免费一级a男人的天堂| 天天躁夜夜躁狠狠久久av| 哪个播放器可以免费观看大片| 九九久久精品国产亚洲av麻豆| 成年女人永久免费观看视频| 国产亚洲精品久久久com| 午夜久久久久精精品| 一级毛片aaaaaa免费看小| 国产成人精品久久久久久| 在线免费观看的www视频| 嫩草影院入口| 18禁动态无遮挡网站| 午夜精品一区二区三区免费看| 国产精品一区二区三区四区免费观看| 亚洲精品aⅴ在线观看| 99热网站在线观看| 国产精品久久视频播放| 成年免费大片在线观看| 国产欧美日韩精品一区二区| 免费av不卡在线播放| 国产一区二区三区av在线| 国产高清视频在线观看网站| 毛片女人毛片| 中文在线观看免费www的网站| 小蜜桃在线观看免费完整版高清| 看免费成人av毛片| 亚洲欧美精品自产自拍| 久久精品久久久久久噜噜老黄 | 精品一区二区免费观看| 日韩三级伦理在线观看| 美女内射精品一级片tv| 久久人人爽人人爽人人片va| 免费电影在线观看免费观看| 少妇丰满av| 九九热线精品视视频播放| 视频中文字幕在线观看| 欧美zozozo另类| 日本三级黄在线观看| 99久久精品一区二区三区| 亚洲成av人片在线播放无| 国内精品美女久久久久久| 久久久亚洲精品成人影院| 深夜a级毛片| 99国产精品一区二区蜜桃av| 日韩 亚洲 欧美在线| 亚洲人成网站高清观看| 亚洲在线观看片| 不卡视频在线观看欧美| 日韩欧美精品v在线| 国产高清视频在线观看网站| 99在线人妻在线中文字幕| 综合色丁香网| 久久韩国三级中文字幕| 午夜福利视频1000在线观看| 国产精品一二三区在线看| 精品国内亚洲2022精品成人| 18+在线观看网站| 久久亚洲国产成人精品v| 久久久久久九九精品二区国产| 欧美日韩一区二区视频在线观看视频在线 | 久久久久久久久久久丰满| 成人欧美大片| 国产一级毛片七仙女欲春2| 日韩欧美 国产精品| 91aial.com中文字幕在线观看| 久久这里只有精品中国| 国产精品精品国产色婷婷| 91精品伊人久久大香线蕉| 国产精品久久视频播放| 欧美成人免费av一区二区三区| 久久人妻av系列| videos熟女内射| .国产精品久久| 两个人的视频大全免费| 长腿黑丝高跟| 亚洲婷婷狠狠爱综合网| 成人二区视频| 久久人妻av系列| 中文资源天堂在线| h日本视频在线播放| 国产精品久久久久久精品电影| 丰满乱子伦码专区| 久久精品国产99精品国产亚洲性色| 亚洲丝袜综合中文字幕| 亚洲丝袜综合中文字幕| 亚洲激情五月婷婷啪啪| 色吧在线观看| 成年女人永久免费观看视频| 人妻系列 视频| 免费人成在线观看视频色| 国产女主播在线喷水免费视频网站 | or卡值多少钱| 国产成人a∨麻豆精品| 中文字幕精品亚洲无线码一区| 在线播放国产精品三级| 国产精品麻豆人妻色哟哟久久 | 狂野欧美激情性xxxx在线观看| 欧美一区二区亚洲| 久久久久精品久久久久真实原创| 亚洲av一区综合| 人妻系列 视频| 国产成人aa在线观看| 又爽又黄a免费视频| 丝袜美腿在线中文| 我的老师免费观看完整版| 久久人人爽人人片av| 午夜精品国产一区二区电影 | 亚洲欧美精品综合久久99| 国产成人a区在线观看| 一卡2卡三卡四卡精品乱码亚洲| 一级黄色大片毛片| 亚洲欧美成人综合另类久久久 | 水蜜桃什么品种好| 内地一区二区视频在线| 色网站视频免费| 国产亚洲5aaaaa淫片| 最近中文字幕高清免费大全6| 免费在线观看成人毛片| 中文字幕av在线有码专区| 91aial.com中文字幕在线观看| 97超视频在线观看视频| 最后的刺客免费高清国语| 国产精品蜜桃在线观看| 亚洲三级黄色毛片| 国产三级中文精品| 麻豆国产97在线/欧美| 干丝袜人妻中文字幕| 成人午夜精彩视频在线观看| 天堂av国产一区二区熟女人妻| 亚洲真实伦在线观看| 青青草视频在线视频观看| 亚洲国产精品成人综合色| 韩国av在线不卡| 国产一区有黄有色的免费视频 | 三级国产精品欧美在线观看| 久久久精品94久久精品| 免费搜索国产男女视频| 亚洲五月天丁香| 免费在线观看成人毛片| 大香蕉久久网| 3wmmmm亚洲av在线观看| 国产不卡一卡二| 日韩高清综合在线| 又爽又黄无遮挡网站| 中国美白少妇内射xxxbb| 亚洲自拍偷在线| 成人特级av手机在线观看| 99久久精品国产国产毛片| 国内揄拍国产精品人妻在线| 久久精品综合一区二区三区| 99热这里只有是精品在线观看| 中文字幕久久专区| 熟妇人妻久久中文字幕3abv| 2022亚洲国产成人精品| 菩萨蛮人人尽说江南好唐韦庄 | 免费av观看视频| 亚洲av免费高清在线观看| 欧美丝袜亚洲另类| 亚洲国产精品专区欧美| 久久99蜜桃精品久久| 天天躁日日操中文字幕| 欧美一级a爱片免费观看看| 日韩欧美 国产精品| 亚洲精品色激情综合| 中文乱码字字幕精品一区二区三区 | 亚洲精品日韩在线中文字幕| 国产av不卡久久| a级毛色黄片| 成人高潮视频无遮挡免费网站| 亚洲精品色激情综合| 亚洲一级一片aⅴ在线观看| 日本免费一区二区三区高清不卡| 五月玫瑰六月丁香| 国产精品久久久久久av不卡| 亚洲av电影不卡..在线观看| kizo精华| 精品人妻偷拍中文字幕| 综合色av麻豆| 非洲黑人性xxxx精品又粗又长| 国产成人精品一,二区| 国产精品久久久久久久电影| 亚洲18禁久久av| 毛片女人毛片| 麻豆成人av视频| 久久99精品国语久久久| 国产亚洲午夜精品一区二区久久 | 尾随美女入室| 老师上课跳d突然被开到最大视频| 国产白丝娇喘喷水9色精品| 日本五十路高清| 狂野欧美激情性xxxx在线观看| 久久韩国三级中文字幕| 日韩视频在线欧美| 中文天堂在线官网| 岛国在线免费视频观看| 亚洲av电影不卡..在线观看| 国产免费福利视频在线观看| 免费电影在线观看免费观看| 男插女下体视频免费在线播放| 亚洲成av人片在线播放无| 神马国产精品三级电影在线观看| 国产av不卡久久| 99在线视频只有这里精品首页| 黄色欧美视频在线观看| 黄色日韩在线| 热99在线观看视频| 能在线免费看毛片的网站| 亚洲欧美成人精品一区二区| 成人毛片60女人毛片免费| 久久99热这里只有精品18| 91久久精品国产一区二区成人| 国产精品久久久久久久久免| 直男gayav资源| 日产精品乱码卡一卡2卡三| 七月丁香在线播放| 午夜精品在线福利| 国产亚洲av片在线观看秒播厂 | 偷拍熟女少妇极品色| 级片在线观看| 边亲边吃奶的免费视频| 永久网站在线| 久久久久久久午夜电影| 精品一区二区三区视频在线| 国产伦在线观看视频一区| 精华霜和精华液先用哪个| 热99在线观看视频| av在线播放精品| 乱码一卡2卡4卡精品| 国产不卡一卡二| 午夜福利视频1000在线观看| 看片在线看免费视频| 国产三级中文精品| 99在线人妻在线中文字幕| 色视频www国产| 看免费成人av毛片| 听说在线观看完整版免费高清| 一卡2卡三卡四卡精品乱码亚洲| 欧美区成人在线视频| 青青草视频在线视频观看| 国产精华一区二区三区| 不卡视频在线观看欧美| 亚洲精品,欧美精品| 狠狠狠狠99中文字幕| 国产伦在线观看视频一区| 久久99热这里只有精品18| 特大巨黑吊av在线直播| 久久久久久久久久久免费av| 国产精品精品国产色婷婷| 欧美丝袜亚洲另类| 免费观看精品视频网站| 亚洲av免费高清在线观看| 一本一本综合久久| 免费无遮挡裸体视频| 爱豆传媒免费全集在线观看| 国产精品一二三区在线看| 国产免费一级a男人的天堂| 最近视频中文字幕2019在线8| 亚洲三级黄色毛片| 国产成人freesex在线| av免费在线看不卡| 晚上一个人看的免费电影| 夜夜看夜夜爽夜夜摸| 久久午夜福利片| 久久鲁丝午夜福利片| 国产私拍福利视频在线观看| 亚洲av电影不卡..在线观看| 汤姆久久久久久久影院中文字幕 | 久久久国产成人免费| 国产日韩欧美在线精品| 久久精品夜夜夜夜夜久久蜜豆| 精品午夜福利在线看| 久久国内精品自在自线图片| 国产真实乱freesex| 国产亚洲精品av在线| 中文字幕熟女人妻在线| 18禁动态无遮挡网站| 全区人妻精品视频| ponron亚洲| 天美传媒精品一区二区| 亚洲成人久久爱视频| 日韩成人伦理影院| 91久久精品国产一区二区三区| 亚洲一区高清亚洲精品| 三级男女做爰猛烈吃奶摸视频| av卡一久久| 黄色欧美视频在线观看| 高清日韩中文字幕在线| 一级毛片我不卡| 国产精品永久免费网站| 国产极品精品免费视频能看的| 国产视频首页在线观看| 亚洲av中文字字幕乱码综合| 国产精品无大码| av卡一久久| 成人一区二区视频在线观看| 亚洲人成网站在线播| 在线播放国产精品三级| 国产亚洲5aaaaa淫片| 免费观看的影片在线观看| 看黄色毛片网站| 一边亲一边摸免费视频| 日韩一区二区三区影片| 欧美不卡视频在线免费观看| 长腿黑丝高跟| 成年免费大片在线观看| 亚洲av中文av极速乱| 男女啪啪激烈高潮av片| av女优亚洲男人天堂| 尾随美女入室| eeuss影院久久| 久久久国产成人免费| 中文资源天堂在线| 激情 狠狠 欧美| 国内精品宾馆在线| 亚洲最大成人中文| 国产一区亚洲一区在线观看| 一级av片app| 午夜福利网站1000一区二区三区| 欧美色视频一区免费| 少妇高潮的动态图| 99国产精品一区二区蜜桃av| 国产麻豆成人av免费视频| 欧美激情国产日韩精品一区| 免费观看人在逋| 精品国内亚洲2022精品成人| 少妇的逼好多水| 超碰97精品在线观看| 人妻夜夜爽99麻豆av|