首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Estimating the dynamics of furbearer populations is challenging because their elusive behavior and low densities make observations difficult. Statistical population reconstruction is a flexible approach to demographic assessment for harvested populations, but the technique has not been applied to furbearers. We extended this approach to furbearers and analyzed 8 yr of age-at-harvest data for American marten (Martes americana) in the Upper Peninsula of Michigan. Marten abundance estimates showed a general downward trend from an estimate of = 1,733.3 animals in 2000 to = 1,163.9 in 2007. The harvest probability of martens increased nearly 5-fold from 0.0542 in 2000 to 0.2637 in 2007, which corresponded to a 5-fold increase in trap-nights. Continued monitoring of martens in the Upper Peninsula, Michigan, and a reassessment of current harvest regulations are necessary given the estimated decreases. Moreover, we do not encourage the use of harvest indices as the sole technique to assess the status and trends of marten and fisher populations. Auxiliary studies in the Upper Peninsula, Michigan, will allow for continued use and improvement in the application of these models. © 2011 The Wildlife Society.  相似文献   

2.
Assessing wildlife management action requires monitoring populations, and abundance often is the parameter monitored. Recent methodological advances have enabled estimation of mean abundance within a habitat using presence–absence or count data obtained via repeated visits to a sample of sites. These methods assume populations are closed and intuitively assume habitats within sites change little during a field season. However, many habitats are highly variable over short periods. We developed a variation of existing occupancy and abundance models that allows for extreme spatio-temporal differences in habitat, and resulting changes in wildlife abundance, among sites and among visits to a site within a field season. We conducted our study in sugarcane habitat within the Everglades Agricultural Area southeast of Lake Okeechobee in south Florida. We counted wintering birds, primarily passerines, within 245 sites usually 5 times at each site during December 2006–March 2007. We estimated occupancy and mean abundance of birds in 6 vegetation states during the sugarcane harvest and allowed these parameters to vary temporally or spatially within a vegetation state. Occupancy and mean abundance of the common yellowthroat (Geothlypis trichas) was affected by structure of sugarcane and uncultivated edge vegetation (occupancy = 1.00 [ = 0.96–1.00] and mean abundance = 7.9 [ = 3.2–19.5] in tall sugarcane with tall edge vegetation versus 0.20 [ = 0.04–0.71] and 0.22 [ = 0.04–1.2], respectively, in short sugarcane with short edge vegetation in one half of the study area). Occupancy and mean abundance of palm warblers (Dendroica palmarum) were constant (occupancy = 1.00, = 0.69–1.00; mean abundance = 18, = 1–270). Our model may enable wildlife managers to assess rigorously effects of future edge habitat management on avian distribution and abundance within agricultural landscapes during winter or the breeding season. The model may also help wildlife managers make similar management decisions involving other dynamic habitats such as wetlands, prairies, and even forested areas if forest management or fires occur during the field season. © 2011 The Wildlife Society.  相似文献   

3.
Environmental factors, such as forest characteristics, have been linked to fawn survival in eastern and southern white-tailed deer (Odocoileus virginianus) populations. In the Great Plains, less is known about how intrinsic and habitat factors influence fawn survival. During 2007–2009, we captured and radiocollared 81 fawns in north-central South Dakota and recorded 23 mortalities, of which 18 died before 1 September. Predation accounted for 52.2% of mortality; remaining mortality included human (hunting, vehicle, and farm accident; 26.1%) and hypothermia (21.7%). Coyotes (Canis latrans) accounted for 83.3% of predation on fawns. We used known-fate analysis in Program MARK to estimate summer (15 May–31 Aug) survival rates and investigated the influence of intrinsic and habitat variables on survival. We developed 2 a priori model sets, including intrinsic variables and a test of annual variation in survival (model set 1) and habitat variables (model set 2). Model set 1 indicated that summer survival varied among years (2007–2009); annual survival rates were 0.94 (SE = 0.06, n = 22), 0.78 (SE = 0.09, n = 27), and 0.54 (SE = 0.10, n = 32), respectively. Model set 2 indicated that survival was further influenced by patch density of cover habitats (Conservation Reserve Program [CRP]-grasslands, forested cover, and wetlands). Mean CRP-grassland and wetland patch density (no. patches/100 ha) were greater (P < 0.001) in home-range areas of surviving fawns ( = 1.81, SE = 0.10, n = 63; = 1.75, SE = 0.14, n = 63, respectively) than in home-range areas of fawns that died ( = 0.16, SE = 0.04, n = 18; = 1.28, SE = 0.10, n = 18, respectively). Mean forested cover patch density was less (P < 0.001) in home-range areas of surviving fawns ( = 0.77, SE = 0.10, n = 63) than in home-range areas of fawns that died ( = 1.49, SE = 0.21, n = 18). Our results indicate that management activities should focus on CRP-grassland and wetland habitats in order to maintain or improve fawn survival in the northern Great Plains, rather than forested cover composed primarily of tree plantings and shelterbelts. © 2012 The Wildlife Society.  相似文献   

