Data from a network of high-frequency (HF) beacons deployed in Peru are used to estimate the regional ionospheric electron density in a volume. Pseudorange, accumulated carrier phase, and signal power measurements for each of the 36 ray paths provided by the network at a 1 min cadence are incorporated in the estimates. Additional data from the Jicamarca incoherent scatter radar, the Jicamarca sounder, and GPS receivers can also be incorporated. The electron density model is estimated as the solution to a global optimization problem that uses ray tracing in the forward model. The electron density is parametrized in terms of B-splines in the horizontal direction and generalized Chapman functions or related functions in the vertical. Variational sensitivity analysis has been added to the method to allow for the utilization of the signal power observable which gives additional information about the morphology of the bottomside F region as well as absorption including absorption in the D and E regions. The goal of the effort is to provide contextual information for improving numerical forecasts of plasma interchange instabilities in the postsunset F region ionosphere associated with equatorial spread F (ESF). Data from two ESF campaigns are presented. In one experiment, the HF data revealed the presence of a large-scale bottomside deformation that seems to have led to instability under otherwise inauspicious conditions. In another experiment, gradual variations in HF signal power were found to be related to the varying shape of the bottomside F layer.