共查询到20条相似文献,搜索用时 31 毫秒
1.
2.
The mean available potential energy released by baroclinic instability into the meso-scale eddy field has to be dissipated in some way and Tandon and Garrett [Tandon, A., Garrett, C., 1996. On a recent parameterization of mesoscale eddies. J. Phys. Oceanogr. 26 (3), 406–416] suggested that this dissipation could ultimately involve irreversible mixing of buoyancy by molecular processes at the small-scale end of the turbulence cascade. We revisit this idea and argue that the presence of dissipation within the thermocline automatically requires that a component of the eddy flux associated with meso-scale eddies must be associated with irreversible mixing of buoyancy within the thermocline. We offer a parameterisation of the implied diapycnal diffusivity based on (i) the dissipation rate for eddy kinetic energy given by the meso-scale eddy closure of Eden and Greatbatch [Eden, C., Greatbatch, R.J., 2008. Towards a meso-scale eddy closure. Ocean Modell. 20, 223–239.] and (ii) a fixed mixing efficiency. The implied eddy-induced diapycnal diffusivity (κ) is implemented in a coarse resolution model of the North Atlantic. In contrast to the vertical diffusivity given by a standard vertical mixing scheme, large lateral inhomogeneities can be found for κ in the interior of the ocean. In general, κ is large, i.e. up to o(10) cm2/s, near the western boundaries and almost vanishing in the interior of the ocean. 相似文献
3.
《Ocean Modelling》2004,6(1):83-100
A series of vertical mixing schemes implemented in a circumpolar coupled ice–ocean model of the BRIOS family is validated against observations of hydrography and sea ice coverage in the Weddell Sea. Assessed parameterizations include the Richardson number-dependent Pacanowski–Philander scheme, the Mellor–Yamada turbulent closure scheme, the K-profile parameterization, a bulk mixed layer model and the ocean penetrative plume scheme (OPPS). Combinations of the Pacanowski–Philander parameterization or the OPPS with a simple diagnostic model depending on the Monin–Obukhov length yield particularly good results. In contrast, experiments using a constant diffusivity and the traditional convective adjustment cannot reproduce the observations. An underestimation of wind-driven mixing in summer leads to an accumulation of salt in the winter water layer, inducing deep convection in the central Weddell Sea and a homogenization of the water column. Large upward heat fluxes in these simulations lead to the formation of unrealistic, large polynyas in the central Weddell Sea after only a few years of integration. Furthermore, spurious open-ocean convection affects the basin-scale circulation and leads to a significant overestimation of meridional overturning rates. We conclude that an adequate parameterization of both wind-induced mixing and buoyancy-driven convection is crucial for realistic simulations of processes in seasonally ice-covered seas. 相似文献
4.
An accurate numerical prediction of the oceanic upper layer velocity is a demanding requirement for many applications at sea and is a function of several near-surface processes that need to be incorporated in a numerical model. Among them, we assess the effects of vertical resolution, different vertical mixing parameterization (the so-called Generic Length Scale –GLS– set of k–ε, k–ω, gen, and the Mellor–Yamada), and surface roughness values on turbulent kinetic energy (k) injection from breaking waves.First, we modified the GLS turbulence closure formulation in the Regional Ocean Modeling System (ROMS) to incorporate the surface flux of turbulent kinetic energy due to wave breaking. Then, we applied the model to idealized test cases, exploring the sensitivity to the above mentioned factors. Last, the model was applied to a realistic situation in the Adriatic Sea driven by numerical meteorological forcings and river discharges. In this case, numerical drifters were released during an intense episode of Bora winds that occurred in mid-February 2003, and their trajectories compared to the displacement of satellite-tracked drifters deployed during the ADRIA02-03 sea-truth campaign.Results indicted that the inclusion of the wave breaking process helps improve the accuracy of the numerical simulations, subject to an increase in the typical value of the surface roughness z0. Specifically, the best performance was obtained using αCH = 56,000 in the Charnok formula, the wave breaking parameterization activated, k–ε as the turbulence closure model. With these options, the relative error with respect to the average distance of the drifter was about 25% (5.5 km/day). The most sensitive factors in the model were found to be the value of αCH enhanced with respect to a standard value, followed by the adoption of wave breaking parameterization and the particular turbulence closure model selected. 相似文献
5.
《Ocean Modelling》2002,4(3-4):291-311
Coupled general circulation models (GCMs) have had weak El Niño/Southern Oscillation variability that has been attributed to a diffuse thermocline in the modeled equatorial Pacific Ocean. Consequently, there have been many attempts to improve the thermocline by developing new or improved ocean vertical mixing schemes. This paper investigates the influence of gradient Richardson Number-based vertical mixing scheme profiles in a tropical Pacific Ocean GCM. It has been common for vertical mixing schemes to be assessed in tropical Pacific Ocean models that have a limited latitudinal domain bounded by zonal walls with sponge layers. However, recent work has shown that warm surface water can accumulate in these models and stop them from achieving the observed sharp equatorial thermocline. The present model employs a parameterized wall heat transport scheme that prevents warm surface water from accumulating. Thus we are able assess the influence of vertical mixing profiles in an ocean model that does not allow warm surface water to accumulate and influence the thermocline.In this paper we evaluate the equatorial performance of three different Richardson number (Ri)-based vertical mixing profiles: an integer power (IP) profile based on the observations of Peters, Gregg and Toole; a form of the Pacanowski and Philander profile modified to have low background mixing; and the Max Planck Institute profile. With the accumulation of warm surface water prevented, each of these profiles is able to achieve a sharp thermocline. When compared with observations, the IP profile achieves a better upwelling velocity distribution. We also examine the influence on equatorial performance of very high mixing coefficients at low Richardson number, and of low background mixing coefficients. 相似文献
6.
We present MLD variability over the North Pacific Ocean in a global ocean general circulation model and impacts of three different vertical mixing schemes on it, based on statistical measures (annual mean difference, root-mean-square difference and correlation coefficient). The constant vertical mixing scheme tends to underestimate MLD over the whole basin. The Pacanowski-Philander scheme tends to overestimate MLD (> 20 m) in the mid- to high latitude during summer, implying that vertical mixing in the mid- to high latitude may not be represented properly by simple internal mixing mechanisms such as stratification or vertical shear of horizontal velocity. On the other hand, the new vertical mixing scheme (Noh et al. 2002) gives the most consistent MLD and its seasonal and spatial variability when compared with observation. These results suggest that parameterization of vertical mixing has significant effects on simulation of the seasonal and spatial variability of MLD over the North Pacific Ocean. 相似文献
7.
S. G. Demyshev 《Physical Oceanography》2010,20(3):184-195
Within the framework of the Mellor–Yamada approach, we realize a numerical scheme for the calculation of the coefficients
of turbulent viscosity and diffusion in the z -system of coordinates for the three-dimensional model of operative prediction of currents in the Black Sea. Some discrete
analogs of the equations for turbulent kinetic energy and turbulence macroscale are studied. Their high sensitivity to the
choice of finite-difference approximations is demonstrated. On the basis of the comparison of the results of prognostic experiments
with the data of observations, we choose the best approximation of the term used to describe the generation pf turbulence
energy. 相似文献
8.
The southern portion of the Brazilian coast is dominated by coastal lagoons formed by sandy barrier spits with small inlets. This coastal configuration is a barrier to the surface flow of freshwater to the sea; thus, we suspect that a significant amount of freshwater flows through the permeable sands, beneath the barrier spits, where it mixes with seawater. We excavated an 18-m-deep well into the barrier spit which separates the Patos Lagoon from the South Atlantic. Using this well, we were able to sample interstitial waters from discrete layers, at 1-m intervals, which were analyzed for salinity, temperature, pH, nutrients (ammonium, nitrate, phosphate, and silicate), uranium, molybdenum, and barium. Similar analyses were made on surface water samples from the Patos Lagoon estuarine mixing zone.Results of well samples show a continuous increase in salinity with depth reaching 18 at the bottom. Ammonium and silicate are high, generally around 100 and 100–150 μM, respectively, throughout the subterranean profile. Phosphate shows a distinct maximum at about 6 m (ca. 25 μM), and nitrate is generally low in all well samples. Uranium and molybdenum exhibit a minimum in the well profile at about the same location where barium exhibits a maximum (greater than 2 μM). When results are compared to the surface lagoon–seawater mixing data, ammonium, phosphate, silicate, and barium in well samples of similar salinity show considerable enrichment, while a comparison of uranium and molybdenum data indicates significant depletion of these metals in most well samples.Based on these and other data, we deduce that the following processes are active: products of remineralization of organic detritus accumulated in lagoon sediments are advected through permeable sediments to the oceans; dissolution of biogenic solids and/or solid silicates mobilizes silicate; phosphate, uranium, and molybdenum are mobilized from phosphate-rich sediment layers; sulfate reducers remove uranium and perhaps molybdenum from solution throughout most of the well profile; barium is desorbed from solids in the subterranean mixing zone. These results demonstrate that freshwater discharged to the ocean through permeable sediments may have a significantly different composition than that discharged at the surface. 相似文献
9.
潮致混合对海洋环流的调整起着重要作用。陆架环流的数值模拟中如果不考虑潮汐作用,往往不能得到与观测相符的垂向温盐结构。本文基于调和分析方法,建立了一套潮致混合参数化方案。该方案通过对垂向混合系数进行调和分析,从而得到随时间变化的潮致混合系数。将该方案用于黄海冷水团数值模拟的结果显示,其能够得到与在数值模式开边界直接加入潮汐强迫相当的冷水团温盐结构。和直接引入潮汐强迫相比,这一潮致混合参数化方案的优势在于,它能够大大节省数值模拟计算机时,因此有望显著提高大规模高分辨率的海洋环流及气候模式的模拟能力和计算效率。 相似文献
10.
Shelf break systems are highly dynamic environments. However little is known about the influence that benthic interactions and water mass mixing may have on vertical distributions of iron in these systems. Dissolved Fe (< 0.4 μm) concentrations were measured in samples from nine vertical profiles across the upper slope (150–2950 m water depth) at the Atlantic Ocean–Celtic Sea shelf break. Dissolved iron concentrations varied between less than 0.2 and 5.4 nM, and the resulting detailed section showed evidence of a range of processes influencing the Fe distributions. The near sea floor data were interpreted in terms of release and removal processes. The concentrations of dissolved Fe present in near seabed waters were consistent with release of Fe from in situ remineralisation of particulate organic matter at two upper slope stations, and possibly release from pore water upon resuspension on shelf. Lateral transport of dissolved iron was evident from elevated Fe concentrations in an intermediate nepheloid layer and its advection along isopycnals. Surface waters at the shelf break also showed evidence of vertical mixing of deeper iron-rich waters. These waters contained macronutrients that sustained primary productivity in these otherwise nutrient-depleted surface waters. The data also suggest some degree of stabilisation of relatively high concentrations of iron, presumably through ligand association or as colloids. This study supports the view that lateral export of dissolved iron to the interior of the ocean from shelf and coastal zones and may have important implications for the global budget of oceanic iron. 相似文献
11.
Anne Levasseur Lei Shi Neil C. Wells Duncan A. Purdie Boris A. Kelly-Gerreyn 《Estuarine, Coastal and Shelf Science》2007,73(3-4):753-767
A three-dimensional hydrodynamic model has been developed to simulate water mass circulation in estuarine systems. This model is based on the primitive equation in Cartesian coordinates with a terrain-following structure, coupled with a Mellor–Yamada 2.5 turbulence scheme. A fractional-step method is applied and the subset of equations is solved with finite volume and finite element methods. A dry–wet process simulates the presence of the tidal flat at low water. River inputs are introduced using a point-source method. The model was applied to a partially mixed, macrotidal, temperate estuary: Southampton Water, UK. The model is validated by comparisons with sea surface elevation, ADCP measurements and salinity data collected in 2001. The mean spring range 2(M2 + S2) and the mean neap range 2(M2 − S2) are modelled with an error relative to observation of 12 and 16%, respectively. The unique tidal regime of the system with the presence of the ‘young flood stand’ corresponding to the slackening conditions occurring at mid flood and ‘double high water’ corresponding to an extension of the slackening conditions at high tide is accurately reproduced in the model. The dynamics of the modelled mean surface and bottom velocity closely match the ADCP measurements during neap tides (rms of the difference is 0.09 and 0.01 m s−1 at the bottom and at the surface, respectively), whereas at spring the difference is greater (rms of the difference is 0.25 and 0.20 m s−1 at bottom and surface, respectively). The spatial and temporal variation of the degree of stratification as indicated by salinity distributions compares well with observations. 相似文献
12.
Peter D. Killworth 《Progress in Oceanography》1976,7(2):59-90
The processes which produce bottom water in the Mediterranean are studied. Observations show three phases. The first, called ‘preconditioning’, is not studied. The second ‘violent mixing’ phase occurs when the cold Mistral begins. Cooling at the surface leads to intense vertical mixing in a narrow chimney-like region. When the wind stops, the ‘sinking and spreading’ phase begins. Interleaving of water masses occurs; horizontal eddying on scales of 15 km is observed; and 600 m columns of water can be lifted up to 500 m.A two-dimensional model is used to explain the mixing phase. Non-penetrative vertical convection explains the observations well. Advection of water from outside the column is small, and slows down the descent of the outer parts of the column as observed. The column does not break up even when it reaches the bottom, or if the strong winds cease.The most efficient mechanism for the spreading phase is baroclinic instability. In conditions of vertical overturning there is a large amount of potential energy available to drive the instability. For an eddy viscosity of 50 m2 s−1, after 10 days sinking the growth time is 3–4 days and the eddies are mainly concentrated at the surface with a major axis length of about 15 km. These tend to stabilize the top few hundred meters of the column. Later, finite amplitude effects will produce slower perturbations in the rest of the column. 相似文献
13.
Barotropic tidal mixing effects in a coupled climate model: Oceanic conditions in the Northern Atlantic 总被引:3,自引:0,他引:3
Impacts of mixing driven by barotropic tides in a coupled climate model are investigated by using an atmosphere–ocean–ice–land coupled climate model, the GFDL CM2.0. We focus on oceanic conditions of the Northern Atlantic. Barotropic tidal mixing effects increase the surface salinity and density in the Northern Atlantic and decrease the RMS error of the model surface salinity and temperature fields related to the observational data. 相似文献
14.
This study considers an important biome in aquatic environments, the subsurface ecosystem that evolves under low mixing conditions, from a theoretical point of view. Employing a conceptual model that involves phytoplankton, a limiting nutrient and sinking detritus, we use a set of key characteristics (thickness, depth, biomass amplitude/productivity) to qualitatively and quantitatively describe subsurface biomass maximum layers (SBMLs) of phytoplankton. These SBMLs are defined by the existence of two community compensation depths in the water column, which confine the layer of net community production; their depth coincides with the upper nutricline. Analysing the results of a large ensemble of simulations with a one-dimensional numerical model, we explore the parameter dependencies to obtain fundamental steady-state relationships that connect primary production, mortality and grazing, remineralization, vertical diffusion and detrital sinking. As a main result, we find that we can distinguish between factors that determine the vertically integrated primary production and others that affect only depth and shape (thickness and biomass amplitude) of this subsurface production layer. A simple relationship is derived analytically, which can be used to estimate the steady-state primary productivity in the subsurface oligotrophic ocean. The fundamental nature of the results provides further insight into the dynamics of these “hidden” ecosystems and their role in marine nutrient cycling. 相似文献
15.
Hydrographic changes in the Labrador Sea, 1960–2005 总被引:1,自引:0,他引:1
The Labrador Sea has exhibited significant temperature and salinity variations over the past five decades. The whole basin was extremely warm and salty between the mid-1960s and early 1970s, and fresh and cold between the late 1980s and mid-1990s. The full column salinity change observed between these periods is equivalent to mixing a 6 m thick freshwater layer into the water column of the early 1970s. The freshening and cooling trends reversed in 1994 starting a new phase of heat and salt accumulation in the Labrador Sea sustained throughout the subsequent years. It took only a decade for the whole water column to lose most of its excessive freshwater, reinstate stratification and accumulate enough salt and heat to approach its record high salt and heat contents observed between the late 1960s and the early 1970s. If the recent tendencies persist, the basin’s storages of salt and heat will fairly soon, likely by 2008, exceed their historic highs.The main process responsible for the net cooling and freshening of the Labrador Sea between 1987 and 1994 was deep winter convection, which during this period progressively developed to its record depths. It was caused by the recurrence of severe winters during these years and in its turn produced the deepest, densest and most voluminous Labrador Sea Water (LSW1987–1994) ever observed. The estimated annual production of this water during the period of 1987–1994 is equivalent to the average volume flux of about 4.5 Sv with some individual annual rates exceeding 7.0 Sv. Once winter convection had lost its strength in the winter of 1994–1995, the deep LSW1987–1994 layer lost “communication” with the mixed layer above, consequently losing its volume, while gaining heat and salt from the intermediate waters outside the Labrador Sea.While the 1000–2000 m layer was steadily becoming warmer and saltier between 1994 and 2005, the upper 1000 m layer experienced another episode of cooling caused by an abrupt increase in the air-sea heat fluxes in the winter of 1999–2000. This change in the atmospheric forcing resulted in fairly intense convective mixing sufficient to produce a new prominent LSW class (LSW2000) penetrating deeper than 1300 m. This layer was steadily sinking or deepening over the years following its production and is presently overlain by even warmer and apparently less dense water mass, implying that LSW2000 is likely to follow the fate of its deeper precursor, LSW1987–1994. The increasing stratification of the intermediate layer implies intensification in the baroclinic component of the boundary currents around the mid-depth perimeter of the Labrador Sea.The near-bottom waters, originating from the Denmark Strait overflow, exhibit strong interannual variability featuring distinct short-term basin-scale events or pulses of anomalously cold and fresh water, separated by warm and salty overflow modifications. Regardless of their sign these anomalies pass through the abyss of the Labrador Sea, first appearing at the Greenland side and then, about a year later, at the Labrador side and in the central Labrador Basin.The Northeast Atlantic Deep Water (2500–3200 m), originating from the Iceland–Scotland Overflow Water, reached its historically freshest state in the 2000–2001 period and has been steadily becoming saltier since then. It is argued that LSW1987–1994 significantly contributed to the freshening, density decrease and volume loss experienced by this water mass between the late 1960s and the mid 1990s via the increased entrainment of freshening LSW, the hydrostatic adjustment to expanding LSW, or both. 相似文献
16.
《Ocean Modelling》2004,6(3-4):245-263
Astronomical data reveals that approximately 3.5 terawatts (TW) of tidal energy is dissipated in the ocean. Tidal models and satellite altimetry suggest that 1 TW of this energy is converted from the barotropic to internal tides in the deep ocean, predominantly around regions of rough topography such as mid-ocean ridges. A global tidal model is used to compute turbulent energy levels associated with the dissipation of internal tides, and the diapycnal mixing supported by this energy flux is computed using a simple parameterization.The mixing parameterization has been incorporated into a coarse resolution numerical model of the global ocean. This parameterization offers an energetically consistent and practical means of improving the representation of ocean mixing processes in climate models. Novel features of this implementation are that the model explicitly accounts for the tidal energy source for mixing, and that the mixing evolves both spatially and temporally with the model state. At equilibrium, the globally averaged diffusivity profile ranges from 0.3 cm2 s−1 at thermocline depths to 7.7 cm2 s−1 in the abyss with a depth average of 0.9 cm2 s−1, in close agreement with inferences from global balances. Water properties are strongly influenced by the combination of weak mixing in the main thermocline and enhanced mixing in the deep ocean. Climatological comparisons show that the parameterized mixing scheme results in a substantial reduction of temperature/salinity bias relative to model solutions with either a uniform vertical diffusivity of 0.9 cm2 s−1 or a horizontally uniform bottom-intensified arctangent mixing profile. This suggests that spatially varying bottom intensified mixing is an essential component of the balances required for the maintenance of the ocean’s abyssal stratification. 相似文献
17.
Experiments with a climate model were conducted under present day and last glacial maximum conditions in order to examine the model’s response to a vertical mixing scheme based on internal tide energy dissipation. The increase in internal tide energy flux caused by a 120 m reduction in sea level had the expected effect on diffusivity values, which were higher under lower sea level conditions. The impact of this vertical diffusivity change on the Atlantic meridional overturning is not straightforward and no clear relationship between diffusivity and overturning is found. There exists a weak positive correlation between overturning and changes to the power consumed by vertical mixing. Most of the climatic response generated by sea level change was not related to alterations in the internal tide energy flux but rather to the direct change in sea level itself. 相似文献
18.
Norman Silverberg Aida Martínez Sergio Aguíiga Jos D. Carriquiry Nancy Romero Evgeni Shumilin Soledad Cota 《Estuarine, Coastal and Shelf Science》2004,59(4):575-587
The vertical flux of particulate matter at 330 m depth in San Lázaro Basin off southern Baja California ranged from 63 to 587 mg m−2 d−1 between August and November 1996. Organic carbon contents were between 5.6 and 14.8%, yielding flux rates of 9–40 mgC m−2 d−1. In December 1997 and January 1998, at the height of the strong El Niño event, the respective fluxes (47–202 mg m−2 d−1 and 3–8 mgC m−2 d−1) were comparable. The February–June 1998 records, however, revealed sharply reduced mass (1–6 mg m−2 d−1) and organic carbon (0.2–0.8 mgC m−2 d−1) fluxes. The organics collected in 1996 were predominantly autochthonous (δ13C=−22‰; C/N=8). The variations in δ15N (8.3–11.0‰) suggest an alternation of new and regenerated production, possibly associated with fluctuations in the intensity of deep mixing during that autumn. The relatively high organic matter fluxes in December 1997 appear to be associated with regenerated production. The average composition from February to June 1998 (δ13C=−23.6‰; 15N=11.7‰; C/N=10.5) indicates degraded material of marine origin. The maximum δ15N value found (14‰) suggests that deeper, denitrified waters were brought to the surface and possibly advected laterally. Regime changes in the waters of the basin occur at 6–10 week intervals, evidenced by concurrent shifts in most of the measured parameters, including fecal pellet types and metal chemistry. The marine snow-dominated detritus collected showed a shift from a mixed diatom-rich-radiolarian-coccolith assemblage in late 1996 to a coccolith-dominated assemblage, including the contents of fecal pellets, during the 1997–1998 El-Niño period. T–S profiles, plankton analysis and chlorophyll contents of the upper water column indicated that the strong phytoplankton bloom, normally associated with seasonal upwelling along the Pacific coast of Baja, did not occur during the spring of 1998. The persistence of oligotrophic conditions during the 1997–1998 El Niño event favored the dominance of nanoplankton and reduced the vertical flux of particles. 相似文献
19.
A coupled wave–tide–surge model has been developed in this study in order to investigate the effect of the interactions among tides, storm surges, and wind waves. The coupled model is based on the synchronous dynamic coupling of a third-generation wave model, WAM cycle 4, and the two-dimensional tide–surge model. The surface stress, which is generated by interactions between wind and wave, is calculated by using the WAM model directly based on an analytical approximation of the results using the quasi-linear theory of wave generation. The changes in bottom friction are created by the interactions between waves and currents and calculated by using simplified bottom boundary layer model. In consequence, the combined wave–current-induced bottom velocity and effective bottom drag coefficient were increased in the shallow waters during the strong storm conditions. 相似文献
20.
在大洋环流模型中对应于惯性内波破碎的垂直混合的一个参数化方案 总被引:1,自引:0,他引:1
Based on the theoretical spectral model of inertial internal wave breaking(fine structure) proposed previously, in which the effects of the horizontal Coriolis frequency component f-tilde on a potential isopycnal are taken into account, a parameterization scheme of vertical mixing in the stably stratified interior below the surface mixed layer in the ocean general circulation model(OGCM) is put forward preliminarily in this paper. Besides turbulence, the impact of sub-mesoscale oceanic processes(including inertial internal wave breaking product) on oceanic interior mixing is emphasized. We suggest that adding the inertial internal wave breaking mixing scheme(F-scheme for short) put forward in this paper to the turbulence mixing scheme of Canuto et al.( T-scheme for short) in the OGCM, except the region from 15°S to 15°N. The numerical results of F-scheme by using WOA09 data and an OGCM(LICOM, LASG/IAP climate system ocean model) over the global ocean are given. A notable improvement in the simulation of salinity and temperature over the global ocean is attained by using T-scheme adding F-scheme, especially in the mid- and high-latitude regions in the simulation of the intermediate water and deep water. We conjecture that the inertial internal wave breaking mixing and inertial forcing of wind might be one of important mechanisms maintaining the ventilation process. The modeling strength of the Atlantic meridional overturning circulation(AMOC) by using T-scheme adding F-scheme may be more reasonable than that by using T-scheme alone, though the physical processes need to be further studied, and the overflow parameterization needs to be incorporated. A shortcoming in F-scheme is that in this paper the error of simulated salinity and temperature by using T-scheme adding F-scheme is larger than that by using T-scheme alone in the subsurface layer. 相似文献