Vertical mixing in the bottom mixed layer near the continental shelf break in the East China Sea

The strength of the vertical mixing in the bottom mixed layer near the continental shelf break in the East China Sea was directly measured with the Micro-Scale Profiler (MSP). It has been shown that there is no significant statistical relation between the turbulent energy dissipation and the degree of the stratificationN ². It seems that the vigorous turbulence occurs not constantly but intermittently in the bottom mixed layer so that a large variation of is found depending on the time. In contrast to , the coefficient of the vertical eddy diffusivityK _z is mostly determined byN such thatK _z is large in the bottom mixed layer and small in the thermocline. Large value ofK _z in the bottom mixed layer is also found in the time series ofK _z estimated in terms of Richardson number calculated from the data obtained with electromagnetic current meters. The value ofK _z more than 10 cm²s^–1 frequently occur in the layer of 20–25 m thick just above the bottom.     

Direct measurements of diapycnal mixing in a fjord reach—Puget Sound's Main Basin

In the spring of 1988, time series of microstructure and ADCP current profiles were collected at four locations in the North Main Basin of Puget Sound, Washington. Depth and time averages of diapycnal diffusivity at the four stations (1.8−67.0×10⁻⁴ m² s⁻¹) were one to three decades above typical open-ocean thermocline levels. The buoyancy frequency-squared N² was near open-ocean levels, but unlike the open-ocean where N²S², finescale shear-squared S² was three to six times N² over significant portions of the water column at two of the stations. The time and space mean of all measurements ( ) is close to inferred vertical eddy diffusivity from a primitive equation model for Puget Sound (K_z=3×10⁻³ m² s⁻¹) (J. Geophys. Res. 96 (1991) 16779). Large time and space variability of K_ρ was found, with differences of inter-station, depth–time means over one decade. A simple scaling argument using the observed K_ρ suggests significant exchange of mass between the layers of the subtidal flow over the basin's residence time. Additionally, measurements show that local mixing may be comparable to volume-weighted sill mixing in modifying the Main Basin's stratification. Both are contrary to the “advective reach” simplification of fjord dynamics. The mixing levels were dominated by the passage of a mid-depth, southward-flowing density intrusion and what we interpret as a strongly advected, non-linear internal tide. These mechanisms elevated profile-averaged K_ρ by more than 10 times background levels, with sustained patches of K_ρ≥1×10⁻² m² s⁻¹. Critical 8-m gradient Richardson numbers (Ri₈<0.25) matching regions of overturns (>20 m) and strong turbulence suggest that shear instabilities dominated the turbulence production, though there was support for double-diffusive convection in the warm core of the density intrusion.     

Data Assimilation Modeling of the Barotropic Tides in the Korea/Tsushima Strait

During 1999–2000, 13 bottom mounted acoustic Doppler current profilers (ADCPs) and 12 wave/tide gauges were deployed along two lines across the Korea/Tsushima Strait, providing long-term measurements of currents and bottom pressure. Tidally analyzed velocity and pressure data from the moorings are used in conjunction with other moored ADCPs, coastal tide gauge measurements, and altimeter measurements in a linear barotropic data assimilation model. The model fits the vertically averaged data to the linear shallow water equations in a least-squares sense by only adjusting the incoming gravity waves along the boundaries. Model predictions are made for the O₁, P₁, K₁, μ₂, N₂, M₂, S₂, and K₂ tides. An extensive analysis of the accuracy of the M₂ surface-height predictions suggests that for broad regions near the mooring lines and in the Jeju Strait the amplitude prediction errors are less than 0.5 cm. Elsewhere, the analysis suggests that errors range from 1 to 4 cm with the exception of small regions where the tides are not well determined by the dataset. The errors in the model predictions are primarily caused by bias error in the model’s physics, numerics, and/or parameterization as opposed to random errors in the observational data. In the model predictions, the highest ranges in sea level height occur for tidal constituents M₂, S₂, K₁, O₁, and N₂, with the highest magnitudes of tidal velocities occurring for M₂, K₁, S₂, and O₁. The tides exhibit a complex structure in which diurnal constituents have higher currents relative to their sea level height ranges than semi-diurnal constituents.     

Formation of the salinity minimum in the Mixed Water Region between the Oyashio and Kuroshio Fronts

