The evolving upper ocean response excited by the passage of hurricane Gilbert (September 14–19, 1988) was investigated using current and temperature observations acquired from the deployment of 79 airborne expendable current profilers (AXCPs) and 51 airborne expendable bathythermographs from the National Oceanic and Atmospheric Administration WP‐3D aircraft in the western Gulf of Mexico. The sea surface temperatures (SSTs), mixed layer depths, and bulk Richardson numbers were objectively analyzed to examine the spatial variability of the upper ocean response to Gilbert. Net decreases of the SSTs of 3°–4°C were observed by the profilers as well as by the airborne infrared thermometer (AIRT) along the flight tracks and advanced very high resolution radiometer (AVHRR) imagery. The AXCPs indicated a marked cooling from 29°C to about 25.5°C on September 17, 1988, which was about 1.2 inertial periods (IP) following storm passage. This pool of cooler water (3.5°) was located further downstream in the hurricane wake by September 19 (2.7 IP following the storm) as a result of the near‐inertial currents in the mixed layer. While there was a bias of about 0.6°C and 1.7°C between the in situ and AVHRR‐derived SSTs, respectively, both the AVHRR images and the objectively analyzed fields indicated a rightward bias in the upper ocean cooling that extended from the storm track to about 4Rmax (where Rmax, the radius of maximum winds, is equal to 50 km). The larger SST offset of 1.7°C was due to the difference between the time of the AVHRR image and the time of the aircraft experiment on September 19. The SSTs derived from the AVHRR images and the AIRT also indicated large gradients between the cold wake and the warm eddy in the central Gulf of Mexico. The mixed layer deepened by about 30–35 m on the right side of the track during the storm and 1.2 IP later, with little evidence of continued deepening afterward. The mixed layer current vectors demonstrate that a strong, near‐inertially rotating current was excited by the passage of Gilbert, with maxima of about 1–1.4 m s−1. The currents, observed during and subsequent to (1.2 IP) the storm, diverged from the storm track, whereas the mixed layer current vectors 2.7 IP after storm passage converged toward the track, with relative maxima of 0.8–1 m s−1. This alternating pattern of convergence and divergence of the mixed layer current was associated with the upwelling and downwelling cycles of the baroclinic response. Considerable current shear existed between the mixed layer and the thermocline currents in the cool wake between the storm track and the 4Rmax. Estimates of the bulk Richardson numbers ranged between 0.2 and 1.0 during Gilbert and at 1.2 IP, which suggests that enhanced current shears were responsible for some of the mixed layer deepening.