首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Transport of larvae by ocean currents is an important dispersal mechanism for many species. The timing and location of spawning can have a large influence on settlement location. Shifts in the known spawning habitat of fish, whether due to climate or the discovery of new spawning stock, can influence the distribution of juveniles and our understanding of connectivity. The globally distributed species; Pomatomus saltatrix, is one such example where a previously unrecognised summer spawning event and a more southern latitudinal extent was recently reported for the southwest Pacific population. Although restrictions are in place to protect the traditional spawning event, the importance of the newly recognised summer spawning event is uncertain. Here, we investigate larval dispersal of P. saltatrix using particle tracking simulations to identify the contributions of the different spawning events to settlement. By modelling dispersal of larvae released in northern and mid‐latitude regions over the Austral spring and summer, we show that the newly recognised mid‐latitude summer spawning event contributes over 50% of the larvae reaching southern latitudes. This is due to a reduced (1–2 days) pelagic larval duration (associated with temperature), resulting in reduced larval mortality, and the seasonal (summer) strengthening of the East Australian Current (EAC) transporting particles ~50 km further south. These findings demonstrate that in dynamic boundary current systems such as the EAC, the final settlement location of larvae that are transported by ocean currents can vary considerably depending on the timing and location of spawning and that multiple spawning events are important for maximum dispersal.  相似文献   

2.
The Eastern Canadian northern shrimp population, representing one of the most important fisheries in the region, decreased dramatically since the mid‐2000s to a historical low in 2017, but changes were not spatially uniform. Applying a biophysical model within Newfoundland and Labrador (NL) management areas, we investigated connectivity processes during the long pelagic larval phase (2–3 months) of Pandalus borealis and key drivers of larval dispersal in different environmental conditions. We selected 3 years representative of contrasting North Atlantic Oscillation (NAO) phases to assess potential larval dispersal patterns of the northern shrimp population in NL, and hierarchically assessed the impact of the timing of release (yearly and daily), release location, and vertical migration behaviour on shrimp larval dispersal. Overall, we found that populations located on the northern Newfoundland and Labrador shelf supplied potential settlers to southern populations because of the dominant Labrador Current. Ocean circulation and current velocities during the NAO positive year differed from other years, generating contrasting settlement spatial patterns. Larval release location and vertical migration behaviour were the two most important influences on the strength of larval supply and settlement patterns. Inclusion of diel and ontogenic swimming behaviour increased settlement success of larvae released from inshore areas, regardless of study years. Our study improves understanding of northern shrimp stock‐recruitment relationships, their sensitivity to changing environmental conditions, and spatially non‐homogeneous population decline for bentho‐pelagic species with a long larval phase, which could potentially help improve management strategies.  相似文献   

3.
Many demersal marine fish species depend on a dispersive larval stage that connects geographically discrete sub‐populations. Understanding connectivity between these sub‐populations is necessary to determine stock structure, which identifies the appropriate spatial scale for fishery management. Such connectivity is poorly understood for King George whiting (Sillaginodes punctatus; Perciformes) in South Australia's gulf system, even though spawning grounds and nursery areas are adequately defined. In response to declines in commercial catches and estimated biomass, this study aimed to determine the most important spawning grounds and nursery areas to recruitment, and the connectivity between them. A biophysical model was seeded with particles according to the distribution and density of eggs throughout the spawning area in 2017 and 2018. Despite inter‐annual differences in the origins of particles, dispersal pathways and predicted settlement areas remained consistent between years. Predicted settlement was generally highest to nursery areas only short distances from regional spawning grounds, consistent with previous hydrodynamic models. However, the model also predicted that spawning in one region could contribute to recruitment in an adjacent region later in the spawning season, which aligned with the breakdown of thermohaline fronts at the entrance of each gulf. The connectivity between spawning grounds and nursery areas predicted by the model is supported by spatio‐temporal patterns in the otolith chemistry of pre‐flexion larvae and settled juveniles. Consequently, the most parsimonious explanation is that the populations of King George whiting in South Australia's gulf system constitute a single, panmictic stock, which has implications for fishery management.  相似文献   