The salinity minimum frequently occurring in the Mixed Water Region between the Oyashio and Kuroshio Fronts seems to originate from the salinity minimum at the density of 26.8σ_θ called the North Pacific Intermediate Water. We examined water exchange of this region with the Oyashio and the Kuroshio Extension using mixing ratio R_K defined as (θ - θ_OY)/(θ_K - θ_OY) × 100, where θ_OY, θ_K, and θ represent potential temperature of the Oyashio and Kuroshio Waters and their mixture on the isopycnal surfaces, respectively. CTD data were obtained by repeated observation from January 1990 to May 1991. R_K increases southward from the Oyashio Front to the Kuroshio Front with the range of −20 to 120%. The gradient of R_K on the isopycnal surfaces is large around the Oyashio Front above the 26.8σ_θ surface, while it is large around the Kuroshio Front below it. This agrees with the average R_K in the Mixed Water Region decreasing greatly with the increase of density at densities less dense than 26.8σ_θ. We calculated thickness and volume transport of the Oyashio between the isopycnal surfaces near the coast of Hokkaido. They increase largely with density at densities less dense than 26.8σ_θ. It is supposed that the salinity minimum in the Mixed Water Region is the upper limit of the water largely influenced by the Oyashio Water. Its density could depend only on the density structure of the Oyashio.     

Sorption of yttrium and rare earth elements by amorphous ferric hydroxide: Influence of pH and ionic strength

The sorption of yttrium and the rare earth elements (YREEs) by amorphous ferric hydroxide at low ionic strength (0.01 M ≤ I ≤ 0.09 M) was investigated over a wide range of pH (3.9 ≤ pH ≤ 7.1). YREE distribution coefficients, defined as _iK_Fe = [MS_i]_T / (M_T[Fe³⁺]_S), where [MS_i]_T is the concentration of YREE sorbed by the precipitate, M_T is the total YREE concentration in solution, and [Fe³⁺]_S is the concentration of precipitated iron, are weakly dependent on ionic strength but strongly dependent on pH. For each YREE, the pH dependence of log _iK_Fe is highly linear over the investigated pH range. The slopes of log _iK_Fe versus pH regressions range between 1.43 ± 0.04 for La and 1.55 ± 0.03 for Lu. Distribution coefficients are well described by an equation of the form _iK_Fe = (_Sβ₁[H⁺]^− 1 + _Sβ₂[H⁺]^− 2) / (_SK₁[H⁺] + 1), where _Sβ_n are stability constants for YREE sorption by surface hydroxyl groups and _SK₁ is a ferric hydroxide surface protonation constant. Best-fit estimates of _Sβ_n for each YREE were obtained with log _SK₁ = 4.76. Distribution coefficient predictions, using this two-site surface complexation model, accurately describe the log _iK_Fe patterns obtained in the present study, as well as distribution coefficient patterns obtained in previous studies at near-neutral pH. Modeled log _iK_Fe results were used to predict YREE sorption patterns appropriate to the open ocean by accounting for YREE solution complexation with the major inorganic YREE ligands in seawater. The predicted log _iK_Fe′ pattern for seawater, while distinctly different from log _iK_Fe observations in synthetic solutions at low ionic strength, is in good agreement with results for natural seawater obtained by others.     

Tidal energy balance in the Seto Inland Sea

The tidal volume transport in the Seto Inland Sea is calculated. The cross-section where the volume transport of the M₂ tide is zero, is located around the western part of Bisan Strait. The tidal energy dissipation of the M₂ tide by friction is 6.30×10¹⁶ ergs s^–1 in the Seto Inland Sea. The quality factorQ for the M₂ tide is 20.2. The total energy dissipation of the M₂, S₂, K₁ and O₁ tides is 7.99×10¹⁶ ergs s^–1.     

珠江河口潮能收支研究

