The importance of altimeter and scatterometer data for ocean prediction

The prediction of ocean circulation using satellite altimeter data is discussed. Three classes of oceanic response to atmospheric forcing are outlined and examined. Storms, surface waves, eddies, and ocean currents were evaluated in terms of forecasting time requirements. Scatterometer and radiometer applications to ocean prediction are briefly reviewed

Sampling strategies and four-dimensional assimilation of altimetric data for ocean monitoring and prediction

Numerical experiments using simulated altimeter data were conducted in order to examine the assimilation of altimeter-derived sea surface heights into numerical ocean circulation models. A reduced-gravity, primitive equation circulation model of the Gulf of Mexico was utilized; the Gulf of Mexico was chosen because of its amenability to modeling and the ability of low vertical-mode models to reproduce the observed dynamical features of the Gulf circulation. The simulated data were obtained by flying an imaginary altimeter over the model ocean and sampling the model sea surface just as real altimeter would observe the true ocean. The data were used to initialize the numerical model and the subsequent forecast was compared to the true numerical solution. Results indicate that for a stationary, circular eddy, approximately three to four tracks (either ascending or descending) across the eddy are sufficient to ensure adequate spatial resolution

Numerical simulations of rotating sunspots

A numerical model of idealized, axisymmetric, rotating sunspots is presented. The model contains a compressible plasma described by the nonlinear MHD equations, with density and temperature gradients simulating the upper layer of the sun’s convection zone. The solution forms a central flux tube in the cylindrical numerical domain, with convection cells pushing the magnetic field to the axis. When the numerical domain is rotated with a constant angular velocity, the umbra rotates as a rigid body while the surrounding convection cells show a swirling, vortical flow. As a result, the azimuthal velocity and magnetic field have their maximum values close to the flux tube, inside the innermost convection cell

Turbulent Compressible Convection with Rotation

The effects of Coriolis forces on compressible convection are studied using three-dimensional numerical simulations carried out within a local modified f-plane model. The physics is simplified by considering a perfect gas occupying a rectilinear domain placed tangentially to a rotating sphere at various latitudes, through which a destabilizing heat flux is driven. The resulting convection is considered for a range of Rayleigh, Taylor, and Prandtl (and thus Rossby) numbers, evaluating conditions where the influence of rotation is both weak and strong. Given the computational demands of these high-resolution simulations, the parameter space is explored sparsely to ascertain the differences between laminar and turbulent rotating convection. The first paper in this series examines the effects of rotation on the flow structure within the convection, its evolution, and some consequences for mixing. Subsequent papers consider the large-scale mean shear flows that are generated by the convection, and the effects of rotation on the convective energetics and transport properties. It is found here that the structure of rotating turbulent convection is similar to earlier nonrotating studies, with a laminar, cellular surface network disguising a fully turbulent interior punctuated by vertically coherent structures. However, the temporal signature of the surface flows is modified by inertial motions to yield new cellular evolution patterns and an overall increase in the mobility of the network. The turbulent convection contains vortex tubes of many scales, including large-scale coherent structures spanning the full vertical extent of the domain involving multiple density scale heights. Remarkably, such structures align with the rotation vector via the influence of Coriolis forces on turbulent motions, in contrast with the zonal tilting of streamlines found in laminar flows. Such novel turbulent mechanisms alter the correlations which drive mean shearing flows and affect the convective transport properties. In contrast to this large-scale anisotropy, small-scale vortex tubes at greater depths are randomly orientated by the rotational mixing of momentum, leading to an increased degree of isotropy on the medium to small scales of motion there. Rotation also influences the thermodynamic mixing properties of the convection. In particular, interaction of the larger coherent vortices causes a loss of correlation between the vertical velocity and the temperature leaving a mean stratification which is not isentropic

Estimating Electric Fields from Vector Magnetogram Sequences

Determining the electric field (E-field) distribution on the Sun's photosphere is essential for quantitative studies of how energy flows from the Sun's photosphere, through the corona, and into the heliosphere. This E-field also provides valuable input for data-driven models of the solar atmosphere and the Sun-Earth system. We show how Faraday's Law can be used with observed vector magnetogram time series to estimate the photospheric E-field, an ill-posed inversion problem. Our method uses a "poloidal-toroidal decomposition" (PTD) of the time derivative of the vector magnetic field. The PTD solutions are not unique; the gradient of a scalar potential can be added to the PTD E-field without affecting consistency with Faraday's Law. We present an iterative technique to determine a potential function consistent with ideal MHD evolution; but this E-field is also not a unique solution to Faraday's Law. Finally, we explore a variational approach that minimizes an energy functional to determine a unique E-field, similar to Longcope's "Minimum Energy Fit". The PTD technique, the iterative technique, and the variational technique are used to estimate E-fields from a pair of synthetic vector magnetograms taken from an MHD simulation; and these E-fields are compared with the simulation's known electric fields. These three techniques are then applied to a pair of vector magnetograms of solar active region NOAA AR8210, to demonstrate the methods with real data.Comment: 41 pages, 10 figure