4.
A major challenge in marine ecology is describing patterns of larval dispersal and population connectivity, as well as their underlying processes. We used a biophysical model to simulate dispersal of eastern oyster, Crassostrea virginica, larvae and connectivity among a network of 10 no‐take reserves in a shallow, wind‐driven estuary to assess the relative importance of spawning location, spawning date, larval behavior, larval mortality, and adult reproductive output to predicted dispersal and connectivity patterns. The location (i.e., natal reserve) and date of spawning relative to physical processes, particularly frequency of wind reversals, were the dominant drivers of dispersal and connectivity patterns. To a lesser extent, larval behavior (i.e., 3D vertical advection and ontogenetic depth regulation) and mortality modified dispersal and connectivity, whereas spatiotemporal variability in adult reproductive output was of minimal importance. Over a 21‐day larval duration, mean dispersal distance of passive surface particles ranged from 5 to 40 km. Reserves were too small (1 km2) relative to mean dispersal distances to promote extensive local retention (median 2%) and spaced too far apart (typically ~50 km) to promote extensive inter‐reserve connectivity (median 2%). Limited connectivity and local retention may preclude the network from being self‐sustainable, thereby limiting its long‐term conservation and management benefits. In reserve systems characterized by limited connectivity, management efforts should focus on increasing connectivity by increasing the number or size of reserves to realize the benefits of improved adult demographics within reserves.  相似文献   

5.
The harvest of bay scallops (Argopecten irradians) from Buzzards Bay, Massachusetts, U.S.A. undergoes large interannual fluctuations, varying by more than an order of magnitude in successive years. To investigate the extent to which these fluctuations may be due to yearly variations in the transport of scallop larvae from spawning areas to suitable juvenile habitat (settlement zones), a high‐resolution hydrodynamic model was used to drive an individual‐based model of scallop larval transport. Model results revealed that scallop spawning in Buzzards Bay occurs during a time when nearshore bay currents were principally directed up‐bay in response to a persistent southwesterly sea breeze. This nearshore flow results in the substantial transport of larvae from lower‐bay spawning areas to settlement zones further up‐bay. Averaged over the entire bay, the spawning‐to‐settlement zone connectivity exhibits little interannual variation. However, connectivities between individual spawning and settlement zones vary by up to an order of magnitude. The model results identified spawning areas that have the greatest probability of transporting larvae to juvenile habitat. Because managers may aim to increase scallop populations either locally or broadly, the high‐connectivity spawning areas were divided into: (i) high larval retention and relatively little larval transport to adjoining settlement areas, (ii) both significant larval retention and transport to more distant settlement areas, and (iii) little larval retention but significant transport to distant settlement areas.  相似文献   

6.
European sardine (Sardina pilchardus) and round sardinella (Sardinella aurita) comprise two‐thirds of total landings of small pelagic fishes in the Canary Current Eastern Boundary Ecosystem (CCEBE). Their spawning habitat is the continental shelf where upwelling is responsible for high productivity. While upwelling intensity is predicted to change through ocean warming, the effects of upwelling intensity on larval fish habitat expansion is not well understood. Larval habitat characteristics of both species were investigated during different upwelling intensity regimes. Three surveys were carried out to sample fish larvae during cold (permanent upwelling) and warm (low upwelling) seasons along the southern coastal upwelling area of the CCEBE (13°–22.5°N). Sardina pilchardus larvae were observed in areas of strong upwelling during both seasons. Larval habitat expansion was restricted from 22.5°N to 17.5°N during cold seasons and to 22.5°N during the warm season. Sardinella aurita larvae were observed from 13°N to 15°N during cold seasons and 16–21°N in the warm season under low upwelling conditions. Generalized additive models predicted upwelling intensity driven larval fish abundance patterns. Observations and modeling revealed species‐specific spawning times and locations, that resulted in a niche partitioning allowing species' co‐existence. Alterations in upwelling intensity may have drastic effects on the spawning behavior, larval survival, and probably recruitment success of a species. The results enable insights into the spawning behavior of major small pelagic fish species in the CCEBE. Understanding biological responses to physical variability are essential in managing marine resources under changing climate conditions.  相似文献   