Tidal energy budget in the Zhujiang(Pearl River) Estuary(ZE) is evaluated by employing high-resolution baroclinic regional ocean modeling system(ROMS). The results obtained via applying the least square method on the model elevations are compared against the tidal harmonic constants at 18 tide stations along the ZE and its adjacent coast. The mean absolute errors between the simulation and the observation of M_2, S_2, K_1 and O_1 are 4.6, 2.8, 3.2 and 2.8 cm in amplitudes and 9.8°, 15.0°, 4.6° and 4.6° in phase-lags, respectively. The comparisons between the simulated and observed sea level heights at 11 tide gauge stations also suggest good model performance. The total tidal energy flux incoming the ZE is estimated to be 343.49 MW in the dry season and larger than 336.18 MW in the wet season, which should due to higher mean sea level height and heavier density in the dry season. M_2, K_1, S_2, O_1 and N_2, the top five barotropic tidal energy flux contributors for the ZE,import 242.23(236.79), 52.97(52.08), 24.49(23.96), 16.22(15.91) and 7.10(6.97) MW energy flux into the ZE in dry(wet) season, successively and respectively. The enhanced turbulent mixing induced by eddies around isolated islands and sharp headlands dominated by bottom friction, interaction between tidal currents and sill topography or constricted narrow waterways together account for the five energy dissipation hotspots, which add up to about 38% of the total energy dissipation inside the ZE.     

The speciation of Fe(II) and Fe(III) in natural waters

The interactions of Fe(II) and Fe(III) with the inorganic anions of natural waters have been examined using the specific interaction and ion pairing models. The specific interaction model as formulated by Pitzer is used to examine the interactions of the major components (Na⁺, Mg²⁺, Ca²⁺, K⁺, Sr²⁺, Cl⁻, SO₄⁻, HCO₃⁻, Br⁻, CO₃²⁻, B(OH)₄⁻, B(OH)₃ and CO₂) of seawater and the ion pairing model is used to account for the strong interaction of Fe(II) and Fe(III) with major and minor ligands (Cl⁻, SO₄²⁻, OH⁻, HCO₃⁻, CO₃²⁻ and HS⁻) in the waters. The model can be used to estimate the activity and speciation of iron in natural waters as a function of composition (major sea salts) and ionic strength (0 to 3 M). The measured stability constants (K_FeX^*) of Fe(II) and Fe(III) have been used to estimate the thermodynamic constants (K_FeX) and the activity coefficient of iron complexes (γ_FeX) with a number of inorganic ligands in NaClO₄ medium at various ionic strengths: In(K_FeX/γ_Feγ_X) = InK_FeX − In(γ_FeX) The activity coefficients for free ions (γ_Fe, γ_x) needed for this extrapolation have been estimated from the Pitzer equations. The activity coefficients of the ion pairs have been used to determine Pitzer parameters (B_FeX, B_FeX⁰, C_FeX^φ) for the iron complexes. These results make it possible to estimate the stability constants for the formation of Fe(II) and Fe(III) complexes over a wide range of ionic strengths and in different media. The model has been used to determine the solubility of Fe(III) in seawater as a function of pH. The results are in good agreement with the measurements of Byrne and Kester and Kuma et al. When the formation of Fe organic complexes is considered, the solubility of Fe(III) in seawater is increased by about 25%.     

Turbulent structure in the upper layer of the western Equatorial Pacific Ocean

The Richardson number dependence of vertical eddy diffusion coefficients in the western Equatorial Pacific Ocean was examined on the basis of a Microstructure Profiler (MSP) observations during the cruise of Natsushima (JAPACS-89). The Richardson numberR _i was estimated by using the mean shear of velocity profile measured by an Acoustic Doppler Current Profiler (ADCP) with the vertical interval of 15 meters within one or two hours of the each MSP cast. The raw data plot of the vertical eddy diffusion coefficientK _p shows a large scatter with increasing tendency belowR _i =0.5. The relation between the mean vertical eddy diffusion coefficientK _p and the Richardson numberR _i , averaged over every 0.025 in theR _i , supports the model of Pacanowski and Philander (1981) in the range ofR _i >0.5, but coincides with the result of Peterset al. (1988) in the range ofR _i <0.5.     

Kinetics of microbially mediated reactions: dissimilatory sulfate reduction in saltmarsh sediments (Sapelo Island, Georgia, USA)