4.
Duckling survival is an important component of mallard (Anas platyrhynchos) recruitment and population growth, yet many factors regulating duckling survival are poorly understood. We investigated factors affecting mallard duckling survival in the drift prairie of northeastern North Dakota, 2006–2007. Mammalian meso-predators were removed by trapping on 4 92.3 km2 study sites and another 4 study sites served as controls. We monitored 169 broods using telemetry and periodic resighting, and we modeled cumulative survival to 30 days of age using the nest survival module in Program MARK. Duckling survival was not affected by predator removal ( , 85% CI: 0.182–0.234; , 85% CI: 0.155–0.211) and was only weakly negatively correlated with duckling density. Duckling survival was higher in 2007 ( , 85% CI: 0.193–0.355) than 2006 ( , 85% CI: 0.084–0.252) and increased with total seasonal and semipermanent wetland area and declined with perennial cover in the surrounding landscape. Broods that hatched earlier in the season (especially in 2006) and ducklings that were heavier at hatch also had higher survival. Our estimates of duckling survival are among the lowest reported for mallards and contradict previous research in Saskatchewan that found predator removal increased duckling survival. However, our results are consistent with other studies suggesting that earlier hatch date, increased wetland availability, and better duckling condition lead to increased survival. Management actions that increase wetland density, improve nest success early in the season, and potentially target brood-specific predators such as mink (Neovison vison) would likely lead to higher duckling survival. © 2011 The Wildlife Society.  相似文献   

5.
In this article, we consider r observations from a non‐homogeneous censored Markov chain, with transition probability matrix P. For the product estimator of P proposed by Aalen and Johansen (1978) and Phelan (1988), we investigate the behavior of Bayesian bootstrap clones to approximate the sampling distribution of , and then construct approximate confidence interval. It is shown that the approximation based on the random‐weighted distribution is first‐order consistent. The performance of the Bayesian bootstrap clones (BBC) is also discussed by small sample simulation. Finally, we illustrate the BBC procedure in the application to the WHO malaria survey data (cf. Singer and Cohen 1970).  相似文献   

6.
Circulating neutrophils isolated from patients 3–4 h after a myocardial infarction produced less $ {\rm O}\frac{ \cdot }{{\rm 2}} $ compared with controls, when stimulated with phorbol myrystate acetate or formyl-methionine-leucine-phenylalanine. Three days after the infraction the $ {\rm O}\frac{ \cdot }{{\rm 2}} $ generation elicited by both stimuli further decreased markedly. Seven and 15 days after infarction the $ {\rm O}\frac{ \cdot }{{\rm 2}} $ stimulated production was only slightly lower than or similar to the control values. The neutrophils of infarcted patients showed an augmented latency period before $ {\rm O}\frac{ \cdot }{{\rm 2}} $ production compared with controls in response to exogenous stimuli, particularly three days after infarction. Electron microscopy revealed that the neutrophils isolated from the infarcted patients displayed signs of cell exhaustion with few alterations of the plasma membranes when stimulated with phorbol ester. In contrast, control neutrophils displayed alterations of the plasma membranes characteristic of active neutrophils. The results of this study indicate that the circulating neutrophils appear exhausted and functionally inhibited immediately after myocardial infarction.  相似文献   