7.
Connectivity between spawning and nursery areas plays a major role in determining the spatial structure of fish populations and the boundaries of stock units. Here, the potential effects of surface current on a red mullet population in the Central Mediterranean were simulated using a physical oceanographic model. Red mullet larvae were represented as Lagrangian drifters released in known spawning areas of the Strait of Sicily (SoS), which represents one of the most productive demersal fishing‐grounds of the Mediterranean. To consider the effect of inter‐annual variability of oceanographic patterns, numerical simulations were performed for the spawning seasons from 1999 to 2012. The main goal was to explore connectivity between population subunits, in terms of spawning and nursery areas, inhabiting the northern (Sicilian‐Maltese) and southern (African) continental shelves of the SoS. The numerical simulations revealed a certain degree of connectivity between the Sicilian–Maltese and the African sides of the SoS. Connectivity is present in both directions, but it is stronger from the Sicilian–Maltese spawning areas to the African nurseries owing to the marine circulation features of the region. However, because the majority of the larvae are transported to areas unsuitable for settlement or outside the SoS, the dispersal process is characterized by a strong loss of potential settlers born in the spawning areas. These results are in agreement with the low genetic heterogeneity reported for this species in the Mediterranean Sea and support the existence of a metapopulation structure of red mullet in the SoS and the adjacent areas.  相似文献   

8.
A simple microcomputer program uses Lagrangian particle tracking to simulate the fates of individual larvae subject to wind- and tide-driven advection fields and diffusion. Space/time advection patterns are provided to the program as inputs; thus, the program can use results from various hydrodynamic models. The program can simultaneously track several groups of larvae that are assigned varying attributes, including (1) body size, to allow simulation of effects of variation in growth rates; (2) spawning locations and times; (3) vertical migration behavior; and (4) settlement strategy (duration of larval period, depth conditions for settlement). The model is intended as an exploratory tool to help identify alternative hypotheses that might explain observed life history patterns and causes of inter-annual variation in recruitment rate. For English sole in the Hecate Strait, B.C., the model leads us to hypothesize that observed spawning locations have resulted from a tradeoff between places that would provide the best feeding opportunities and places that would minimize risk of advective export to unfavorable rearing habitats.  相似文献   

9.
The pelagic larval duration (PLD) period of fish can influence dispersal, recruitment, and population connectivity, thereby potentially informing best strategies for fisheries management. Computer models were used to simulate the dispersal of larvae of three species, representing a range of PLDs, from the Pacific island of Guam and neighboring islands for a 9‐year period (2004–2012) to gain insight into the best management strategies for these species. The species included two springtime spawners with shorter and longer PLDs, scribbled rabbitfish (Siganus spinus; 33 days) and yellowfin goatfish (Mulloidichthys flavolineatus; ~90 days maximum), and a fall spawner with a similarly long PLD, bluespine unicornfish (Naso unicornis; ~94 days maximum). An ocean circulation model coupled with a particle dispersal model provided simulated numbers of larvae settling at each island in relation to the island where they were spawned. Graph analysis was used to examine generational connections between islands. For S. spinus, self‐seeding was the dominant means of replenishment at Guam. Local management actions to maintain adequate spawning stock may be a primary control on long‐term sustainability for that fishery. In contrast, N. unicornis and M. flavolineatus populations at Guam were reliant on outside sources for 92%–98% of larval supply. For them, identifying and negotiating the preservation of upstream spawning potential in the Marshall Islands and Federated States of Micronesia will be needed. Guam played a relatively minor role in generational connectivity across the region. Shortest paths spanning the region often did not pass through Guam, or there were equally short paths through other islands.  相似文献   

10.
The Strait of Georgia (SoG), between Vancouver Island and mainland British Columbia, is a larval rearing ground for both hake and herring stocks, which are commercially important. Year‐to‐year variability in larval retention within the strait is examined by simulating drift tracks of larvae for these species using an ocean circulation model and a particle‐tracking model. Larvae with different vertical swimming behaviors were tracked in the springs of 2007, 2008, and 2009. Since herring larvae mostly stay near the surface, their distribution is heavily influenced by the wind. Strong winds to the north soon after the hatching period tend to wash herring larvae out of SoG and winds to the south help retain herring larvae inside the Strait. In 2007, the model indicates a massive wind‐driven export of herring larvae which may have led to the observed failure of herring production. In contrast, hake larvae reside deeper in the water column (50–200 m). Their distribution is less sensitive to surface forcing but is shaped by a deep gyre with cross‐strait currents. This study also suggests that the northern and southern SoG are weakly connected for herring larvae dispersal, which makes both regions potentially important to recruitment.  相似文献   