A sediment disk reactor was tested in once flow-through mode to retrieve kinetic parameters for the Monod rate law that describes sulfate reduction. The experimental method was compared with a previously described procedure by the authors where a sediment plug-flow reactor was operated in a recirculation mode. In recirculation mode, accumulation of metabolic byproducts in certain cases may result in negative feedback, thus preventing accurate determination of kinetic information. The method described in this article provides an alternative to the recirculation sediment plug-flow-through reactor technique for retrieving kinetic parameters of microbially mediated reactions in aquatic sediments.For sulfate reduction in a saltmarsh site, a maximum estimate of the half-saturation concentration, K_s, of 204±26 μM and a maximum reaction rate, R_m, of 2846±129 nmol cm_{(wetsediment)}³ d⁻¹ was determined. The K_s value obtained was consistent with the one estimated previously (K_s=240±20 μM) from a different site within the same saltmarsh mud flat using a recirculating reactor. From the R_m value and reduction rates determined using ³⁵SO₄²⁻ incubation experiments, we infer that sulfate reduction is limited in the field. Substrate availability is not the main contributor for the limitation, however. Competition from other microbes, such as iron reducers affects the activity of sulfate reducers in the suboxic to anoxic zones, whereas aerobes compete in the oxic zone. High sulfide concentration in the pore water may also have acted as a toxin to the sulfate reducers in the field.     

The inception of sheet flow in oscillatory flow

A simple model is developed to study the inception of sheet flow in oscillatory flow based on the available experimental data. The inception of sheet flow in oscillatory flow is well defined by the simple model: A/d=KA²ω/ν+B, where A is the semi-excursion of wave orbital motion near the bed, d is the grain size, ω is the angular frequency, ν is the kinematic viscosity of water, and K and B are the coefficients and dependent on sediment properties only. The inception velocity of sheet flow derived from the model is shown to be the function of grain size d, oscillatory period T and specific sediment density s. For a given sediment, the inception velocity is found to increase sharply initially with T and then approach a constant at T>6.0 s. The present model is quite simple and gives good agreement with the available experimental data.     

Eddy Diffusion Coefficients of Suspended Particulate Matter: Effects of Wind Energy Transfer and Stratification

Gross sedimentation rates (GSR) have been measured using sediment traps placed at nine different levels above the bed (0·3, 0·5, 0·8, 1·0, 2·0, 4·0, 6·0, 8·0 and 10·0 m). The sediment traps were deployed for 1·25 years and recovered 28 times during the study period. Low average GSR values of 5·5 g m^-2 day^-1 were obtained at 10·0 m, and high average GSR values of 114·8 g m^-2 day^-1 were obtained at 0·3 m. An expression for the eddy diffusion coefficient of suspended particulate matter (K_s), based on the measured GSR is given. The expression has been used for modelling of K_s at the different trap levels above the bed. High values (≈42 cm² s^-1) of K_s were obtained at the upper traps, whereas low values (≈2 cm² s^-1) were obtained near the bed. Comparison between level of turbulent energy in terms of shear stress at the boundaries of the water column, i.e. from the wind and the bed flow, showed that wind energy exceeded that of the bed flow by a factor 16. At 5·0 m K_s was positively correlated (r=0·66) to the eddy diffusion coefficient of momentum (K_m) derived from the wind energy transfer to the water, giving an average β of 0·5 for K_s =βK_m. The density difference between surface and bottom waters has been designated a parameter of stratification, and is discussed in relation to variations of K_s and K_m .     

Solubility of calcite in the ocean

The apparent solubility product K′_sp of calcite in seawater was measured as a function of temperature, salinity, and pressure using potentiometric saturometry techniques. The temperature effect was hardly discernible experimentally. The value of K′_sp at 25°C was 4.59·10⁻⁷ mole²/(kg seawater)² at 35‰S, 5.34·10⁻⁷ at 43‰S, and 3.24·10⁻⁷ at 27‰S. The apparent partial molal volume was found to be −34.4 cm³ at 25°C and −42.3 cm³ at 2°C from a linear fit of log(K′_sp ^P/K′_sp ¹). These results were used in conjunction with field data to calculate the degree of saturation in the oceans and showed undersaturation at shallower depths than previously reported.     

Relationships between suspended particulate material, light attenuation and Secchi depth in UK marine waters

