Enhancing Optical Gradient Forces with Metamaterials

We demonstrate how the optical gradient force between two waveguides can be enhanced using transformation optics. A thin layer of double-negative or single-negative metamaterial can shrink the interwaveguide distance perceived by light, resulting in a more than tenfold enhancement of the optical force. This process is remarkably robust to the dissipative loss normally observed in metamaterials. Our results provide an alternative way to boost optical gradient forces in nanophotonic actuation systems and may be combined with existing resonator-based enhancement methods to produce optical forces with an unprecedented amplitude.Comment: 5 pages, 4 figures; supplemental information available from AP

Improving the Ocean Initialization of Coupled Hurricane–Ocean Models Using Feature-Based Data Assimilation

Coupled hurricane–ocean forecast models require proper initialization of the ocean thermal structure. Here, a new feature-based (F-B) ocean initialization procedure in the GFDL/University of Rhode Island (URI) coupled hurricane prediction system is presented to account for spatial and temporal variability of mesoscale oceanic features in the Gulf of Mexico, including the Loop Current (LC), Loop Current eddies [i.e., warm-core rings (WCRs)], and cold-core rings (CCRs). Using only near-real-time satellite altimetry for the “SHA-assimilated” case, the LC, a single WCR, and a single CCR are assimilated into NAVOCEANO’s Global Digitized Environmental Model (GDEM) ocean temperature and salinity climatology along with satellite-derived daily sea surface temperature (SST) data from 15 September 2005 to produce a more realistic three-dimensional temperature field valid on the model initialization date (15 September 2005). For the “fully assimilated” case, both near-real-time altimetry and real-time in situ airborne XBT (AXBT) temperature profiles are assimilated into GDEM along with SST to produce the three-dimensional temperature field. Vertical profiles from the resulting SHA-assimilated and fully assimilated temperature fields are compared to 18 real-time AXBT temperature profiles, the ocean climatology (GDEM), and an alternative data-assimilated product [the daily North and Equatorial Atlantic Ocean Prediction System Best Estimate (RSMAS HYCOM), which uses an Optimal Interpolation (OI) based assimilation technique] to determine the relative accuracy of the F-B initialization procedure presented here. Also, the tropical cyclone heat potential (TCHP) from each of these profiles is calculated by integrating the oceanic heat content from the surface to the depth of the 26°C isotherm. Assuming the AXBT profiles are truth, the TCHP rms error for the F-B SHA-assimilated case, the F-B fully assimilated case, the GDEM ocean climatology, and the RSMAS HYCOM product is 12, 10, 45, and 26 kJ cm−2, respectively

Broadband metasurfaces enabling arbitrarily large delay-bandwidth products

Metasurfaces allow for advanced manipulation of optical signals by imposing phase discontinuities across flat interfaces. Unfortunately, these phase shifts remain restricted to values between 0 and 2 pi, limiting the delay-bandwidth product of such sheets. Here, we develop an analytical tool to design metasurfaces that mimic three-dimensional materials of arbitrary thickness. In this way, we demonstrate how large phase discontinuities can be realized by combining several subwavelength Lorentzian resonances in the unit cell of the surface. Our methods open up the temporal response of metasurfaces and may lead to the construction of metasurfaces with a plethora of new optical functions

A Numerical Investigation of the Local Ocean Response to Westerly Wind Burst Forcing in the Western Equatorial Pacific

Numerical simulations of the local equatorial ocean response to idealized westerly wind burst (WWB) forcing are described. In particular, the authors examine the development and evolution of the subsurface westward jet (SSWJ) that has been observed to accompany these wind events. This westward current is interpreted as the signature of equatorial waves that accompany the downwelling and upwelling that occurs along the edges of the wind forcing region. Some important features of the SSWJ include maximum intensity toward the eastern edge of the forcing region, a time lag between the wind forcing and peak SSWJ development, and an eastward spreading of the SSWJ with time. The effect of wind burst zonal profile, magnitude, duration, and fetch on the SSWJ are explored. The response of an initially resting ocean to WWB forcing is compared with that for model oceans that are spun up with annual-mean surface fluxes and monthly varying fluxes. It is demonstrated that the gross features of the response for the spun up simulations can be well approximated by adding the background zonal current structure prior to the introduction of the wind burst to the initially resting ocean current response to the WWB. This result suggests that the zonal current structure that is present prior to the commencement of WWB forcing plays a key role in determining whether or not a SSWJ will develop

Design of a Multiply Nested Primitive Equation Ocean Model

A new multiply nested primitive equation ocean model is presented. The model employs a two-way interactive nesting technique successfully applied for many years in the Geophysical Fluid Dynamics Laboratory–NOAA hurricane prediction model. The formulation of the mesh nesting algorithm allows flexibility in deciding the number of meshes and the ratio of grid resolutions between adjacent meshes. Other advanced features include realistic coastline geometry and spatially variable grid spacing. The results of various idealized experiments indicate good performance of the nesting technique. The most important feature of the model is the ability to combine large-scale and regional-scale predictions. The model is tested as a general circulation model (GCM) in a 3-yr spinup experiment of the large-scale circulation in the tropical Pacific Ocean. It demonstrates skill comparable to that of other recently developed GCMs. The resulting large-scale fields are then used in the nested configuration as initial conditions for simulations of the ocean response to a westerly wind burst and a tropical cyclone. Significant improvements over a coarse, single-mesh model have been achieved in resolving finescale features of the wind-induced current and temperature fields. These results highlight the importance of model resolution for realistic simulations of mesoscale ocean variability