11.
We assessed by numerical modeling the coastal fish larval dispersion along the southern coast of Mallorca (Balearic Islands, NW Mediterranean) with the objective of determining the factors that contribute to successful recruitment. We assumed that fish larvae dispersal is mainly regulated by physical transport. Currents are mainly wind driven in this area; therefore, changes in wind forcing have a first‐order impact on larval transport. The synoptic wind patterns were systematically analyzed based on self‐organizing map analysis. The wind fields were clustered using a neural network pattern recognition approach into two modes, producing opposite along‐shelf flow. The seasonal changes between spring and summer in the dominance of either mode modulate the along‐shelf circulation, producing flow shifts under some circumstances. This variability in the wind regime was consistent throughout the 10 years analyzed (2000–2009). Using the Princeton Ocean Model (POM) and a particle‐tracking algorithm, we analyzed the effect of wind‐forced currents in the connectivity among near‐shore habitats. We show that, at the spatial scale considered, the coastal morphology and stochastic wind forcing favor local recruitment (mean of 30% self‐recruitment). Maximum transport distances of 20–30 km were typically associated with particles left to drift for 21 days. The implications for the performance of the four marine protected areas near SW Mallorca Island are discussed. Our results suggest that, although wind episodes determine the fate of short‐time spawning events, on a seasonal basis, regular larval supply to coastal zones is ensured by wind stochasticity.  相似文献   

12.
The poleward flowing East Australian Current (EAC) drives sporadic upwelling, entrains coastal water and forms the western Tasman Front (wTF), creating a mosaic of water types and larval transport routes along south eastern Australia. The spatial distribution, otolith chemistry and growth rates of larval sardine (Sardinops sagax) were examined to infer spawning location and larval transport. A gradient of increasing larval size from north to south along the shelf was not detected but was evident between the shelf and offshore in the wTF. Here larvae were larger and older. Based on the occurrence of newly hatched larvae, spawning by S. sagax between southern Queensland and mid New South Wales (NSW) was more extensive than previously reported. The otolith chemistry from two wTF larval size classes differed, implying different origins. The otolith chemistry of wTF post‐flexion larvae was similar to larvae from northern NSW, whereas wTF flexion larvae were similar to larvae observed nearby from mid‐NSW. Two possible larval transport routes, direct and indirect, are inferred from otolith chemistry, current velocities and a previously published particle tracking study. Either larvae from northern NSW were advected south and entrained with younger larvae directly into the wTF, or larvae from a range of shelf regions were advected around the southern edge of an anticyclonic eddy, to join younger larvae directly entrained into the wTF. Based on the co‐occurrence of larval ages and sizes in the wTF and their advection routes, the wTF appears to be an important larval retention zone.  相似文献   

13.
To study the transport of plaice (Pleuronectes platessa L.) eggs and larvae in the eastern Irish Sea, we constructed a 3D‐baroclinic physical model and coupled it to a particle‐tracking scheme that allowed aspects of larval behaviour to be simulated. Starting positions for eggs were based upon data from a series of ichthyoplankton surveys and final positions were compared with results of settled plaice distributions from two beam trawl surveys conducted on beaches around the eastern Irish Sea. If simulated larval behaviour was limited to passive drift or horizontal swimming, the particles diffused away from the spawning areas but failed to reach nursery grounds in significant numbers (85–90% remaining offshore). In contrast, switching on circatidal vertical swimming significantly increased the numbers of larvae reaching the coast (only 23–30% remained offshore). Particles tended to accumulate in bays and estuaries and this pattern compared well with the distribution of settled plaice from the field surveys. Studies in the southern North Sea (where spawning and nursery grounds are widely separated) have also demonstrated the importance of selective tidal stream transport for successful recruitment of settling plaice to nursery grounds. Although our understanding of the ontogeny of this behaviour is still poor, the model results presented suggest that this aspect of behaviour is a key factor influencing plaice settlement success.  相似文献   