Measurements of sub-surface light attenuation (K_d), Secchi depth and suspended particulate material (SPM) were made at 382 locations in transitional, coastal and offshore waters around the United Kingdom (hereafter UK) between August 2004 and December 2005. Data were analysed statistically in relation to a marine water typology characterised by differences in tidal range, mixing and salinity. There was a strong statistically significant linear relationship between SPM and K_d for the full data set. We show that slightly better results are obtained by fitting separate models to data from transitional waters and coastal and offshore waters combined. These linear models were used to predict K_d from SPM. Using a statistic (D) to quantify the error of prediction of K_d from SPM, we found an overall prediction error rate of 23.1%. Statistically significant linear relationships were also evident between the log of Secchi depth and the log of K_d in waters around the UK. Again, statistically significant improvements were obtained by fitting separate models to estuarine and combined coastal/offshore data – however, the prediction error was improved only marginally, from 31.6% to 29.7%. Prediction was poor in transitional waters (D = 39.5%) but relatively good in coastal/offshore waters (D = 26.9%).SPM data were extracted from long term monitoring data sites held by the UK Environment Agency. The appropriate linear models (estuarine or combined coastal/offshore) were applied to the SPM data to obtain representative K_d values from estuarine, coastal and offshore sites. Estuarine waters typically had higher concentrations of SPM (8.2–73.8 mg l⁻¹) compared to coastal waters (3.0–24.1 mg l⁻¹) and offshore waters (9.3 mg l⁻¹). The higher SPM values in estuarine waters corresponded to higher values of K_d (0.8–5.6 m⁻¹). Water types that were identified by large tidal ranges and exposure typically had the highest K_d ranges in both estuarine and coastal waters. In terms of susceptibility to eutrophication, large macrotidal, well mixed estuarine waters, such as the Thames embayment and the Humber estuary were identified at least risk from eutrophic conditions due to light-limiting conditions of the water type.     

Kinetic and equilibrium studies of copper-dissolved organic matter complexation in water column of the stratified Krka River estuary (Croatia)

An interaction of dissolved natural organic matter (DNOM) with copper ions in the water column of the stratified Krka River estuary (Croatia) was studied. The experimental methodology was based on the differential pulse anodic stripping voltammetric (DPASV) determination of labile copper species by titrating the sample using increments of copper additions uniformly distributed on the logarithmic scale. A classical at-equilibrium approach (determination of copper complexing capacity, CuCC) and a kinetic approach (tracing of equilibrium reconstitution) of copper complexation were considered and compared. A model of discrete distribution of organic ligands forming inert copper complexes was applied. For both approaches, a home-written fitting program was used for the determination of apparent stability constants (K_i^equ), total ligands concentration (L_iT) and association/dissociation rate constants (k_i¹,k_i^- 1).A non-conservative behaviour of dissolved organic matter (DOC) and total copper concentration in a water column was registered. An enhanced biological activity at the freshwater–seawater interface (FSI) triggered an increase of total copper concentration and total ligand concentration in this water layer. The copper complexation in fresh water of Krka River was characterised by one type of binding ligands, while in most of the estuarine and marine samples two classes of ligands were identified. The distribution of apparent stability constants (log K₁^equ: 11.2–13.0, log K₂^equ:8.8–10.0) showed increasing trend towards higher salinities, indicating stronger copper complexation by autochthonous seawater organic matter.Copper complexation parameters (ligand concentrations and apparent stability constants) obtained by at-equilibrium model are in very good accordance with those of kinetic model. Calculated association rate constants (k₁¹:6.1–20 × 10³ (M s)^− 1, k₂¹: 1.3–6.3 × 10³ (M s)^− 1) indicate that copper complexation by DNOM takes place relatively slowly. The time needed to achieve a new pseudo-equilibrium induced by an increase of copper concentration (which is common for Krka River estuary during summer period due to the nautical traffic), is estimated to be from 2 to 4 h.It is found that in such oligotrophic environment (dissolved organic carbon content under 83 µM_C, i.e. 1 mg_CL^− 1) an increase of the total copper concentration above 12 nM could enhance a free copper concentration exceeding the level considered as potentially toxic for microorganisms (10 pM).     

Tidal Currents in the Tsushima Straits Estimated from ADCP Data by Ferryboat

Tidal currents in the Tsushima Straits have been analyzed using measurements obtained since February 1997 by an acoustic Doppler current profiler (ADCP) mounted on the ferryboat Camellia. Tidal current constituents (M ₂, S ₂, K ₁, O ₁) are dominant among the ten tidal current constituents (Q ₁, O ₁, P ₁, K ₁, N ₂, M ₂, S ₂, K ₂, MSf, Mf), and generally 1.4–2.1 times stronger at the western channel of the straits than those at the eastern channel. The ratio between amplitude of M ₂, S ₂, K ₁ and O ₁ averaged along the ferryboat track is 1:0.45:0.59:0.51. The major axis directions of tidal current ellipses are generally SW to NE, exceptionally in the vicinity of the Tsushima Islands. Approaching the Tsushima Islands from the Korean Peninsula side, the major axis gradually rotates clockwise. At the western channel, the M ₂ and K ₁ constituents change the rotation direction of current vectors from clockwise to counterclockwise at about 90–130 m depth. The contributions of the tidal currents to the mean kinetic energy and the mean eddy kinetic energy along the ferryboat track are, on average, 0.56 and 0.71, respectively. This suggests that tidal current activities are generally more dominant than the mean current activities and much more dominant than eddy activities. The only region where the eddy activities are comparable to the tidal current activities is located on the east side of the Tsushima Islands. This revised version was published online in July 2006 with corrections to the Cover Date.     

