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 ab...