14.
Determining the extent of dispersal in exploited marine fishes is essential for understanding their population dynamics and optimising management. The West Australian dhufish, Glaucosoma hebraicum, is a highly sought‐after, large and long‐lived reef‐dwelling species, endemic to south‐western Australia. Stock assessments indicate that this indicator species is overexploited. The fishery is managed using a zone‐based system, which implicitly assumes a high degree of demographic independence among zones. While tagging studies indicate limited movement of adult G. hebraicum, there is no understanding of the spatial scale of dispersal of its larvae and thus the true extent of interdependence of management zones. We analysed 13 microsatellite DNA loci to characterise the extent of gene flow, and conducted particle tracking simulations to model larval transport in this species. Genetic data demonstrated that some local recruitment was likely, but that on a broad scale gene flow between the management zones was extensive, and the entire fishery represents a single genetic stock. Hydrodynamic modelling predicted that the majority of dhufish larvae recruit from within the management zone where they are spawned, and that inter‐annual variation in current velocities has limited effect on the extent of larval transport. Because management zones are likely to be largely independent in terms of both larval and adult recruitment, heavy localised fishing pressure has the potential to reduce the abundance and reproductive capacity of this species in highly fished areas, but it should have limited impact on neutral genetic diversity.  相似文献   

15.
The connectivity and three‐dimensional (3D) dispersion of the larvae of giant red shrimp, potentially released from known spawning areas along the Sardinia slope in the western Mediterranean Sea, were assessed using Lagrangian simulations forced by a 3D submesoscale permitting a regional ocean model. Biophysical simulations using the hydrodynamic conditions of the year 2012 were run to track propagules released from known spawning areas during the spawning period (May to September). Passive transport (PT) and vertical migration (VM) scenarios were tested, each with two possible pelagic larval durations (PLDs) of 21 or 42 days. Dispersion of propagules in the PT and VM scenarios differed in terms of travelled distance, export out of the domain (larger for VM), and depth distribution (shallower and bimodal for VM due to the larger variability of encountered currents). Connectivity patterns were investigated among eight release areas, and four predetermined Eco‐Regions, and results showed strong connectivity among the North‐Western (NW), Western (W), and Southern (S) regions of Sardinia, whereas the Eastern region was more segregated. Differences in connectivity patterns among scenarios were related mainly to the tendency of greater retention of propagules in the release area for the PT scenarios. This finding, together with existing hypotheses of vertical migration likely occurring during first egg‐larval phases, suggest that the VM scenarios are the most probable of the two hypotheses tested. Strong connectivity between the W and S sides of Sardinia and the relative isolation of the E side could have significant implications for the protection of this important resource.  相似文献   

16.
For many marine fish species, recruitment is strongly related to larval survival and dispersal to nursery areas. Simulating larval drift should help assessing the sensitivity of recruitment variability to early life history. An individual‐based model (IBM) coupled to a hydrodynamic model was used to simulate common sole larval supply from spawning areas to coastal and estuarine nursery grounds at the population scale in the eastern Channel on a 14‐yr time series, from 1991 to 2004. The IBM allowed each particle released to be transported by currents from the hydrodynamic model, to grow with temperature, to migrate vertically giving stage development, and possibly to die according to drift duration, representing the life history from spawning to metamorphosis. Despite sensitivity to the larval mortality rate, the model provided realistic simulations of cohort decline and spatio‐temporal variability of larval supply. The model outputs were analysed to explore the effects of hydrodynamics and life history on the interannual variability of settled sole larvae in coastal nurseries. Different hypotheses of the spawning spatial distribution were also tested, comparing homogeneous egg distribution to observation and potential larval survival (PLS) maps. The sensitivity analyses demonstrated that larval supply is more sensitive to the life history along larval drift than to the phenology and volume of spawning, providing explanations for the lack of significant stock–recruitment relationship. Nevertheless, larval supply is sensitive to spawning distribution. Results also suggested a very low connectivity between supposed different sub‐populations in the eastern Channel.  相似文献   