7.
In this paper it is shown that if N= \documentclass{article}\pagestyle{empty}\begin{document}$ \mathop \sum \limits_{i = 1}^{S_h} $\end{document} cihNih, where cih are some non-negative integer numbers and Nih are such incidence matrices that Ah = \documentclass{article}\pagestyle{empty}\begin{document}$ \mathop \sum \limits_{i = 1}^{S_h} $\end{document} i Nih is a balanced matrix defined by SHAH (1959), for h = 1, 2,…, p, then a block design with an incidence matrix Ñ = [N, N,…,N] is an equi-replicated balanced block design. Here the balance of a block design is defined in terms of the matrix M0 introduced by CALI?SKI (1971).  相似文献   

8.
Let x1x2x3 … ≤xr be the r smallest observations out of n observations from a location-scale family with density $ \frac{1}{\sigma}f\left({\frac{{x - \mu}}{\sigma}} \right) $ where μ and σ are the location and the scale parameters respectively. The goal is to construct a prediction interval of the form $ \left({\hat \mu + k_1 \hat \sigma,\,\hat \mu + k_2 \hat \sigma} \right) $ for a location-scale invariant function, T(Y) = T(Y1, …, Ym), of m future observations from the same distribution. Given any invariant estimators $ \hat \mu $ and $ \hat \sigma $, we have developed a general procedure for how to compute the values of k1 and k2. The two attractive features of the procedure are that it does not require any distributional knowledge of the joint distribution of the estimators beyond their first two raw moments and $ \hat \mu $ and $ \hat \sigma $ can be any invariant estimators of μ and σ. Examples with real data have been given and extensive simulation study showing the performance of the procedure is also offered.  相似文献   

9.
G. Weill  J. Sturm 《Biopolymers》1975,14(12):2537-2553
Theoretical expressions are derived for the change in the polarized components of the fluorescence, resulting from the orientation of a rigid molecule bearing a chromophore with arbitrary angles for the absorption and transition moments \documentclass{article}\pagestyle{empty}\begin{document}$ \vec \mu _a $\end{document} and \documentclass{article}\pagestyle{empty}\begin{document}$ \vec \mu _e $\end{document} with respect to the molecular axis. The break in the symmetry relation HV = VH is related to the tilt angle between \documentclass{article}\pagestyle{empty}\begin{document}$ \vec \mu _a $\end{document} and \documentclass{article}\pagestyle{empty}\begin{document}$ \vec \mu _e $\end{document}. The theory is applied to a sonicated DNA–2-hydroxy-4,4′-diamidinostilbene complex, in the blue and red emission bands of this peculiar dye. Simultaneous measurements of linear dichroism and fluorescence lead to the determination of an angle of 47° between a fluorescent bound dye and the DNA axis, with no difference for the blue- and red-emitting species, but confirm the presence of nonfluorescent bound dye in a more perpendicular arrangement.  相似文献   

10.
The relationship between cytosolic concentrations of Ca2+ (Ca) and Na+ (Na) were studied in preparations of rat submandibular and pancreatic acini loaded with the Ca2+-sensitive dye Fura-2 or the Na+-sensitive dye SBFI. Pancreatic acini showed no changes in Na during either transient or persistent changes in Ca. Increases in Ca produced by exposure of submandibular gland acini to carbachol, a muscarinic cholinergic agonist, were followed by an increase in Na after a delay of 5–10 s. When Ca2+ stores were mobilized without Ca2+ influx Na also increased, but in acini loaded with BAPTA, a nonfluorescent Ca2+ chelator, the transient increase in Ca2+ caused by mobilization of stored Ca2+ was virtually abolished, as was the increase in Na. In the presence of ionomycin, increases in Ca were followed by increases in Na. Ca2+-dependent increases in Na were abolished in Na+-free buffer and by the presence of furosemide, a blocker of Na+-K+-2Cl cotransport. In other studies, extracellular ATP (ATPo) produced an increase in Ca and Na. The steady-state increase in Ca was reduced by increasing extracellular Na+ concentrations (Na) in dose-dependent fashion (IC50 = 16.4 ± 4.7 mM Na+). Likewise, increasing Na reduced ATPo-stimulated 45Ca2+ uptake at steady state (IC50 = 15.8 ± 9.2 mM Na+). Changing Na had no effect on carbachol-stimulated increases in Ca. We conclude that, in rat submandibular gland acini, ATPo promotes an increase in Ca and Na via a common influx pathway and that, under physiologic conditions, Na+ significantly limits the ATPo-stimulated increase in Ca. In the presence of carbachol, however, Na rises in Ca-dependent fashion in submandibular gland acini via stimulation of Na+-K+-2Cl cotransport. © 1996 Wiley-Liss, Inc.  相似文献   

