Multiscale asymptotics are used to derive three systems of equations connecting the planetary geostrophic (PG) equations for gyre-scale flow to a quasigeostrophic (QG) equation set for mesoscale eddies. Pedlosky (1984), following similar analysis, found eddy buoyancy fluxes to have only a small effect on the large-scale flow; however, numerical simulations disagree. While the impact of eddies is relatively small in most regions, in keeping with Pedlosky's result, eddies have a significant effect on the mean flow in the vicinity of strong, narrow currents. First, the multiple-scales analysis of Pedlosky is reviewed and amplified. Novel results of this analysis include new multiple-scales models connecting large-scale PG equations to sets of QG eddy equations. However, only introducing anisotropic scaling of the large-scale coordinates allows us to derive a model with strong two-way coupling between the QG eddies and the PG mean flow. This finding reconciles the analysis with simulations, viz, that strong two-way coupling is observed in the vicinity of anisotropic features of the mean flow like boundary currents and jets. The relevant coupling terms are shown to be eddy buoyancy fluxes. Using the Gent-McWilliams parameterization to approximate these fluxes allows solution of the PG equations with closed tracer fluxes in a closed domain, which is not possible without mesoscale eddy (or other small-scale) effects. The boundary layer width is comparable to an eddy mixing length when the typical eddy velocity is taken to be the long Rossby wave phase speed, which is the same result found by Fox-Kemper and Ferrari (2009) in a reduced gravity layer. (C) 2011 Elsevier B.V. All rights reserved.