17.
Winter flounder, Pseudopleuronectes americanus, from southern New England and the Gulf of Maine were historically considered to be obligate estuarine spawners. However, recent experiments and observations document that winter flounder in the Gulf of Maine also utilize coastal waters for spawning. An individual‐based modeling approach was used to investigate the transport of winter flounder larvae from three hypothesized coastal spawning grounds in the Gulf of Maine. Transport success rates were greatest for larvae released from Ipswich Bay, intermediate for Stellwagen Bank and least successful for those released from Jeffreys Ledge. There was substantial interannual variability in larval transport and geographic patterns of potential connectivity. Furthermore, the date of spawning had an important influence on transport success. Model results suggest that certain coastal spawning grounds used by winter flounder may serve as an important source of larvae to estuaries and nearshore nursery areas. The potential influx of coastal spawned larvae could have implications for the resilience, productivity and gene flow in local populations. Model results provide further support for the conclusion that winter flounder in the Gulf of Maine may not be solely dependent upon estuaries for spawning. Results also suggest that coastal spawning groups should be considered explicitly in the management of winter flounder, and protected under Essential Fish Habitat regulations.  相似文献   

18.
A bio‐physical semi‐Langrangian model was developed to follow the drift and abundance of lobster larvae from hatching to settlement as post‐larvae. The geographic domain encompasses Northumberland Strait and the areas surrounding Prince Edward Island in the southern Gulf of St. Lawrence. The model was run for ten larval seasons, 1983–87 and 1997–2001, representing two periods of egg production. The model was forced with tides, winds, heat fluxes, and marine and freshwater fluxes at open boundaries. Biological inputs were location, date, and density of larvae at hatching, development rates, a time window and a minimum bottom temperature required for successful settlement, and two scenarios of daily mortality. Net drift was west to east but stronger on the north than the south side of Prince Edward Island. The hatch was greater in the second 5‐yr period, but the spatial patterns of settlement were similar. For both mortality scenarios, the same five of the 24 larval source areas were important in providing settlers. Horizontal shear of larval distribution indicates fishing communities were dependent on the hatch from fishing grounds of multiple upstream communities. Variation in daily and annual post‐larval settlement was greatest north of Prince Edward Island. From the 24 source areas there was an eightfold range in the fraction of larvae surviving to settlement, with advection into deep water an important cause of mortality. Four to five possible fishery management areas were identified using multidimensional scaling to group sink areas having shared source areas. We hypothesize leaky discontinuity in connectivity between these groups.  相似文献   

19.
In the marine environment, pelagic dispersal is important for determining the distribution and abundance of populations, as well as providing connections among populations. Estimates of larval dispersal from spawning grounds are important to determining temporal and spatial patterns in recruitment that may have significant influences on the dynamics of the population. We present a case study of the dispersal of Centropristis striata (black sea bass) larvae on the southeast U.S. continental shelf. We use a coupled larval behavior – 3D circulation model to compare the effects of the timing and location of spawning against that of larval vertical migration on larval dispersal. Using the results of field data on larval vertical distributions, we compare the dispersal of virtual 'larvae' which have ontogenetic changes in vertical behavior with that of particles fixed near the surface and near the bottom. Larvae were released at potential spawning sites four times throughout the spawning season (February through May) for 3 yr (2002–04) and tracked for the assumed larval duration (from 27 to 37 days including the egg stage). Results indicate that adult behavior, in the form of spawning time and location, may be more important than larval vertical behavior in determining larval dispersal on the inner- and mid- continental shelves of this region.  相似文献   

20.
The reproductive biology of Centroberyx gerrardi (Günther) was investigated across ~2,000 km of its southern Australian distribution, encompassing different jurisdictions and varying environmental features. Greater gonad mass and prevalence of spawning fish, along with lower ratios of lengths at maturity:maximum lengths and ages at maturity:maximum ages, were identified at the western‐most (Capes) and eastern‐most (Great Australian Bight; GAB) regions. Across the study region, spawning peaks in summer/autumn, when water temperatures are warmest. Regional differences in potential “reproductive output,” while not consistent with the eastward decline in mean monthly water temperature, may instead be related to summer upwelling in the Capes and GAB, driving greater oceanic productivity prior to peak spawning, supporting larval survival. In autumn, the prevailing southward and eastward flowing, downwelling Leeuwin Current (LC) strengthens, providing a dispersal mechanism along the west and south Australian coasts, but limiting upwelling effects. Predicted changes in environmental conditions and their potential impacts on C. gerrardi are discussed, in particular how these factors may affect recruitment to stocks and fisheries, requiring a better understanding of source‐sink relationships for this species. As environmental changes occur, management strategies to sustain fish resources must adapt to spatially variable and changing reproductive output and be collaborative across jurisdictions.  相似文献   

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

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