Turbulence investigation in a tidal coastal region

A high-frequency (1.2 MHz) four-beam Acoustic Doppler Current Profiler (ADCP) moored on the sea bottom was used for the direct measurements of the turbulence parameters in the shallow (20 m) coastal zone of the eastern English Channel. The measurements were as long as four tidal cycles during the period of the spring tide development. The measurements in the ocean and estimates showed that the Reynolds stress variability coincided with the semidiurnal tide. Their maximum values during the flood phase were approximately 1.5 Pa, while, during the ebb phase, they reached −1.2 Pa. The variations of the turbulence’s kinetic energy (TKE) and the rate of its production (P) coincided with the period of the tidal harmonic M4. Their maximum values were found during the flood phase near the bottom, and they were approximately equal to 0.03 m²/s² and 0.8 W/m³, respectively. These values decreased rapidly with the distance from the bottom. During the periods of low stagnant water, the values of TKE and P in the water column decreased to the minimum values (2 × 10⁻³ m²/s² and 3 × 10⁻⁵ W/m³, respectively), which coincided with the moment of the current’s reversal flow. The results demonstrated the dominating role of the tidal motion, which controls the structure and intensity of the turbulence in the bottom layer, and revealed the characteristic asymmetry of its distribution related to the nonlinear character of the tidal cycle.     

A regeneration model to estimate addition fluxes and removal constants of trace metals in seawater

A model is presented for flux evaluation of trace metals in situ. For the practically feasible simplified model, mere estimation of metals in various dissolved states (e.g. Cr(III) and Cr(VI)) and in particulate phase, in a single vertical profile at least, together with temperature, salinity, O₂, titration alkalinity (TA), Ca²⁺, NO₃⁻, HPO₄²⁻ and H₄SiO₄ parameters enable evaluation of the amount of metal removed in situ, to a first approximation. From the evaluated net metal removed (M_R), estimation could be made for removal constants, i.e. removal times (t), residence times (τ) and removal rate constants (ψ), not only for the average ocean but also for the oxygen-minimum zone. Besides, M_R can be used to calculate t, τ and ψ of various species of the metal concerned, that are not yet available. The present model evaluates t, τ and ψ for total chromium from a single station, which are consistent with the reported values, and shows that chromium is removed with fast-settling particulate matter. This study also reveals that t, τ and ψ values of 0.001 years 7 years and 0.150 y⁻¹, respectively, for Cr(III) in the oxygen-minimum zone and 0.006 y, 211 y and 0.005 y⁻¹ in the marine environment, respectively.     

On the importance of resolving the western boundary layer in wind-driven ocean general circulation models

The western boundary layer (WBL) plays a fundamental role in basin-scale wind-driven ocean circulations. In idealized ocean models with flat bottom topography, this layer is required not only to balance the interior Sverdrup transport to close the gyre circulation, but also to dissipate the vorticity imposed by the wind-stress curl. The width of the WBL in Munk-type models is estimated to be δ_M(A_H/β)^1/3, where A_H and β are horizontal eddy viscosity and the meridional derivative of the Coriolis parameter respectively. For commonly used values of A_H, the boundary-layer width δ_M ranges from 30 to less than 200 km in the mid-latitude ocean. This scale is often poorly resolved in large-scale climate models.This paper intends to demonstrate some consequences when the western boundary layer is not adequately resolved. It is found that coarse resolution models reach equilibrium states by distorting some important dynamics in order to dissipate wind-imposed vorticity. In three-dimensional models, for instance, very strong spurious upwelling and downwelling can occur along the WBL. In models of two-dimensional flow, however, spurious recirculations may develop near the boundary. These false features can be removed when the boundary layer is better resolved. We propose a method in which a spatially varying A_H is used to broaden the WBL without affecting mixing in the interior. The method improves the model results considerably.