11.
Box turtles are being extirpated from much of their former range and remaining populations often live in association with anthropogenically altered habitats. This is particularly evident at the northern distributional limit of eastern box turtles (Terrapene carolina carolina) and is an important factor to consider during the winter months when their ability to respond to microclimatic change is limited. Using temperature dataloggers, we studied the hibernal microclimate of box turtles and associated habitat following timber harvests. We monitored the body temperatures of 38 eastern box turtles and collected detailed air and soil profile temperatures of 12 box turtle hibernacula, 6 clearcuts, and 6 adjacent forested areas during the hibernal season (winter 2009–2010). We partitioned the hibernal season into 2 biologically significant thermal periods: hibernation and emergence. The mean hibernation body temperature averaged (3.28° C, SE = 0.09) and corresponded to an average depth of 10 cm. Clearcuts were consistently colder ( = 1.91° C) than forests ( = 2.68° C) and hibernacula ( = 2.77° C) during hibernation, but became the warmest areas during emergence ( = 9.96° C). We found that in the average clearcut, turtles could burrow to approximately 20 cm to attain the average hibernation body temperature or to approximately 15 cm to attain a body temperature no different than those overwintering on colder, northeast-facing slopes in the forest ( = 2.83° C). Alternatively, we found that southwest-facing slopes were warmer and if turtles chose to overwinter only in clearcuts on those slopes, they could remain shallower. All but 1 turtle overwintered in forested areas; however, our study suggests that some timber harvested areas offer various microhabitats exploitable by hibernating box turtles based on soil profile temperatures, slope aspect, and depth of hibernation. © 2012 The Wildlife Society.  相似文献   

12.
We describe conditions which lead to complete helix formation of poly(I) in the presence of NH. Binding of NH is shown to be specific in the presence of Li+, which does not by itself support helix formation under these conditions. The NH–poly(I) complex is characterized by uv, CD, and ir spectroscopy. The CD spectrum is strikingly different from those of the Na+ or K+ complexes, the first extremum being changed from negative for the metal ions to positive for NH. A stereospecific model is proposed for the NH–poly(I) helix in which the N of NH is located on the axis of the four-stranded helix, midway between planar tetramers formed by the bases. The model is consistent with the tetrahedral symmetry of NH, the requirement for four acceptable hydrogen bonds, the observed stability of the helix, and the accepted geometry of the backbone.  相似文献   

13.
The use of methane for production of bacterial protein   总被引:3,自引:0,他引:3  
Single cell protein production was studied in a mixed bacterial culture grown in methane using batch and continuous culture techniques. Overall productivity was found to be higher in the continuous culture which gave a maximum productivity value p = 0.15 g/l/h. Methane and oxygen were consumed in the relation 1 : 1.7. Yield coefficients for methane, oxygen, and ammonium chloride were Y = 0.90, Y = 0.26, and Y = 0.14. The crude protein content of the biomass was 71%.  相似文献   

14.
We studied the interactions of the substitution-inert inversion-labile complexes Fe(bipy) and Fe(phen) [and the inversion-stable complex Ru(bipy)] with DNA. The association of these complexes to DNA is mainly electrostatic, and Fe(phen) shows a more effective binding to DNA than the two bipyridyl complexes, possibly owing to a different binding mode. The interactions are enantioselective, leading to a Pfeiffer shift in the diastereomeric inversion equilibria and an excess of the Δ-enantiomer of Fe(phen) and Fe(bipy), which is directly monitorable through CD. The partition constants for the inversion equilibrium range from 1.3 to 2.0 for Fe(bipy) and Fe(phen), depending on ionic conditions. From flow LD information about the orientation of the complexes on DNA was obtained: it is consistent with a fit of the Δ-enantiomer in the major groove of the right-handed DNA helix. The mechanisms of interaction are discussed against equilibrium, spectroscopic, and kinetic data.  相似文献   