On the Currents and Transports Connected With the Atlantic Meridional Overturning Circulation in the Subpolar North Atlantic

Results from an interannually forced, 0.08 degrees eddy-resolving simulation based on the Hybrid Coordinate Ocean Model, in conjunction with a small but well-determined transport database, are used to investigate the currents and transports associated with the Atlantic meridional overturning circulation (AMOC) in the subpolar North Atlantic (SPNA). The model results yield a consistent warming in the western SPNA since the early 1990s, along with mean transports similar to those observed for the trans-basin AMOC across the World Ocean Circulation Experiment hydrographic section AR19 (16.4 Sv) and boundary currents at the exit of the Labrador Sea near 53 degrees N (39.0 Sv) and east of the Grand Banks near 43 degrees N (15.9 Sv). Over a 34 year integration, the model-determined AMOC across the AR19 section and the western boundary current near 53 degrees N both exhibit no systematic trend but some long-term (interannual and longer) variabilities, including a decadal transport variation of 3-4 Sv from relatively high in the 1990s to low in the 2000s. The decadal variability of the model boundary current transport near 53 degrees N lags the observed winter time North Atlantic Oscillation index by about 2 years and leads the model AMOC across the AR19 section by about 1 year. The model results also show that the long-term variabilities are low compared to those on shorter time scales. Thus, rapid sampling of the current over long time intervals is required to filter out high-frequency variabilities in order to determine the lower frequency variabilities of interest. Citation: Xu, X., H. E. Hurlburt, W. J. Schmitz Jr., R. Zantopp, J. Fischer, and P. J. Hogan (2013), On the currents and transports connected with the atlantic meridional overturning circulation in the subpolar North Atlantic, J. Geophys. Res. Oceans, 118, 502-516, doi:10.1002/jgrc.20065

Global ocean modeling and state estimation in support of climate research

During the last decade it has become obvious that the ocean circulation shows vigorous variability on a wide range of time and space scales and that the concept of a "sluggish" and slowly varying circulation is rather elusive. Increasing emphasis has to be put, therefore, on observing the rapidly changing ocean state on time scales ranging from weeks to decades and beyond, and on understanding the ocean's response to changing atmospheric forcing conditions. As outlined in various strategy and implementation documents (e.g., the implementation plans of WOCE, AMS, CLIVAR, and GODAE) a combination of the global ocean data sets with a state-of-the-art numerical circulation model is required to interpret the various diverse data sets and to produce the best possible estimates of the time-varying ocean circulation. The mechanism of ocean state estimates is a powerful tool for such a "synthesis" of observations, obtained on very complex space-time pattern, into one dynamically consistent picture of the global time-evolving ocean circulation. This process has much in common with ongoing analysis and reanalysis activities in the atmospheric community. But because the ocean is, and will remain for the foreseeable future, substantially under-sampled, the burden put on the modeling and estimations components is substantially larger than in the atmosphere. Moreover, the smaller dynamical eddy scales which need to be properly parameterized or resolved in ocean model simulations, put stringent requirements on computational resources for ongoing and participated climate research

An Interface Region Imaging Spectrograph first view on Solar Spicules

Solar spicules have eluded modelers and observers for decades. Since the discovery of the more energetic type II, spicules have become a heated topic but their contribution to the energy balance of the low solar atmosphere remains unknown. Here we give a first glimpse of what quiet Sun spicules look like when observed with NASA's recently launched Interface Region Imaging Spectrograph (IRIS). Using IRIS spectra and filtergrams that sample the chromosphere and transition region we compare the properties and evolution of spicules as observed in a coordinated campaign with Hinode and the Atmospheric Imaging Assembly. Our IRIS observations allow us to follow the thermal evolution of type II spicules and finally confirm that the fading of Ca II H spicules appears to be caused by rapid heating to higher temperatures. The IRIS spicules do not fade but continue evolving, reaching higher and falling back down after 500-800 s. Ca II H type II spicules are thus the initial stages of violent and hotter events that mostly remain invisible in Ca II H filtergrams. These events have very different properties from type I spicules, which show lower velocities and no fading from chromospheric passbands. The IRIS spectra of spicules show the same signature as their proposed disk counterparts, reinforcing earlier work. Spectroheliograms from spectral rasters also confirm that quiet Sun spicules originate in bushes from the magnetic network. Our results suggest that type II spicules are indeed the site of vigorous heating (to at least transition region temperatures) along extensive parts of the upward moving spicular plasma.Comment: 6 pages, 4 figures, accepted for publication in ApJ Letters. For associated movies, see http://folk.uio.no/tiago/iris_spic