New parameterization of air-sea exchange coefficients and its impact on intensity prediction under major tropical cyclones

Understanding and quantifying air-sea exchanges of enthalpy and momentum fluxes are crucial for the advanced prediction of tropical cyclone (TC) intensity. Here, we present a new parameterization of air-sea fluxes at extreme wind speeds from 40 m s−1 to 75 m s−1, which covers the range of major TCs. Our approach assumes that the TC can reach its maximum potential intensity (MPI) if there are no influences of external forces such as vertical wind shear or other environmental constraints.This method can estimate the ratio of the enthalpy and momentum exchange coefficient (Ck/Cd) under the most intense TCs without direct flux measurements. The estimation showed that Ck/Cd increases with wind speed at extreme winds above 40 m s−1. Two types of surface layer schemes of the Hurricane Weather and Research Forecast (HWRF) were designed based on the wind speed dependency of the Ck/Cd found at high winds: (i) an increase of Ck/Cd based on decreasing Cd (Cd_DC) and (ii) an increase of Ck/Cd based on increasing Ck (Ck_IC). The modified surface layer schemes were compared to the original HWRF scheme (using nearly fixed Cd and Ck at extreme winds; CTRL) through idealized experiments and real-case predictions. The idealized experiments showed that Cd_DC reduced frictional dissipation in the air-sea interface as well as significantly reduced sea surface cooling, making the TC stronger than other schemes. As a result, Cd_DC reduced the mean absolute error and negative bias by 15.0% (21.0%) and 19.1% (32.0%), respectively, for all lead times of Hurricane Irma in 2017 (Typhoon Mangkhut in 2018) compared to CTRL. This result suggests that new parameterization of Ck/Cd with decreasing Cd at high winds can help improve TC intensity prediction, which currently suffers from underestimating the intensity of the strongest TCs

Nonbreaking wave‐induced mixing in upper ocean during tropical cyclones using coupled hurricane‐ocean‐wave modeling

The effects of turbulence generated by nonbreaking waves have been investigated by testing and evaluating a new nonbreaking wave parameterization in a coupled hurricane‐ocean‐wave model. The MPI version of the Princeton Ocean Model (POM) with hurricane forcing is coupled with the WAVEWATCH‐III (WW3) surface wave model. Hurricane Ivan is chosen as the test case due to its extreme intensity and availability of field data during its passage. The model results are validated against field observations of wave heights and sea surface temperatures (SSTs) from the National Data Buoy Centre (NDBC) during Hurricane Ivan and against limited in situ current and bottom temperature data. A series of numerical experiments is set up to examine the influence of the nonbreaking wave parameterization on the mixing of upper ocean. The SST response from the modeling experiments indicates that the nonbreaking wave‐induced mixing leads to significant cooling of the SST and deepening of the mixed layer. It was found that the nondimensional constant b1 in the nonbreaking wave parameterization has different impacts on the weak and the strong sides of the storm track. A constant value of b1 leads to improved predictions on the strong side of the storm while a steepness‐dependent b1 provides a better agreement with in situ observations on the weak side. A separate simulation of the intense tropical cyclone Olwyn in north‐west Australia revealed the same trend for b1 on the strong side of the tropical cyclone

Tunable terahertz frequency comb generation using time-dependent graphene sheets

We investigate the interaction between electromagnetic pulses and two-dimensional current sheets whose conductivity is controlled as a function of time by the generation of photocarriers, and we discuss its applicability to tunable frequency comb generation. To this aim, we develop an analytical model that permits the calculation of the scattered waves off a thin sheet with time-dependent, dispersive sheet conductivity. We evaluate the transmitted spectrum as a function of the dispersive behavior and the modulation frequency of the number of photocarriers. We conclude that such active materials, e.g., time-dependent graphene sheets, open up the possibility to manipulate the frequency of incident pulses and, hence, could lead to highly tunable, miniaturized frequency comb generation

Description and Analysis of the Ocean Component of NOAA’s Operational Hurricane Weather Research and Forecasting (HWRF) Model

The Princeton Ocean Model for Tropical Cyclones (POM-TC), a version of the three-dimensional primitive equation numerical ocean model known as the Princeton Ocean Model, was the ocean component of NOAA’s operational Hurricane Weather Research and Forecast Model (HWRF) from 2007 to 2013. The coupled HWRF–POM-TC system facilitates accurate tropical cyclone intensity forecasts through proper simulation of the evolving SST field under simulated tropical cyclones. In this study, the 2013 operational version of HWRF is used to analyze the POM-TC ocean temperature response in retrospective HWRF–POM-TC forecasts of Atlantic Hurricanes Earl (2010), Igor (2010), Irene (2011), Isaac (2012), and Leslie (2012) against remotely sensed and in situ SST and subsurface ocean temperature observations. The model generally underestimates the hurricane-induced upper-ocean cooling, particularly far from the storm track, as well as the upwelling and downwelling oscillation in the cold wake, compared with observations. Nonetheless, the timing of the model SST cooling is generally accurate (after accounting for along-track timing errors), and the ocean model’s vertical temperature structure is generally in good agreement with observed temperature profiles from airborne expendable bathythermographs