15.
The ethanol-inhibitory behaviour of the yeast Saccharomyces cerevisiae Sc 5 was found to be characterized by a continual-linear relation between the specific ethanol formation rate and the ethanol concentration. Therefore the simple equation could be applied for it. It is shown that this model is correct only then, if all of the process parameters are in their optimum. Out of the optimum temperature range the characteristics of the function ν = f(P) change in such a way that in regard to the ethanol concentration P twc linear relations exist for each suboptimum temperature: and a non-linear equation is current for each superoptimum temperature: where bT is also a function of the temperature and always less than 1. Taking as a basis these equations the specific ethanol formation rate of the used strain can be calculated for the whole biokinetic P/T-sphere of ethanol production.  相似文献   

16.
Models are important tools that can help managers and researchers understand the population dynamics of a species and how different habitat or population management scenarios impact that species. We used radio-telemetry data from northern bobwhites (Colinus virginianus) in southern Texas from 2000 to 2005 to develop a stochastic simulation model for bobwhite populations. Our model is based on difference equations, with stochastic variables drawn from normal and Weibull distributions. We simulated bobwhite populations to 100 yr and evaluated our model by comparing results with independent estimates of 4 population parameters (spring and fall density, finite rate of increase in the fall population [λ], and winter juv:ad age ratios). Using a quasi-extinction criterion of ≤40 birds (density = ≤0.05 birds/ha), probability of persistence to 100 yr was 88.3% (106 of 120 simulations) for the spring population and 96.7% (116 of 120 simulations) for the fall population. Using a less restrictive quasi-extinction criteria (≤14 birds), probability of persistence was 93.3% (112 of 120 simulations) for the spring population and 98.3% (118 of 120 simulations) for the fall population. Simulated population parameters were similar to independent estimates for 4 of 4 population parameters. Winter age ratios differed between our model ( juv:ad, n = 120, SE = 0.32) and empirical age ratios from harvested bobwhites on our study area ( juv:ad, n = 25, SE = 0.24). However, when we corrected harvest age ratios for bias in juvenile harvest ( juv:ad, n = 25, SE = 0.32) simulated and empirical estimates were similar. Our model appears to be a reliable predictor of bobwhite populations in the southern Texas. Our simulation results indicate that bobwhite hunters and managers can expect excellent bobwhite hunting (fall populations ≥2.2 birds per ha) in about one of 10 yr. © 2011 The Wildlife Society  相似文献   

17.
Spring harvest rates of male wild turkeys (Meleagris gallapavo) influence the number and proportion of adult males in the population and turkey population models have treated harvest as additive to other sources of mortality. Therefore, hunting regulations and their effect on spring harvest rates have direct implications for hunter satisfaction. We used tag recovery models to estimate survival rates, investigate spatial, temporal, and demographic variability in harvest rates, and assess how harvest rates may be related to management strategies and landscape characteristics. We banded 3,266 male wild turkeys throughout New York, Ohio, and Pennsylvania during 2006–2009. We found little evidence that harvest rates varied by year or management zone. The proportion of the landscape that was forested within 6.5 km of the capture location was negatively related to harvest rates; however, even though the proportion forested ranged from 0.008 to 0.96 across our study area, this corresponded to differences in harvest rates of only 2–5%. Annual survival was approximately twice as high for juveniles as adults . In turn, spring harvest rates for adult turkeys were greater for adults than juveniles . We estimated the population of male turkeys in New York and Pennsylvania ranged from 104,000 to 132,000 in all years and ranged from 63,000 to 75,000 in Ohio. Because of greater harvest rates for adult males, the proportion of adult males in the population was less than in the harvest and ranged from 0.40 to 0.81 among all states and years. The high harvest rates observed for adults may be offset by greater recruitment of juveniles into the adult age class the following year such that these states can sustain high harvest rates yet still maintain a relative high proportion of adult males in the harvest and population. © 2011 The Wildlife Society.  相似文献   

18.
Phorbol ester treatment of granulocytes triggers release of superoxide (O) and a concomitant burst of DNA strand breaks. The relationship between the amount of O and the number of DNA breaks has not previously been explored. To quantify the relatively large amount of O generated over a 40-min period by 1 × 106 granulocytes/mL, a discontinuous “10-min pulse” method employing cytochrome c was used; 140 nmol O per 1 × 106 cells was detected. DNA strand breaks were quantified by fluorimetric analysis of DNA unwinding (FADU). To vary the level of O released by cells, inhibitors of the respiratory burst were used. Sodium fluoride (1–10 mM) and staurosporine (2–10 nM) both inhibited O production. In both cases, however, inhibition of strand breakage was considerably more pronounced than inhibition of O. Zinc chloride (50–200 μM) inhibited both O and DNA breaks, approximately equally. Dinophysistoxin-1 (okadaic acid) inhibited O production more effectively than it inhibited DNA breaks. O dismutes to H2O2, a reactive oxygen species known to cause DNA breaks. The addition of catalase to remove extracellular H2O2 had no effect on DNA breakage. Using pulse field gel electrophoresis, few double-stranded breaks were detected compared to the number detected by FADU, indicating that about 95% of breaks were single-stranded. The level of DNA breaks is not directly related to the amount of extracellular O or H2O2 in PMA-stimulated granulocytes. We conclude that either an intracellular pool of these reactive oxygen species is involved in breakage or that the metabolic inhibitors are affecting a novel strand break pathway. J. Cell. Biochem. 66:219–228, 1997. © 1997 Wiley-Liss Inc.  相似文献   

19.
Four commonly used formulae for measuring percentage similarity (PS) of biological communities were tested for their usefulness in relating to two plankton community properties, species proportional differences and total density differences. The formula best combining species proportionality and total density in the expression of PS is new: where min (xi,yi) is the lesser percentage (doubly standardized) of a species in two samples X and Y and where 2 q, 2xi and 2yi are the total quantities of all species in samples 8,X and Y, where \documentclass{article}\pagestyle{empty}\begin{document}$ \sum\limits_i {z_i } ,\,\sum\limits_i {x_i } \,and\sum\limits_i {y_i } $\end{document} are the total quantities of all species in samples Z, X and Y, respectively. Sample 2 contains the highest density of all species in the set; \documentclass{article}\pagestyle{empty}\begin{document}$ \sum\limits_i {z_i \, > \,(\sum\limits_i {x_i ,\,} \sum\limits_i {y_i } )} $\end{document}. The new expression of PS is simple to use and has the additional advantage of offering the analyst an unlimited choice of weighting factors or importance values for proportionality of species content and total density. The method has been applied to data from Gravenhurst Bay (Ontario) and effectively demonstrates the consequences of phosphorus loading reductions for phytoplankton communities.  相似文献   

20.
Wetlands in the Upper Mississippi River and Great Lakes Region (UMRGLR) must annually sustain populations of migrating waterfowl from the mid-continent of North America. We used multi-stage sampling to estimate plant and invertebrate food biomasses (kg/ha) for ducks in 3 wetland habitat types at 6 stop-over locations in the UMRGLR during 2006 and 2007. Total biomass was greatest in palustrine emergent (PEM; = 208 kg/ha, SE = 23, median = 120), followed by palustrine forested (PF; = 87 kg/ha, SE = 7; median = 43), and lacustrine–riverine (LR; = 52 kg/ha, SE = 7; median = 27) wetlands. Ducks that foraged in forested and LR wetlands encountered the least food abundance during spring in the UMRGLR. Our estimates of food abundance were the lowest reported among other landscape scale surveys from mid-continent North America. About 1 in every 5 PEM wetlands and over half of our PF and LR wetlands that we sampled contained <50 kg/ha of food, suggesting many had little or no forage value to ducks during spring. Biomass of plant foods generally exceeded invertebrate biomass in all habitat types, although invertebrate biomass estimates exceeded plant biomass in 8 of 29 sites when considered by wetland type and year. Total food biomass estimates varied widely ( = 6–425 kg/ha) between years and among habitats; thus, using global arithmetic means to estimate food abundance for conservation planning obscures fine scale temporal and spatial variation that may be necessary for management on local and sub-regional levels. Distributions of food biomass estimates were right-skewed, causing us to question whether arithmetic means realistically represent levels of food abundance that all ducks encounter during spring migration. Alternative measures of central tendency (e.g., median) may be more biologically realistic, particularly if spring-migrating ducks are not distributed in an ideal-free manner with respect to food abundance. Future research should determine how ducks distribute themselves in relation to variation in food abundance in space and time during spring migration to strengthen the biological approach to conservation planning in non-breeding Joint Venture areas of the North American Waterfowl Management Plan. © 2011 The Wildlife Society.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号